Researcher profile

A. Piccolroaz

A. Piccolroaz contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
31works
0followers
7topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

31 published item(s)

preprint2022arXiv

A way to hypo-elastic artificial materials without a strain potential and displaying flutter instability

Cauchy-elastic solids include hyper-elasticity and a subset of elastic materials for which the stress does not follow from a scalar strain potential. More in general, hypo-elastic materials are only defined incrementally and comprise Cauchy-elasticity. Infringement of the hyper-elastic 'dogma' is so far unattempted and normally believed to be impossible, as it apparently violates thermodynamics, because energy may be produced in closed strain cycles. Contrary to this belief, we show that non-hyper-elastic behavior is possible and we indicate the way to a practical realization of this new concept. In particular, a design paradigm is established for artificial materials where follower forces, so far ignored in homogenization schemes, are introduced as loads prestressing an elastic two-dimensional grid made up of linear elastic rods (reacting to elongation, flexure and shear). A rigorous application of Floquet-Bloch wave asymptotics yields an unsymmetric acoustic tensor governing the \textit{incremental} dynamics of the effective material. The latter is therefore the incremental response of a hypo-elastic solid, which does not follow from a strain potential and thus does not belong to hyper-elasticity. Through the externally applied follower forces (which could be originated via interaction with a fluid, or gas, or by application of Coulomb friction, or non-holonomic constraints), the artificial material may adsorb/release energy from/to the environment and therefore produce energy in a closed strain loop, without violating any rule of thermodynamics. The solid is also shown to display flutter, a material instability corresponding to a Hopf bifurcation, which was advocated as possible in plastic solids, but never experimentally found and so far believed to be impossible in elasticity.

preprint2022arXiv

Tensile material instabilities in elastic beam lattices lead to a bounded stability domain

Homogenization of the incremental response of grids made up of preloaded elastic rods leads to homogeneous effective continua which may suffer macroscopic instability, occurring at the same time in both the grid and the effective continuum. This instability corresponds to the loss of ellipticity in the effective material and the formation of localized responses as, for instance, shear bands. Using lattice models of elastic rods, loss of ellipticity has always been found to occur for stress states involving compression of the rods, as usually these structural elements buckle only under compression. In this way, the locus of material stability for the effective solid is unbounded in tension, i.e. the material is always stable for a tensile prestress. A rigorous application of homogenization theory is proposed to show that the inclusion of sliders (constraints imposing axial and rotational continuity, but allowing shear jumps) in the grid of rods leads to loss of ellipticity in tension, so that the locus for material instability becomes bounded. This result explains (i.) how to design elastic materials passible of localization of deformation and shear banding for all radial stress paths; (ii.) how for all these paths a material may fail by developing strain localization and without involving cracking.

preprint2020arXiv

Bifurcation of elastic solids with sliding interfaces

Lubricated sliding contact between soft solids is an interesting topic in biomechanics and for the design of small-scale engineering devices. As a model of this mechanical set-up, two elastic nonlinear solids are considered jointed through a frictionless and bilateral surface, so that continuity of the normal component of the Cauchy traction holds across the surface, but the tangential component is null. Moreover, the displacement can develop only in a way that the bodies in contact do neither detach, nor overlap. Surprisingly, this finite strain problem has not been correctly formulated until now, so that this formulation is the objective of the present paper. The incremental equations are shown to be non-trivial and different from previously (and erroneously) employed conditions. In particular, an exclusion condition for bifurcation is derived to show that previous formulations based on frictionless contact or 'spring-type' interfacial conditions are not able to predict bifurcations in tension, while experiments -- one of which, ad hoc designed, is reported -- show that these bifurcations are a reality and become possible when the correct sliding interface model is used. The presented results introduce a methodology for the determination of bifurcations and instabilities occurring during lubricated sliding between soft bodies in contact.

preprint2020arXiv

Incremental constitutive tensors and strain localization for prestressed elastic lattices: Part I -- quasi-static response

