Interaction of an Ellipsoidal Inclusion with an Elliptic Crack in an Elastic Medium

An approach is considered to how to allow for the interaction between an ellipsoidal heterogeneity (inclusion) and an elliptic crack in an elastic medium. Using the superposition of perturbed stress states, the boundary conditions are satisfied on the ellipsoidal surface by the method of equivalent inclusion and on the crack surface by the least-squares method. A numerical analysis is carried out. Typical mechanical effects are revealed. In the calculations, the stress state near the ellipsoidal heterogeneity is approximated by a polynomial of the second degree in Cartesian coordinates, whereas the load on the crack surface is simulated by a polynomial of the fourth degree in Cartesian coordinates. In particular cases, the results are in good agreement with the data obtained by other authors     

Mechanical response of multicrystalline thin films in mesoscale field dislocation mechanics

A continuum model of plasticity, Phenomenological Mesoscopic Field Dislocation Mechanics (PMFDM), is used to study the effect of surface passivation, grain orientation, grain boundary constraints, and film thickness on the mechanical response of multicrystalline thin films. The numerical experiments presented in this paper show that a surface passivation layer on thin films introduces thickness dependence of the mechanical response. However, the effect of passivation decreases in films with impenetrable grain boundaries. The orientation of individual grains of the multicrystal also has a significant effect on the mechanical response. Our results are in qualitative agreement with experimental observations. A primary contribution of this work is the implementation of a jump condition that enables the modeling of important limits of grain boundary constraints to plastic flow, independent of ad-hoc constitutive assumptions and interface conditions.     

Thickness effects in polycrystalline thin films: Surface constraint versus interior constraint

The uniaxial tension behavior of polycrystalline thin films, in which all grain boundaries (GBs) are penetrable by dislocations, is investigated by two-dimensional discrete dislocation dynamics (DDD) method with a penetrable dislocation-GB interaction model. In order to study thickness effect on the tensile strength of thin films with and without surface treatment, three types of thin films are comparatively considered, including the thin films without surface treatment, with surface passivation layers (SPLs) of nanometer thickness and with surface grain refinement zones (SGRZs) consisting of nano-sized grains. Our results show that thickness effects and their underlying dislocation mechanisms are quite distinct among different types of thin films. The thicker thin films without surface treatment are stronger than the thinner ones; however, opposite thickness effects are captured in the thin films with SPLs or SGRZs. Moreover, the underlying dislocation mechanisms of the same thickness effects of thin films with SPLs and SGRZs are different. In the thin films with SPLs, the thickness effect is caused by the sharp increase of dislocation density near the film-passivation interface, while it is mainly due to the sharp decrease of dislocation density within the refined surface grains of the thin films with SGRZs. No matter in what type of thin films, thickness effect gradually disappears when the number of grains in the thickness direction is large enough. Our analysis reveals that general mechanism of those thickness effects lies in the competition between the exterior surface-constraint and interior GB-constraint on gliding dislocations.     

Contact analysis in the presence of an ellipsoidal inhomogeneity within a half space

Many materials contain inhomogeneities or inclusions that may greatly affect their mechanical properties. Such inhomogeneities are for example encountered in the case of composite materials or materials containing precipitates. This paper presents an analysis of contact pressure and subsurface stress field for contact problems in the presence of anisotropic elastic inhomogeneities of ellipsoidal shape. Accounting for any orientation and material properties of the inhomogeneities are the major novelties of this work. The semi-analytical method proposed to solve the contact problem is based on Eshelby’s formalism and uses 2D and 3D Fast Fourier Transforms to speed up the computation. The time and memory necessary are greatly reduced in comparison with the classical finite element method. The model can be seen as an enrichment technique where the enrichment fields from the heterogeneous solution are superimposed to the homogeneous problem. The definition of complex geometries made by combination of inclusions can easily be achieved. A parametric analysis on the effect of elastic properties and geometrical features of the inhomogeneity (size, depth and orientation) is proposed. The model allows to obtain the contact pressure distribution – disturbed by the presence of inhomogeneities – as well as subsurface and matrix/inhomogeneity interface stresses. It is shown that the presence of an inclusion below the contact surface affects significantly the contact pressure and subsurfaces stress distributions when located at a depth lower than 0.7 times the contact radius. The anisotropy directions and material data are also key elements that strongly affect the elastic contact solution. In the case of normal contact between a spherical indenter and an elastic half space containing a single inhomogeneity whose center is located straight below the contact center, the normal stress at the inhomogeneity/matrix interface is mostly compressive. Finally when the axes of the ellipsoidal inclusion do not coincide with the contact problem axes, the pressure distribution is not symmetrical.     

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.     

