Differential nonhomogeneous models for elastic randomly cracked solids

A differential scheme is developed to approximate the elastic behaviour of randomly cracked solids, accounting for possible locally nonhomogeneous distribution of the cracks. Certain inclusion or matrix spaces within the solids are modelled as forbidden regions for the cracks. At small to intermediate values of the crack density and proportion of forbidden regions, the effective elastic moduli of the models do not differ much from each other, but the differences become profound at higher values of those parameters: the effective moduli can be very small and large (toward those of the uncracked solids) depending upon the particular nonhomogeneous arrangements of the cracks.     

Effective elastic moduli of inhomogeneous solids by embedded cell model

An embedded cell model is presented to obtain the effective elastic moduli for three-dimensional two-phase composites which is an exact analytic formula without any simplified approximation and can be expressed in an explicit form. For the different cells such as spherical inclusions and cracks surrounded by sphere and oblate ellipsoidal matrix, the effective elastic moduli are evaluated and the results are compared with those from various micromechanics models. These results show that the present model is direct, simple and efficient to deal with three-dimensional two-phase composites. The project supported by the National Natural Science Foundation of China (No. 19704100) and the National Natural Science Foundation of Chinese Academy of Sciences (No. KJ951-1-201)     

Elastic interaction of partially debonded circular inclusions. II. Application to fibrous composite

A complete analytical solution has been obtained of the elasticity problem for a plane containing periodically distributed, partially debonded circular inclusions, regarded as the representative unit cell model of fibrous composite with interface damage. The displacement solution is written in terms of periodic complex potentials and extends the approach recently developed by Kushch et al. (2010) to the cell type models. By analytical averaging the local strain and stress fields, the exact formulas for the effective transverse elastic moduli have been derived. A series of the test problems have been solved to check an accuracy and numerical efficiency of the method. An effect of interface crack density on the effective elastic moduli of periodic and random structure FRC with interface damage has been evaluated. The developed approach provides a detailed analysis of the progressive debonding phenomenon including the interface cracks cluster formation, overall stiffness reduction and damage-induced anisotropy of the effective elastic moduli of composite.     

含微裂纹材料的损伤理论

本文从含微裂纹材料的变形能出发引出了裂纹的方位张量。在考虑裂纹受压闭合与滑动摩擦的基础上,给出了损伤张量、损伤应变及有效弹性常数。文中给出了损伤机构离散化的方法,并对方位密度给出了演化方程。最后给出一个单向拉压的应力应变关系例子,并揭示了裂纹扩展时的应力突跌现象。     

Electromechanical response of (3–0, 3–1) particulate,fibrous, and porous piezoelectric composites with anisotropic constituents: A model based on the homogenization method

An analytical framework based on the homogenization method has been developed to predict the effective electromechanical properties of periodic, particulate and porous, piezoelectric composites with anisotropic constituents. Expressions are provided for the effective moduli tensors of n-phase composites based on the respective strain and electric field concentration tensors. By taking into account the shape and distribution of the inclusion and by invoking a simple numerical procedure, solutions for the electromechanical properties of a general anisotropic inclusion in an anisotropic matrix are obtained. While analytical forms are provided for predicting the electroelastic moduli of composites with spherical and cylindrical inclusions, numerical evaluation of integrals over the composite microstructure is required in order to obtain the corresponding expressions for a general ellipsoidal particle in a piezoelectric matrix. The electroelastic moduli of piezoelectric composites predicted by the analytical model developed in the present study demonstrate excellent agreement with results obtained from three-dimensional finite-element models for several piezoelectric systems that exhibit varying degrees of elastic anisotropy.     

Dynamic bulk and shear moduli due to grain-scale local fluid flow in fluid-saturated cracked poroelastic rocks: Theoretical model

