From coherent to incoherent mismatched interfaces: A generalized continuum formulation of surface stresses

The equilibrium of coherent and incoherent mismatched interfaces is reformulated in the context of continuum mechanics based on the Gibbs dividing surface concept. Two surface stresses are introduced: a coherent surface stress and an incoherent surface stress, as well as a transverse excess strain. The coherent surface stress and the transverse excess strain represent the thermodynamic driving forces of stretching the interface while the incoherent surface stress represents the driving force of stretching one crystal while holding the other fixed and thereby altering the structure of the interface. These three quantities fully characterize the elastic behavior of coherent and incoherent interfaces as a function of the in-plane strain, the transverse stress and the mismatch strain. The isotropic case is developed in detail and particular attention is paid to the case of interfacial thermo-elasticity. This exercise provides an insight on the physical significance of the interfacial elastic constants introduced in the formulation and illustrates the obvious coupling between the interface structure and its associated thermodynamics quantities. Finally, an example based on atomistic simulations of Cu/Cu₂O interfaces is given to demonstrate the relevance of the generalized interfacial formulation and to emphasize the dependence of the interfacial thermodynamic quantities on the incoherency strain with an actual material system.     

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

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

双周期分布圆形弹性夹杂平面热弹性问题

研究了含双周期分布圆形弹性夹杂的无限弹性平面在均匀拉伸和均匀温变下的弹性响应问题．运用Isida的区域单元法和复势函数的级数展开技术,将问题转化为线性方程组的求解．数值结果表明：相邻夹杂问距过大或过小都不利于减小界面应力,当相邻夹杂中心间距与夹杂半径之比为2．2～2．8时,界面剪切应力与环向应力的极大值最小;比值为2．5～3．5时,界面最大径向应力值最小;并且该比值范围不随两相材料弹性模量之比和热膨胀系数之比而变化．     

An Asymptotic Study¶of the Debonding of Thin Films

We examine the asymptotic behavior of a bilayer thin film using the notion of Γ-convergence. We allow for debonding at the interface, but penalize it using an interfacial energy; thus the functional we consider consists of the elastic energy of the two layers and the interfacial energy with penalized debonding. We show that the asymptotic theory or Γ-limit depends on the particular form of the interfacial energy, and derive detailed results for both the cohesive and the brittle interface.     

含界面相Voronoi单元的等效弹性常数的计算

采用含界面相Voronoi单元有限元法,根据广义胡克定律,计算了在给定边界条件下,颗粒增强复合材料的等效弹性常数。建立了含多个随机分布的椭圆形夹杂及界面相的VCFEM模型,分析了夹杂体分比,界面相厚度和界面相弹性模量等因素对颗粒增强复合材料等效弹性常数的影响,并利用普通有限元方法对比验证。结果表明,当界面相弹性模量小于基体与夹杂时,材料的等效弹性模量会随着界面相厚度的增大而减小,随着夹杂体分比的增大而减小,并且界面过薄时,材料的等效弹性模量会随着夹杂体分比的增大而增大;当界面相弹性模量大于基体或夹杂时,材料的等效弹性模量会随着夹杂体分比和界面相厚度的增大而增大。而界面相的厚度和弹性模量对材料的等效泊松比的影响较小,材料的等效泊松比主要受夹杂体分比的影响,与其呈反比关系。     

涂层式裂纹监测系统中基体裂纹穿越行为研究

建立了智能涂层的两相模型与三相模型,基于能量准则分别用这两种模型研究了基体裂纹达到涂层界面后的穿越/偏转行为. 用有限元法分析了相对裂纹扩展长度、弹性错配参数及界面层厚度对偏转裂纹与穿越裂纹能量释放率之比的影响,结果表明当基体裂纹到达驱动层与基体界面时,能量释放率之比不仅与基体和驱动层之间的弹性错配相关,而且当驱动层较薄时对驱动层与传感层之间的弹性错配亦有较强的依赖性. 此外,随着驱动层厚度的增加,能量释放率之比对驱动层与传感层之间的弹性错配的依赖性逐渐降低. 通过与实验结果相比,建立的模型能够较好的解释基体裂纹在界面的扩展行为,可用于智能涂层裂纹传感器的优化设计.      

