Research Article  Open Access
KeunHyeok Yang, Yongjei Lee, JuHyun Mun, "A StressStrain Model for Unconfined Concrete in Compression considering the Size Effect", Advances in Materials Science and Engineering, vol. 2019, Article ID 2498916, 13 pages, 2019. https://doi.org/10.1155/2019/2498916
A StressStrain Model for Unconfined Concrete in Compression considering the Size Effect
Abstract
In this study, a stressstrain model for unconfined concrete with the consideration of the size effect was proposed. The compressive strength model that is based on the function of specimen width and aspect ratio was used for determining the maximum stress. In addition, in stressstrain relationship, a strain at the maximum stress was formulated as a function of compressive strength considering the size effect using the nonlinear regression analysis of data records compiled from a wide variety of specimens. The descending branch after the maximum stress was formulated with the consideration of the effect of decreasing area of fracture energy with the increase in equivalent diameter and aspect ratio of the specimen in the compression damage zone (CDZ) model. The key parameter for the slope of the descending branch was formulated as a function of equivalent diameter and aspect ratio of the specimen, concrete density, and compressive strength of concrete. Consequently, a rational stressstrain model for unconfined concrete was proposed. This model reflects trends that the maximum stress and strain at the peak stress decrease and the slope of the descending branch increases, when the equivalent diameter and aspect ratio of the specimen increase. The proposed model agrees well with the test results, irrespective of the compressive strength of concrete, concrete type, equivalent diameter, and aspect ratio of the specimen.
1. Introduction
The stressstrain relationship for unconfined concrete is a fundamental material property for the design and analysis of structural elements [1–3]. Generally, in stressstrain relationship, the ascending and descending branches are dependent on concrete type and compressive strength, as well as the maximum diameter of aggregate, specimen width or diameter, and aspect ratio in the descending branch from the crack propagation in the fracture zone [4–7]. The slope of the ascending branch commonly increases with the increase in compressive strength of concrete and the decrease in specimen width or diameter and aspect ratio, while that of the descending branch increases with concrete compressive strength, specimen width or diameter, and aspect ratio. To study this trend, researchers [5, 8–10] proposed concrete compressive strength models with the size effect in various approaches based on the fracture energy theory. Bažant and Planas [8] reported that concrete compressive strength was considerably affected by the specimen width or diameter, indicating that it decreased by 10%, when the specimen width or diameter increased twice. Sim et al. [5] emphasized that the size effect on compressive strength for lightweight concrete (LWC) was more notable than that for normal weight concrete (NWC). In particular, because cracks at the failure zone for LWC pass through lightweight aggregate particles, the crack band zone is more localized in LWC than NWC. Hence, the size effect of concrete on the peak stress and descending branch behavior that is directly related to crack propagations in the failure zone could be more notable in LWC than those in NWC [5]. However, for the size effect of concrete on the descending branch from the crack propagation in a localized crack band zone, only few studies have been conducted. In particular, very little literature on the size effect in LWC is available. Furthermore, the existing proposed models [6, 7] for the stressstrain relationship regarding the size effect of concrete on the descending branch are typically determined from NWC test data records, rather than from LWC, with limited ranges of variables.
Markeset and Hillerborg [6] generalized the compression damage zone (CDZ) model to consider the size effect of concrete on the descending branch using a function of strain dissipated by the shear band including the fracture energy. However, the descending branch behavior in Markeset and Hillerborg’s model are drawn from the limited ranges of variables. In addition, this model is limited for a practical equation because the specific information about the strain of the starting point for the softening behavior is not available. Samani and Attard [7] considered the size effect on the descending branch in stressstrain relationship proposed by Attard and Setunge’s model [11], using the CDZ model [6]. As the Samani and Attard’s [7] model has an identical strain model of shear band as that of Markeset and Hillerborg’s [6], information regarding the material property factors is necessary to predict the stressstrain relationship. In addition, the descending branch of these models does not fully consider the effect of aggregate property on crack propagation and the localized fracture zone. For example, in LWC, the contribution of stress transfer at the crack plane for aggregate interlocking action is little [5] because most of the cracks at the failure plane pass through lightweight aggregate particles. In addition, the strength and elastic modulus of aggregates such as magnetite in heavyweight concrete (HWC) are typically higher than those in NWC, which can cause a wider fracture zone by crack propagation. Hence, the size effect of LWC and HWC could be different from that of NWC because the size effect on concrete depends on the areas of failure and crack propagation in the fracture zone. However, as the existing models [6, 7] incorporating the size effect of concrete on the descending branch were derived from limited NWC test data records, a limited number of studies for the size effect of concrete using lightweight and heavyweight aggregates are available.
The objective of this study is to propose a model for the stressstrain curve considering the size effect of various concrete types. In this model, the basic formula for the stressstrain curve and the key parameter that determines the slopes of the ascending and descending branches established by Yang et al. [1] was used to generate a complete nonlinear curve. The concrete compressive strength model proposed by Sim et al. [5] that considers the size effect was used for the peak stress. The strain at the peak stress was generalized with a simple equation using the regression analysis of data records compiled from specimens with various ranges of equivalent diameter and aspect ratio of specimen, compressive strength, and density. In the softening behavior of the descending branch in the stressstrain relationship, the size effect and fracture energy were considered using Markeset and Hillerborg’s CDZ model [6]. The key parameter for the softening behavior was determined from the secant modulus joining the origin and , where is the compressive strength of concrete considering the size effect. The strain model at to determine the key parameter was generalized with various ranges of variables using nonlinear regression analysis. Finally, the key parameter that determines the slope of the descending branch was formulated as a function of equivalent diameter and aspect ratio of specimen, concrete density, and concrete compressive strength, using parametric numerical analysis. The accuracy of the proposed model was evaluated using a normalized rootmeansquare error obtained from the comparisons of predicted curves with the test results.
2. Database
To formulate the material properties, Yang et al. [1] compiled 3295 data records for the elastic modulus of concrete (), 415 data records for strain at the peak stress (), and 96 data records for strain at 50% of the peak stress () in the descending branch. In the compiled data records, the numbers available for the concrete compressive strength () and density () were 3295 and 3274, respectively, and varied from 8.4 MPa to 170 MPa and from 1200 kg/m^{3} to 4500 kg/m^{3}, respectively. To apply the effect of concrete type to the empirical formulations, the data records were divided into LWC, NWC, and HWC according to . The concrete density () varied from 1200 kg/m^{3} to 2000 kg/m^{3} for LWC, from 2000 kg/m^{3} to 2500 kg/m^{3} for NWC, and from 2500 kg/m^{3} to 4500 kg/m^{3} for HWC. The data records compiled by Yang et al. [1] were all measured from a cylinder of diameter 100 mm and height 200 mm. Hence, the data records compiled by Yang et al. [1] lack the test results for different specimen size. To compensate for this, additional test results for the size effect of concrete were compiled. The compiled additional test results [4, 5, 12–15] for the stressstrain relationship were 38 data records for LWC and 26 data records for NWC, as shown in Table 1. To consider various specimen shapes, an equivalent diameter was introduced, where and are the sectional width and depth of prism specimen, respectively, assuming that the area of the prism specimens is identical to that of cylinder specimen. For example, for prism specimens with a of 300 mm and of 400 mm, corresponds to 398.9 mm. In the data records, , aspect ratio (), , and were varied from 50 mm to 800 mm, 0.5 to 8, 4 mm to 20 mm, 17.1 MPa to 90.2 MPa, and 1500 kg/m^{3} to 2464 kg/m^{3}, respectively. The specimen equivalent diameter () and were varied from 100 mm to 350 mm and 1 to 2 for LWC and 50 mm to 800 mm and 0.5 to 8 for NWC, respectively.

