Local-buckling-induced elastic interaction between inclusions in a free-standing film

The local-buckling-induced elastic interaction between two circular inclusions in a free-standing film is reported using numerical simulation. The simulation relies on a continuum model based on the modified Föppl-von Kármán plate theory for a film with arbitrarily distributed eigenstrain and eigencurvature. It is shown that due to the overlapping of the nonlinear local buckling the elastic interaction between the two inclusions with the same eigencurvature is repulsive, while the interaction between them with the opposite eigencurvature is attractive. The interaction strength in both cases decays with their mutual distance. In addition, the inclusion with positive/negative eigenstrain above critical values can trigger an axisymmetric/non-axisymmetric buckling, respectively, and the buckling induced elastic interaction between the two inclusions with eigenstrain shows a nonmonotonic behavior.     

Prediction of elastic properties of heterogeneous materials with complex microstructures

The phase-field microelasticity (PFM) is adapted into a homogenization process to predict all the effective elastic constants of three-dimensional heterogeneous materials with complex microstructures. Comparison between the PFM approach and the Hashin-Shtrikman variational approach is also given. Using 3D images of two-phase heterogeneous media with regular and irregular microstructures, results indicate that the PFM approach can accurately take into account the effects of both elastic anisotropy and inhomogeneity of materials with arbitrary microstructure geometry, such as complex porous media with suspended inclusions.     

平面弹性中双圆柱夹杂问题的格林函数

通过复势法并结合共形映射、解析开拓、奇点分析、Cauchy积分公式、圆环域上的Laurent级数展开等技术的联合运用,给出了含有双圆柱异相夹杂的无限大弹性基体这种三相复合系统在受到面内集中力和刃型位错作用时的级数形式的精确解答。该奇点可以作用于基体上或两夹杂内部,当一刃型位错作用于基体上时,也推导了该位错与双圆柱夹杂的交互作用能以及作用于位错上的广义。所给出的算例验证了该弹性解答的正确性并直观显示出奇点与双圆柱夹杂的相互作用。     

Interaction between a dislocation and a core–shell nanowire with interface effects

The problem of a screw dislocation interacting with a core–shell nanowire (coated nanowire) containing interface effects (interface stresses) is first investigated. The interaction energy and the interaction force are calculated. The interaction force and the equilibrium position of the dislocation are examined for variable parameters (interface stress and material mismatch). The influence of the core–shell nanowire and the interface stresses on the interaction between two screw dislocations is also considered. The results show that the impact of the interface stresses on the motion and the equilibrium position of the dislocation near the core–shell nanowire is very significant when the radius of the nanowire is reduced to nanometer scale.     

On stress analysis for a penny-shaped crack interacting with inclusions and voids

An analytical approach to calculate the stress of an arbitrary located penny-shaped crack interacting with inclusions and voids is presented. First, the interaction between a penny-shaped crack and two spherical inclusions is analyzed by considering the three-dimensional problem of an infinite solid, composed of an elastic matrix, a penny-shaped crack and two spherical inclusions, under tension. Based on Eshelby’s equivalent inclusion method, superposition theory of elasticity and an approximation according to the Saint–Venant principle, the interaction between the crack and the inclusions is systematically analyzed. The stress intensity factor for the crack is evaluated to investigate the effect of the existence of inclusions and the crack–inclusions interaction on the crack propagation. To validate the current framework, the present predictions are compared with a noninteracting solution, an interacting solution for one spherical inclusion, and other theoretical approximations. Finally, the proposed analytical approach is extended to study the interaction of a crack with two voids and the interaction of a crack with an inclusion and a void.     