Indentation of a spherical cavity in an elastic body by a rigid spherical inclusion: influence of non-classical interface conditions

This paper examines the indentation of an elastic body by a rigid spherical inclusion. In contrast to conventional treatments where the contact between a rigid inclusion and the elastic medium is regarded as being perfectly bonded, we examine the influence of non-classical interface conditions including frictionless bilateral contact, separation and Coulomb friction on the load–displacement behaviour of the spherical rigid inclusion. Both analytical methods and boundary element techniques are used to examine the inclusion/elastic medium interaction problems. This paper also provides a comprehensive review of non-classical interface conditions between inclusions and the surrounding elastic media.     

Shape control of thin rigid inclusions and cracks in elastic bodies

The paper is concerned with a control of thin rigid inclusion and crack shapes in elastic bodies. It is assumed that rigid inclusions are delaminated; thus, cracks are located on the boundary of inclusions as well as outside of inclusions. We provide the problem formulations and analyze the shape sensitivity with respect to geometrical perturbations in the frame of free boundary models. Inequality type boundary conditions are considered at the crack faces to guarantee a mutual non-penetration between crack faces. Inclusion and crack shapes are considered as control functions. The cost functional, which is based on the Griffith rupture criterion, characterizes the energy release rate and provides the shape sensitivity with respect to a change of the geometry of the structure. We prove an existence of optimal shapes in the problems considered.     

Structural interfaces in linear elasticity. Part II: Effective properties and neutrality

The model of structural interfaces developed in Part I of this paper allows us to analytically attack and solve different problems of stress concentration and composites. In particular, (i) new formulae are given for effective properties of composite materials containing dilute suspensions of (randomly oriented) reinforced elliptical voids or inclusions; (ii) a new definition is proposed for inclusion neutrality (to account for the fact that the matrix is always ‘overstressed’, and thus non-neutral in a classical sense, at the contacts with the interfacial structure), which is shown to provide interesting stress optimality conditions. More generally, it is shown that the incorporation of an interfacial structure at the contact between two elastic solids exhibits properties that cannot be obtained using the more conventional approach of the zero-thickness, linear interface. For instance: contrary to the zero-thickness interface, both bulk and shear effective moduli can be optimized for a structural interface; effective properties higher that those possible with a perfect interface can be attained with a structural interface; and neutrality holds with a structural interface for a substantially broader range of parameters than for a zero-thickness interface.     

Interfacial excess energy, excess stress and excess strain in elastic solids: Planar interfaces

In this paper, interfacial excess energy and interfacial excess stress for coherent interfaces in an elastic solid are reformulated within the framework of continuum mechanics. It is shown that the well-known Shuttleworth relationship between the interfacial excess energy and interfacial excess stress is valid only when the interface is free of transverse stresses. To account for the transverse stress, a new relationship is derived between the interfacial excess energy and interfacial excess stress. Dually, the concept of transverse interfacial excess strain is also introduced, and the complementary Shuttleworth equation is derived that relates the interfacial excess energy to the newly introduced transverse interfacial excess strain. This new formulation of interfacial excess stress and excess strain naturally leads to the definition of an in-plane interfacial stiffness tensor, a transverse interfacial compliance tensor and a coupling tensor that accounts for the Poisson's effect of the interface. These tensors fully describe the elastic behavior of a coherent interface upon deformation.     

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

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

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.     

A continuum description of partially coherent interfaces

