Source author record

Gennady Mishuris

Gennady Mishuris appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

24 published item(s)

preprint2022arXiv

A temporal-averaging based approach to toughness homogenisation in heterogeneous material

A new approach to defining the effective fracture toughness for heterogeneous materials is proposed. This temporal averaging approach is process-dependent, incorporating the crack velocity and material toughness. The effectiveness of the new technique is investigated in the context of hydraulic fracture through heterogeneous rock with a periodic material toughness. The plane strain model is considered without fluid leak-off, to more easily investigate different regimes (toughness/viscosity). Numerical simulations are used to examine the effectiveness of the new homogenisation strategy, with comparison against the recently-proposed maximum toughness strategy. Simulations are conducted using an extremely effective (in house-built) time-space adaptive solver. The regimes in which each strategy is effective are determined.

preprint2022arXiv

Integral equation method for a Robin-type traction problem in a periodic domain

In this note, we consider a Robin-type traction problem for a linearly elastic body occupying an infinite periodically perforated domain. After proving the uniqueness of the solution we use periodic elastic layer potentials to show that the solution can be written as the sum of a single layer potential, a constant function and a linear function of the space variable. The density of the periodic single layer potential and the constant are identified as the unique solutions of a certain integral equation.

preprint2016arXiv

Analysis of dynamic failure of the discrete chain structure with non-local interactions

In the present work the steady-state crack propagation in a chain of oscillators with non-local interactions is considered. The interactions are modeled as linear springs while the crack is presented by the absence of extra springs. The problem is reduced to the Wiener-Hopf type and solution is presented in terms of inverse Fourier transform. It is shown that the non-local interactions may change the structure of the problem solution well-known from the classical local interactions formulation. In particular, it may change the range of the region of stable crack motion. The conclusions of the analysis are supported by numerical results. Namely, the observed phenomenon is partially clarified by evaluation of the structure profiles on the crack line ahead.

preprint2016arXiv

Energy Release Rate in hydraulic fracture: can we neglect an impact of the hydraulically induced shear stress?

A novel hydraulic fracture (HF) formulation is introduced which accounts for the hydraulically induced shear stress at the crack faces. It utilizes a general form of the elasticity operator alongside a revised fracture propagation condition based on the critical value of the energy release rate. It is shown that the revised formulation describes the underlying physics of HF in a more accurate way and is in agreement with the asymptotic behaviour of the linear elastic fracture mechanics. A number of numerical simulations by means of the universal HF algorithm previously developed in Wrobel & Mishuris (2015) are performed in order to: i) compare the modified HF formulation with its classic counterpart and ii) investigate the peculiarities of the former. Computational advantages of the revised HF model are demonstrated. Asymptotic estimations of the main solution elements are provided for the cases of small and large toughness. The modified formulation opens new ways to analyse the physical phenomenon of HF and also improves the reliability and efficiency of its numerical simulations.

preprint2016arXiv

Three-dimensional contact of TITH cartilage layers, a closed-form solution

Inhomogeneity and anisotropy play a crucial role in attributing articular cartilage its properties. The frictionless contact model constructed here consists in two thin biphasic transversely isotropic transversely homogeneous (TITH) cartilage layers firmly attached onto rigid substrates and shaped as elliptic paraboloids of different radii. Using asymptotic techniques, a solution to the deformation problem of such material has been recently obtained extending previous ones referred to homogeneous materials. The layer itself is thin in comparison with the size of the contact area and the observed time is shorter than the hydrogel characteristic time. The emerging three-dimensional contact problem is solved in closed-form and numerical benchmarks for constant and oscillating loads are given. The results are shown in terms of contact pressure and approach of the bones. The latter is derived to be directly proportional to the contact area. Existing experimental data are reinterpreted in view of the current model formulation. Comparisons are made with existing solutions for homogeneous biphasic materials in order to underline the functional importance of inhomogeneity in spreading the contact pressure distribution across the contact area. Particular attention is paid to the applicability of the retrieved formulas for interpreting measurements of in vivo experiments. Future directions are also prospected.

preprint2015arXiv

An asymptotic model for a thin bonded elastic layer coated with an elastic membrane

