Search

The paper studies systems of rigid bodies with randomly generated geometry interconnected by normal and tangential bonds. The stiffness of these bonds determines themacroscopic elastic modulus while the macroscopic Poisson's ratio of the system is determined solely by the normal/tangential stiffness ratio. Discrete models with no directional bias have the same probability of element orientation for any direction and therefore the same mechanical properties in a statistical sense at any point and direction. However, the layers of elements in the vicinity of the boundary exhibit biased orientation, preferring elements parallel with the boundary. As a consequence, when strain occurs in this direction, the boundary layer becomes stiffer than the interior for the normal/tangential stiffness ratio larger than one, and vice versa. Nonlinear constitutive laws are typically such that the straining of an element in shear results in higher strength and ductility than straining in tension. Since the boundary layer tends, due to the bias in the elemental orientation, to involve more tension than shear at the contacts, it also becomes weaker and less ductile. The paper documents these observations and compares them to the results of theoretical analysis. [ABSTRACT FROM AUTHOR]

Discrete models of solids have been motivated, in large part, by the discontinuous and heterogeneous nature of material structure and its breakdown under loading. The capabilities of discrete models have evolved over the past several decades, offering novel means for investigating material structure–property relationships. However, lack of understanding of both the utilities and disadvantages of discrete models limits their further development and applications. This paper reviews relevant features of discrete approaches applied to modeling the mechanical behavior of geomaterials, concrete materials in particular. The discrete models are classified according to their form and abilities to represent elastic and fracture behaviors in the presence of large-scale material heterogeneity. Discretization of the material domain plays a large role in this respect. Emphasis is placed on particle-based lattice models. The relative merits of various strategies for introducing reinforcing components, which are essential for many applications, are outlined. Recent advances are highlighted, including the use of discrete models for coupled, multi-field analysis. The merits of discrete approaches are summarized in the conclusions. • Discrete methods for simulating concrete fracture are critiqued. • Model capabilities depend on the choice of physical or non-physical forms of discretization. • The topics of stress oscillation, Poisson's ratio and volumetric effects are examined. • Dual-lattice models are suited to modeling hygro-mechanical couplings in cracked media. • Application examples demonstrate the utilities and expanding use of discrete methods. [ABSTRACT FROM AUTHOR]

• Mesoscale discrete model of concrete enhanced by spatially random material parameters. • Notched and unnotched three point bending and uniaxial tension are analyzed. • Fracture process zone statistically evaluated for various random field parameters. • Nontrivial dependence of peak loads on correlation length and variance is reported. • Local averaging and weakest-link are identified as the underlying mechanisms. Numerical simulations of concrete fracture performed with a probabilistic mesoscale discrete model are presented. The model represents a substantial part of material randomness by assigning random locations to the largest aggregates. The remaining part of randomness is introduced by causing material parameters to fluctuate randomly via a homogeneous random field. An extensive numerical study performed with the model considers prisms loaded in uniaxial tension with both fixed and rotating platens, and also beams with and without a notch loaded in three point bending. The results show the nontrivial effect of (i) autocorrelation length and (ii) variance of the random field on the fracture behavior of the model. Statistics of the peak load are presented as well as the size and shape of the fracture process zone at the moment when the maximum load is attained. Local averaging within the fracture process zone and weakest-link are identified as underlying mechanisms explaining the reported results. The companion paper, Part II (Vořechovský and Eliáš, 2020), introduces an analytical model capable of predicting the distribution of the peak load obtained with the probabilistic discrete model via the simple estimation of extremes of a random field obtained as moving average of local strength. [ABSTRACT FROM AUTHOR]