Grain-scale local fluid flow is an important loss mechanism for attenuating waves in cracked fluid-saturated poroelastic rocks. In this study, a dynamic elastic modulus model is developed to quantify local flow effect on wave attenuation and velocity dispersion in porous isotropic rocks. The Eshelby transform technique, inclusion-based effective medium model (the Mori–Tanaka scheme), fluid dynamics and mass conservation principle are combined to analyze pore-fluid pressure relaxation and its influences on overall elastic properties. The derivation gives fully analytic, frequency-dependent effective bulk and shear moduli of a fluid-saturated porous rock. It is shown that the derived bulk and shear moduli rigorously satisfy the Biot-Gassmann relationship of poroelasticity in the low-frequency limit, while they are consistent with isolated-pore effective medium theory in the high-frequency limit. In particular, a simplified model is proposed to quantify the squirt-flow dispersion for frequencies lower than stiff-pore relaxation frequency. The main advantage of the proposed model over previous models is its ability to predict the dispersion due to squirt flow between pores and cracks with distributed aspect ratio instead of flow in a simply conceptual double-porosity structure. Independent input parameters include pore aspect ratio distribution, fluid bulk modulus and viscosity, and bulk and shear moduli of the solid grain. Physical assumptions made in this model include (1) pores are inter-connected and (2) crack thickness is smaller than the viscous skin depth. This study is restricted to linear elastic, well-consolidated granular rocks.     

Overall moduli of solids with microcracks: Load-induced anisotropy

For a linearly elastic brittle solid containing microcracks that may be closed or may undergo frictional sliding, a general method is developed for estimating the overall instantaneous moduli which depend on the loading conditions. When the cracks are all open and when they are randomly distributed, then the overall response is isotropic. The moduli for this case have been obtained by Budiansky and O'C onnell (1976). On the other hand, when some cracks close, and when some closed cracks undergo frictional sliding, then the overall response becomes anisotropic and dependent on the loading conditions, as well as on the loading path. The self-consistent method is used to estimate the overall moduli. The effects of crack closure and loadinduced anisotropy are included. Several illustrative examples are worked out, showing the important influence of the load path on the overall response when crack closure and frictional sliding are involved.     

Elastic moduli of a cracked solid

Calculations on the basis of the self-consistent method are made for the elastic moduli of bodies containing randomly distributed flat cracks, with or without fluid in their interiors. General concepts are outlined for arbitrary cracks and explicit derivations together with numerical results are given for elliptic cracks. Parameters are identified which adapt the elliptic-crack results to arbitrary convex crack shapes. Finally, some geometrical relations involving randomly distributed cracks and their traces on cross-sections are presented.     

Rigorous bounds on the effective moduli of composites and inhomogeneous bodies with negative-stiffness phases

We review the theoretical bounds on the effective properties of linear elastic inhomogeneous solids (including composite materials) in the presence of constituents having non-positive-definite elastic moduli (so-called negative-stiffness phases). Using arguments of Hill and Koiter, we show that for statically stable bodies the classical displacement-based variational principles for Dirichlet and Neumann boundary problems hold but that the dual variational principle for traction boundary problems does not apply. We illustrate our findings by the example of a coated spherical inclusion whose stability conditions are obtained from the variational principles. We further show that the classical Voigt upper bound on the linear elastic moduli in multi-phase inhomogeneous bodies and composites applies and that it imposes a stability condition: overall stability requires that the effective moduli do not surpass the Voigt upper bound. This particularly implies that, while the geometric constraints among constituents in a composite can stabilize negative-stiffness phases, the stabilization is insufficient to allow for extreme overall static elastic moduli (exceeding those of the constituents). Stronger bounds on the effective elastic moduli of isotropic composites can be obtained from the Hashin–Shtrikman variational inequalities, which are also shown to hold in the presence of negative stiffness.     

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.     

Bounds for effective elastic moduli of inhomogeneous solid bodies