A lattice of elastic rods organized in a parallelepiped geometry can be axially loaded up to an arbitrary amount without distortion and then be subject to incremental displacements. Using quasi-static homogenization theory, this lattice can be made equivalent to a prestressed elastic solid subject to incremental deformation, in such a way to obtain extremely localized mechanical responses. These responses can be analyzed with reference to a mechanical model which can, in principle, be realized, so that features such as for instance shear bands inclination, or emergence of a single shear band, or competition between micro (occurring in the lattice but not in the equivalent solid) and macro (present in both the lattice and the equivalent continuum) instabilities become all designable features. The analysis of localizations is performed using a Green's function-based perturbative approach to highlight the correspondence between micromechanics of the composite and homogenized response of the equivalent solid. The presented results, limited to quasi-static behaviour, provide a new understanding of strain localization in a continuum and open new possibilities for the realization and experimentation of materials exhibiting these extreme mechanical behaviours. Dynamic homogenization and vibrational localization are deferred to Part II of this study.

preprint2019arXiv

Redirection of a crack driven by viscous fluid taking into account plastic effects in the process zone

In this paper the problem of redirection of a crack driven by viscous fluid under mixed mode loading is considered. The loading includes the classical Modes I - III and the hydraulically induced tangential traction on the fracture walls. The effect of the plastic deformation of the near tip zone is accounted for. Different criteria to determine the fracture deflection angle are examined and compared.

preprint2016arXiv

Computational modelling and experimental validation of industrial forming processes by cold pressing of aluminum silicate powder

A novel approach to the modelling and simulation of the industrial compaction process of ceramic powders is proposed, based on a combination of: (i) continuum mechanics modelling of the constitutive response of the material; (ii) finite element discretization and computer implementation of the mechanical model; (iii) parametric identification by a multi-objective optimization of simulated experimental tests. The capabilities of the proposed approach are highlighted through computer simulations of realistic industrial compaction processes, namely, the forming of an axisymmetric tablet and of a three-dimensional ceramic tile. The presented methods and the pointed out results pave the way for the introduction of so-called virtual prototyping into the industrial practice of ceramic forming processes.

preprint2016arXiv

Overall thermomechanical properties of layered materials for energy devices applications

This paper is concerned with the analysis of effective thermomechanical properties of multi- layered materials of interest for solid oxide fuel cells (SOFC) and lithium ions batteries fabrication. The recently developed asymptotic homogenization procedure is applied in order to express the overall thermoelastic constants of the first order equivalent continuum in terms of microfluctuations functions, and these functions are obtained by the solution of the corresponding recursive cell problems. The effects of thermal stresses on periodic multi-layered thermoelastic composite reproducing the characteristics of solid oxide fuel cells (SOFC-like) are studied assuming periodic body forces and heat sources, and the solution derived by means of the asymptotic homogenization approach is compared with the results obtained by finite elements analysis of the associate heterogeneous material.

preprint2015arXiv

Boundary integral formulation for interfacial cracks in thermodiffusive bimaterials

An original boundary integral formulation is proposed for the problem of a semi-infinite crack at the interface between two dissimilar elastic materials in the presence of heat flows and mass diffusion. Symmetric and skew-symmetric weight function matrices are used together with a generalized Betti's reciprocity theorem in order to derive a system of integral equations that relate the applied loading, the temperature and mass concentration fields, the heat and mass fluxes on the fracture surfaces and the resulting crack opening. The obtained integral identities can have many relevant applications, such as for the modelling of crack and damage processes at the interface between different components in electrochemical energy devices characterized by multi-layered structures (solid oxide fuel cells and lithium ions batteries).

preprint2015arXiv

General transmission conditions for thin elasto-plastic pressure-dependent interphase between dissimilar materials

We consider a thin soft adhesive interphase between dissimilar elastic media. The material of the intermediate layer is modelled by elasto-plastic pressure-sensitive constitutive law. An asymptotic procedure, together with a novel formulation of the deformation theory of plasticity for pressure-sensitive materials, is used in order to derive nonlinear transmission conditions for the corresponding imperfect zero-thickness interface. A FEM analysis of the original three-phase structure is performed to validate the transmission conditions for the simplified bimaterial structure.

preprint2015arXiv

Multiscale asymptotic homogenization analysis of thermo-diffusive composite materials

In this paper an asymptotic homogenization method for the analysis of composite materials with periodic microstructure in presence of thermodiffusion is described. Appropriate down-scaling relations correlating the microscopic fields to the macroscopic displacements, temperature and mass concentration are introduced. The effects of the material inhomogeneities are described by perturbation functions derived from the solution of recursive cell problems. Exact expressions for the overall elastic and thermodiffusive constants of the equivalent first order thermodiffusive continuum are derived. The proposed approach is applied to the case of a two-dimensional bi-phase orthotropic layered material, where the effective elastic and thermodiffusive properties can be determined analytically. Considering this illustrative example and assuming periodic body forces, heat and mass sources acting on the medium, the solution performed by the first order homogenization approach is compared with the numerical results obtained by the heterogeneous model.