Interface models for thin interfacial layers

There have already been several interface models for the analyses of thin interfacial layers in bonded materials. To distinguish their corresponding advantages or limitations, a comparative study is carried out, and a new constitutive-based interface model is proposed. Through numerical examinations, the limitations of typical models are clarified. It is found that the new interface model is an efficient and accurate model, by which both the traction and the displacement jumps across the modelled interface with the thickness of zero are allowed, and the stresses within the interfacial layer can also be analyzed.     

具有光滑界面的椭球夹杂的轴对称弹性场

本文在旋转椭球坐标系下,利用Papkovich—Neuber位移通解求解了具有光滑界面椭球夹杂由于均匀的特征应变引起的轴对称弹性场,与理想界面不同,在夹杂与基体界面不能经受剪应力而可自由滑动的情况下,解答只能是无穷级数形式,因此文中给出了数值算例。     

Propagation and scattering of ultrasonic waves in polycrystals with arbitrary crystallite and macroscopic texture symmetries

A general ultrasonic attenuation model for a polycrystal with arbitrary macroscopic texture and triclinic ellipsoidal grains is described with proper accounting for the anisotropic Green’s function for the reference medium. The texture and the ellipsoidal grain frames in the model are independent and the wave propagation direction is arbitrary. The attenuation coefficients are obtained in the Born approximation accompanied by the Rayleigh and stochastic asymptotes. The scattering model displays statistical anisotropy due to two independent factors: (1) shape of the oriented grains and (2) preferred crystallographic orientation of the grains leading to macroscopic anisotropy of the homogenized reference medium. The model is applicable to most single phase polycrystalline materials that may occur as a result of thermomechanical manufacturing processes leading to different macrotextures and elongated-shaped grains. It predicts the strength of ultrasonic scattering and its dependence on frequency and propagation direction as a function of grain shape, grain crystallographic symmetry and macroscopic texture parameters and provides the texture-induced dependence of macroscopic ultrasonic velocity on propagation angle. It considers proper wave polarizations due to macroscopic anisotropy and scattering-induced transformations of waves with different polarizations. Competing effects of grain shape and texture on the attenuation are observed. In contrast to the macroscopically isotropic case, where in the stochastic regime the attenuation is highest in the direction of the longest ellipsoidal axis of the grain, the wave attenuation in the elongation direction may be suppressed or amplified by the texture with different effects on the quasilongitudinal and quasitransverse waves. The frequency behavior is also interestingly affected by texture: a hump in the total attenuation coefficient is found for the fast quasitransverse wave which is purely the result of macroscopic anisotropy and the existence of two quasitransverse waves; this hump is not observed in the macroscopically isotropic case. Striking differences of the texture effect on the directional dependences of the attenuation coefficients are found at low versus high frequencies.     

Stress intensity factor for an edge crack in two bonded dissimilar materials under convective cooling

The transient thermal stress crack problem for two bonded dissimilar materials subjected to a convective cooling on the surface containing an edge crack perpendicular to the interface is considered. The problem is solved using the principle of superposition and the uncoupled quasi-static thermoelasticity. The crack problem is formulated by applying the transient thermal stresses obtained from the uncracked medium with opposite sign on the crack surfaces to be the only external loads. Fourier integral transform is used to solve the perturbation problem resulting in a singular integral equation of Cauchy type in which the derivative of the crack surface displacement is the unknown function. The numerical results of the stress intensity factors are calculated for both the edge crack and the crack terminating at the interface using two different composite materials and illustrated as a function of time, crack length, coefficient of heat transfer, and the thickness ratio.     

A microscopic damage model considering the change of void shape and application in the void closing

A microscopic damage model of ellipsoidal body containing ellipsoidal void for nonlinear matrix materials is developed under a particular coordinate. The change of void shape is considered in this model. The viscous restrained equation obtained from the model is affected by stress ∑_(ij), void volume fraction f, material strain rate exponent m as well as the void shape. Gurson's equation is modified from the numerical solution. The modified equation is suitable for the case of nonlinear matrix materials and changeable voids. Lastly, the model is used to analyze the closing process of voids.     

Micromechanics-based elastic damage modeling of particulate composites with weakened interfaces