The deformation problem for a transversely isotropic elastic layer bonded to a rigid substrate and coated with a very thin elastic layer made of another transversely isotropic material is considered. The leading-order asymptotic models (for compressible and incompressible layers) are constructed based on the simplifying assumptions that the generalized plane stress conditions apply to the coating layer, and the flexural stiffness of the coating layer is negligible compared to its tensile stiffness.

preprint2015arXiv

An asymptotic model for the deformation of a transversely isotropic, transversely homogeneous biphasic cartilage layer

In the present paper, an asymptotic model is constructed for the short-time deformation of an articular cartilage layer modeled as transversely isotropic, transversely homogeneous (TITH) biphasic material. It is assumed that the layer thickness is relatively small compared with the characteristic size of the normal surface load applied to the upper surface of the cartilage layer, while the bottom surface is assumed to be firmly attached to a rigid impermeable substrate. In view of applications to articular contact problems it is assumed that the interstitial fluid is not allowed to escape through the articular surface.

preprint2015arXiv

Analysis of residual stresses in thermoelastic multilayer cylinders

The topic faced here is the modeling of an axisymmetric multilayer structure. An exact analytical formulation is proposed in the framework of the plane strain problem for a cylinder encircled by annular layers. Isotropic linear thermoelastic materials constitute the body. Perfect and imperfect contacts between the layers are made available for the analysis. The derived exact solution is compared to finite element simulations. Numerical applications are shown in order to study the dependency of the residual stress distribution on the constituent material properties during a cooling process. The role of residual stresses in brittle materials, particularly ceramics, is discussed.

preprint2015arXiv

Factorization of a class of matrix-functions with stable partial indices

A new effective method for factorization of a class of nonrational $n\times n$ matrix-functions with \emph{stable partial indices} is proposed. The method is a generalization of the one recently proposed by the authors which was valid for the canonical factorization only. The class of considered matrices is motivated by problems originated from applications. The properties and details of the asymptotic procedure are illustrated by examples. The efficiency of the procedure is highlighted by numerical results.

preprint2015arXiv

Particle velocity based universal algorithm for numerical simulation of hydraulic fractures

In the paper, we propose a new effective mathematical formulation and resulting universal numerical algorithm capable of tackling various HF models in the framework of a unified approach. The presented numerical scheme is not limited to any particular elasticity model or crack propagation regime. Its basic assumptions are: i) proper choice of independent and dependent variables (with the direct utilization of a new one - the reduced particle velocity), ii) tracing the fracture front by use of the speed equation which can be integrated in a closed form and sets an explicit relation between the crack propagation speed and the coefficients in the asymptotic expansion of the crack opening, iii) proper regularization techniques, iv) improved temporal approximation, v) modular algorithm architecture. The application of the new dependent variable, the reduced particle velocity, instead of the usual fluid flow rate, facilitates the computation of the crack propagation speed from the local relation based on the speed equation. As a result, the position of the crack front is accurately determined from an explicit formula derived from the speed equation. The underlying ideas employed in the algorithm are combined together producing a robust and efficient numerical scheme. Its performance is demonstrated using classical examples of 1D models for hydraulic fracturing: PKN and KGD under various fracture propagation regimes. Solution accuracy is verified against dedicated analytical benchmarks and other solutions available in the literature. Most of the ideas developed here, can be directly extended to more general 2D and 3D cases.

preprint2015arXiv

Universal hydrofracturing algorithm for shear-thinning fluids: particle velocity based simulation

A universal particle velocity based algorithm for simulating hydraulic fractures, previously proposed for Newtonian fluids, is extended to the class of shear-thinning fluids. The scheme is not limited to any particular elasticity operator or crack propagation regime. The computations are based on two dependent variables: the crack opening and the reduced particle velocity. The application of the latter facilitates utilization of the local condition of Stefan type (speed equation) to trace the fracture front. The condition is given in a general explicit form which relates the crack propagation speed (and the crack length) to the solution tip asymptotics. The utilization of a modular structure, and the adaptive character of its basic blocks, result in a flexible numerical scheme. The computational accuracy of the proposed algorithm is validated against a number of analytical benchmark solutions.

preprint2014arXiv

Brittle fracture in a periodic structure with internal potential energy. Spontaneous crack propagation