preprint2015arXiv

Plastically-driven variation of elastic stiffness in green bodies during powder compaction: Part I Experiments and elastoplastic coupling

Cold compaction of ceramic powders is driven by plastic strain, during which the elastic stiffness of the material progressively increases from values typical of granular matter to those representative of a fully dense solid. This increase of stiffness strongly affects the mechanical behaviour of the green body and is crucial in the modelling of forming processes for ceramics. A protocol for ultrasonic experimental investigation (via P and S waves transmission) is proposed to quantify the elastic constants (Young modulus and Poisson's ratio) as functions of the forming pressure. Experimental results performed in uniaxial strain allow for the introduction of laws that describe the variation of the elastic constants during densification. These laws are motivated in terms of elastoplastic coupling through the simulation of an isostatic pressure compaction process of alumina powder. A micromechanical explanation of the stiffening of elastic properties during densification is deferred to Part II of this study.

preprint2015arXiv

Plastically-driven variation of elastic stiffness in green bodies during powder compaction. Part II: Micromechanical modelling

A micromechanical approach is set-up to analyse the increase in elastic stiffness related to development of plastic deformation (the elastoplastic coupling concept) occurring during the compaction of a ceramic powder. Numerical simulations on cubic (square for 2D) and hexagonal packings of elastoplastic cylinders and spheres validate both the variation of the elastic modulus with the forming pressure and the linear dependence of it on the relative density as experimentally found in Part~I of this study, while the dependence of the Poisson's ratio on the green's density is only qualitatively explained.

preprint2015arXiv

Strain localization and shear banding in ductile materials

A model of a shear band as a zero-thickness nonlinear interface is proposed and tested using finite element simulations. An imperfection approach is used in this model where a shear band, that is assumed to lie in a ductile matrix material (obeying von Mises plasticity with linear hardening), is present from the beginning of loading and is considered to be a zone in which yielding occurs before the rest of the matrix. This approach is contrasted with a perturbative approach, developed for a J$_2$-deformation theory material, in which the shear band is modelled to emerge at a certain stage of a uniform deformation. Both approaches concur in showing that the shear bands (differently from cracks) propagate rectilinearly under shear loading and that a strong stress concentration should be expected to be present at the tip of the shear band, two key features in the understanding of failure mechanisms of ductile materials.

preprint2014arXiv

A simple and robust elastoplastic constitutive model for concrete

An elasto-plastic model for concrete, based on a recently-proposed yield surface and simple hardening laws, is formulated, implemented, numerically tested and validated against available test results. The yield surface is smooth and particularly suited to represent the behaviour of rock-like materials, such as concrete, mortar, ceramic and rock. A new class of isotropic hardening laws is proposed, which can be given both an incremental and the corresponding finite form. These laws describe a smooth transition from linear elastic to plastic behaviour, incorporating linear and nonlinear hardening, and may approach the perfectly plastic limit in the latter case. The reliability of the model is demonstrated by its capability of correctly describing the results yielded by a number of well documented triaxial tests on concrete subjected to various confinement levels. Thanks to its simplicity, the model turns out to be very robust and well suited to be used in complex design situations, as those involving dynamic loads.

preprint2014arXiv

Boundary integral formulation for cracks at imperfect interfaces

We consider an infinite bi-material plane containing a semi-infinite crack situated on a soft imperfect interface. The crack is 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, we derive the corresponding boundary integral formulation, relating physical quantities. The boundary integral equations derived in this paper in the imperfect interface setting show a weak singularity, in contrast to the perfect interface case, where the kernel is of the Cauchy type. We further present three alternative variants of the boundary integral equations which offer computationally favourable alternatives for certain sets of parameters.

preprint2014arXiv

Dispersion and localization in structured Rayleigh beams

This paper brings a comparative analysis between dynamic models of couple-stress elastic materials and structured Rayleigh beams on a Winkler foundation. Although physical phenomena have different physical origins, the underlying equations appear to be similar, and hence mathematical models have a lot in common. In the present work, our main focus is on the analysis of dispersive waves, band-gaps and localized waveforms in structured Rayleigh beams. The Rayleigh beam theory includes the effects of rotational inertia which are neglected in the Euler--Bernoulli beam theory. This makes the approach applicable to higher frequency regimes. Special attention is given to waves in pre-stressed Rayleigh beams on elastic foundations.

preprint2014arXiv

Effective elastic properties of planar SOFCs: A non-local dynamic homogenization approach

The focus of the article is on the analysis of effective elastic properties of planar Solid Oxide Fuell Cell (SOFC) devices. An ideal periodic multi-layered composite (SOFC-like) reproducing the overall properties of multi-layer SOFC devices is defined. Adopting a non-local dynamic homogenization method, explicit expressions for overall elastic moduli and inertial terms of this material are derived in terms of micro-fluctuation functions. These micro-fluctuation function are then obtained solving the cell problems by means of finite element techniques. The effects of the temperature variation on overall elastic and inertial properties of the fuel cells are studied. Dispersion relations for acoustic waves in SOFC-like multilayered materials are derived as functions of the overall constants, and the results obtained by the proposed computational homogenization approach are compared with those provided by rigorous Floquet-Boch theory. Finally, the influence of the temperature and of the elastic properties variation on the Bloch spectrum is investigated.

preprint2014arXiv

Elastoplastic coupling to model cold ceramic powder compaction

The simulation of industrial processes involving cold compaction of powders allows for the optimization of the production of both traditional and advanced ceramics. The capabilities of a constitutive model previously proposed by the authors are explored to simulate simple forming processes, both in the small and in the large strain formulation. The model is based on the concept of elastoplastic coupling providing a relation between density changes and variation of elastic properties and has been tailored to describe the transition between a granular ceramic powder and a dense green body. Finite element simulations have been compared with experiments on an alumina ready-to-press powder and an aluminum silicate spray-dried granulate. The simulations show that it is possible to take into account friction at the die wall and to predict the state of residual stress, density distribution and elastic properties in the green body at the end of the forming process.

preprint2014arXiv

Experimental investigation of the elastoplastic response of aluminum silicate spray dried powder during cold compaction

Mechanical experiments have been designed and performed to investigate the elasto-plastic behaviour of green bodies formed from an aluminum silicate spray dried powder used for tiles production. Experiments have been executed on samples obtained from cold compaction into a cylindrical mould and include: uniaxial strain, equi-biaxial flexure and high-pressure triaxial compression/extension tests. Two types of powders have been used to realize the green body samples, differing in the values of water content, which have been taken equal to those usually employed in the industrial forming of traditional ceramics. Yielding of the green body during compaction has been characterized in terms of yield surface shape, failure envelope, and evolution of cohesion and void ratio with the forming pressure, confirming the validity of previously proposed constitutive models for dense materials obtained through cold compaction of granulates.

preprint2014arXiv

Finite-strain formulation and FE implementation of a constitutive model for powder compaction

A finite-strain formulation is developed, implemented and tested for a constitutive model capable of describing the transition from granular to fully dense state during cold forming of ceramic powder. This constitutive model (as well as many others employed for geomaterials) embodies a number of features, such as pressure-sensitive yielding, complex hardening rules and elastoplastic coupling, posing considerable problems in a finite-strain formulation and numerical implementation. A number of strategies are proposed to overcome the related problems, in particular, a neo-Hookean type of modification to the elastic potential and the adoption of the second Piola-Kirchhoff stress referred to the intermediate configuration to describe yielding. An incremental scheme compatible with the formulation for elastoplastic coupling at finite strain is also developed, and the corresponding constitutive update problem is solved by applying a return mapping algorithm.

preprint2014arXiv

Implicit yield function formulation for granular and rock-like materials

The constitutive modelling of granular, porous and quasi-brittle materials is based on yield (or damage) functions, which may exhibit features (for instance, lack of convexity, or branches where the values go to infinity, or false elastic domains) preventing the use of efficient return-mapping integration schemes. This problem is solved by proposing a general construction strategy to define an implicitly defined convex yield function starting from any convex yield surface. Based on this implicit definition of the yield function, a return-mapping integration scheme is implemented and tested for elastic-plastic (or -damaging) rate equations. The scheme is general and, although it introduces a numerical cost when compared to situations where the scheme is not needed, is demonstrated to perform correctly and accurately.

preprint2014arXiv

Integration algorithms of elastoplasticity for ceramic powder compaction

Inelastic deformation of ceramic powders (and of a broad class of rock-like and granular materials), can be described with the yield function proposed by Bigoni and Piccolroaz (2004, Yield criteria for quasibrittle and frictional materials. Int. J. Solids and Structures, 41, 2855-2878). This yield function is not defined outside the yield locus, so that 'gradient-based' integration algorithms of elastoplasticity cannot be directly employed. Therefore, we propose two ad hoc algorithms: (i.) an explicit integration scheme based on a forward Euler technique with a 'centre-of-mass' return correction and (ii.) an implicit integration scheme based on a 'cutoff-substepping' return algorithm. Iso-error maps and comparisons of the results provided by the two algorithms with two exact solutions (the compaction of a ceramic powder against a rigid spherical cup and the expansion of a thick spherical shell made up of a green body), show that both the proposed algorithms perform correctly and accurately.

preprint2014arXiv

Remarks on the energy release rate for an antiplane moving crack in couple stress elasticity

This paper is concerned with the steady-state propagation of an antiplane semi-infinite crack in couple stress elastic materials. A distributed loading applied at the crack faces and moving with the same velocity of the crack tip is considered, and the influence of the loading profile variations and microstructural effects on the dynamic energy release rate is investigated. The behaviour of both energy release rate and maximum total shear stress when the crack tip speed approaches the critical speed (either that of the shear waves or that of the localised surface waves) is studied. The limit case corresponding to vanishing characteristic scale lengths is addressed both numerically and analytically by means of a comparison with classical elasticity results.

preprint2014arXiv

Steady-State Propagation of a Mode II Crack in Couple Stress Elasticity

The present work deals with the problem of a semi-infinite crack steadily propagating in an elastic body subject to plane-strain shear loading. It is assumed that the mechanical response of the body is governed by the theory of couple-stress elasticity including also micro-rotational inertial effects. This theory introduces characteristic material lengths in order to describe the pertinent scale effects that emerge from the underlying microstructure and has proved to be very effective for modeling complex microstructured materials. It is assumed that the crack propagates at a constant sub-Rayleigh speed. An exact full field solution is then obtained based on integral transforms and the Wiener-Hopf technique. Numerical results are presented illustrating the dependence of the stress intensity factor and the energy release rate upon the propagation velocity and the characteristic material lengths in couple-stress elasticity. The present analysis confirms and extends previous results within the context of couple-stress elasticity concerning stationary cracks by including inertial and micro-inertial effects.

preprint2013arXiv

Integral identities for a semi-infinite interfacial crack in anisotropic elastic bimaterials

The focus of the article is on the analysis of a semi-infinite crack at the interface between two dissimilar anisotropic elastic materials, loaded by a general asymmetrical system of forces acting on the crack faces. Recently derived symmetric and skew-symmetric weight function matrices are introduced for both plane strain and antiplane shear cracks, and used together with the fundamental reciprocal identity (Betti formula) in order to formulate the elastic fracture problem in terms of singular integral equations relating the applied loading and the resulting crack opening. The proposed compact formulation can be used to solve many problems in linear elastic fracture mechanics (for example various classic crack problems in homogeneous and heterogeneous anisotropic media, as piezoceramics or composite materials). This formulation is also fundamental in many multifield theories, where the elastic problem is coupled with other concurrent physical phenomena.

preprint2013arXiv

On fracture criteria for dynamic crack propagation in elastic materials with couple stresses

The focus of the article is on fracture criteria for dynamic crack propagation in elastic materials with microstructures. Steady-state propagation of a Mode III semi-infinite crack subject to loading applied on the crack surfaces is considered. The micropolar behavior of the material is described by the theory of couple-stress elasticity developed by Koiter. This constitutive model includes the characteristic lengths in bending and torsion, and thus it is able to account for the underlying microstructures of the material. Both translational and micro-rotational inertial terms are included in the balance equations, and the behavior of the solution near to the crack tip is investigated by means of an asymptotic analysis. The asymptotic fields are used to evaluate the dynamic J-integral for a couple-stress material, and the energy release rate is derived by the corresponding conservation law. The propagation stability is studied according to the energy-based Griffith criterion and the obtained results are compared to those derived by the application of the maximum total shear stress criterion.

preprint2013arXiv

Perturbation analysis for an imperfect interface crack problem using weight function techniques

We analyse a problem of anti-plane shear in a bi-material plane containing a semi-infinite crack situated on a soft imperfect interface. The plane also contains a small thin inclusion (for instance an ellipse with high eccentricity) whose influence on the propagation of the main crack we investigate. An important element of our approach is the derivation of a new weight function (a special solution to a homogeneous boundary value problem) in the imperfect interface setting. The weight function is derived using Fourier transform and Wiener-Hopf techniques and allows us to obtain an expression for an important constant (which may be used in a fracture criterion) that describes the leading order of tractions near the crack tip for the unperturbed problem. We present computations that demonstrate how this constant varies depending on the extent of interface imperfection and contrast in material stiffness. We then perform perturbation analysis to derive an expression for the change in the leading order of tractions near the tip of the main crack induced by the presence of the small defect, whose sign can be interpreted as the inclusion's presence having an amplifying or shielding effect on the propagation of the main crack.

preprint2012arXiv

Steady-state propagation of a Mode III crack in couple stress elastic materials

This paper is concerned with the problem of a semi-infinite crack steadily propagating in an elastic solid with microstructures subject to antiplane loading applied on the crack surfaces. The loading is moving with the same constant velocity as that of the crack tip. We assume subsonic regime, that is the crack velocity is smaller than the shear wave velocity. The material behaviour is described by the indeterminate theory of couple stress elasticity developed by Koiter. This constitutive model includes the characteristic lengths in bending and torsion and thus it is able to account for the underlying microstructure of the material as well as for the strong size effects arising at small scales and observed when the representative scale of the deformation field becomes comparable with the length scale of the microstructure, such as the grain size in a polycrystalline or granular aggregate. The present analysis confirms and extends earlier results on the static case by including the effects of crack velocity and rotational inertia. By adopting the criterion of maximum total shear stress, we discuss the effects of microstructural parameters on the stability of crack propagation.

preprint2011arXiv

Mode III crack propagation in a bimaterial plane driven by a channel of small line defects

We consider the quasi-static propagation of a Mode III crack along the interface in a bimaterial plane containing a finite array of small line defects (microcracks and rigid line inclusions). The microdefects are arranged to form a channel around the interface that can facilitate (or prevent) the crack propagation. The two dissimilar elastic materials are assumed to be weakly bonded, so that there is no kinking of the main crack from the straight path. On the basis of asymptotic formulae obtained by the authors, the propagation is analysed as a perturbation problem and the incremental crack advance is analytically derived at each position of the crack tip along the interface relative to the position of the defects. Numerical examples are provided showing potential applications of the proposed approach in the analysis of failure of composite materials. Extension to the case of infinite number of defects is discussed.

preprint2010arXiv

Perturbation of Mode III interfacial cracks

We consider the perturbation problem of a Mode III interfacial crack. The perturbation is of geometrical type and can be both perturbation of the crack faces and perturbation of the interface, which can deviate from the initial straight line configuration. Asymptotic formulae are derived for the first-order perturbation of the stress intensity factor. It is shown that, due to the unsymmetrical nature of the problem, the Mode III skew-symmetric weight function derived in Piccolroaz et al. (2009) is essential for the derivation of the correct asymptotic formulae. To illustrate the method, we present the numerical results for different geometrical perturbations of a half-plane interfacial crack in an infinite bimaterial structure. Discussion on the extension of the method to finite bodies is also presented.

preprint2009arXiv

Symmetric and skew-symmetric weight functions in perturbation models of 2D interfacial cracks

In this paper we address the vector problem of a 2D half-plane interfacial crack loaded by a general asymmetric distribution of forces acting on its faces. It is shown that the general integral formula for the evaluation of stress intensity factors, as well as high-order terms, requires both symmetric and skew-symmetric weight function matrices. The symmetric weight function matrix is obtained via the solution of a Wiener-Hopf functional equation, whereas the derivation of the skew-symmetric weight function matrix requires the construction of the corresponding full field singular solution. The weight function matrices are then used in the perturbation analysis of a crack advancing quasi-statically along the interface between two dissimilar media. A general and rigorous asymptotic procedure is developed to compute the perturbations of stress intensity factors as well as high-order terms.