Elastic Inclusion Research Articles

The transient dynamic wave process in unidirectional composites

The propagation of nonstationary waves in unidirectional composite material is simulated by finite-element analysis. The composite material is explicitly formulated, i.e., discrete, stiff, elastic inclusions uniformly distributed in a viscoelastic matrix are modeled. An original procedure rooted in the linear Boltzmann–Volterra equations accounts for the viscoelastic effects in the matrix behavior. The chosen approach involves finite-element simulations of the transient wave process in two indicative configurations to study the effects attributed to wave scattering by the inclusions. The transient wave propagation in the unidirectional composite material and a viscoelastic specimen without inclusions are compared. Peculiarities in the numerical results due to the finite-element algorithm employed are decoupled from the output that reveals purely mechanical phenomena. For this purpose, the transient wave propagation in a one-dimensional elastic continuum, for which the exact analytical solution is well known, is effectively obtained by finite-element analysis. Furthermore, based on the unidirectional composite a homogenized model is defined by applying the mixture rule. An additional comparison between the composite containing discrete inclusions and the analogous homogenized material provides a basis to assess the output of models obtained by implementing a self-consistent scheme in a broader context.

Multi-functional periodically heterogeneous structures for energy harvesting and vibration attenuation-effects of piezoelectricity and shunting circuits

Abstract We explore a kind of metamaterial plate structures intended for simultaneous energy harvesting and vibration control. These structures are designed using a periodically perforated piezoelectric plate (the matrix) with elastic inclusions situated in the holes and serving for the resonators. The design options comprise two- and three-phase configurations related to the mechanical connection between the matrix and inclusions. By introducing a singularity—the focal spot created as a defect in the perfectly periodic structure and using the theory of super-cell, an enhanced piezoelectric energy harvester is obtained. It is observed that such a meta-structure serves as a dual-purpose system: efficiently capturing vibrational energy at a focal spot while maintaining the overall vibration attenuation throughout the structure. The band gap analysis based on the Bloch’s wave decomposition theory shows that by concentrating energy and halting vibration propagation, approximately 10 times energy harvesting enhancement and a remarkable 100 dB reduction in vibrations are achieved simultaneously. Besides the passive response of these meta-structures, we consider its extension by an external electric circuit (EC). Such modified configurations enable to exploit ‘actively’ the piezoelectric plate property to transmit the mechanical response between two, or more distant locations. Due to nonlocal interactions introduced by means the controllable EC, we consider optimization of the EC impedance to reduce the vibrations at a selected location of the whole structure without any external energy supply. The computational study discovers perspectives and benefits of designing such active self-powered meta-structures.

Elastic coupled phase theory based on the Cosserat equations: Propagation of coherent waves

We develop a new coupled phase theory (CPT) in order to model the propagation of elastic coherent waves in homogeneous solid media containing randomly distributed spherical elastic inclusions. To this end, the Buyevich’s theory previously developed for fluid media (Y. Buyevich et I. Shchelchkova, Progress in Aerospace Sciences 18, 1978.) has been adapted and applied to the Cosserat equations. A key point is to introduce the Independent Scattering Approximation (ISA). We show that the coherent wavenumbers depend explicitly on the external force applied to spheres, the torque characterizing their rotation and the stresslet which is the result of the resistance of the particles to the straining motion. Numerical results are in good agreement with experimental data.

Scattering of elastic waves by a cluster of cylindrical inclusions: multiple scattering analysis and equivalent scatterer modeling

Abstract The two-dimensional scattering of elastic waves by a cluster of aligned circular cylindrical elastic inclusions embedded in a circular cylindrical domain of an infinitely extended elastic medium is analysed numerically. The analysis is based on a semi-analytical method where the scattered waves are expressed as the eigenfunction expansions of the time-harmonic displacement potentials and the expansion coefficients are determined by a numerical collocation technique. The wave fields within and around the cluster are demonstrated for different spatial arrangements of inclusions in the cluster. It is also examined if the scattering characteristics of the cluster can be simulated by a single equivalent scatterer which has the effective homogenized properties of the matrix/inclusion composite. As a result, it is shown that for sufficiently small inclusion radius-wavelength ratios, the scattering characteristics of the cluster are relatively insensitive to the spatial inclusion arrangements therein, and can be reasonably simulated by those of the equivalent scatterer.