THE INTERACTION PROBLEM BETWEEN THE ELASTIC LINE INCLUSIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＯｆａｌｌｔｈｅｆｉｂｅｒ＿ｒｅｉｎｆｏｒｃｅｄｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｓ,ｔｈｅｓｈｏｒｔ＿ｆｉｂｅｒｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｎｏｔｏｎｌｙｓｔｒｅｎｇｔｈｅｎｓｔｈｅｍａｔｒｉｘｂｕｔａｖｏｉｄｓｄｅｆｅｃｔｉｏｎｓｏｆｔｈｅｌｏｎｇ＿ｆｉｂｅｒｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌａｓｗｅｌｌ.Ｔｈｅｍｉｃｒｏ＿ｍｅｃｈａｎｉｃｓａｂｏｕｔｉｔｓｕｃｈａｓｆｒａｃｔｕｒｅ ,ｆａｔｉｇｕｅａｎｄｄａｍａｇｅｉｓｖｅｒｙｃｏｍｐｌｅｘ .Ｉｎｔｈｅｐｒｅｖｉ…     

双周期圆柱形夹杂纵向剪切问题的精确解

研究无限介质中矩形排列双周期圆柱形夹杂的纵向剪切问题．利用Eshelby等效夹杂理论并结合双周期与双准周期解析函数工具，为这类考虑夹杂相互影响的问题提供了一个严格又实用的分析方法，求得了问题的全场级数解．作为退化情形得到单夹杂问题的经典解答，双周期孔洞、双周期刚性夹杂及单行(列)周期弹性夹杂等问题也可作为特殊情况被解决．数值结果揭示了这类非均匀材料力学性质随微结构参数变化的规律．     

Intensity of singular stress at the fiber end in a hexagonal array of fibers

To evaluate the mechanical strength of fiber-reinforced composites it is necessary to consider singular stresses at the end of fibers because they cause crack initiation, propagation, and final failure. The singular stress field is controlled by generalized stress intensity factor (GSIF) defined at the fiber end. In this study, periodic and zigzag arrays of cylindrical inclusions under longitudinal tension are considered in comparison with the results for a single fiber. The unit cell region is approximated as an axi-symmetric cell; then, the body force method is applied, which requires the stress and displacement fields due to ring forces in infinite bodies having the same elastic constants as those of the matrix and inclusions. The given problem is solved on the superposition of two auxiliary problems under different boundary conditions. To obtain the GSIF accurately, the unknown body force densities are expressed as piecewise smooth functions using fundamental densities and power series. Here, the fundamental densities are chosen to represent the symmetric stress singularity, and the skew-symmetric stress singularity. The GSIFs are systematically calculated with varying the elastic modulus ratio and spacing of fibers. The effects of volume fraction and spacing of fibers are discussed in fiber-reinforced plastics.     

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

Bauschinger and size effects in thin-film plasticity due to defect-energy of geometrical necessary dislocations

The Bauschinger and size effects in the thinfilm plasticity theory arising from the defect-energy of geometrically necessary dislocations (GNDs) are analytically investigated in this paper. Firstly, this defect-energy is deduced based on the elastic interactions of coupling dislocations (or pile-ups) moving on the closed neighboring slip plane. This energy is a quadratic function of the GNDs density, and includes an elastic interaction coefficient and an energetic length scale L. By incorporating it into the work- conjugate strain gradient plasticity theory of Gurtin, an energetic stress associated with this defect energy is obtained, which just plays the role of back stress in the kinematic hardening model. Then this back-stress hardening model is used to investigate the Bauschinger and size effects in the tension problem of single crystal Al films with passivation layers. The tension stress in the film shows a reverse dependence on the film thickness h. By comparing it with discrete-dislocation simulation results, the length scale L is determined, which is just several slip plane spacing, and accords well with our physical interpretation for the defect- energy. The Bauschinger effect after unloading is analyzed by combining this back-stress hardening model with a friction model. The effects of film thickness and pre-strain on the reversed plastic strain after unloading are quantified and qualitatively compared with experiment results.     

Bounds on the overall properties of composites with ellipsoidal inclusions

