On the representation formula for well‐ordered elastic composites: a convergence of measure approach

The aim of this paper is to derive an integral representation formula for the effective elasticity tensor for a two‐component well‐ordered composite of elastic materials without using a third reference medium and without assuming the completeness of the eigenspace of the operator ? defined in (2.16) in (J. Mech. Phys. Solids 1984; 32 (1):41–62). As shown in (J. Mech. Phys. Solids 1984; 32 (1):41–62) and (Math. Meth. Appl. Sci. 2006; 29 (6):655–664), this integral representation formula implies a relation which links the effective elastic moduli to the N‐point correlation functions of the microstructure. Such relation not only facilitates a powerful scheme for systematic incorporation of microstructural information into bounds on the effective elastic moduli but also provides a theoretical foundation for de‐homogenization. The analysis presented in this paper can be generalized to an n‐component composite of elastic materials. The relations developed here can be applied to the de‐homogenization for a special class of linear viscoelastic composites. The results will be presented in another paper. Copyright © 2006 John Wiley & Sons, Ltd.     

Thermoelastic homogenization of periodic composites using an eigenstrain-based micromechanical model

This paper presents a study on effective thermoelastic properties of composite materials with periodic microstructures. The overall elastic moduli and coefficients of thermal expansion of such materials are evaluated by a micromechanical model based on the Eshelby equivalent inclusion approach. The model employs Fourier series in the representation of the periodic strain and displacement fields involved in the homogenization procedures and uses the Levin's formula for determining the effective coefficients of thermal expansion. Two main objectives can be highlighted in the work. The first of them is the implementation and application of an efficient strategy for computation of the average eigenstrain vector which represents a crucial task required by the thermoelastic homogenization model. The second objective consists in a detailed investigation on the behavior of the model, considering the convergence of results and efficiency of the strategy used to obtain the approximate solution of the elastic homogenization problem. Analyses on the complexity of the eigenstrain fields in function of the inclusion volume fractions and contrasts between the elastic moduli of the constituent phases are also included in the investigation. Comparisons with results provided by other micromechanical methods and experimental data demonstrate the very good performance of the presented model.     

Boundary element method for calculation of effective elastic moduli in 3D linear elasticity

A boundary integral formulation of the linear elasticity problem for a multi‐component composite is given. The fast BEM solver based on the adaptive cross approximation is then obtained by the data‐sparse representation of the resulting Galerkin matrices. The solver is used to obtain effective elastic moduli of fibre and particle reinforced composites in three dimensions by means of the strain energy equivalence principle. Copyright © 2009 John Wiley & Sons, Ltd.     

Stieltjes representation of the 3D Bruggeman effective medium and Padé approximation

The paper deals with Bruggeman effective medium approximation (EMA) which is often used to model effective complex permittivity of a two-phase composite. We derive the Stieltjes integral representation of the 3D Bruggeman effective medium and use constrained Padé approximation method introduced in [39] to numerically reconstruct the spectral density function in this representation from the effective complex permittivity known in a range of frequencies. The problem of reconstruction of the Stieltjes integral representation arises in inverse homogenization problem where information about the spectral function recovered from the effective properties of the composite, is used to characterize its geometric structure. We present two different proofs of the Stieltjes analytical representation for the effective complex permittivity in the 3D Bruggeman effective medium model: one proof is based on direct calculation, the other one is the derivation of the representation using Stieltjes inversion formula. We show that the continuous spectral density in the integral representation for the Bruggeman EMA model can be efficiently approximated by a rational function. A rational approximation of the spectral density is obtained from the solution of a constrained minimization problem followed by the partial fractions decomposition. We show results of numerical rational approximation of Bruggeman continuous spectral density and use these results for estimation of fractions of components in a composite from simulated effective permittivity of the medium. The volume fractions of the constituents in the composite calculated from the recovered spectral function show good agreement between theoretical and predicted values.     

Asymptotic evaluation of effective complex moduli of fibre‐reinforced viscoelastic composite materials