Elastic inclusion effects on deformation behavior of indented elastic–plastic solids

Elastic inclusion effects on deformation behavior of indented elastic–plastic solids

Mechanical response and failure mechanism of circular inclusion embedded in brittle materials under dynamic impact

Inclusions are prevalent in both natural and synthetic materials. Gaining a comprehensive understanding of their dynamic response and interaction with the surrounding matrix is essential in fundamental mechanics and engineering applications. This study aims to achieve such understanding through theoretical analysis, experimental investigations, and numerical simulations, focusing on the dynamic destruction of inclusions and the underlying mechanisms. Firstly, the circumferential, radial, and shear stresses around elastic heterogeneous inclusions were derived by the wave function expansion and Duhamel integration methods. Subsequently, a set of experiments on sandstone specimens containing three types of inclusions (plaster, epoxy resin, cement) were conducted utilizing the Split Hopkinson Pressure Bar (SHPB) system in conjunction with high-speed photography and Digital Image Correlation (DIC) techniques to obtain the surface deformation of inclusion specimens. Eventually, a series of Finite Element Method (FEM) numerical simulations adopted suitable materials were also carried out to investigate the fracture process of inclusion and matrix under dynamic impact. The results demonstrate that the stress distributions and fracture mode of matrix and inclusion are highly dependent on the physical and mechanical properties of the inclusions and the surrounding matrix, specifically, density ρ, elastic modulus E, Poisson's ratio v, and strength. With an increase in the E and v, there is a reduction in the concentration of circumferential stress, while the radial and shear stresses experience an increase. The experimental and numerical results corroborated the theoretical findings and indicated that the localized dynamic stress concentrations induced by wave scattering around the inclusions directly dominated both local and overall specimen failure. Due to the dissimilarities in elastic parameters and strength between the inclusion and the matrix, the stress state of the inclusion and the interface is not homogeneous. Under dynamic loading, the weaker inclusion experiences tensile cracking at both ends of the loading, and as the mechanical properties (ρ, E, v, and strength) of the inclusion rise, a transition from tensile- hybrid-shear failure occurs within the inclusion. The findings of this study help us understand the dynamic failure mechanisms of dissimilar inclusions in brittle material.

Thin inclusion at the junction of two elastic bodies: non-coercive case.

This article addresses an analysis of the non-coercive boundary value problem describing an equilibrium state of two contacting elastic bodies connected by a thin elastic inclusion. Nonlinear conditions of inequality type are imposed at the joint boundary of the bodies providing a mutual non-penetration. As for conditions at the external boundary, they are Neumann type and imply the non-coercivity of the problem. Assuming that external forces satisfy suitable conditions, a solution existence of the problem analysed is proved. Passages to limits are justified as the rigidity parameters of the inclusion and the elastic body tend to infinity.This article is part of the theme issue 'Non-smooth variational problems with applications in mechanics'.

Reconstruction of elastic inclusions in layered medium

We consider the recovery of small inclusions in a two-layered and an arbitrary multi-layered medium for the elastic equation, respectively. We use layer potential techniques and asymptotic analysis to obtain asymptotic expansions of the perturbed elastic field in a two-layered and an arbitrary multi-layered medium, respectively. Furthermore, we show the uniqueness of the recovery of the locations and Lamé constants of small inclusions through a single measurement.

Ultrasound Shear Elastography With Expanded Bandwidth (USEWEB): A Novel Method for 2D Shear Phase Velocity Imaging of Soft Tissues.