3. Model Generalization
3.1. Basic Approach
The stressstrain relationship for unconfined concrete in compression is a parabola with ascending and descending branches, and a vertex at the peak stress [1–3]. This shape can be generalized using the following equation [1]:where is the normalized stress, is the normalized strain, and are the concrete stress and strain at some point in stressstrain curve, respectively, and is the key parameter that determines the slopes of the ascending and descending branches. The ascending branch can be determined from , which is defined as the slope of the line joining the origin and 40% of the peak stress [16]. In addition, the descending branch can be determined from the secant modulus joining the origin and [1]. In accordance with Yang et al.’s model [1], the equation for the key parameter that determines the slopes of the ascending and descending branches can be expressed in the following forms:where and .
Bažant [9] proposed the crack band zone by a crack width with microcrack propagation for concrete failure. In addition, Bažant [9] idealized the crack band width as a function of , assuming that microcracks propagated the interfaces between aggregates and pastes. Sim et al. [5] proposed a smaller area of the crack band zone in LWC than NWC as shown in Figure 1, based on the crack band theory [9]. This model includes the effect of reduced area of the crack band zone caused by the cracks at the failure zone passing through lightweight aggregate particles and also considers the size effect on concrete due to the decrease in as increases, where is the crack length. The above previous models clearly revealed that the size effect exists in compressive strength of concrete due to different propagations of longitudinal splitting cracks at different specimen sizes, although further elaborated analytical approach would be needed to account for the effect of splitting cracks on the size effect in different concrete types. Sim et al. [5] derived the equation for the compressive strength of concrete () considering the size effect from the energy balance up to the peak stress in the crack band zone idealizing the propagation of longitudinal splitting cracks, in that the strain energy for concrete deformation dissipated by the crack band zone equals the total energy consumed by the band of the axial microsplitting cracks. The proposed model is as follows:where is , is the number of microcracks in the band, is , is , is , is , and are the modification factors to account for the volume of the crack band zone, is the strainsoftening modulus, and , , , and are the experimental constants. In equation (3), indicates the compressive strength of concrete measured in a reference specimen with of 150 mm and of 2. From equation (3), this indicates that is considerably affected by the functions of , , and [5, 8–10]. Sim et al. [5] determined the functions of , , , , , and in equation (3) from 1509 data records with LWC and NWC test results and proposed as follows:
In the model of Sim et al. [5], equation (4) considers a function of that reflects the trend that the size effect is more notable in LWC than that in NWC. The stressstrain relationship is shown in Figure 2, using equation (4) for the peak stress. In Figure 2, the key parameter can be determined with the following equations:where , , is the strain at the peak stress considering the size effect, and is the strain at after the peak stress. In equation (1), the key parameter that determines the ascending and descending branches requires information about the functions of , , and , as expressed in equations (5) and (6).
3.2. Determination of in Ascending Branch
The size effect on concrete is based on the crack band theory indicating the crack width with the propagation of microcracks, as shown in Figure 3. It implies that the behavior of the points of “O” and “A” without any cracks in the graph is not affected by the size effect. In addition, Taylor and Broms [17] reported that the bond cracks can be definitely observed between 38 and 42% of the peak stress at the ascending branch of the stressstrain curve. As a result, the strain at the ascending branch begins to increase rapidly with the propagation of the bond cracks, showing the nonlinear curve, as pointed out by Neville [16]. According to fracture mechanics, the size effect on concrete is primarily caused by the crack propagation. Thus, the present study considered that representing the quasilinear relationship at the ascending branch is marginally affected by the size effect. Yang et al. [1] formulated the empirical formulation for on the basis of the regression analysis of 2680 data records for NWC, 370 data records for LWC with ranging from 1200 to 2000 kg/m^{3}, and 245 data records for HWC with ranging from 2500 to 4450 kg/m^{3}. Noguchi et al. [18] reported that can be expressed as a function of and and included correction factors to account for the effects of aggregate type and presence of supplementary cementitious materials (SCMs). In Yang et al.’ model [1], the various unusual aggregates such as artificial lightweight aggregates and heavyweight magnetite particles were implicitly considered in the use of the parameter of . In addition, the correction factor for SCMs can be implicitly included when the empirical constants are obtained from regression analysis, resulting in negligible errors. Yang et al.’s model reflected the trend that a lower increasing rate in than that in was considered by using a power function of to represent the nonlinear relationship between the two parameters. Consequently, Yang et al.’s model [1] was used for as follows:
(a)
(b)
Meanwhile, as shown in Figure 2, equation (4) can be expressed as follows from the relation of and :
Consequently, can be expressed as follows:where is the secant modulus joining the origin and the peak stress. Equation (9) shows that is fully affected by the functions of and . However, data records or predicted model for is not available in the literature. Hence, in this study, can be proposed as follows using a certain relation with , as shown in Figure 2:where is a coefficient to account for the relation of and , which can be determined from the test results. From the data records compiled in this study, equation (10) can be proposed as follows (Figure 4):
The key parameter that determines the slope of the ascending branch can be solved by substituting equations (7) and (11) into equation (5). The key parameter was calculated using the Newton–Raphson method, identical to Yang et al.’s model [1]. The values determined for different concrete properties need to be formulated as a simple equation for practical application of the proposed stressstrain relationship of concrete. Thus, the present research involved a parametric study to generalize under the comprehensive ranges of parameters as follows: between 10 MPa and 180 MPa; between 1400 kg/m^{3} and 4000 kg/m^{3}; between 50 mm and 500 mm; and between 0.5 and 5. Note that the geometrical conditions of specimens were considered in the parametric study because the given parameters in equations (4) and (11) are affected by the equivalent width and aspect ratio of the specimen. From the regression analysis using the solutions obtained from the Newton–Raphson method, the key parameter to account for the slope of the ascending branch can be simply formulated as follows (Figure 5):
Equation (12) shows that the slope of the ascending branch increases with the increase in or decrease in . Consequently, the slope of the ascending branch includes the size effect of concrete with the generalization of a function of , considering the parameters and .
3.3. Determination of in Descending Branch
The compressive failure behavior of unconfined concrete is commonly characterized by the mode I (pure tensile) and mode III (sliding displacement due to diagonal shear) in estimating the fracture zone. As pointed out by Markeset and Hillerborg, microsplitting cracks (mode I) are formulated due to Poisson’s effect under pure compressive stresses, whereas sliding displacement along diagonal tensile cracks (mode III) occurs at 45degree slope relative to the principal normal stresses when the maximum shearing stresses reach the shear capacity of concrete. Hence, Markeset and Hillerborg assumed that longitudinal crack zones of concrete relate to the tensile fracture energy zone, whereas diagonal shear crack zones after the peak stress identify the shear fracture energy zone (Figure 6). In fact, most of cracks in concrete are caused by tensile stresses rather than compressive stresses because the tensile resistance of concrete is highly lower than the compressive resistance. Overall, the compressive failure of concrete is commonly caused by the tensile fracture due to Poisson’s effect. The descending branch behavior after the peak stress is determined from the localized deformation developed in the damaged or failure zone, while the undamaged zone elastically unloads [6]. Hence, the undamaged zone is generated only when is greater than the damaged zone height (). For compressive behavior, Markeset and Hillerborg [6] idealized the CDZ model generated from the longitudinal microsplitting cracks and localized diagonal tensile shear band crack in the damaged zone (Figure 6). In CDZ model, microsplitting cracks can be idealized as a crack band zone with several microsplitting cracks because their propagation requires energy release. Thus, the size effect of concrete on compression is expected as demonstrated by lots of previous researches. Bažant and Planas also idealized the fracture energy zone of concrete under compression, considering the tensile and diagonal shear cracks to consider the size effect. According to Markeset and Hillerborg [6], the total strain in the softening behavior is the sum of the strain during the unloading region after the peak stress in the undamaged zone, the strain while microsplitting crack occurs in the longitudinal direction, and the strain caused by the diagonal shear band ():where is the height of the region propagated by the longitudinal microsplitting cracks. is the localized deformation assumed between 0.4 and 0.7 for NWC and less than 0.3 for LWC. The longitudinal microsplitting cracks in unconfined concrete under compression commonly develop at approximately 75∼90% of the peak stress. From this finding, in the CDZ model, the amount of energy () released in the unloading zones can be obtained using the following equation:where is the fracture energy and is the factor based on the material properties. is the factor to account for average spacing between longitudinal spitting microcracks. Assuming that the strain () in the region propagated by the longitudinal microsplitting cracks is proportional to the tensile fracture energy (), it can be proposed as follows:
As expressed in equation (15), the CDZ model proposed by Markeset and Hillerborg [6] includes a function of , in the descending branch. However, , , and require calibration according to various concrete types because they are based on the material properties, which are too demanding for a practical application. Furthermore, because is proposed only for a of 16 mm, the use of a practical equation is limited for other specimens with larger aggregate. Hence, to obtain information about these factors, a comprehensive test is required with various influencing parameters including concrete type, , , and . To improve Markeset and Hillerborg’s model [6], the key parameter by Yang et al. [1] was applied to the descending branch behavior. The peak stress from equation (4) by Sim et al. [5] and from equations (13a) and (13b) are used to produce the following equation:where , , and are known values. In equations (16a) and (16b), the first term in the right is moved to the left side and can be arranged as follows:where the values of and are experimental constants in predicting the softening in the CDZ model. On the basis of test results, Markeset and Hillerborg assumed the value of as 3.0 for NWC and 1.0 for LWC. However, the value of for HWC is still unknown because of the lack of test data. Markeset and Hillerborg also introduced the factor to account for the average spacing of the longitudinal microsplitting cracks due to the primary tensile stresses. They assumed the value of as 1.25 for the maximum aggregate size of 16 mm. However, there is no further information on the value of for different aggregate sizes although the spacing of the longitudinal microsplitting cracks can be significantly affected by the aggregate size due to the aggregate interlock action. In addition, the value of depends on the equivalent width and aspect ratio of the specimen because the energy release at the crack band zone is affected by the spacing of the longitudinal microsplitting cracks.
The present study conducted the regression analysis of test data on to simply generalize the righthand side of equations (17a) and (17b) including the factors and . Test results of according to , , and are shown in Figure 7. In Figure 7(a), nonlinearly decreased with the increase in , indicating that it decreased by approximately 18% for specimens with of 2, and approximately 31% for specimens with of 1, when increased by 3 times. almost linearly decreased with increased, irrespective of . ranged between 0.0026 and 0.0049 for specimens with of 2 and 0.0018 for specimens with of 5.5, indicating that it decreased by 43%, when increased by 3 times. In addition, increased with the increase in , and its increasing rate according to was almost constant. ranged between 0.002 and 0.0038 for LWC ( less than 2000 kg/m^{3}) and 0.0029 to 0.0044 for HWC ( more than 2500 kg/m^{3}). These imply that the descending branch behavior in the stressstrain relationship for unconfined concrete is considerably affected by the functions of , , and . Based on this analysis, was generalized as functions of , , , and (Figure 8), using regression analysis from the test results [4, 5, 12–15, 19–29] for 45 data records for LWC, 91 data records for NWC, and 24 data records for HWC:where the CEBFIP model [30] for that includes functions of and was expressed as follows:where is 0.025 N/mm, 0.03 N/mm, and 0.05 N/mm for of 8 mm, 16 mm, and 32 mm, respectively. Overall, the present model can predict the softening performance of concrete with different parameters including the compressive strength and density of concrete, equivalent width and aspect ratio of specimens, and aggregate sizes even though the values for and are not determined from test specimens. The value of in the descending branch can be solved using equations (11) and (18). The solution of in the descending branch was also calculated using the Newton–Raphson method as in the ascending branch. Finally, the key parameter was formulated using the analytical parametric study. In the analytical parametric study, , , , , and were selected from 10 MPa to 180 MPa, 50 mm to 500 mm, 0.5 to 5, 4 mm to 25 mm, and 1400 kg/m^{3} to 4000 kg/m^{3}, respectively. From the analytical results, statistical optimization was performed to generalize the key parameter that determines the slope of the descending branch as follows (Figure 9):
(a)
(b)
(c)
Finally, the stressstrain relationship for unconfined concrete can be proposed as follows:where is given by equation (11), is given by equation (4), and key parameter is given by equation (12) or (20). The proposed stressstrain relationship for unconfined concrete can consider the size effect on concrete in the ascending and descending branches, using the power functions of the key parameters and .
4. Comparisons with Test Results
The test results compiled from the available literatures [4, 5, 12–15, 19–29] were compared with predictions of this study and the existing models [1, 6, 7, 11]. The existing models for the strainstress relationship proposed by Markeset and Hillerborg [6] and Samani and Attard [7] were selected as summarized in Table 2. Figure 10 shows comparisons of the predicted and measured stressstrain curves [4, 19–25]. The comparative analysis focused on the effect of , , , and on the stressstrain curve. Table 3 summarizes the normalized rootmeansquare error (NRMSE) obtained from the comparisons of test results with predictions. In Table 3, and are the mean and standard deviation of the NRMSE, respectively. It is noteworthy that the comparisons of test results with the predictions of Markeset and Hillerborg [6] were conducted only in the descending branch because Markeset and Hillerborg’s model [6] provides the equations only for the descending branch behavior.
 