In connection with the extensive use of various kinds of inhomogeneous materials (glass, carbon and boron reinforced plastics, cermets, concrete, reinforced materials, etc.) in technology, there arises a need to calculate the elastic properties of such systems. Here in each case it is necessary to work out specific methods for finding both elastic fields and effective moduli. Since, as a rule, such methods do not take into account the character of distribution of inhomogeneities in space, which is reflected on the form of the central moment functions [1], they can be referred to a single class and, consequently, can be obtained by a common method [2], In the given paper, by means of the method of solution of stochastic problems for microinhomogeneous solid bodies proposed in the work of the author [2], we find elastic fields and effective moduli in an arbitrary approximation. Depending on the choice of parameters, the latter form bounds within which there lie the exact values of the effective moduli. It is shown that the conditions used earlier for finding these parameters [3] are not the best ones. The effective elastic moduli of an inhomogeneous medium are calculated, and bounds, narrower than the bounds formed in [3], are found for them.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhniki, No. 5, pp. 144–150, September–October, 1973.     

Simple algebraic approximations for the effective elastic moduli of cubic arrays of spheres

Recently, Cohen and Bergman (Phys. Rev. B 68 (2003a) 24104) applied the method of elastostatic resonances to the three-dimensional problem of nonoverlapping spherical isotropic inclusions arranged in a cubic array in order to calculate the effective elastic moduli. The leading order in this systematic perturbation expansion, which is related to the Clausius-Mossotti approximation of electrostatics, was obtained in the form of simple algebraic expressions for the elastic moduli. Explicit expressions were derived for the case of a simple cubic array of spheres, and comparison was made with some accurate results. Here, we present explicit expressions for the effective elastic moduli of base-centered and face-centered cubic arrays as well, and make a comparison with other estimates and with accurate numerical results. The simple algebraic expressions provide accurate results at low volume fractions of the inclusions and are good estimates at moderate volume fractions even when the contrast is high.     

Anisotropic effective higher-order response of heterogeneous Cauchy elastic materials

The homogenization results obtained by Bacca et al. (2013a), to identify the effective second-gradient elastic materials from heterogeneous Cauchy elastic solids, are extended here to the case of phases having non-isotropic tensors of inertia. It is shown that the nonlocal constitutive tensor for the homogenized material depends on both the inertia properties of the RVE and the difference between the effective and the matrix local elastic tensors. Results show that: (i) orthotropic nonlocal effects follow from homogenization of a dilute distribution of aligned elliptical holes and, in the limit case, of cracks; (ii) even under the dilute assumption and isotropic local effective behaviour, homogenization may lead to effective nonlocal orthotropic properties.     

On collinear microcrack–inclusion arrays interaction and related problems

Investigated is a crack problem for an array of collinear microcracks in composite matrix. Inclusions are situated in between the neighbouring microcracks tips and exhibit different elastic properties than matrix. The problem is solved using the technique of distributed dislocations. A developed approximate fundamental solution for a single dislocation lying in a general point between inclusions is employed in the distribution of continuously distributed dislocation to cracks modelling. Stress intensity factor is calculated for various cracks/inclusions geometries and elastic moduli mismatches. Stability and/or instability of the straight microcrack paths is investigated for slowly growing microcracks with inclusions located in between the neighbouring microcracks tips. Applications to periodic microcrack tunnelling and microcracks weakening ahead of the main crack are discussed.     

含柔性涂层的颗粒增强复合材料弹性模量估计

采用线弹簧型弱界面模型来模拟柔性涂层,研究柔性涂层对复合材料宏观弹性模量的影响。首先利用Mori-Tanaka方法和弱界面球形夹杂问题的弹性解,获得单夹杂内部的平均应力和平均应变,进而求得具有柔性涂层的复合材料的宏观弹性模量,并研究界面柔度对复合材料弹性模量的影响。     

Effective elastic properties of matrix composites with transversely-isotropic phases

The present work addresses the problem of calculation of the macroscopic effective elastic properties of composites containing transversely isotropic phases. As a first step, the contribution of a single inhomogeneity to the effective elastic properties is quantified. Relevant stiffness and compliance contribution tensors are derived for spheroidal inhomogeneities. The limiting cases of spherical, penny-shaped and cylindrical shapes are discussed in detail. The property contribution tensors are used to derive the effective elastic moduli of composite materials formed by transversely isotropic phases in two approximations: non-interaction approximation and effective field method. The results are compared with elastic moduli of quasi-random composites.     