Ultrasound shear wave elastography (SWE) is a noninvasive approach for evaluating mechanical properties of soft tissues. In SWE either group velocity measured in the time-domain or phase velocity measured in the frequency-domain can be reported. Frequency-domain methods have the advantage over time-domain methods in providing a response for a specific frequency, while time-domain methods average the wave velocity over the entire frequency band. Current frequency-domain approaches struggle to reconstruct SWE images over full frequency bandwidth. This is especially important in the case of viscoelastic tissues, where tissue viscoelasticity is often studied by analyzing the shear wave phase velocity dispersion. For characterizing cancerous lesions, it has been shown that considerable biases can occur with group velocity-based measurements. However, using phase velocities at higher frequencies can provide more accurate evaluations. In this paper, we propose a new method called Ultrasound Shear Elastography with Expanded Bandwidth (USEWEB) used for two-dimensional (2D) shear wave phase velocity imaging. We tested the USEWEB method on data from homogeneous tissue-mimicking liver fibrosis phantoms, custom-made viscoelastic phantom measurements, phantoms with cylindrical inclusions experiments, and in vivo renal transplants scanned with a clinical scanner. We compared results from the USEWEB method with a Local Phase Velocity Imaging (LPVI) approach over a wide frequency range, i.e., up to 200-2000 Hz. Tests carried out revealed that the USEWEB approach provides 2D phase velocity images with a coefficient of variation below 5% over a wider frequency band for smaller processing window size in comparison to LPVI, especially in viscoelastic materials. In addition, USEWEB can produce correct phase velocity images for much higher frequencies, up to 1800 Hz, compared to LPVI, which can be used to characterize viscoelastic materials and elastic inclusions.

Elastic prolate spheroidal inclusion in an infinite elastic solid—an exact analysis of the inclusion stress by an engineering treatment of the problem based on the corresponding cavity solutions

The paper demonstrates the analysis of the stress and strain states of an inhomogeneous structure with an axial symmetrical spheroidal inclusion in an infinite solid under a remote uniform tension load derived from the conditions of the theory of elasticity. With respect to the inhomogeneous structure, the deformations produce constraints that require a complete 3-dimensional analysis in the z, r-coordinate system. The solution generates the stress state of the inclusion and at the interface of the matrix. Spheroidal inclusions in an uniform outer tension stress field deform self-similarly to a rather elongated spheroid. With respect to the compatibility condition and depending on the different elastic moduli and Poisson’s ratios in the inclusion and matrix, the solution procedure leads to a system of two linear equations with the magnitudes of the two tractions as unknowns. In this way, the solution is shortened to a twofold statically indeterminate system. The analysis performed leads to an exact solution of the general spheroidal problem in a formulation of stress concentration factors considering the different stiffness parameters of inclusion and matrix.

Probing active nematics with in situ microfabricated elastic inclusions

In this work, we report a direct measurement of the forces exerted by a tubulin/kinesin active nematic gel as well as its complete rheological characterization, including the quantification of its shear viscosity, η, and its activity parameter, α. For this, we develop a method that allows us to rapidly photo-polymerize compliant elastic inclusions in the continuously remodeling active system. Moreover, we quantitatively settle long-standing theoretical predictions, such as a postulated relationship encoding the intrinsic time scale of the active nematic in terms of η and α. In parallel, we infer a value for the nematic elasticity constant, K, by combining our measurements with the theorized scaling of the active length scale. On top of the microrheology capabilities, we demonstrate strategies for defect encapsulation, quantification of defect mechanics, and defect interactions, enabled by the versatility of the microfabrication strategy that allows to combine elastic motifs of different shapes and stiffnesses that are fabricated in situ.

First evaluation of stiff-in-soft host–inclusion systems: experimental synthesis of zircon inclusions in quartz crystals

Quartz crystals with zircon inclusions were synthesized using a piston-cylinder apparatus to experimentally evaluate the use of inclusions in “soft” host minerals for elastic thermobarometry. Synthesized zircon inclusion strains and, therefore, pressures (Pinc) were measured using Raman spectroscopy and then compared with the expected inclusion strains and pressures calculated from elastic models. Measured inclusion strains and inclusion pressures are systematically more tensile than the expected values and, thus, re-calculated entrapment pressures are overestimated. These discrepancies are not caused by analytical biases or assumptions in the elastic models and strain calculations. Analysis shows that inclusion strain discrepancies progressively decrease with decreasing experimental temperature in the α-quartz field. This behavior is consistent with inelastic deformation of the host–inclusion pairs induced by the development of large differential stresses during experimental cooling. Therefore, inclusion strains are more reliable for inclusions trapped at lower temperature conditions in the α-quartz field where there is less inelastic deformation of the host–inclusion systems. On the other hand, entrapment isomekes of zircon inclusions entrapped in the β-quartz stability field plot along the α–β quartz phase boundary, suggesting that the inclusion strains were mechanically reset at the phase boundary during experimental cooling and decompression. Therefore, inclusions contained in soft host minerals can be used for elastic thermobarometry and inclusions contained in β-quartz may provide constraints on the P–T at which the host–inclusion system crossed the phase boundary during exhumation.