A new model is put forward to bound the effective elastic moduli of composites with ellipsoidal inclusions. In the present paper, transition layer for each ellipsoidal inclusion is introduced to make the trial displacement field for the upper bound and the trial stress field for the lower bound satisfy the continuous interface conditions which are absolutely necessary for the application of variational principles. According to the principles of minimum potential energy and minimum complementary energy, the upper and lower bounds on the effective elastic moduli of composites with ellipsoidal inclusions are rigorously derived. The effects of the distribution and geometric parameters of ellipsoidal inclusions on the bounds of the effective elastic moduli are analyzed in details. The present upper and lower bounds are still finite when the bulk and shear moduli of ellipsoidal inclusions tend to infinity and zero, respectively. It should be mentioned that the present method is simple and needs not calculate the complex integrals of multi-point correlation functions. Meanwhile, the present paper provides an entirely different way to bound the effective elastic moduli of composites with ellipsoidal inclusions, which can be developed to obtain a series of bounds by taking different trial displacement and stress fields.     

Elastic wave scattering and dynamic stress in composite with fiber

IntroductionTheproblemofscatteringofelasticwavesinsolidstructureshasnotonlytheoreticalsignificancebutalsowideoutlookinengineering .Thisproblemhasabsorbedmanypeople’sattentionsinmanyfieldssuchasaeronautics,compositemechanics,civilengineeringandearthquakeengineering[1～ 3 ] .Thepropagatingvelocityanddirectionofelasticwaveareinvariableasitpropagatesinaninfiniteuniformmedium .Butscatteringofelasticwavescanoccurifthereexistsincontinuitysuchasinclusion ,crackandcavityinelasticmedium .Fiber_reinfor…     

多圆形刚性夹杂的反平面问题

构造任意分布且相互影响的多个圆形刚性夹杂模型的复应力函数，采用复变函数方法，达到满足各个夹杂的边界条件，利用坐标变换和围线积分将求解方程组化为线性代数方程组，推导出了圆形刚性夹杂任意分布的界面应力解析表达式，算例对多夹杂与单夹杂两种模型的界面应力最大值进行了对比，同时还给出了界面应力最大值随夹杂间距的变化规律，求出了刚性夹杂的合理间距。本文发展的分析方法为研究夹杂材料的细观机理探索了一条有效的分析途径。     

Analysis of a transversely isotropic rod containing a single cylindrical inclusion with axisymmetric eigenstrains

This paper studies a transversely isotropic rod containing a single cylindrical inclusion with axisymmetric eigenstrains. The analytical elastic solution is obtained for the displacements, stresses and elastic strain energy of the rod. The effects of microstructural parameters and its evolution on the elastic stress and strain fields as well as the strain energy of the rod are quantitatively demonstrated through examples.     

Stress intensity factors for ring-shaped crack surrounded by spherical inclusions

Interaction of a ring-shaped crack with inhomogeneities such as inclusions is analyzed for the resulting three-dimensional stress field. Considered for the composite solid with a given volume fraction of inclusions are the two cases of (a) spherical voids and (b) spherical inclusions with elastic moduli different from the matrix. A ring-shaped crack is initiated at the equator of one of the voids or inclusions. A three-phase model is used to examine the interaction between the crack and surrounding inhomogeneities. Finite element method is then applied to calculate the stress intensity factor for different configurations. The effects of volume fraction of inhomogeneities, relative size of crack to inclusions, and material constants on crack behavior are discussed.     

Interaction between dislocation and two circular cylindrical inclusions

In this paper, the elastic field of the infinite homogeneous medium with two circular cylindrical inclusions under the action of a screw dislocation is investigated and the corresponding analytical solution is obtained. Here, the conformal mapping and the theorem of analytical continuation are used. From the results obtained, it can be seen that the elastic field depends on the shear moduli of individual phases, the geometric parameters of the system, and the position and relative slip of the screw dislocation. In addition, the corresponding specific cases are also considered in this paper when two circular cylindrical inclusions are tangent to each other and they are holes and/or rigid inclusions. Finally, numerical results are illustrated to show the interaction between the screw dislocation and two circular cylindrical inclusions.     

Fracture saturation and critical thickness in layered materials