Theory of the elastic moduli of a heterogeneous medium

It is shown that the theory of random functions permits the expansion of the effective tensor X~jkl for the elastic moduli with respect to correlation functions and that it leads in the second approximation in the Voigt-Reuss scheme to values that lie to one side of the Xijkl, while in the third approximation it brackets the latter. The analysis is used to refine the Hashin limits to the elastic moduli for a mechanical mixture of isotrcpic components and polycrystalline aggregates of cubic structure.There are two methods for calculating the effective elastic moduli of heterogeneous solids: virial expansion [2] (as a power series in the concentration of one of the components) and the method of correlation functions [2] (expansion with respect to relative fluctuation of the elastic moduli). Identical results should be obtained in the two cases if all terms are incorporated, but great mathematical difficulties restrict one to the lowest approximations. The first approximation in the virial method gives better results when the concentration of one component is low, while the method of correlation functions gives better results when the fluctuations in the elastic moduli are small and the concentrations are similar.Methods have been developed for determining the upper and lower bounds in both approaches, and various schemes of averaging are used for this purpose in the correlation-function method. The upper bound is established by renormalizing the equation of equilibrium, while the lower one is found by renormalizing the equation of incompatibility. The range of the bracketing can be reduced by means of higher approximations. The range can be reduced in the limit to zero, which implies passing from an approximate effective tensor to the true one, which relates the means in stress and strain over the material. Here we show that the two methods of renormalization give identical results when all terms of the series are summed.If the tensor has a Gaussian distribution, the moment functions of odd order are zero, while the even ones are expressed via combinations of the binary functions [3]. However, a mechanical mixture of several components is not Gaussian, and the odd moments are not zero. Splitting of the higher-order correlation functions is possible also for mechanical mixtures having determinate phase interfaces, but this involves various simplifying assumptions. A derivation is given for a moment of arbitrary order, which allows one to formulate the conditions under which such splitting is possible. The results are used in calculating the exact value of the effective bulk modulus for a medium with a homogeneous shear modulus.We are indebted to V. V. Bolotin for a discussion.     

Mechanics of low density materials

Mechanics analyses are used to derive the effective elastic moduli for low density materials. Both open cell and closed cell geometric models are employed in the case of isotropic media. The five independent effective moduli are derived for a low density transversely isotropic medium. Compressive strength, as defined by elastic stability, is also derived for open cell and closed cell isotropic materials. The theoretical results are compared with some experimental results, and also are assessed with respect to previous work.     

确定裂纹体等效弹性模量的边界元方法

采用边界元方法计算含有序分布裂纹的裂纹体在压缩载荷作用下的等效弹性模量，利用一种能适当考虑裂纹有间相作用的自洽理论，建立了相应的迭代格式，通过算例研究了裂纹方向，裂纹面间摩擦系数对裂纹体等效弹性模量的影响。     

The boundaries of the effective elastic moduli for inhomogeneous solids

The calculation of the effective elastic moduli of inhomogeneous solids, which connect the stresses and strains averaged for the material, is accompanied by certain mathematical difficulties owing to correlation relationships of arbitrary orders. Neglect of correlation relationships leads to average elastic moduli, where averaging according to Voigt and Reuss establishes boundaries containing the effective elastic moduli [1]. Approximate values of the latter can be found by taking into account the correlation relationships of the second order in both calculation schemes [2, 3]. Another method of evaluating the true moduli consists of narrowing the boundaries of Voigt and Reuss on the basis of model representations [4-6]. The approximate effective elastic moduli for a series of polycrystals with various common-angle values are presented in [7]. An analysis of the effect of the correlation relationships between the grains of a mechanical mixture of isotropic components on the effective elastic moduli is carried out in [8], although in all the papers just mentioned the use of correlative corrections to narrow the range of elastic moduli is not investigated. Below it is shown that the calculation of the correlation corrections in the second approximation allows the range for the effective moduli to be narrowed.