Mean field homogenization schemes using averages of stress, strain and strain energy density and applications

The average of stress (AS) and average of strain (AN) were often used in the development of conventional mean-field homogenization (MFH) schemes, but their correlation with the average of strain energy density (AE) has rarely been discussed. In this work, we examined the correlations of AE with AS and AN, and found that AE can be reduced to the product of AS or AN under uniform boundary stress or linear surface displacement associated with a uniform strain field. By introducing a reference matrix and using any two among AE, AS and AN, three MFH schemes were developed: NS (using AN and AS); NE (using AN and AE); and SE (using AS and AE). Following the concept of the Mori-Tanaka's (MT) method, three MFH schemes obtained above were reduced to MT-NS, MT-NE and MT-SE, respectively. Effective elastic properties of composites consisting of isotropic elastic matrices and spherical isotropic elastic inclusions were predicted using these schemes and compared with experimental results. The comparison revealed that MT-SE and MT-NS can reasonably replicate the experimental results in the case of moderate particle volume fractions and ratios between the Young's moduli of matrix and the reinforcement particles, demonstrating the validity of the MFH schemes. This work is significant because it not only clarifies the correlation of AE with AN and AS, but also provides more options for the development of MFH schemes, which can enrich the theory and methods for the analysis of the effective properties of composites.

Plane Problem of the Theory of Elasticity on the Identification of Nodal Points of a Quadrature Inclusion

The problem of detection and identification of an elastic inclusion in an isotropic, linearly elastic plane is considered. It is assumed that constant stresses are set at infinity. It is also assumed that on some closed curve containing an inclusion inside, the acting forces and displacements are known quantities. For the case of quadrature domain occupied by the inclusion, a method for identifying its nodal points has been developed. The developed method is based on the application of the principle of reciprocity. Numerical examples are given.

A coupled acoustic/elastic discontinuous Galerkin finite element method: Application to ultrasonic imaging of 3D-printed synthetic materials

We present the derivation of upwind numerical fluxes for the space discontinuous Galerkin (dG) finite element method applied to the numerical modeling of wave propagation in multidimensional coupled acoustic/elastic media. The space dG method is formulated using the first-order velocity-pressure and velocity-stress systems for acoustic and elastic media, respectively. After eigenanalysis of the first-order hyperbolic systems highlighting the eigenmodes of wave propagation, the upwind numerical fluxes on the acoustic/acoustic and acoustic/elastic interface are obtained in terms of exact solutions of relevant Riemann problems. Thanks to the proposed approach, explicit closed-form expressions of the upwind numerical fluxes are obtained on the acoustic/elastic interface for the general case of multidimensional anisotropic heterogeneous solid media coupled with acoustic fluids. The developed numerical fluxes are validated by analytical/numerical comparison considering the example of an acoustic domain with a circular elastic inclusion. Finally, the coupled solver is used to perform a multiparametric study on the microstructure's echogenicity in a 3D-printed synthetic material under ultrasonic imaging.

On Numerical Solving of Junction Problem for the Thin Rigid and Elastic Inclusions in Elastic Body

On Numerical Solving of Junction Problem for the Thin Rigid and Elastic Inclusions in Elastic Body

Evaluation of the elastic field in phase‐field crystal simulations

AbstractThe phase‐field crystal model (PFC) describes crystal structures at diffusive timescales through a periodic order parameter representing the atomic density. One of its main features is that it naturally incorporates elastic and plastic deformation. To correctly interpret numerical simulation results or devise extensions related to the elasticity description, it is important to have direct access to the elastic field. In this work, we discuss its evaluation in classical PFC models based on the Swift–Hohenberg energy functional. We consider approaches where the stress field can be derived from the microscopic density field (i.e., the order parameter) and a simple novel numerical routine is proposed. By numerical simulations, we demonstrate that it overcomes some limitations of currently used methods. Moreover, we shed light on the elasticity description conveyed by classical PFC models, characterizing a residual stress effect present at equilibrium. We show explicitly and discuss the evaluation of the elastic fields in prototypical representative cases involving an elastic inclusion, a grain boundary, and dislocations.

