A hyperelastic constitutive model is developed for textile composite reinforcement at large strain. A potential is proposed, which is the addition of two tension and one shear energies. The proposed potential is a function of the right Cauchy Green and structural tensor invariants whose choice corresponds to textile composite reinforcement mechanical behavior which exhibits weak elongations in the fiber directions and large angular variations in the fabric plane. The model is implemented in a Vumat user routine of ABAQUS/Explicit. Some elementary tests are performed in order to identify the model and verify its validity. It is then used to simulate the hemispherical punch forming of balanced and unbalanced fabrics. A correct agreement is obtained with experimental forming processes.

International audienc

This paper describes implementation of anisotropic damage mechanics in the material point method. The approach was based on previously proposed, fourth-rank anisotropic damage tenors. For implementation, it was convenient to recast the stress update using a new damage strain partitioning tensor. This new tensor simplifies numerical implementation (a detailed algorithm is provided) and clarifies the connection between cracking strain and an implied physical crack with crack opening displacements. By using 2 softening laws and 3 damage parameters corresponding to 1 normal and 2 shear cracking strains, damage evolution can be directly connected to mixed tensile and shear fracture mechanics. Several examples illustrate interesting properties of robust anisotropic damage mechanics such as modeling of necking, multiple cracking in coatings, and compression failure. Direct comparisons between explicit crack modeling and damage mechanics in the same material point method code show that damage mechanics can quantitatively reproduce many features of explicit crack modeling. A caveat is that strengths and energies assigned to damage mechanics materials must be changed from measured material properties to apparent properties before damage mechanics can agree with fracture mechanics.

We have developed a new geomechanical workflow to study the mechanics of hydraulic fracturing in naturally fractured unconventional reservoirs. This workflow used the material point method (MPM) for computational mechanics and an equivalent fracture model derived from continuous fracture modeling to represent natural fractures (NFs). We first used the workflow to test the effect of different stress anisotropies on the propagation path of a single NF intersected by a hydraulic fracture. In these elementary studies, increasing the stress anisotropy was found to decrease the curving of a propagating NF, and this could be used to explain the observed trends in the microseismic data. The workflow was applied to Marcellus and Eagle Ford wells, where multiple geomechanical results were validated with microseismic data and tracer tests. Application of the workflow to a Marcellus well provides a strain field that correlates well with microseismicity, and a maximum energy release rate, or [Formula: see text] integral at each completion stage, which appeared to correlate to the production log and could be used to quantify the impact of skipping the completion stages. On the first of two Eagle Ford wells considered, the MPM workflow provided a horizontal differential stress map that showed significant variability imparted by NFs perturbing the regional stress field. Additionally, a map of the strain distribution after stimulating the well showed the same features as the interpreted microseismic data: three distinct regions of microseismic character, supported by tracer tests and explained by the MPM differential stress map. Finally, the workflow was able to estimate, in the second well with no microseismic data, its main performance characteristics as validated by tracer tests. The field-validated MPM geomechanical workflow is a powerful tool for completion optimization in the presence of NFs, which affect in multiple ways the final outcome of hydraulic fracturing.

Transverse compression of wood is a process that induces large deformations. The process is dominated by elastic and plastic cell wall buckling. This work reports a numerical study of the transverse compression and densification of wood using a large-deformation, elastic-plastic constitutive law. The model is isotropic, formulated within the framework of hyperelasticity, and implemented in explicit material point method (MPM) software. The model was first validated for modeling of cellular materials by compression of an isotropic cellular model specimen. Next, it was used to model compression of wood by first validating use of isotropic, transverse plane properties for tangential compression of hardwood, and then by investigating both tangential and radial compression of softwood. Importantly, the discretization of wood specimens used MPM methods to reproduce accurately the complex morphology of wood anatomy for different species. The simulations have reproduced observations of stress-strain response during wood compression including details of inhomogeneous deformation caused by variations in wood anatomy.