Eshebly tensors for a finite spherical domain with an axisymmetric inclusion

In recent papers the finite Eshelby tensors for a concentrically placed spherical inclusion in a finite spherical domain have been computed and applied to numerous micromechanical problems. The present work is the extension of the computation of finite Eshelby tensors to general inclusions that are axisymmetric with respect to enclosing spherical domain. The problem of finding the finite Eshelby tensors is transformed into the integral equation. It is shown in the paper that the integral equation has a unique solution. Existence of the solution is proved by exploiting the symmetry of the problem which induce invariant subspaces of the integral equation. In the particular case for a excentrically placed spherical inclusion the problem is explicitly solved. Using computer algebra the solution is found in a closed form up to the second order.     

Eshelby tensors for an ellipsoidal inclusion in a microstretch material

Eshelby tensors for an ellipsoidal inclusion in a microstretch material are derived in analytical form, involving only one-dimensional integral. As micropolar Eshelby tensor, the microstretch Eshelby tensors are not uniform inside of the ellipsoidal inclusion. However, different from micropolar Eshelby tensor, it is found that when the size of inclusion is large compared to the characteristic length of microstretch material, the microstretch Eshelby tensor cannot be reduced to the corresponding classical one. The reason for this is analyzed in details. It is found that under a pure hydrostatic loading, the bulk modulus of a microstretch material is not the same as the one in the corresponding classical material. A modified bulk modulus for the microstretch material is proposed, the microstretch Eshelby tensor is shown to be reduced to the modified classical Eshelby tensor at large size limit of inclusion. The fully analytical expressions of microstretch Eshelby tensors for a cylindrical inclusion are also derived.     

三维立方准晶椭球夹杂的Eshelby张量

准晶体颗粒复合智能材料具有优异的物化性能和应用前景。不同于传统的各向同性材料,三维立方准晶体材料包含声子场,相位子场,及声子-相位子耦合场。为更好地研究准晶体颗粒夹杂问题,揭示准晶体材料夹杂问题的物理现象,本文利用本征应变公式和柯西留数定理,考虑椭球体夹杂,获得了三维立方准晶材料夹杂问题的Eshelby张量,并给出了统一的表达式。进而,当三维立方准晶夹杂形状为球形、棒状、扁平状和带状时,获得了封闭形式的三维立方准晶Eshelby张量表达式。同时,给出了椭球体长径比变化时Eshelby张量的变化规律,这对研究准晶体颗粒夹杂问题具有重要的理论意义。     

Electroelastic behavior modeling of piezoelectric composite materials containing spatially oriented reinforcements

In this work, a modeling of electroelastic composite materials is proposed. The extension of the heterogeneous inclusion problem of Eshelby for elastic to electroelastic behavior is formulated in terms of four interaction tensors related to Eshelby’s electroelastic tensors. Analytical formulations of interaction tensors are presented for ellipsoidal inclusions. These tensors are basically used to derive the self-consistent model, Mori–Tanaka and dilute approaches. Numerical solutions are based on numerical computations of these tensors for various types of inclusions. Using the obtained results, effective electroelastic moduli of piezoelectric multiphase composites are investigated by an iterative procedure in the context of self-consistent scheme. Generalised Mori–Tanaka’s model and dilute approach are re-formulated and the three models are deeply analysed. Concentration tensors corresponding to each model are presented and relationships of effective coefficients are given. Numerical results of effective electroelastic moduli are presented for various types of piezoelectric inclusions and for various orientations and compared to existing experimental and theoretical ones.     

Solutions of inhomogeneity problems with graded shells and application to core-shell nanoparticles and composites

