Evolution of bond fractures in a randomly distributed fiber network

Fracture in a planar randomly ordered fiber network subjected to approximately homogenous macroscopic stress and strain field is considered. A theory describing material degradation on a macroscopic scale is derived via Griffith’s energy balance for an internal fractured area in the network assuming the active fracture process on the microscopic level is fiber–fiber bond breakage. Attention is confined to a purely mechanical theory assuming isothermal processes and the theory relies on equations commonly used in theories of statistical physics. In the theory, a bond breaking driving force is stated to be equal to the elastic strain energy density of a non-fractured network. A debond fraction can be coupled to a linearly decrease of the network’s macroscopic stiffness. The rate of the fracture processes is determined by the network’s inherent properties (bond and fiber density, bond strength, etc.). During the loading process, until onset of localization, the bond breaks occur at randomly distributed locations spread over the fiber network and the theory estimate material degradation on a macroscopic level. When localization takes place, the fracture process changes from a two-dimensional randomly distributed process to a one-dimensional process and other theories have to be included to describe post-localization behavior. An approximately in-plane isotropic low-density paper is used in tensile experiments while monitoring acoustic emission activity to evaluate the theory. The experimentally obtained results support the theory surprisingly well.     

The influence of particle fluctuations on the average rotation in an idealized granular material

Granular materials are a simple example of a Cosserat continuum in that the average particle rotations may differ from the rotation of the average deformation. In the absence of couple stress, this difference insures that the stress is symmetric. This has been shown in theories that assume that the displacement at particle contacts is given by the average deformation and spin. Here, we indicate how the difference between the average rotation of the particles and the average rotation of the deformation can be determined when fluctuations in particle displacements and rotations satisfy local force and moment equilibria in a random aggregate of identical spheres. The predictions based on this model are in better agreement with numerical simulation than that given by the simple average strain theory.     

GPU-based simulations of fracture in idealized brick and mortar composites

Stiff ceramic platelets (or bricks) that are aligned and bonded to a second ductile phase with low volume fraction (mortar) are a promising pathway to produce stiff, high-toughness composites. For certain ranges of constituent properties, including those of some synthetic analogs to nacre, one can demonstrate that the deformation is dominated by relative brick motions. This paper describes simulations of fracture that explicitly track the motions of individual rigid bricks in an idealized microstructure; cohesive tractions acting between the bricks introduce elastic, plastic and rupture behaviors. Results are presented for the stresses and damage near macroscopic cracks with different brick orientations relative to the loading orientation. The anisotropic macroscopic initiation toughness is computed for small-scale yielding conditions and is shown to be independent of specimen geometry and loading configuration. The results are shown to be in agreement with previously published experiments on synthetic nacre.     

Adaptive coupling between damage mechanics and peridynamics: A route for objective simulation of material degradation up to complete failure

The objective (mesh-independent) simulation of evolving discontinuities, such as cracks, remains a challenge. Current techniques are highly complex or involve intractable computational costs, making simulations up to complete failure difficult. We propose a framework as a new route toward solving this problem that adaptively couples local-continuum damage mechanics with peridynamics to objectively simulate all the steps that lead to material failure: damage nucleation, crack formation and propagation. Local-continuum damage mechanics successfully describes the degradation related to dispersed microdefects before the formation of a macrocrack. However, when damage localizes, it suffers spurious mesh dependency, making the simulation of macrocracks challenging. On the other hand, the peridynamic theory is promising for the simulation of fractures, as it naturally allows discontinuities in the displacement field. Here, we present a hybrid local-continuum damage/peridynamic model. Local-continuum damage mechanics is used to describe “volume” damage before localization. Once localization is detected at a point, the remaining part of the energy is dissipated through an adaptive peridynamic model capable of the transition to a “surface” degradation, typically a crack. We believe that this framework, which actually mimics the real physical process of crack formation, is the first bridge between continuum damage theories and peridynamics. Two-dimensional numerical examples are used to illustrate that an objective simulation of material failure can be achieved by this method.     

The transverse permeability of disordered fiber arrays: a statistical correlation in terms of the mean nearest interfiber spacing

In the porous media literature, unidirectional fibrous systems are broadly categorized as ordered or disordered. The former class, easily tractable for analysis purposes but limited in its relation to reality, involves square, hexagonal and various staggered arrays. The latter class involves everything else. While the dimensionless hydraulic permeability of ordered fibrous media is known to be a deterministic function of their porosity ϕ, the parameters affecting the permeability of disordered fiber arrays are not very well understood. The objective of this study is to computationally investigate flow across many unidirectional arrays of randomly placed fibers and derive a correlation between K and some measure of their microstructure. In the process, we explain the wide scatter in permeability values observed computationally as well as experimentally. This task is achieved using a parallel implementation of the Boundary Element Method (BEM). Over 600 simulations are carried out in two-dimensional geometries consisting of 576 fiber cross-sections placed within a square unit cell by a Monte Carlo procedure. The porosity varies from 0.45 to 0.90. The computed permeabilities are compared with earlier theoretical results and experimental data. Analysis of the computational results reveals that the permeability of disordered arrays with ϕ < 0.7 is reduced as the non-uniformity of the fiber distribution increases. This reduction can be substantial at low porosities. The key finding of this study is a direct correlation between K and the mean nearest inter-fiber spacing , the latter depending on the microstructure of the fibrous medium.     