Spontaneous brittle fracture is studied based on the recently introduced model (Mishuris and Slepyan, Brittle fracture in a periodic structure with internal potential energy. Proc. Roy. Soc. A, in press). A periodic structure is considered, where only the prospective crack-path layer is specified as a discrete set of alternating initially stretched and compressed bonds. A bridged crack destroying initially stretched bonds may propagate under a certain level of the internal energy without external sources. The general analytical solution with the crack speed $-$ energy relation is presented in terms of the crack-related dynamic Green's function. For the anisotropic two-line chain and lattice considered earlier in quasi-statics, the dynamic problem is examined in detail. The crack speed is found to grow unboundedly as the energy approaches its upper limit. It is revealed that the spontaneous fracture can occur in the form of a pure bridged, partially bridged or fully open crack depending on the internal energy level. Generally, the steady-state mode of the crack propagation is found to be realised, whereas an irregular growth, clustering and the crack speed oscillations are detected in a vicinity of the lower bound of the energy.

preprint2014arXiv

Exact solution to a nonlinear heat conduction problem in doubly periodic 2D composite materials

An analytic solution to a stationary heat conduction problem in 2D unbounded doubly periodic composite materials with temperature dependent conductivities of their components is given. Corresponding nonlinear boundary value problem is reduced a Laplace equation with nonlinear transmission conditions. For special relationships between the conductivity coefficients of the matrix and inclusions, the problem is transformed to fully linear boundary value problem for doubly periodic analytic functions. This allows to reconstruct the solution of the originally nonlinear composite and to find its effective properties. The results are illustrated by numerical examples.

preprint2014arXiv

Forerunning mode transition in a continuous waveguide

We have discovered a new, forerunning mode transition as the periodic transition wave propagating in a uniform continuous waveguide. The latter is represented by an elastic beam separating from the elastic foundation under the action of sinusoidal waves. The critical displacement is the separation criterion. We show that the steady-state separation mode, where the separation front speed is independent of the wave amplitude, exists only in a bounded speed-dependent range of the wave amplitude. As the latter exceeds the critical value the steady-state mode is replaced by a more complicated regime with porous-like local separations emerging at a distance ahead of the main transition front. The crucial feature of this simple model is that the wave group speed is greater than the phase speed. The analytical solution allows us to show in detail how the steady-state mode transforms into the forerunning one. The established forerunning mode studied numerically manifests itself as a periodic process. As the incident wave amplitude grows the period decreases, while the separation front speed averaged over the period increases to the group velocity of the wave. In addition, the complete set of relations is presented for the waves excited by the oscillating force moving along the free beam, including the resonant case corresponding to a certain relation between the load's speed and frequency.

preprint2013arXiv

Efficient pseudo-spectral solvers for the PKN model of hydrofracturing

In the paper, a novel algorithm employing pseudo-spectral approach is developed for the PKN model of hydrofracturing. The respective solvers based on this approach compute both the solution and its temporal derivative. In comparison with conventional solvers, they demonstrate excellent cost effectiveness in terms of balance between the accuracy of computations and densities of the temporal and spatial meshes. Various leak-off regimes are considered.

preprint2013arXiv

Improved algorithm for analytical solution of the heat conduction problem in doubly periodic 2D composite materials

We consider a boundary value problem in unbounded 2D doubly periodic composite with circular inclusions having arbitrary constant conductivities. By introducing complex potentials, the boundary value problem for the Laplace equation is transformed to a special R-linear BVP for doubly periodic analytic functions. This problem is solved with use of the method of functional equations. The R-linear BVP is transformed to a system of functional equations. A new improved algorithm for solution of the system is proposed. It allows one not only to compute the average property but to reconstruct the solution components (temperature and flux) at an arbitrary point of the composite. Several computational examples are discussed in details demonstrating high efficiency of the method. Indirect estimate of the algorithm accuracy has been also provided.

preprint2012arXiv

Integral identities for a semi-infinite interfacial crack in 2D and 3D elasticity

The paper is concerned with the problem of a semi-infinite crack at the interface between two dissimilar elastic half-spaces, loaded by a general asymmetrical system of forces distributed along the crack faces. On the basis of the weight function approach and the fundamental reciprocal identity (Betti formula), we formulate the elasticity problem in terms of singular integral equations relating the applied loading and the resulting crack opening. Such formulation is fundamental in the theory of elasticity and extensively used to solve several problems in linear elastic fracture mechanics (for instance various classic crack problems in homogeneous and heterogeneous media). This formulation is also crucial in important recent multiphysics applications, where the elastic problem is coupled with other concurrent physical phenomena. A paradigmatic example is hydraulic fracturing, where the elasticity equations are coupled with fluid dynamics.