Investigation of microstructure evolution accounting for crystal plasticity in the multiphase‐field method

AbstractRegarding microstructured materials, a quantitative prediction of phase transformation processes is highly desirable for a wide range of applications. With respect to polycrstalline materials, the plastic material behavior is commonly investigated using a crystal plasticity (CP) theory, since it accounts for the underlying microstructure, that is, slip systems of the crystal lattice. In classical continuum mechanics, grain boundaries (GBs) are commonly modeled as material singular surfaces. However, the tracking of moving GBs, present during phase transformation processes, is numerically challenging and costly. This can be circumvented by the use of a multiphase‐field method (MPFM), which provides a numerically highly efficient method for the treatment of moving interfaces, considered as diffuse interfaces of finite thickness. In this work, the microstructural evolution is investigated within the MPFM accounting for CP. The implementation of the constitutive material behavior within the diffuse interface region accounts for phase‐specific plastic fields and the jump condition approach. To improve the understanding of the impact of plastic deformation on the phase evolution, a single inclusion problem is analyzed. Within a plastically deformed matrix, the shape evolution of a purely elastic inclusion with a different Young's modulus, referred to as inhomogeneity, is investigated. It is shown, how the anisotropic plastic behavior affects the phase evolution. The resulting equilibrium shapes are illustrated and examined.

Modelling the time‐dependent mechanical behaviour of clay rocks based on meso‐ and micro‐structural viscous properties

AbstractClay rocks are multiphase porous media having a complex structure and behaviour characterised by heterogeneity, damage and viscosity, existing on a wide range of scales. The mesoscopic scale of mineral inclusions embedded in a clay matrix has an important role in the mechanisms of deformation under mechanical loading by cracking and creeping. This study introduces a micromechanical approach to model the time‐dependent mechanical behaviour of clay rocks. A heterogeneous clay rock is represented at the mesoscopic scale as a composite material consisting of rigid elastic mineral inclusions (quartz, calcite and pyrite) embedded in a clay matrix. To describe the damageable rock behaviour and its failure modes at the small scale, interfaces between different mineral phases and within the clay matrix are considered. Viscous effects are incorporated inside the clay aggregates, with intergranular microfractures propagating in the clay matrix, in order to investigate their contribution to the creep behaviour of clay rock at the macroscale. The mesostructure of the clay rock is represented in digital 2D Representative Elementary Areas (REAs). The overall mesoscale behaviour of the clay rock under mechanical solicitation is numerically obtained from the REA by computational homogenisation within a two‐scale finite element squared framework. Then, the model is validated at mesoscale against experimental data. The variability of the material response and the time evolution of the mineral interfacial damage state are investigated in relation to the small‐scale properties and failure, while considering mesostructure variability. The results can give some valuable insights into creep behaviour of the clay rock from a small‐scale perspective.

Modelling the creep behaviour and induced failure of clay rock from microscale viscosity to large-scale time-dependant gallery convergences using a multiscale numerical approach

The large-scale creep behaviour of Callovo-Oxfordian (COx) clay rock is modelled from small-scale viscous mechanisms of the rock using a multiscale numerical approach in the context of radioactive waste repositories. At the mesoscale, the saturated non-homogeneous rock is represented in digital 2D Representative Elementary Areas (REAs) as a cracked composite material consisting of rigid elastic mineral inclusions (quartz, calcite, and pyrite) embedded in a clay matrix. The modelling of these materials is considered within double-scale finite element framework, in which the homogenised responses of the REA to the enforced global kinematics serve as a numerical constitutive relation at the macroscale. To reproduce the rock creeping, two viscous mechanisms have been introduced to consider the creep of the clay matrix: either the viscoplasticity of the clay aggregates or the viscoelasticity of their contacts, or both. Firstly, laboratory biaxial creep tests of clay rock samples are simulated. A three-stage creep stage is reproduced and the creep failure process is discussed. Then, large-scale underground engineering structures are modelled to reproduce its time-dependent behaviour. It is found that the developed multiscale model is able to provide some valuable insights into the large-scale creep behaviour of clay rocks through the morphological and material small-scale characterization of REA.