This paper first presents the Eshelby tensors and stress concentration tensors for a spherical inhomogeneity with a graded shell embedded in an alien infinite matrix. The solution is then specialized to inhomogeneous inclusions in finite spherical domains with fixed displacement or traction-free boundary conditions. The Eshelby tensors in the infinite and finite domains and the stress concentration tensors are especially useful for solving many problems in mechanics and materials science. This is demonstrated on two examples. In the first example, the strain distributions in core-shell nanoparticles with eigenstrains induced by lattice mismatches are calculated using the Eshelby tensors in the finite domains. In the second example, the Eshelby and stress concentration tensors in the three-phase configuration are used to formulate the generalized self-consistent prediction of the effective moduli of composites containing spherical particles within the framework of the equivalent inclusion method. The advantage of this micromechanical scheme is that, whilst its predictions are almost identical to the classical generalized self-consistent method and the third-order approximation, the expressions for the effective moduli have simple closed forms.     

Magnetoelectric Green's functions and their application to the inclusion and inhomogeneity problems

Explicit expressions of magnetoelectric Green's functions are obtained for a transversely isotropic medium exhibiting coupling between the static electric and magnetic fields utilizing the contour integral representation. Four Green's functions exist which represent the coupled static electric and magnetic response to a unit point electric or magnetic charge. The Green's functions are applied to analyze the inclusion and inhomogeneity problems in an infinite magnetoelectric medium, and explicit, closed form expressions are obtained for the Eshelby type tensors. The magnetoelectric Eshelby's tensors can be readily used in the solution of numerous problems in the mechanics and physics of magnetoelectric solids.     

Stationarity of the strain energy density for some classes of anisotropic solids

Homogeneous, anisotropic and linearly elastic solids, subjected to a given state of strain (or stress), are considered. The problem dealt with consists in finding the mutual orientations of the principal directions of strain to the material symmetry axes in order to make the strain energy density stationary. Such relative orientations are described through three Euler’s angles. When the stationarity problem is formulated for the generally anisotropic solid, it is shown that the necessary condition for stationarity demands for coaxiality of the stress and the strain tensors. From this feature, a procedure which leads to closed form solutions is proposed. To this end, tetragonal and cubic symmetry classes, together with transverse isotropy, are carefully dealt with, and for each case all the sets of Euler’s angles corresponding to critical points of the energy density are found and discussed. For these symmetries, three material parameters are then defined, which play a crucial role in ordering the energy values corresponding to each solution.     

New explicit expressions of the Hill polarization tensor for general anisotropic elastic solids

Except for particular cases, the classical expressions of the Eshelby or Hill polarization tensors, depend, respectively, on a simple or double integral for a fully anisotropic two-dimensional or three-dimensional elastic body. When the body is two-dimensional, we take advantage of Cauchy’s theory of residues to derive a new explicit expression which depends on the two pairs of complex conjugate roots of a quartic equation. If the body exhibits orthotropic symmetry, these roots are explicitly given as a function of the independent components of the elasticity tensor. Similarly, the double integral is reduced to a simple one when the body is three-dimensional. The corresponding integrand depends on the three pairs of complex conjugate roots of a sextic equation which reduces to a cubic one for orthotropic symmetry. This new expression improves significantly the computation times when the degree of anisotropy is high. For both two and three-dimensional bodies, degenerate cases are also studied to yield valid expressions in any events.     

On the generalized virial theorem and Eshelby tensors

The formal relationships between the scalar and tensorial virials and Eshelby tensors have been presently investigated. The key idea is to evaluate the Eshelby stress from discrete or atomistic simulations for a structured body, conceived as a numerical homogenization method to reconstitute the macroscopic continuum behavior in multiscale modelling approaches. Extending first the writing of the scalar virial to a material format, it is shown that the average of the elaborated scalar material virial is the trace of the (material) Eshelby stress. The spatial and material virials are further related to eachother in the framework of hyperelasticity, and a tensorial extension of the material virial is provided. Interpretation of those results from the microscopic point of view shows that Eshelby stress may be identified and calculated at the discrete level from the average of the virial tensor. Consideration of the material version of the virial theorem further leads to express Eshelby stress versus the average of the internal tensorial material virial and of the kinetic energy. The average scalar virial is further identified to the grand potential in a thermodynamic context. A definition of the material scalar virial for a second order continuum is lastly proposed, based on the identification of a second order Eshelby stress and in line with the second order Cauchy–Born rule.     