On some connections between equivalent single material and mixture theory models for fiber reinforced hyperelastic materials

There are two approaches that can be used to model the large strain mechanical response of material systems in which elastic fibers are embedded in an elastic matrix. In the first approach, a fiber reinforced material undergoing large deformation is homogenized in the sense that it is assumed to act as an equivalent single material that is transversely isotropic and hyperelastic. Both constituents then share a common reference configuration, which is typically assumed to be a natural or stress-free configuration for the equivalent single material. The stress depends on a single deformation gradient defined with respect to the natural configuration.In the second approach, the fiber/matrix system is treated as a mixture, with the matrix and the fibrous constituents having their own reference configurations and material symmetries. The total stress depends on the deformation gradients and material symmetries for both constituents, defined with respect to their reference configurations.Under appropriate assumptions, the constitutive theory developed using mixture theory can coincide with the constitutive theory assuming an equivalent single material that is transversely isotropic and hyperelastic. This paper explores the connection between the two approaches by considering the various reference configurations and material symmetries.     

On the interplay between material flaws and dynamic necking

In this paper we investigate the interplay between material defects and flow localization in elastoplastic bars subjected to dynamic tension. For that task, we have developed a 1D finite difference scheme within a large deformation framework in which the material is modelled using rate-dependent J₂ plasticity. A perturbation of the initial yield stress is introduced in each node of the finite difference mesh to model localized material flaws. Numerical computations are carried out within a wide spectrum of strain rates ranging from 500 s⁻¹ to 2500 s⁻¹. On the one hand, our calculations reveal the effect of the material defects in the necking process. On the other hand, our results show that the necking inception, instead of being a random type process, is the deterministic result of the interplay between the mechanical behaviour of the material and the boundary conditions. This conclusion agrees with the experimental evidence reported by Rittel et al. [1] and Rotbaum et al. [2].     

On the coupling mechanism between nonlinear solitary waves and slender beams

This paper describes the coupling mechanism between highly nonlinear solitary waves propagating along a granular system and slender beams in contact with the granular medium. Nonlinear solitary waves are compact non-dispersive waves that can form and travel in nonlinear systems such as one-dimensional chains of particles, where they are conventionally generated by the mechanical impact of a striker. These waves have a constant spatial wavelength and their speed, amplitude, and duration can be tuned by modifying the particles’ material or size, or the velocity of the striker. In the study presented in this article we investigated numerically the interaction between solitary waves propagating along a chain of granular particles and a slender beam. Some of the numerical findings were validated experimentally. We found that the geometric and mechanical properties of the beam or thermal stress applied to the beam alter certain features of the solitary waves. In the future, these findings may be used to develop a novel sensing system for the nondestructive assessment of beams.     

Crack shielding and material deterioration in damaged materials: an antiplane shear fracture problem

Summary A simple damage evolution model is proposed for a quasibrittle material in the framework of continuum damage mechanics. The model is used to obtain a closed form solution for a mode-III stationary crack under small scale damage conditions. It is found that the crack tip stress intensity factor is reduced, i.e., the crack is shielded by the damage. However, this shielding effect is completely offset by the material deterioration caused by the damage. It also holds for steady state crack growth. When the most effective shielding is reached for the stationary crack, the zone dominated by the stress intensity factor shrinks to the crack tip. The results for the antiplane shear problem should shed some light on the in- plane fracture problem. Received 4 August 1997; accepted for publication 7 October 1997     

A network theory for the thermal conductivity of an amorphous polymeric material

A model to relate the thermal conductivity tensor to the deformation of an amorphous polymeric material above the glass transition temperature is presented. The basis of the model is formed by the transient network theory for polymer melts. With this theory it is possible to calculate the average orientation of the macromolecular segments as a function of the history of the deformation. Combined with an expression which relates the thermal conductivity to the orientation of the molecules, this provides us with the information needed to calculate the heat conduction tensor. Despite the fact that the simplest possible network model is chosen, there is good agreement with the sparse, experimental results.     

Experimental measurements and computational modeling of the flow field in an idealized human oropharynx

The velocity field in the central sagittal plane of an idealized representation of the human oropharynx (HOP) during steady inspiration, simulating oral inhalation through an inhaler mouthpiece, was measured experimentally using endoscopic particle image velocimetry (PIV). Measurements were made at three flow rates: 15, 30, and 90 L/min, which correspond to a wide range of physiological conditions. Extensive tests were performed to verify the veracity of the PIV data. The flow was also modeled computationally using Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) methods. The PIV data clearly indicate the complex nature of HOP flow, with three-dimensionality and several regions of separation and recirculation evident. Comparison of the experimental and computational results shows that, although the RANS CFD reproduces the basic features of the flow, it does not adequately capture the increased viscous effects at lower Reynolds numbers. The results demonstrate the need for more development and validation of CFD modeling, in particular RANS methods, in these flows.     