preprint2012arXiv

On modeling hydraulic fracture in proper variables: stiffness, accuracy, sensitivity

The problem of hydraulic fracture propagation is considered by using its recently suggested modified formulation in terms of the particle velocity, the opening in the proper degree, appropriate spatial coordinates and $\varepsilon$-regularization. We show that the formulation may serve for significant increasing the efficiency of numerical tracing the fracture propagation. Its advantages are illustrated by re-visiting the Nordgren problem. It is shown that the modified formulation facilitates (i) possibility to have various stiffness of differential equations resulting after spatial discretization, (ii) obtaining highly accurate and stable numerical results with moderate computational effort, and (iii) sensitivity analysis. The exposition is extensively illustrated by numerical examples.

preprint2012arXiv

Propagation of Slepyan's crack in a non-uniform elastic lattice

We model and derive the solution for the problem of a Mode I semi-infinite crack propagating in a discrete triangular lattice with bonds having a contrast in stiffness in the principal lattice directions. The corresponding Green's kernel is found and from this wave dispersion dependencies are obtained in explicit form. An equation of the Wiener-Hopf type is also derived and solved along the crack face, in order to compute the stress intensity factor for the semi-infinite crack. The crack stability is analysed via the evaluation of the energy release rate for different contrasts in stiffness of the bonds.

preprint2011arXiv

Interaction of an interfacial crack with linear micro-defects under out-of-plane shear loading

The interaction of an interfacial crack with small impurities is analysed on the basis of an asymptotic formula derived by the authors. The interaction between the main crack and the defects (e.g. small cracks or inclusions) is described asymptotically by analysing the dipole fields and the corresponding dipole matrices of the defects in question. The method is generic, and it serves interfacial cracks with general distributed loading on the crack faces, taking into account possible asymmetry in the boundary conditions, and in a particular configuration for a crack in a homogeneous medium results agree with those obtained earlier by Gong (1995). Shielding and amplification effects of the defects on the propagation of the main crack along the interface are investigated. Numerical computations based on the explicit analytical formulae show potential applications in the design of composite and fiber reinforced materials.

preprint2011arXiv

Mode III interfacial crack in the presence of couple stress elastic materials

In this paper we are concerned with the problem of a crack lying at the interface between dissimilar materials with microstructure undergoing antiplane deformations. The micropolar behaviour of the materials is described by the theory of couple stress elasticity developed by Koiter (1964). This constitutive model includes the characteristic lengths in bending and torsion and thus it is able to account for the underlying microstructure of the two materials. We perform an asymptotic analysis to investigate the behaviour of the solution near the crack tip. It turns out that the stress singularity at the crack tip is strongly influenced by the microstructural parameters and it may or may not show oscillatory behaviour depending on the ratio between the characteristic lengths.

preprint2011arXiv

Perturbation analysis of Mode III interfacial cracks advancing in a dilute heterogeneous material

The paper addresses the problem of a Mode III interfacial crack advancing quasi-statically in a heterogeneous composite material, that is a two-phase material containing elastic inclusions, both soft and stiff, and defects, such as microcracks, rigid line inclusions and voids. It is assumed that the bonding between dissimilar elastic materials is weak so that the interface is a preferential path for the crack. The perturbation analysis is made possible by means of the fundamental solutions (symmetric and skew-symmetric weight functions) derived in Piccolroaz et al. (2009). We derive the dipole matrices of the defects in question and use the corresponding dipole fields to evaluate effective tractions along the crack faces and interface to describe the interaction between the main interfacial crack and the defects. For a stable propagation of the crack, the perturbation of the stress intensity factor induced by the defects is then balanced by the elongation of the crack along the interface, thus giving an explicit asymptotic formula for the calculation of the crack advance. The method is general and applicable to interfacial cracks with general distributed loading on the crack faces, taking into account possible asymmetry in the boundary conditions. The analytical results are used to analyse the shielding and amplification effects of various types of defects in different configurations. Numerical computations based on the explicit analytical formulae allows for the analysis of crack propagation and arrest.