, , , and are in MPa; , and are in mm; is in N/mm; is in kg/m^{3}. 
(a)
(b)
(c)

Markeset and Hillerborg [6] idealized the CDZ model considering fracture energy and proposed the descending branch behavior in the stressstrain relationship. In the CDZ model, the total strain is a combination of the strain () in the region where the undamaged zone elastically unloads after the peak stress, the strain () in the region propagated by the longitudinal microsplitting cracks, and the strain () by the diagonal tensile band crack. In the model, is based on the assumption that inelastic deformation in the damaged zone determines the descending branch behavior. To consider inelasticity in the descending branch, introduces as that expressed in equation (8); nevertheless, the descending branch is predicted as a virtually linear curve (Figure 10(a)). However, the shapes of the descending branch measured in the existing test results are primarily curved rather than linear. Hence, the accuracy of Markeset and Hillerborg’s model [6] according to the concrete type fluctuates with large deviations. The values of obtained by Markeset and Hillerborg’s model [6] are 0.39 for LWC, 0.29 for NWC, and 0.29 for HWC. In particular, Markeset and Hillerborg’s model [6] underestimates the compressive stress of HWC. This is because this model does not consider the decreasing effect of the slope for HWC in the descending branch because the factor of to determine according to concrete type is 3 and 1, only for NWC and LWC, respectively, without considering HWC.
Samani and Attard [7] applied the size effect to the equations for the descending branch proposed by Attard and Setunge’s model [11]. As summarized in Table 2, the descending branch behavior proposed by Samani and Attard [7] is also based on the CDZ model and its strain equation is similar to that of Markeset and Hillerborg’ s model [6]. In addition, Markeset and Hillerborg’s model [6] for composed of functions of , , and is used without modification. The descending branch behavior in Samani and Attard’s model [7], however, is different from that in Markeset and Hillerborg [6], indicating that it is predicted as a curve with an inflection point. As shown in Figure 10(b), the ascending branch behavior composed of functions and in Samani and Attard’s model [7] is identical to that in Attard and Setunge’s model [11]. In the descending branch behavior of Samani and Attard’s model [7], the increasing effect of the decreasing slope is explained well with the increase in , as shown in Figure 10(b). However, Samani and Attard’s model [7] underestimates the compressive stress for HWC. This is because, in this model, the factor related with the material property does not consider the decreasing effect of the descending slope for HWC. This implies that the factor in Samani and Attard’s model [7] requires calibration using additional test results. In addition, because Samani and Attard’s model [7] does not consider the size effect on the peak stress, it overestimates the compressive strength of concrete for LWC. The overestimation increased with the increase in . The values of obtained by Samani and Attard’s model [7] are 0.41 for LWC, 0.23 for NWC, and 0.73 for HWC.
The proposed model in this study shows better agreement with test results, irrespective of , , concrete type, and . The values of and are 0.24 and 0.16 for LWC, 0.19 and 0.12 for NWC, and 0.10 and 0.01 for HWC, respectively. The results are lower than those of the models of Markeset and Hillerborg [6] and Samani and Attard [7]. The proposed values of and are 0.19 and 0.13, respectively, which are the lowest among other models. Based on the CDZ model, a rational stressstrain model for unconfined concrete considering the size effect is proposed, using the key parameter formulated by functions of , , , and . Note that most tests to investigate the stressstrain curves of concrete in compression were conducted using standard cylindrical specimens of 100 × 200 mm or 150 × 300 mm. Moreover, very few specimens with the equivalent diameter exceeding 200 mm are available in the literatures because of the capacity limitation of the testing machine. Thus, the proposed models need to be further examined in LWC and HWC specimens with a larger size.
5. Conclusions
From the proposed stressstrain relationship model for various unconfined concrete types considering the size effect based on the CDZ model, the following conclusions were derived:(1)Although the concrete commonly has microcracks, the elastic modulus of concrete typically defined as where bond cracks occurred in was not affected by the size effect, whereas the strain at the peak stress was affected by the size effect because of propagation of cracks.(2) that closely related to the slopes of descending branches in stressstrain relationship decreased averagely by 25% and 43%, respectively, when equivalent diameter and aspect ratio increased by 3 times. The corresponding values for lightweight concrete (LWC) were lower than those for normal weight concrete (NWC) and heavyweight concrete (HWC).(3)The key parameter determining the slope of the ascending and descending branches could be proposed as an exponential function of and , respectively.(4)The proposed model of the stressstrain relationship for unconfined concrete showed good agreements with the test results, irrespective of equivalent diameter and aspect ratio of specimen, concrete density, and compressive strength.
Notations
:  Equivalent diameter 
:  Maximum size of aggregate 
:  Height of specimen 
:  Aspect ratio of specimen 
:  Reference height of specimen 
:  Damage zone height of specimen 
:  Secant modulus at the peak stress 
:  Elastic modulus of concrete 
:  Strainsoftening modulus 
:  Stress at stressstrain curve 
:  Compressive strength of concrete measured in the standard specimen 
:  Inflection stress in descending branch 
:  Compressive strength concrete considering the size effect 
:  Fracture energy 
:  Factor relating to material property 
:  Conversion coefficient 
:  Number of microcracks in the band 
:  Amount of energy released in the unloading zones 
, , , and :  Experimental constants 
, :  Modification functions to account for the volume of the crack band zone 
:  Key parameter that determines the slope of ascending and descending branches 
:  Factor relating to height of specimen 
:  Mean of normalized rootmeansquare error 
:  Standard deviation of normalized rootmeansquare error 
:  Strain induced from elastically unloading in the undamaged zone 
:  Strain at after the peak stress 
:  Stain at the peak stress 
:  A strain at stressstrain curve 
:  Strain in the damaged zone relating to longitudinal microsplitting cracking 
:  Total strain occurred in compression damage zone (CDZ) 
:  Strain at inflection stress in descending branch 
:  Strain at the peak stress considering the size effect 
:  Strain at after the peak stress 
:  Factor relating to type of coarse aggregate 
:  Concrete density 
:  Coefficient to account for the relation of and 
:  Localized deformation 
:  Strain of diagonal tensile shear band. 
Data Availability
The data records used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This research was supported by the research grant from Kyonggi University through the Korea Agency for Infrastructure Technology Advancement funded by the Ministry of Land, Infrastructure and Transport of the Korean Government (Project no. 19TBIPC12647003).
References
 K. H. Yang, J. H. Mun, M. S. Cho, and T. H.K. Kang, “Stressstrain model for various unconfined concretes in compression,” ACI Structural Journal, vol. 111, no. 4, pp. 819–826, 2014. View at: Publisher Site  Google Scholar
 T. H. Almusallam and S. H. Alsayed, “Stressstrain relationship of normal, highstrength and lightweight concrete,” Magazine of Concrete Research, vol. 47, no. 170, pp. 39–44, 1995. View at: Publisher Site  Google Scholar
 B.I. Bae, H.K. Choi, B.S. Lee, and C.H. Bang, “Compressive behavior and mechanical characteristics and their application to stressstrain relationship of steel fiberreinforced reactive powder concrete,” Advances in Materials Science and Engineering, vol. 2016, Article ID 6465218, 11 pages, 2016. View at: Publisher Site  Google Scholar
 D. C. Jansen and S. P. Shah, “Effect of length on compressive strain softening of concrete,” Journal of Engineering Mechanics, vol. 123, no. 1, pp. 25–35, 1997. View at: Publisher Site  Google Scholar
 J.I. Sim, K.H. Yang, H.Y. Kim, and B.J. Choi, “Size and shape effects on compressive strength of lightweight concrete,” Construction and Building Materials, vol. 38, pp. 854–864, 2013. View at: Publisher Site  Google Scholar
 G. Markeset and A. Hillerborg, “Softening of concrete in compression—localization and size effects,” Cement and Concrete Research, vol. 25, no. 4, pp. 702–708, 1995. View at: Publisher Site  Google Scholar
 A. K. Samani and M. M. Attard, “A stress–strain model for uniaxial and confined concrete under compression,” Engineering Structures, vol. 41, pp. 335–349, 2012. View at: Publisher Site  Google Scholar
 Z. P. Bažant and J. Planas, Fracture and Size Effect in Concrete and Other Quasibrittle Materials, CRC Press, New York, NY, USA, 1998.
 Z. P. Bažant, “Size effect in blunt fracture: concrete, rock, metal,” Journal of Engineering Mechanics, vol. 110, no. 4, pp. 518–535, 1984. View at: Google Scholar
 J.K. Kim and S. H. Eo, “Size effect in concrete specimens with dissimilar initial cracks,” Magazine of Concrete Research, vol. 42, no. 153, pp. 233–238, 1990. View at: Publisher Site  Google Scholar
 M. Attard and S. Setunge, “Stressstrain relationship of confined and unconfined concrete,” ACI Materials Journal, vol. 93, no. 5, pp. 432–442, 1996. View at: Google Scholar
 R. Vonk, Softening of Concrete Loaded in Compression, Eindhoven University of Technology, Eindhoven, Netherlands, 1992, Ph.D. thesis.
 K. H. Lee, K. H. Yang, J. H. Mun, and S. J. Kwon, “Mechanical properties of concrete made from different expanded lightweight aggregates,” ACI Materials Journal, 2018, In press. View at: Publisher Site  Google Scholar
 J.S. Mun, J.H. Mun, K.H. Yang, and H. Lee, “Effect of substituting normalweight coarse aggregate on the workability and mechanical properties of heavyweight magnetite concrete,” Journal of the Korea Concrete Institute, vol. 25, no. 4, pp. 439–446, 2013, in Korean. View at: Publisher Site  Google Scholar
 G. Muciaccia, G. Rosati, and G. Di Luzio, “Compressive failure and size effect in plain concrete cylindrical specimens,” Construction and Building Materials, vol. 137, pp. 185–194, 2017. View at: Publisher Site  Google Scholar
 A. Neville, Properties of Concrete, Pearson Education Limited, Harlow, UK, 5th edition, 2011.
 M. A. Taylor and B. B. Broms, “Shear bond strength between coarse aggregate and cement paste or mortar,” ACI Proceedings, vol. 61, no. 8, pp. 939–958, 1964. View at: Google Scholar
 T. Noguchi, F. Tomosawa, K. M. Nemati, B. M. Chiaia, and A. P. Fantilli, “A practical equation for elastic modulus of concrete,” ACI Structural Journal, vol. 106, no. 5, pp. 690–696, 2009. View at: Google Scholar
 P. T. Wang, S. P. Shah, and A. E. Naaman, “Stressstrain curves of normal and lightweight concrete in compression,” ACI Journal Proceedings, vol. 75, no. 11, pp. 603–611, 1978. View at: Google Scholar
 K. H. Yang and J. I. Sim, “Modeling of the mechanical properties of structural lightweight concrete based on size effects,” Technical Report, Department of Plant∙Architectural Engineering, Kyonggi University, Suwon, Republic of Korea, 2011, in Korean. View at: Google Scholar
 H. C. Mertol, S. J. Kim, A. Mirmiran, S. Rizkalla, and P. Zia, “Behavior and design of HSC members subjected to axial compression and flexure,” in Proceedings of the 7th International Symposium on Utilization of HighStrength/High–Performance Concrete (SP228), H. G. Russell, Ed., pp. 395–420, American Concrete Institute, Washington, DC, USA, 2005. View at: Google Scholar
 M. H. Zhang and O. E. Gjørv, “Mechanical properties of high strength lightweight concrete,” ACI Materials Journal, vol. 88, no. 3, pp. 240–247, 1991. View at: Publisher Site  Google Scholar
 K. H. Yang, J. S. Mun, and H. Lee, “Workability and mechanical properties of heavyweight magnetite concrete,” ACI Materials Journal, vol. 111, no. 3, pp. 273–282, 2014. View at: Publisher Site  Google Scholar
 T. H. Wee, M. S. Chin, and M. A. Mansur, “Stressstrain relationship of highstrength concrete in compression,” Journal of Materials in Civil Engineering, vol. 8, no. 2, pp. 70–76, 1996. View at: Publisher Site  Google Scholar
 A. H. Nilson, “High–strength concrete: an overview of Cornell research,” in Proceedings of Symposium on Utilization of High–Strength Concrete, pp. 27–37, Stavanger, Norway, June 1987. View at: Google Scholar
 S.T. Yi, J.K. Kim, and T.K. Oh, “Effect of strength and age on the stressstrain curves of concrete specimens,” Cement and Concrete Research, vol. 33, no. 8, pp. 1235–1244, 2003. View at: Publisher Site  Google Scholar
 T. C. Liu, Stress–Strain Response and Fracture of Concrete in Biaxial Compression, Cornell University, Ithaca, NY, USA, 1971.
 A. Tomaszewicz, “Betongens arbeidsdiagram,” SINTEF, Trondheim, Norway, 1984, SINTEF Report No. STF 65A84065. View at: Google Scholar
 J. H. Mun, J. S. Mun, and K. H. Yang, “Stress–strain relationship of heavyweight concrete using magnetite aggregate,” Journal of Architectural Institute of Korea, vol. 29, no. 8, pp. 85–92, 2013, in Korean. View at: Google Scholar
 Comité EuroInternational du Beton (CEBFIP), Structural Concrete: Textbook on Behaviour, Design and Performance, International Federation for Structural Concrete (fib), Lausanne, Switzerland, 1999.
Copyright
Copyright © 2019 KeunHyeok Yang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.