On the relationship between critical tensile stress and fracture toughness in mild steel

An analysis is presented which relates the critical value of tensile stress (σ_f) for unstable cleavage fracture to the fracture toughness (K_Ic) for a high-nitrogen mild steel under plane strain conditions. The correlation is based on (i) the model for cleavage cracking developed by E. Smith and (ii) accurate plastic-elastic solutions for the stress distributions ahead of a sharp crack derived by J. R. Rice and co-workers. Unstable fracture is found to be consistent with the attainment of a stress intensification close to the tip such that the maximum principal stress σ_yy exceeds σ_f over a characteristic distance, determined as twice the grain size. The model is seen to predict the experimentally determined variation of K_Ic with temperature over the range -150 to -75°C from a knowledge of the yield stress and hardening properties. It is further shown that the onset of fibrous fracture ahead of the tip can be deduced from the position of the maximum achievable stress intensification. The relationship between the model for fracture ahead of a sharp crack, and that ahead of a rounded notch, is discussed in detail.     

The dynamic growth of a void in a plastic material and an application to fracture

The problem of the dynamic growth of a spherical void in applied stress and strain-rate fields is solved approximately using a variational method. The solution is applied to a model of fracture by void growth from brittle or incoherent inclusions. It is shown that the speed of fracture is limited by inertia near the crack tip. Some numerical results are obtained from data for high strength steels.     

On the influence of material non-linearities in geometric modeling of kink band instabilities in unidirectional fiber composites

The axial compressive failure of aligned fiber composites triggered by kink band instabilities is the topic of investigation herein. Particular emphasis is put on the accurate prediction of the post-failure regime, where fiber composites are known to exhibit substantial post-failure strength. In this regard, a previous analytical model, based on geometric relationships and energy principles, is enhanced by consistently taking into account material non-linearities. Therefore, a non-linear constitutive law is introduced herein based on a newly developed exponential formulation. This non-linear constitutive law is subsequently implemented into the stress–strain response in interlaminar shearing as well as the compression response. The model enhancements are validated against published experimental data yielding excellent comparisons. Furthermore, the relevance of modeling non-linear material behavior in interlaminar dilation and bending is assessed using a bilinear constitutive law. However, implementing non-linear material behavior does not yield any improvements and can therefore be neglected.     

Mode III fracture of an arbitrary oriented crack in two dimensional functionally graded material

In this paper, a two dimensional functionally graded material (2D-FGM) under an anti-plane load with an internal crack is considered. The crack is oriented in an arbitrary direction. The material properties are assumed to vary exponentially in two planar directions. The problem is analyzed and solved by two different methods namely Fourier integral transforms with singular integral equation technique, and then by the finite element method. The effects of crack orientation, material non-homogeneity, and other parameters on the value of stress intensity factor (SIF) are studied. Finally, the obtained results for Mode III stress intensity factor of different methods are compared.     

A mesomechanical fracture model for an orthotropic material with different tensile and compressive strengths

The paper addresses a fracture problem for an orthotropic cracked plate made of a material with different tensile and compressive strengths and subjected to biaxial loading. The problem is solved using a micromechanical fracture model proposed earlier by the authors. It is assumed that the fracture of the material in the fracture process zones at the crack front is described by the Gol’denblat–Kopnov failure criterion. Strength curves for an orthotropic cracked plate with different strength and fracture-toughness parameters are plotted     

On kinematic,thermodynamic, and kinetic coupling of a damage theory for polycrystalline material

The paper proposes a new consistent formulation of polycrystalline finite-strain elasto-plasticity coupling kinematics and thermodynamics with damage using an extended multiplicative decomposition of the deformation gradient that accounts for temperature effects. The macroscopic deformation gradient comprises four terms: thermal deformation associated with the thermal expansion, the deviatoric plastic deformation attributed to the history of dislocation glide/movement, the volumetric deformation gradient associated with dissipative volume change of the material, and the elastic or recoverable deformation associated with the lattice rotation/stretch. Such a macroscopic decomposition of the deformation gradient is physically motivated by the mechanisms underlying lattice deformation, plastic flow, and evolution of damage in polycrystalline materials. It is shown that prescribing plasticity and damage evolution equations in their physical intermediate configurations leads to physically justified evolution equations in the current configuration. In the past, these equations have been modified in order to represent experimentally observed behavior with regard to damage evolution, whereas in this paper, these modifications appear naturally through mappings by the multiplicative decomposition of the deformation gradient. The prescribed kinematics captures precisely the damage deformation (of any rank) and does not require introducing a fictitious undamaged configuration or mechanically equivalent of the real damaged configuration as used in the past.