On the stress singularities generated by anisotropic eigenstrains and the hydrostatic stress due to annular inhomogeneities

The problems of singularity formation and hydrostatic stress created by an inhomogeneity with eigenstrain in an incompressible isotropic hyperelastic material are considered. For both a spherical ball and a cylindrical bar with a radially symmetric distribution of finite possibly anisotropic eigenstrains, we show that the anisotropy of these eigenstrains at the center (the center of the sphere or the axis of the cylinder) controls the stress singularity. If they are equal at the center no stress singularity develops but if they are not equal then stress always develops a logarithmic singularity. In both cases, the energy density and strains are everywhere finite. As a related problem, we consider annular inclusions for which the eigenstrains vanish in a core around the center. We show that even for an anisotropic distribution of eigenstrains, the stress inside the core is always hydrostatic. We show how these general results are connected to recent claims on similar problems in the limit of small eigenstrains.     

Variational formulation of the equivalent eigenstrain method with an application to a problem with radial eigenstrains

In this paper a variational formulation of the equivalent eigenstrain method is established. A functional of the Hashin–Shtrikman type is proposed such that the solution of the equivalent eigenstrain equation is a unique minimizer of the functional. Moreover, it is also shown that the equivalent eigenstrain equation is the Euler–Lagrange equation of the potential energy of the inclusions. An approximate solution of the equivalent eigenstrain equation is then found as a minimizer of the functional on a finite dimensional span of basic eigenstrains. Special attention is paid to possible symmetries of the problem. The variational formulation is illustrated by determination of effective linear elastic properties. In particular, material with a simple cubic microstructure is considered in detail. A solution for the polynomial radial basic eigenstrains approximation is found. In particular, for the homogeneous eigenstrain approximation, the effective moduli are derived in an exact closed form.     

The axisymmetrical elastic field of an ellipsoidal inclusion with a slipping interface

This paper is concerned with the axisymmetrical elastic fields caused by an ellipsoidal inclusion with a slipping interface which undergoes a uniform eigenstrain. The problem is solved under a revolving ellipsoidal coordinate with the aid of Papkovich-Neuber general dipacement formula. In contrast to the perfectly bonded interface, when the interface between the inclusion and the matrix cannot sustain shear stress, and is free to slip, the solution cannot be expressed in closed form and involves infinite series. Therefore, the results are illustrated by numerical examples.     

弱界面异质球形夹杂本征应变问题的解析解

本文研究了具有线弹簧弱界面的异质球形夹杂的本征应变问题，所采用的线弹簧界面模型既能界面的切线方向滑动，又能考虑界面的法线方向张开，根据叠加原理、原问题的弹性场可分成三部分；二部分由真实均匀本征应变所引起，另一部分由等效的非均匀本征应变所引起，后一部分则由虚拟的Ｓｏｍｉｇｌｉａｎａ位错场所产生。本文求得了等效非均匀本征应变和虚拟位错场的Ｂｕｒｇｅｒ矢量的解析表达式，进而确定的问题的弹性场。     

The principle of equivalent eigenstrain for inhomogeneous inclusion problems

In this paper, based on the principle of virtual work, we formulate the equivalent eigenstrain approach for inhomogeneous inclusions. It allows calculating the elastic deformation of an arbitrarily connected and shaped inhomogeneous inclusion, by replacing it with an equivalent homogeneous inclusion problem, whose eigenstrain distribution is determined by an integral equation. The equivalent homogeneous inclusion problem has an explicit solution in terms of a definite integral. The approach allows solving the problems about inclusions of arbitrary shape, multiple inclusion problems, and lends itself to residual stress analysis in non-uniform, heterogeneous media. The fundamental formulation introduced here will find application in the mechanics of composites, inclusions, phase transformation analysis, plasticity, fracture mechanics, etc.     

Boundary Effects in the Eigenstrain Method

We present a comprehensive study of the effects of internal boundaries on the accuracy of residual stress values obtained from the eigenstrain method. In the experimental part of this effort, a composite specimen, consisting of an aluminum cylinder sandwiched between steel cylinders of the same diameter, was uniformly heated under axial displacement constraint. During the experiment, the sample temperature and the reaction stresses in the load frame in response to changes in sample temperature were monitored. In addition, the local (elastic) lattice strain distribution within the specimen was measured using neutron diffraction. The eigenstrain method, utilizing finite element modeling, was then used to predict the stress field existing within the sample in response to the constraint imposed by the load frame against axial thermal expansion. Our comparison of the computed and measured stress distributions showed that, while the eigenstrain method predicted acceptable stress values away from the cylinder interfaces, its predictions did not match experimentally measured values near them. These observations indicate that the eigenstrain method is not valid for sample geometries with this type of internal boundaries.     