Inclusions and inhomogeneities in strain gradient elasticity with couple stresses and related problems

At small length scales and/or in presence of large field gradients, the implicit long wavelength assumption of classical elasticity breaks down. Postulating a form of second gradient elasticity with couple stresses as a suitable phenomenological model for small-scale elastic phenomena, we herein extend Eshelby’s classical formulation for inclusions and inhomogeneities. While the modified size-dependent Eshelby’s tensor and hence the complete elastic state of inclusions containing transformation strains or eigenstrains is explicitly derived, the corresponding inhomogeneity problem leads to integrals equations which do not appear to have closed-form solutions. To that end, Eshelby’s equivalent inclusion method is extended to the present framework in form of a perturbation series that then can be used to approximate the elastic state of inhomogeneities. The approximate scheme for inhomogeneities also serves as the basis for establishing expressions for the effective properties of composites in second gradient elasticity with couple stresses. The present work is expected to find application towards nano-inclusions and certain types of composites in addition to being the basis for subsequent non-linear homogenization schemes.     

General solution for Eshelby’s problem of 2D arbitrarily shaped piezoelectric inclusions

Eshelby’s problem of piezoelectric inclusions arises sometimes in exploiting the electromechanical coupling effect in piezoelectric media. For example, it intervenes in the nanostructure design of strained semiconductor devices involving strain-induced quantum dot (QD) and quantum wire (QWR) growth. Using the extended Stroh formalism, the present work gives a general analytical solution for Eshelby’s problem of two-dimensional arbitrarily shaped piezoelectric inclusions. The key step toward obtaining this general solution is the derivation of a simple and compact boundary integral expression for the eigenfunctions in the extended Stroh formalism applied to Eshelby’s problem. The simplicity and compactness of the boundary integral expression derived make it much less difficult to analytically tackle Eshelby’s piezoelectric problem for a large variety of non-elliptical inclusions. In the present work, explicit analytical solutions are obtained and detailed for all polygonal inclusions and for the inclusions characterized by Jordan’s curves and Laurent’s polynomials. By considering the piezoelectric material GaAs (110), the analytical solutions provided are illustrated numerically to verify the coincidence between different expressions, and to clarify the jump across the boundary of the inclusion and the singularity around the corner of the inclusion.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations(BIE)and solved with the newly developed boundary point method(BPM).The model is closely derived from the concept of the equivalent inclusion of Eshelby tensors.Eigenstrains are iteratively determined for each short.fiber embedded in the matrix with various properties via the Eshelby tensors,which can be readily obtained beforehand either through analytical or numerical means.As unknown variables appear only on the boundary of the solution domain,the solution scale of the inhomogeneity problem with the model is greatly reduced.This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM.The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element(RVE),showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Solutions to the Pólya–Szegö Conjecture and the Weak Eshelby Conjecture

Eshelby showed that if an inclusion is of elliptic or ellipsoidal shape then for any uniform elastic loading the field inside the inclusion is uniform. He then conjectured that the converse is true, that is, that if the field inside an inclusion is uniform for all uniform loadings, then the inclusion is of elliptic or ellipsoidal shape. We call this the weak Eshelby conjecture. In this paper we prove this conjecture in three dimensions. In two dimensions, a stronger conjecture, which we call the strong Eshelby conjecture, has been proved: if the field inside an inclusion is uniform for a single uniform loading, then the inclusion is of elliptic shape. We give an alternative proof of Eshelby’s conjecture in two dimensions using a hodographic transformation. As a consequence of the weak Eshelby’s conjecture, we prove in two and three dimensions a conjecture of Pólya and Szegö on the isoperimetric inequalities for the polarization tensors (PTs). The Pólya–Szegö conjecture asserts that the inclusion whose electrical PT has the minimal trace takes the shape of a disk or a ball.     

The Eshelby stress tensor,angular momentum tensor and dilatation flux in gradient elasticity