A micromechanical framework is proposed to predict the effective elastic behavior and weakened interface evolution of particulate composites. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials 14, 269–281] is adopted to model spherical particles having imperfect interfaces in the composites and is incorporated into the micromechanical framework. Based on the Eshelby’s micromechanics, the effective elastic moduli of three-phase particulate composites are derived. A damage model is subsequently considered in accordance with the Weibull’s probabilistic function to characterize the varying probability of evolution of weakened interface between the inclusion and the matrix. The proposed micromechanical elastic damage model is applied to the uniaxial, biaxial and triaxial tensile loadings to predict the various stress–strain responses. Comparisons between the present predictions with other numerical and analytical predictions and available experimental data are conducted to assess the potential of the present framework.     

Stability against localization of softening into ellipsoids and bands: Parameter study

Strain-localization instabilities due to strain-softening which result from distributed damage such as cracking in heterogeneous brittle materials are analysed. Attention is restricted to the stability problem of equilibrium states. This problem is not equivalent to bifurcation of the equilibrium path, which may occur before stability of equilibrium is lost. The continuum is local but is enhanced by the localization limiter used in the crack band model, consisting of a lower bound on the minimum dimension of the strain-localization region, which is regarded as a material property. Presented are derivations of the critical state conditions for localization of initially uniform strain into ellipsoidal domains within an infinite continuum and into a planar band within a layer of finite thickness. These derivations are simpler than the previous Ba ant's derivations of the general stability conditions for these localizations. A numerical parameter study of the critical states is made for a broad range of material properties as well as various initial stress states and relative sizes of the strain-softening region. The material is described by Drucker-Pragcr plasticity with strain-softening that is caused by yield limit degradation. The Hatter the ellipsoidal domain, or the larger the size of the body (layer thickness), the smaller is found to be the strain-softening slope magnitude at which the critical state is reached. A softening Drucker-Prager material is found to be stable against planar-band localizations in infinite continuum for a certain range of softening material parameters.     

关于土—结相互作用界面力学行为的数值模拟

根据土 -结 (桩 )相互作用典型的界面力学行为特点——两者相互作用直至剪切破坏的滑动面一般都是在土体一侧内部发生的 ,滑动层的厚薄与土性及结构物表面粗糙度等有关等 ,提出了把二者所谓相互作用的界面在概念上做广义化的处理 ,即泛化为一有限的区域 ,认为结构与土体相互作用的界面应是一个广义的相互作用影响范围而非绝对地仅只限于两者的相互接触界面 ,两者间相互作用的界面行为可近似地通过这有限薄层范围内土体材料的本构行为来模拟 ,从而形成了本文关于土 -结相互作用界面行为的模拟关键实质上可以归结为土体的非线性本构行为及其数值分析途径的研究等观点 ,仅供同行参考     

考虑粘结面厚度的局部粘结单元法

传统无厚度粘结单元法CFEM (Cohesive finite element method)在模拟脆性材料断裂方面具有很强的优势,但也存在很大问题.一是单元尺寸增大,收敛性变差;二是单元尺寸变小,模型刚度发生折减.为了克服这两个问题,发展了考虑厚度的局部粘结单元法,即在裂纹可能扩展区插入具有一定厚度的粘结面单元.粘结面单元采用拓展虚内键本构(Augmented virtual internal bond)描述.由于考虑了厚度,粘结面交叉处会形成多边形空缺.为了弥补这一空缺,将其看作多边形键元胞,采用离散虚内键模型(Discretized virtual internal bond)对其建模,保证了模型的几何完整性.模拟结果表明,本文方法有效,克服了传统CFEM方法的刚度折减问题,提高了计算稳定性和收敛性.     

Differential geometrical method in elastic composite with imperfect interfaces

I.IntroductionWhethertheinterfacesofcompositematerialsareperfectornotwillaffectitsmacromechanicaloreffectivepropertiesimportantly.Butsofar,almostallofthestudiesontheeffectivepropertiesofcompositematerialsarebasedontheassumptionthattheinterfacesareperfectl"2].Infact,thisisnotappropriateforallinterfaces[31.Thusthestudiesonmechanicalpropertyofcompositematerialswithimperfaceintert'acehavebeenconsideredrecentlyinsomeliteratures.Hashin16]hasextendedtheelasticextremumprinciplesofminimumpotentialandm…     

An elastoplastic multi-level damage model for ductile matrix composites considering evolutionary weakened interface