We propose an asymptotic approach for the evaluation of effective complex moduli of viscoelastic fibre‐reinforced composite materials. Our method is based on the homogenization technique. We start with a non‐trivial expansion of the input plane‐strain boundary value problem by ratios of visco‐elastic constants. This allows to simplify the governing equations to forms analogous to the complex transport problem. Then we apply the asymptotic homogenization method, coming from the original problem on multi‐connected domain to the cell problem, defined on a unit cell of the periodic structure. For the analytical solution of the cell problem we apply the boundary perturbation technique, the asymptotic expansion by a distance between two neighbouring fibres and the method of two‐point Padé approximants. As results we derive uniform analytical representations for effective complex moduli, valid for all values of the components volume fractions and properties.     

非完美界面弹性复合材料中的微分几何方法

首次用微分几何方法计算了含一般旋转椭球体嵌入相的非完美界面弹性复合材料的有效模量·用内蕴几何量表出了能量泛函中的全部界面积分项，由此得到了这种统一嵌入相模型的复合材料有效模量的上下界限·在三种极限情况，即球、盘和针状嵌入相下，本文的结果将退化到Ｈａｓｈｉｎ（１９９２）的结果·     

Numerical calculation of effective elastic moduli for incompressible porous material

Incompressible elastic material with a periodic system of pores is considered. Processes are studied with a typical length which is much more than the typical diameter of pores and the typical distance between pores. Porous material behaves as a certain “effective” material without pores in such processes. The method of calculation of effective moduli based on mathematical homogenization theory is described. The estimates for the effective moduli are proved. The results of numerical calculations of effective moduli for materials with spherical and cubic pores are presented. The dependence of the effective moduli on the volume fracture of pores is investigated. The explicit formulae for effective coefficients are deduced. Comparison with the effective moduli for compressible materials is performed. Presented at the Ninth International Conference on the Mechanics of Composite Materials, Riga, October, 1995. Lomonosov Moscow State University, Department of Mathematics and Mechanics, Moscow. Published in Mekhanika Kompozitnykh Materialov, No. 5, pp. 579–587, September–October, 1996.     

Rheological characteristics of orthotropically reinforced polymeric materials

It is proposed to calculate the rheological characteristics of orthotropically reinforced polymeric materials using a model in accordance with which the composite is assumed to consist of layered blocks, the planes of the layers having one of three possible orientations with given probabilities. The effective rheological characteristics of each block are exactly determined, while those of the composite as a whole are determined in the Voight or Reuss approximations. The first scheme is convenient for calcuating the shear modulus operators, the second for calculating the Young's modulus and Poisson's ratio operators. A numerical computer analysis has been made for the Young's moduli. It is shown that for the anisotropic composites in question each of the technical elastic moduli is characterized by two resolvent Q* operators of real or complex argument. The conditions under which one of the Q* operators can be neglected are determined. The integral representation of the kernel of the Q* operator of complex argument is given in the case where the latter is fractional-exponential with fractional index =–1/2 and the corresponding curves are presented.Moscow. Translated from Mekhanika Polimerov, No. 2, pp. 276–283, March–April, 1972.     

Inverse homogenization with diagonal Padé approximants

The paper formulates inverse homogenization problem as a problem of recovery of Markov function using diagonal Padé approximants. Inverse homogenization or de-homogenization problem is a problem of deriving information about the micro-geometry of composite material from its effective properties. The approach is based on reconstruction of the spectral measure in the analytic Stieltjes representation of the effective tensor of two-component composite. This representation relates the n-point correlation functions of the microstructure to the moments of the spectral measure, which contains all information about the microgeometry. The problem of identification of the spectral function from effective measurements in an interval of frequency has a unique solution. The problem is formulated as an optimization problem which results in diagonal Padé approximation and exact formulas for the moments of the measure. The reconstructed spectral function can be used to evaluate geometric parameters of the structure and to compute other effective parameters of the same composite; this gives solution to the problem of coupling of different effective properties of a two-component composite material with random microstructure. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On the functional properties of the confluent hypergeometric zeta-function

A zeta-function associated with Kummer’s confluent hypergeometric function is introduced as a classical Dirichlet series. An integral representation, a transformation formula, and relation formulas between contiguous functions and one generalization of Ramanujan’s formula are given. The inverse Laplace transform of confluent hypergeometric functions is essentially used to derive the integral representation.     