A continuum model of a two-phase crystal-crystal system is constructed in which the structure of the interface between the phases is determined by energy minimization, rather than by being specified a priori. The interfacial structure is parameterized by a variable? corresponding to the jump in the surface deformation gradient (or strain) at the interface, so that coherence is defined locally by the condition? = 0. The energy of the system is taken to be the sum of the bulk and interfacial energies, where the interfacial energy densityf ^xs depends on?. In order to explore how the equilibrium interfacial structure depends on the functionf ^xs (?), a model system consisting of an elastic film on a rigid substrate is studied, and the interfacial energy density is taken to be nonconvex with a sharp minimum associated with coherence. In this case, it can be shown that the energy of the system is driven to its infimum by separating the interface into coherent and incoherent regions, which may be viewed as a continuum analog to a partially coherent interface. Further, this solution only appears above a certain critical thickness of the film, in agreement with misfit dislocation models of partially coherent interfaces.     

An anisotropic micromechanics model for predicting the rafting direction in Ni-based single crystal superalloys

An anisotropic micromechanics model based on the equivalent inclusion method is developed to investigate the rafting direction of Ni-based single crystal superalloys. The micromechanical model considers actual cubic structure and orthogonal anisotropy properties. The von Mises stress, elastic strain energy density, and hydrostatic pressure in different inclusions of micromechanical model are calculated when applying a tensile or compressive loading along the [001] direction. The calculated results can successfully predict the rafting direction for alloys exhibiting a positive or a negative mismatch, which are in agreement with pervious experimental and theoretical studies. Moreover, the elastic constant differences and mismatch degree of the matrix and precipitate phases and their influences on the rafting direction are carefully discussed.     

Timoshenko inclusions in elastic bodies crossing an external boundary at zero angle

The paper concerns an analysis of equilibrium problems for 2D elastic bodies with a thin Timoshenko inclusion crossing an external boundary at zero angle. The inclusion is assumed to be delaminated, thus forming a crack between the inclusion and the body. We consider elastic inclusions as well as rigid inclusions. To prevent a mutual penetration between the crack faces, inequality type boundary conditions are imposed at the crack faces.Theorems of existence and uniqueness are established. Passages to limits are investigated as a rigidity parameter of the elastic inclusion going to infinity.     

THE INTERACTION PROBLEM BETWEEN THE ELASTIC LINE INCLUSIONS

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

The low frequency scatter of SH waves by a rough interface in an elastic plate

The study of the reflection and transmission of low frequency SH waves incident upon a rough interface in an elastic plate is undertaken by employing a theory of acoustic wave scattering from rough surfaces originally due to Biot and subsequently generalised to the case of elastic media. In this theory the interface is replaced by a distribution of voids/asperities whose individual size is small compared to the excitation wavelength. We plot the absolute values of the reflection and transmission coefficients versus frequency when a single symmetric SH plate mode is used as the input excitation. The different types of inclusions are used to simulate the rough surface are the hollow, fluid filled and aluminum spheres. Lastly, the loss of energy due to scattering is also estimated for the different inclusion distributions considered.     

Edge dislocation interacting with an interfacial crack along a circular inhomogeneity

The elastic interaction of an edge dislocation, which is located either outside or inside a circular inhomogeneity, with an interfacial crack is dealt with. Using Riemann–Schwarz’s symmetry principle integrated with the analysis of singularity of the complex potentials, the closed form solutions for the elastic fields in the matrix and inhomogeneity regions are derived explicitly. The image force on the dislocation is then determined by using the Peach–Keohler formula. The influence of the crack geometry and material mismatch on the dislocation force is evaluated and discussed when the dislocation is located in the matrix. It is shown that the interfacial crack has significant effect on the equilibrium position of the edge dislocation near a circular interface. The results also reveal a strong dependency of the dislocation force on the mismatch of the shear moduli and Poisson’s ratios between the matrix and inhomogeneity.     

Asymptotic behavior of the energy functional for a three-dimensional body with a rigid inclusion and a crack

A model of a three-dimensional elastic body containing a rigid inclusion and a crack located on the interface between the inclusion and the body is considered. Natural boundary conditions are imposed on the crack. A derivative of the energy functional with respect to the perturbation parameter is derived for an arbitrary, rather smooth perturbation of the domain, in particular, the Griffith formula is obtained.