The Eshelby (static energy momentum) stress tensor, the angular momentum tensor and the dilatation flux are derived for anisotropic linear gradient elasticity in non-homogeneous materials. The divergence of these tensors gives the configurational forces, moments and work terms in gradient elasticity. There are several types of configurational forces, acting on the dislocation density and its gradient, on the inhomogeneities, proportional to the distortion, and linear and quadratic in the distortion gradient, and on the body force.     

局部裂纹损伤简支梁的曲率模态特性

将裂缝损伤简化成矩形凹槽,采用delta函数表示简支梁的裂纹损伤位置,得到了全梁范围内截面转动惯量和单位长度质量的表达式,建立了局部裂缝损伤简支梁的横向自由振动方程.利用摄动方法给出了裂纹摄动项的一般表达式,根据摄动项和完整梁都同时满足边界条件的特点,将一阶和二阶摄动项都表示成完整梁模态的线性组合,结合delta函数的性质,最终获得了受损简支梁的特征值和模态振型的解析表达式.最后,通过数值计算得到结构模态参数,对比了一阶摄动和二阶摄动对计算结果的影响,分析了不同阶固有频率和模态曲率的变动量,为简支梁的损伤监控和检测提供了理论依据.     

The energy of a growing elastic surface

The potential energy of the elastic surface of an elastic body which is growing by the coherent addition of material is derived. Several equivalent expressions are presented for the energy required to add a single atom, also known as the chemical potential. The simplest involves the Eshelby stress tensors for the bulk medium and for the surface. Dual Lagrangian/Eulerian expressions are obtained which are formally similar to each other. The analysis employs two distinct types of variations to derive the governing bulk and surface equations for an accreting elastic solid. The total energy of the system is assumed to comprise bulk and surface energies, while the presence of an external medium can be taken into account through an applied surface forcing. A detailed account is given of the various formulations possible in material and current coordinates, using four types of bulk and surface stresses: the Piola-Kirchhoff stress, the Cauchy stress, the Eshelby stress and a fourth, called the nominal energy-momentum stress. It is shown that inhomogeneity surface forces arise naturally if the surface energy density is allowed to be position dependent.     

On the derivative of the square root of a tensor and Guo's rate theorems

Formulas for the derivative of the square root of a positive definite, symmetric, second-order tensor are derived and used to obtain expressions for the material time derivatives of the right and left stretch tensors.     

Geometric shapes of the interface surface of bicomponent flows between two concentric rotating cylinders

In this paper, the shape problem of interface of bicomponent flows between two concentric rotating cylinders is investigated. With tensor analysis, the problem is reduced to an energy functional isoperimetric problem when neglecting the effects of the dissipative energy caused by viscosity. We derive the associated Euler-Lagrangian equation, which is a nonlinear elliptic boundary value problem of the second order. Moreover, by considering the effects of the dissipative energy, we propose another total energy functional to characterize the geometric shape of the interface, and obtain the corresponding Euler-Lagrangian equation, which is also a nonlinear elliptic boundary value problem of the second order. Thus, the problem of the geometric shape is converted into a nonlinear boundary value problem of the second order in both cases.     

On configurational compatibility and multiscale energy momentum tensors

In this work the continuum theory of defects has been revised through the development of kinematic defect potentials. These defect potentials and their corresponding variational principles provide a basis for constructing a new class of conservation laws associated with the compatibility conditions of continua. These conservation laws represent configurational compatibility conditions which are independent of the constitutive behavior of the continuum. They lead to the development of a new concept termed configurational compatibility, dual to the concept of configurational force. The contour integral of the corresponding conserved quantity is path-independent, if the domain encompassed by the integral is defect-free. It is shown that the Peach-Koehler force can be recovered as one of these invariant integrals. Based on the proposed defect potentials and their corresponding defect energies, two-field multiscale mixed variational principles can be employed to construct multiscale energy momentum tensors. An application is outlined in the form of a mode III elasto-plastic crack problem for which the new configurational quantities are calculated.