Some estimates for the torsional rigidity of composite rods

A well‐known problem in elasticity consists in placing two linearly elastic materials (of different shear moduli) in a given plane domain Ω, so as to maximize the torsional rigidity of the resulting rod; moreover, the proportion of these materials is prescribed. Such a problem may not have a classical solution as the optimal design may contain homogenization regions, where the two materials are mixed in a microscopic scale. Then, the optimal torsional rigidity becomes difficult to compute. In this paper we give some different theoretical upper and lower bounds for the optimal torsional rigidity, and we compare them on some significant domains. (© 2007 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Some Riemann boundary value problems in Clifford analysis

In this paper, we mainly study the R_m (m>0) Riemann boundary value problems for functions with values in a Clifford algebra C?(V_{3, 3}). We prove a generalized Liouville‐type theorem for harmonic functions and biharmonic functions by combining the growth behaviour estimates with the series expansions for k‐monogenic functions. We obtain the result under only one growth condition at infinity by using the integral representation formulas for harmonic functions and biharmonic functions. By using the Plemelj formula and the integral representation formulas, a more generalized Liouville theorem for harmonic functions and biharmonic functions are presented. Combining the Plemelj formula and the integral representation formulas with the above generalized Liouville theorem, we prove that the R_m (m>0) Riemann boundary value problems for monogenic functions, harmonic functions and biharmonic functions are solvable. Explicit representation formulas of the solutions are given. Copyright © 2009 John Wiley & Sons, Ltd.     

Homogenization of Unidirectional and Arbitrarily Oriented Fiber-Reinforced Materials by the Method of Conditional Moments

In this paper the method of conditional moments is developed for the case of a two–component matrix composite with randomly distributed unidirectional and arbitrarily oriented ellipsoidal inclusions. The algorithm for determination of the effective elastic properties of composites from the given elastic constants of the components and geometrical parameters and orientation of inclusions is discussed. It is assumed that the components of the composite show orthotropic symmetry of thermoelastic properties. As a numerical example arbolite (straw particle inclusions in a cement matrix) is considered. The dependencies of Young's moduli, Poisson's ratios and shear moduli from the concentration of inclusions and for certain orientations of the inclusions are predicted and discussed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Exact relations for effective tensors of composites: Necessary conditions and sufficient conditions

Typically the elastic and electrical properties of composite materials are strongly microstructure dependent. So it comes as a nice surprise to come across exact formulae for effective moduli that are universally valid no matter what the microstructure. Such exact formulae provide useful benchmarks for testing numerical and actual experimental data and for evaluating the merit of various approximation schemes. They can also be regarded as fundamental invariances existing in a given physical context. Classic examples include Hill's formulae for the effective bulk modulus of a two‐phase mixture when the phases have equal shear moduli, Levin's formulae linking the effective thermal expansion coefficient and effective bulk modulus of two‐phase mixtures, and Dykhne's result for the effective conductivity of an isotropic two‐dimensional polycrystalline material. Here we present a systematic theory of exact relations embracing the known exact relations and establishing new ones. The search for exact relations is reduced to a search for matrix subspaces having a structure of special Jordan algebras. One of many new exact relations is for the effective shear modulus of a class of three‐dimensional polycrystalline materials. We present complete lists of exact relations for three‐dimensional thermoelectricity and for three‐dimensional thermopiezoelectric composites that include all exact relations for elasticity, thermoelasticity, and piezoelectricity as particular cases. © 2000 John Wiley & Sons, Inc.     

Determination of the elastic moduli of reinforced materials

A method is proposed for determining the elastic moduli of reinforced materials at arbitrary concentrations of the components. The method is based on the successive symmetrization of the frequency function of the moduli of the component elements of a two-component material with subsequent application of the correlation theory [1]. Algorithms are constructed for determining the elastic moduli of layered, fibrous, and granular materials and typical calculations are presented.Mekhanika Polimerov, Vol. 4, No. 1, pp. 78–85, 1968     

Effective transverse elastic moduli of composites at non-dilute concentration of a random field of aligned fibers

We consider a transversal loading of a linearly elastic isotropic media containing the identical isotropic aligned circular fibers at non-dilute concentration c. By the use of solution obtained by the Kolosov–Muskhelishvili complex potential method for two interacting circles subjected to three different applied stresses at infinity, and exact integral representations for both the stress and strain distributions in a microinhomogeneous medium, one estimates the effective moduli of the composite accurately to order c². Received: March 4, 2003; revised: August 8, 2003     

COMPUTATION OF THE RELAXATION OSCILLATION PERIOD IN LOTKA‐VOLTERRA SYSTEMS

ABSTRACT. This work surveys techniques of Grasman and Veling [1973], Vasil'eva and Belyanin [1988] and Shih [1996] for computing the relaxation oscillation period of singularly perturbed Lotka‐Volterra systems. Grasman and Veling [1973] used an implicit function theorem to derive an asymptotic formula for the period; Vasil'eva and Belyanin [1988] employed a method of matched asymptotic expansions to obtain an approximation to the period; Shih [1996] obtained two (exact) integral representations for the period in terms of two inverse functions W(–k, x) of xexp(x). These results are compared numerically and asymptotically. In particular, the integral representation of the period in Shih [1996] is computed numerically using a Gauss‐Tschebyscheff integration rule of the first kind, and is further investigated asymptotically by virtue of the asymptotics of W(–k, x), Laplace's method, and a method of consequent representation. Computational results indicate that the Gauss‐Tschebyscheff approximation of the period in Shih [1996] is uniformly accurate for a wide range of the singular parameter (? in the paper).     

Well‐conditioned boundary integral formulations for high‐frequency elastic scattering problems in three dimensions

We construct and analyze a family of well‐conditioned boundary integral equations for the Krylov iterative solution of three‐dimensional elastic scattering problems by a bounded rigid obstacle. We develop a new potential theory using a rewriting of the Somigliana integral representation formula. From these results, we generalize to linear elasticity the well‐known Brakhage–Werner and combined field integral equation formulations. We use a suitable approximation of the Dirichlet‐to‐Neumann map as a regularizing operator in the proposed boundary integral equations. The construction of the approximate Dirichlet‐to‐Neumann map is inspired by the on‐surface radiation conditions method. We prove that the associated integral equations are uniquely solvable and possess very interesting spectral properties. Promising analytical and numerical investigations, in terms of spherical harmonics, with the elastic sphere are provided. Copyright © 2014 John Wiley & Sons, Ltd.     

横向剪切对双模量复合材料叠层矩形板非线性弯曲的影响

本文探讨了动力松弛(DR)法在双模量复合材料叠层矩形板非线性弯曲问题中的应用。在分析中分别采用叠层板大挠度经典理论和计及大转动(在Von Karman意义上)的复合材料叠层板剪切变形理论。我们发现,对于考虑横向剪切变形的非线性弯曲问题,如何计算虚拟密度以控制数值计算的稳定性,仍然需要进一步研究。本文提出了一种虚拟密度的计算方法,从而保证了本课题数值计算的稳定性。文中介绍了用DR法求解双模量复合材料叠层板非线性弯曲的主要步骤,给出了由轻度双模量材料(Born-Epxy(B-E))和高度双模量材料(Aramid-Rubber(A-R)和Polyester-Rubber(P-R))的两层正交叠层简支矩形板在正弦分布载荷及均布载荷作用下的非线性弯曲特性的数值结果。将所得结果和小挠度分析结果及普通复合材料的结果作了比较,并分析了横向剪切变形对无量纲中心挠度的影响。     

The transmutation operator method for efficient solution of the inverse Sturm‐Liouville problem on a half‐line

The inverse Sturm‐Liouville problem on a half‐line is considered. With the aid of a Fourier‐Legendre series representation of the transmutation integral kernel and the Gel'fand‐Levitan equation, the numerical solution of the problem is reduced to a system of linear algebraic equations. The potential q is recovered from the first coefficient of the Fourier‐Legendre series. The resulting numerical method is direct and simple. The results of the numerical experiments are presented.