Development of plasticity and damage in vibrating structural elements performing guided rigid-body motions

Summary  A numerical algorithm for studying the development of plastic and damaged zones in a vibrating structural element with a large, guided rigid-body motion is presented. Beam-type elements vibrating in the small-strain regime are considered. A machine element performing rotatory motions, similar to an element of a slider-crank mechanism, is treated as a benchmark problem. Microstructural changes in the deforming material are described by the mesolevel variables of plastic strain and damage, which are consistently included into a macroscopic analysis of the overall beam motion. The method is based on an eigenstrain formulation, considering plastic strain and damage to contribute to an eigenstrain loading of a linear elastic background structure. Rigid-body coordinates are incorporated into this beam-type structural formulation, and an implicit numerical scheme is presented for iterative computation of the eigenstrains from the mesolevel constitutive behavior. Owing to the eigenstrain formulation, any of the existing constitutive models with internal variables could in principle be implemented. Linear elastic/perfectly plastic behavior is exemplarily treated in a numerical study, where plastic strain is connected to the Kachanov damage parameter by a simple damage law. Inelastic effects like plastic shakedown and damage-induced low-cycle rupture are shown to occur in the examplary problems. Received 1 September 1999; accepted for publication 9 March 2000     

压电半平面梯形夹杂的梯度特征应变问题

针对边界处自由和绝缘以及固定和绝缘两种不同的条件,分别计算分析了均匀和梯度特征应变下梯形夹杂内部和外部诱发产生的弹性场和电场,并且讨论了梯形夹杂角点处的奇异性．最后,计算了平均梯度特征应变为零时梯形夹杂内部产生的平均弹性场和电场．所得结果揭示了基体不同边界条件对诱导场的影响．     

Two-dimensional polynomial eigenstrain formulation of boundary integral equation with numerical verification

The low-order polynomial-distributed eigenstrain formulation of the boundary integral equation (BIE) and the corresponding definition of the Eshelby tensors are proposed for the elliptical inhomogeneities in two-dimensional elastic media. Taking the results of the traditional subdomain boundary element method (BEM) as the control, the effectiveness of the present algorithm is verified for the elastic media with a single elliptical inhomogeneity. With the present computational model and algorithm, significant improvements are achieved in terms of the efficiency as compared with the traditional BEM and the accuracy as compared with the constant eigenstrain formulation of the BIE.     

滑动界面的球形夹杂问题

滑动界面对多相介质力学性能的影响日益受到重视．但已有的解析结果往往假定界面是自由滑动的．即假设界面上的剪应力为零，这与大多数实际情况并不相符．本文假定界面上剪应力不为零并满足线弹簧型界面条件，在这一前提下，首次获得了球形夹杂本征应变问题的解析解．     

Inclusion of an arbitrary polygon with graded eigenstrain in an anisotropic piezoelectric half plane

This paper presents an exact closed-form solution for the Eshelby problem of a polygonal inclusion with graded eigenstrains in an anisotropic piezoelectric half plane with traction-free on its surface. Using the line-source Green’s function, the line integral is carried out analytically for the linear eigenstrain case, with the final expression involving only elementary functions. The solutions are applied to the semiconductor quantum wire (QWR) of square, triangular, and rectangular shapes, with results clearly illustrating various influencing factors on the induced fields. The exact closed-form solution should be useful to the analysis of nanoscale QWR structures where large strain and electric fields could be induced by the non-uniform misfit strain.     

Oscillations of a Circular Cylinder in a Linearly Stratified Fluid

The plane problem of the motions of a three-layer fluid initiated by the oscillations of a circular cylinder is solved in the linear formulation in the Boussinesq approximation. The cylinder is completely immersed in the linearly stratified middle layer, and the upper and lower layers are homogeneous and bounded by rigid horizontal walls. The fluid is assumed to be ideal and incompressible. The added mass and damping coefficients are calculated as functions of the oscillation frequency of the cylinder and the layer thicknesses.     

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.     

Analytical solution for a strained reinforcement layer bonded to lip-shaped crack under remote mode III uniform load and concentrated load

In this paper, the analytical solution of stress field for a strained reinforcement layer bonded to a lip-shaped crack under a remote mode III uniform load and a concentrated load is obtained explicitly in the series form by using the technical of conformal mapping and the method of analytic continuation. The effects of material combinations, bond of interface and geometric configurations on interfaciai stresses generated by eigenstrain, remote load and concentrated load are studied. The results show that the stress concentration and interfaciai stresses can be reduced by rational material combinations and geometric configurations designs for different load forms.     