Opening-mode fractures in layered materials are commonly found in a layer with uniform spacing that is nearly proportional to the thickness of the fractured layer. However, when fracture spacing reduces to a certain value, fracture density is saturated and no new fracture forms. If a loading condition is fixed, there exists a critical thickness of the layer, below which no fracture forms. This paper presents a three-layer model, containing a weak layer between two stronger layers, to interpret the fracture saturation and critical thickness of layered materials. Using elastic governing equations and a weak form stress boundary condition, a closed-form solution of elastic fields in the weak layer is derived and the energy release rate for opening-mode fracture is obtained. Interestingly, the normal stress between such fractures undergoes a transition from tensile to compressive with increasing applied tensile loading, which causes fracture saturation. Explicit expressions of critical fracture-spacing-to-layer-thickness ratio and critical thickness are derived for fracture saturation and fracture free conditions, respectively. Comparison with the existing numerical simulation results demonstrates the capability of this model. This explicit, analytical solution is useful to structural design and geosciences.     

纳米线涂层中失配位错偶极子分析

研究了纳米线环形晶体薄膜涂层中失配位错偶极子与纳米线的干涉效应,并考虑纳米尺度应力效应及纳米线晶格失配的影响。运用弹性复势方法,分别获得了涂层和纳米线区域复势函数的精确解析解。利用求得的应力场和Peach-Koehler公式,得到了作用在螺型位错偶极子上像力和失配应力的精确表达式。算例结果表明,涂层纳米线界面应力和失配应力对涂层中失配位错偶极子的作用影响很大,由于界面应力和失配应力的存在,可以改变涂层内位错偶极子与纳米线干涉的引斥规律。与宏观尺度下的线型材料相比,纳米线材料由于界面效应的影响,位错偶极子在涂层中平衡位置的数量增加,涂层中更容易产生失配位错偶极子。考虑纳米线晶格失配的影响后,位错偶极子在涂层中平衡位置的数量增加。由此可知,减小纳米线晶格失配的影响,可以控制失配位错偶极子的产生。     

Interaction of Anisotropic Ellipsoidal Inclusions in an Isotropic Medium under Mechanical Loads and Uniform Heating

Eshelby's equivalent-inclusion method is extended to a finite number of arbitrarily oriented anisotropic ellipsoidal inclusions in an elastic isotropic matrix under polynomial mechanical loading and heating. The interaction of two identical and two different triaxial ellipsoidal inclusions in an elastic medium is studied as numerical examples. In special cases, the results are compared with those obtained by other authors     

Mean Green operators and Eshelby tensors for hemispherical inclusions and hemisphere interactions in spheres. Application to bi-material spherical inclusions in isotropic spaces

This paper mainly presents an exact expression for the mean shape function of a hemispherical inclusion, from which are obtained analytical forms for the mean Green operator (GO) and Eshelby tensor of this hemi-sphere as well as for the related mean pair interaction Green operator (IGO) between the two hemi-spheres of a sphere, in media with isotropic (elastic or dielectric) properties. We secondly address the problem of bi-material inclusions, in the sense of a two-phase compact set of two or a few elementary domains, a particular inclusion pattern case for which we give an estimate of the mean stress and strain in each phase accounting for interactions. This estimate results from knowing the mean GO (or Eshelby tensor) for each pattern element plus the mean IGO between element pairs, what is rarely fulfilled analytically. The here solved case for bi-material spherical inclusions made of two different hemispherical elements adds to the recently made available solution for bi-material cylindrical inclusions made of piled coaxial finite cylinders. The obtained mean stress estimates are exemplified able to satisfactorily match with FEM calculations up to highly contrasted bi-material inclusions. Other types of bi-material spherical inclusions are mentioned for which the mean GOs for the sub-domains and their pair IGO can be obtained without calculation, owing to particular symmetries of the phase arrangement. Mean GOs and IGOs are also useful in certain homogenization frameworks yielding overall property estimates for inclusion-reinforced matrices. Further discussions and specific applications will be presented in forthcoming papers.