2.1. Parameter Selection
Concrete is an artificial mixed material with cement as the main cementing material. Sand, stone, and water are mixed in certain ratios, and the mixture gradually solidifies and hardens after molding injection, vibration, curing, and other procedures. It is a kind of non-uniform quasi-brittle material with complex mechanical properties. In this paper, we set the parameters of the model according to the difference in bond strength between aggregate and mortar. The particle contact model (
Figure 1a) and particle bond model (
Figure 1b) form the foundation of the PFC2D model. The stick model can only relay contact force, and the parallel bond model in 2D/3D represents the line/surface bonding rather than that at a point. This can better reflect the capability of concrete mortar to resist torque as well as the interfacial transition zone between particles. Therefore, in the concrete model in this paper, different contact models are set between aggregate and aggregate, and between mortar and mortar, so as to more truly reflect the actual performance of concrete. The parallel bonding model was used to simulate the mechanical behaviors of mortar and the aggregate, and the linear contact bonding model was used to simulate the contact-related behavior of coarse aggregate.
The strength and deformation characteristics of concrete were preliminarily characterized by using the three macroscopic parameters of peak strength, elastic modulus, and Poisson’s ratio, and were measured using uniaxial compression tests. The peak strength σu, elastic modulus E and Poisson’s ratio ν were then used as test indices (objects), and the influence of the test factors (mesoscopic parameters of concrete) on them was examined. In the process of uniaxial compression in the particle discrete element model, σu is taken as the stress corresponding to the apex of the stress-strain curve, E is the tangent slope of the curve when the stress reaches half the peak stress for the first time, and ν is the absolute value of the ratio of lateral strain to axial strain at this time.
In simulations of the uniaxial compression test, the geometric parameters (L,W,n, Rmax/Rmin), load parameters (, , , ), parameters of the particles (,
Ec, kn/ks, μ) and key binding parameters (,
, ) affect the results. Because the parameters of the particles (including particles of cement mortar (ball) and aggregate particles (clump)) and bonds are the focus of this study, the geometric and load parameters were controlled in advance. The compression of members in engineering is mostly prismatic, because of which prismatic specimens can better reflect the compressive capacity of concrete than cubic specimens. A model size of 200 × 100 mm2 was thus set. Porosity (n) was set to 0.12, the radius ratio of the particles of cement mortar (Rmax/Rmin) was 2.5, particle density (ρball) was 2000 kg/m3, and aggregate density was 2650 kg/m3. We preset the average radius of particles of cement mortar (
) to 0.5 mm. The bonding radius factor (
) is normally set to one, and remains constant during the calculation; thus, it was not be considered. To reduce the number of unknown parameters, the stiffness parameters of the particles of cement mortar, their bonds, and the upper and lower loading plates were set to be consistent, that is, kn/ks =
=
, Ec =
=
. Parameters of the particles of aggregate were set to be consistent with those of contact stiffness, i.e., (kn/ks)g = ( and (EC)g = ()g given a wall loading rate ( ) of 0.02 m/s, the friction coefficient of the wall (μw) was zero.
Finally, the three macroscopic parameters σu, E and ν were determined as test indicators. Parameters of the particles of cement mortar and their bonds (Ec, kn/ks, , and μ), the aggregate ((Ec)g, (kn/ks)g, μg), and the interfacial transition zone and its bonds ((EC)j, (kn/ks)j, ()j, ()j, μj), a total of 13 mesoscopic parameters, were unknown. They were chosen as the experimental factors. y = Ec, kn/ks, , μ, (Ec)g, (kn/ks)g, μg, (Ec)j, (kn/ks)j, ()j, ()j, μj.
The relationship is as follows:
As concrete is a complex three-phase composite material with obvious anisotropy, the microscopic parameters of the aggregate, cement mortar, and the weak interfacial transition zone all need to be considered. The interfacial transition zone features both the aggregate particles and the cement paste, where the distribution of the cement particles is affected by the aggregate surface at the parts in contact. Due to the presence of many original cracks in the interior, most of the damage under external load starts from the interfacial transition zone, because of which it is the weakest part of concrete. The strength of the zone affects the overall mechanical properties and failure characteristics of concrete [
9] as well as the fluctuation in its stress–strain curve. Therefore, it must be considered to enable the model to represent empirical scenarios involving concrete. The aggregate is the “skeleton” of concrete that plays an important load-bearing role. In the PFC2D software, the thickness of the interfacial transition zone cannot be expressed because the bond between particles is set through contact. In this paper, the influence of the thickness of the interfacial transition zone on the mechanical properties of concrete was ignored.
The parameters of the model are shown in
Table 1.
2.3. Aggregate Generation
Past work has used the CT scanning technique to simulate the aggregate, and the scanned images are subjected to threshold processing using third-party software to obtain the boundary information to establish the model. However, this often reflects only one or more conditions while ignoring the random aggregate model. This paper uses a random algorithm for generating a convex polygonal aggregate.
In this algorithm, the circle is used as the base to generate the aggregate to ensure that the model is convex polygonal. The number of points on the circle determines the number of sides of a convex polygon. To ensure the randomness of the number of edges of the generated aggregate, a minimum value and a maximum value were set for the number of points on the circle. The range of the number of edges of the convex polygon was determined by calculating the difference between them, which was then multiplied by a random number. To round up the results, the resulting integer was set as the number of points on the circle, which was also the number of sides of the convex polygon:
where
N is the number of points on the boundary of the circle,
Nmin is the minimum number of points generated.
Nmax is the maximum number of points generated, range is a pseudo-random number uniformly distributed between zero and one, and
round is a rounding of the entire function.
Assuming that one of the circles has radius
R and coordinates (
X,
Y), the point coordinates (
X1,
Y1) on each circle are calculated by Equation (8):
where
θ is the anticlockwise angle between the line of each point and the center of the circle and has a positive value along the
X-axis.
Given the randomness to ensure the positional angle, the normal vector of the convex polygon contains information on internal damage that is used to determine the scope of the convex polygons generated after the outline. The aim is to build a set of algorithms to calculate the convex polygons of each edge point with respect to its internal normal vector.
Based on the proposed algorithm, multiple aggregate models were built. The fewer edges the aggregate model had, the sharper was its shape. The model was used to simulate an elongated aggregate. The more edges the aggregate model had, the smoother the shape of the aggregate was, and the more round it tended to be. This model was used to simulate the circular aggregate.
To simulate aggregate of any shape for the numerical analysis of concrete, multiple particles can be combined in the particle flow program to form a “block.” The block formed by a clump participates in the cyclic calculation as a rigid body. During the calculation process, the distance and contact force between the internal particles do not change with the steps of calculation, and the deformation of the block occurs only along the boundary. Previous studies have shown that the compressive strength of the aggregate is twice as high as that of cement mortar, and it is not damaged in the process of concrete compression. Therefore, the unbreakable unit “clump’’ is used as the object of generation of aggregate in the simulation. In the calculation cycle, contact between the particles inside the “block’’ is ignored, which reduces calculation time.
The range of sizes of the coarse aggregate used was 5~15 mm, and the Walraven Equation [
11] (Equation (9)) was used to calculate its gradation. The percentage (
Pk) of the volume of coarse aggregate in the concrete specimen was 75%. The average gradation curve obtained by using a slice of the non-circular aggregate concrete specimen and subjecting it to image processing using the “equivalent particle size” (
Dequivalent particle size = 2(
Sarregate area/
π)
0.5) conformed to a 3D gradation curve and a 2D conversion curve. Therefore, the circular area of the aggregate was used as the equivalent area of the convex polygonal aggregate. Aggregate sizes of 5~7 mm, 7~9 mm, 9~11 mm, 11~13 mm, and 13~15 mm was used, and the results are shown in
Table 2:
where
Pk is the percentage by volume of coarse aggregate in the concrete specimen,
PC (
D <
D0) is the probability that the particle size
D in the section is less than the size of the sieve
D0, and
Dmax is the maximum aggregate particle size.
The number of particles within the range of each particle size is calculated according to the following equation:
where
ni is the number of particles in a certain range,
pi is the probability that the particle size in the section is smaller than the sieve in the Walraven Equation, and
A and
Ai are the sectional area and aggregate area, respectively.
We set-up four to six polygons sides in MATLAB to generate a random aggregate model and imported it into the PFC2D model for use as a model of the concrete aggregate.
The generation model is shown in
Figure 2.