Non-uniform eigenstrain induced stress field in an elliptic inhomogeneity embedded in orthotropic media with complex roots

This paper deals with an elastic orthotropic inhomogeneity problem due to non-uniform eigenstrains. The specific form of the distribution of eigenstrains is assumed to be a linear function in Cartesian coordinates of the points of the inhomogeneity. Based on the polynomial conservation theorem, the induced stress field inside the inhomogeneity which is also linear, is determined by the evaluation of 10 unknown real coefficients. These coefficients are derived analytically based on the principle of minimum potential energy of the elastic inhomogeneity/matrix system together with the complex function method and conformal transformation. The resulting stress field in the inhomogeneity is verified using the continuity conditions for the normal and shear stresses on the boundary. In addition, the present analytic solution can be reduced to known results for the case of uniform eigenstrain.     

Temperature distributions in a periodically stratified layer with slant lamination

The paper deals with the two-dimensional stationary temperature distribution in a composite layer. The nonhomogenous body is assumed to be composed of periodically repeated two-layered laminae. The layering is inclined with an arbitrary angle to the boundary planes. The lower and upper boundary planes are assumed to be kept at given temperatures. The considered problem is solved within the framework of the homogenized model with microlocal parameters, where the continuity thermal conditions on interfaces are satisfied.     

Transient heat conduction in two-dimensional functionally graded hollow cylinder with finite length

In this paper a thick hollow cylinder with finite length made of two dimensional functionally graded material (2D-FGM) subjected to transient thermal boundary conditions is considered. The volume fraction distribution of materials, geometry and thermal boundary conditions are assumed to be axisymmetric but not uniform along the axial direction. The finite element method with graded material properties within each element is used to model the structure and the Crank–Nicolson finite difference method is implemented to solve time dependent equations of the heat transfer problem. Two-dimensional heat conduction in the cylinder is considered and variation of temperature with time as well as temperature distribution through the cylinder are investigated. Effects of variation of material distribution in two radial and axial directions on the temperature distribution and time response are studied. The achieved results show that using two-dimensional FGM leads to a more flexible design so that transient temperature, maximum amplitude and uniformity of temperature distributions can be modified to achieve required specifications by selecting a suitable material distribution profile in two directions.     

Determination of the thermal stress wave propagation in orthotropic hollow cylinder based on classical theory of thermoelasticity

The thermoelasticity problem in a thick-walled orthotropic hollow cylinder is solved analytically using finite Hankel transform and Laplace transform. Time-dependent thermal and mechanical boundary conditions are applied on the inner and the outer surfaces of the cylinder. For solving the energy equation, the temperature itself is considered as boundary condition to be applied on both the inner and the outer surfaces of the orthotropic cylinder. Two different cases are assumed for solving the equation of motion: traction–traction problem (tractions are prescribed on both the inner and the outer surfaces) and traction–displacement (traction is prescribed on the inner surface and displacement is prescribed on the outer surface of the hollow orthotropic cylinder). Due to considering uncoupled theory, after obtaining temperature distribution, the dynamical structural problem is solved and closed-form relations are derived for radial displacement, radial and hoop stress. As a case study, exponentially decaying temperature with respect to time is prescribed on the inner surface of the cylinder and the temperature of the outer surface is considered to be zero. Owing to solving dynamical problem, the stress wave propagation and its reflections were observed after plotting the results in both cases.     

Axisymmetric micromechanics of elasticity-perfectly plastic fibrous composites under uniaxial tension loading

The uniaxial response of a continuous fiber elastic-perfectly plastic composite is modelled herein as a two-element composite cylinder. An axisymmetric analytical micromechanics solution is obtained for the rate-independent elastic-plastic response of the two-element composite cylinder subjected to tensile loading in the fiber direction for the case wherein the core fiber is assumed to be a transversely isotropic elastic-plastic material obeying Tsai-Hill's yield criterion, with yielding simulating fiber failure. The matrix is assumed to be an isotropic elastic-plastic material obeying Tresca's yield criterion. It is found that there are three different circumstances that depend on the fiber and matrix properties: (1) fiber yield, followed by matrix yielding; (2) complete matrix yield, followed by fiber yielding; and (3) partial matrix yield, followed by fiber yielding, followed by complete matrix yield. The order in which these phenomena occur is shown to have a pronounced effect on the predicted uniaxial effective composite response.