An elastoplastic multi-level damage model considering evolutionary weakened interface is developed in this work to predict the effective elastoplastic behavior and multi-level damage evolution in particle reinforced ductile matrix composites (PRDMCs). The elastoplastic multi-level damage model is micromechanically derived on the basis of the ensemble-volume averaging procedure and the first-order effects of eigenstrains. The Eshelby’s tensor for an ellipsoidal inclusion with slightly weakened interface [Qu, J., 1993a. Eshelby tensor for an elastic inclusion with slightly weakened interfaces. Journal of Applied Mechanics 60 (4), 1048–1050; Qu, J., 1993b. The effect of slightly weakened interfaces on the overall elastic properties of composite materials. Mechanics of Materials, 14, 269–281] is adopted to model particles having mildly or severely weakened interface, and a multi-level damage model [Lee, H.K., Pyo, S.H., in press. Multi-level modeling of effective elastic behavior and progressive weakened interface in particulate composites. Composites Science and Technology] in accordance with the Weibull’s probabilistic function is employed to describe the sequential, progressive weakened interface in the composites. Numerical examples corresponding to uniaxial, biaxial and triaxial tension loadings are solved to illustrate the potential of the proposed micromechanical framework. A series of parametric analysis are carried out to investigate the influence of model parameters on the progression of weakened interface in the composites. Furthermore, the present prediction is compared with available experimental data in the literature to verify the proposed elastoplastic multi-level damage model.     

Stress state of a piezoelectric elastic medium with an arbitrarily oriented triaxial ellipsoidal inclusion

We determine the electrostressed state of a piezoceramic medium with an arbitrarily oriented triaxial ellipsoidal inclusion under homogeneous mechanical and electric loads. Use is made of Eshelby’s equivalent inclusion method generalized to the case of a piezoelectric medium. Solving the problem for a spheroidal cavity with the axis of revolution aligned with the polarization axis demonstrates the high efficiency of the approach. A numerical analysis is carried out. The stress distribution along the surface of the arbitrarily oriented triaxial ellipsoidal inclusion is studied     

基于构型力内变量的界面损伤问题研究

本文基于材料构型力的基本理论和损伤力学中含有内变量的热力学框架,提出了新的损伤变量定义方式,为研究界面损伤问题提供了一种新思路。首先,基于双相弹性体的能量分析,给出界面材料构型力表达式,通过构型力的离散化方法,实现了其在有限元中的数值计算。其次,定义构型力为界面损伤内变量,进而提出一种新的损伤演化模型,并采用刚度劣化的方法,对该界面损伤模型进行数值实现。最后,通过对复合材料界面损伤问题（有裂纹或无裂纹）进行数值模拟,分析了其界面损伤发展趋势,探讨了此模型的合理性和优越性。基于构型力内变量的界面损伤模型,可为研究复合材料的界面损伤失效问题提供一种普适性的方法。     

Cahn–Hilliard modeling of particles suspended in two‐phase flows

In this paper, we present a model for the dynamics of particles suspended in two‐phase flows by coupling the Cahn–Hilliard theory with the extended finite element method (XFEM). In the Cahn–Hilliard model the interface is considered to have a small but finite thickness, which circumvents explicit tracking of the interface. For the direct numerical simulation of particle‐suspended flows, we incorporate an XFEM, in which the particle domain is decoupled from the fluid domain. To cope with the movement of the particles, a temporary ALE scheme is used for the mapping of field variables at the previous time levels onto the computational mesh at the current time level. By combining the Cahn–Hilliard model with the XFEM, the particle motion at an interface can be simulated on a fixed Eulerian mesh without any need of re‐meshing. The model is general, but to demonstrate and validate the technique, here the dynamics of a single particle at a fluid–fluid interface is studied. First, we apply a small disturbance on a particle resting at an interface between two fluids, and investigate the particle movement towards its equilibrium position. In particular, we are interested in the effect of interfacial thickness, surface tension, particle size and viscosity ratio of two fluids on the particle movement towards its equilibrium position. Finally, we show the movement of a particle passing through multiple layers of fluids. Copyright © 2011 John Wiley & Sons, Ltd.     

Zero-moment ellipsoidal domes

This paper is concerned with the problem of designing ellipsoidal domes so that they are in a membrane state of stress. The investigation is based on a method proposed by Vekua, in which an adjustment of the gravity of the shell is used to ensure a membrane state. The theory is specialized to ellipsoidal shells and some numerical results of the gravity distribution (i.e. thickness of the shell), and the corresponding stresses and linear elastic displacements are presented.However, it should be noted that the solutions obtained in this way are not unique.