Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions

Size: px
Start display at page:

Download "Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions"

Transcription

1 Arch Appl Mech DOI 1.17/s ORIGINAL R. Rodríguez-Ramos R. Guinovart-Díaz J. Bravo-Castillero F. J. Sabina H. Berger S. Kari U. Gabbert Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions Received: 31 March 28 / Accepted: 9 May 28 Springer-Verlag 28 Abstract In the present work, unified formulae for the overall elastic bounds for multiphase transversely isotropic composites with different geometrical types of inclusions embedded in a matrix are calculated, including the spherical and long or short continuous cylindrical fiber cases. The influence of the different geometrical configurations of the inclusions on the composites is studied. The transversely isotropic effective bounds are obtained by applying the variational formulation for anisotropic composites developed by Willis, which relies on expressions for the static transversely isotropic Green s function. Some numerical calculations and comparisons with the effective coefficients derived from the self-consistent approach, asymptotic homogenization method, and finite element method (FEM) are shown for different aspect ratio values, exhibiting good agreement. Keywords Particle-reinforced composites Variational bounds Transversely isotropic components Green s function 1 Introduction A classical problem in solid mechanics is the determination of the effective elastic properties of a composite material made up of a statistically isotropic random distribution of isotropic and elastinclusions embedded in a continuous isotropic and elastic matrix. Various models have been developed to compute the two independent elastic constants of the composite,as in the Mori Tanakaestimates [2,25], the self-consistent scheme [16,17], and the differential method [22]. They are based on the mean-field approximation, which assumes that the stress and strain fields in the matrix and in the inclusion are adequately represented by their volume-averaged values, and differ in the way they account for the elastinteraction between the inclusions. In other cases, rigorous bounds for the elastic properties of isotropic two-phase composites were obtained from variational principles. Besides the Voigt Reuss bounds, which are usually very wide, Hashin and Shtrikman [12] derived much tighter bounds by employing variational principles which involve polarization fields. The two-point correlation bounds are the best possible ones given only the volume fraction information. They were improved R. Rodríguez-Ramos R. Guinovart-Díaz J. Bravo-Castillero Facultad de Matemática y Computación, Universidad de La Habana, San Lázaro y L, Vedado, CP 14 Habana 4, Cuba F. J. Sabina Instituto de Investigaciones en Matemáticas Aplicadas y en Sistemas, Universidad Nacional Autónoma de México, Apartado Postal 2-726, Delegación de Alvaro Obregón, 1 Mexico D.F., Mexico H. Berger (B) S. Kari U. Gabbert Institute of Mechanics, Otto-von-Guericke-University of Magdeburg, Universitaetsplatz 2, 3916 Magdeburg, Germany Tel.: Harald.Berger@Masch-Bau.Uni-Magdeburg.DE

2 R. Rodríguez-Ramos et al. by three-point bounds [4,24,35], which incorporate information about the phase arrangement through the statistical correlation parameters mentioned above. Derivation of the effective coefficients for anisotropic elastic composites using the asymptotic homogenization method (AHM) and self-consistent methods (SCM) has been studied by Pobedrya [29] and Willis [39], and more recently by Levin et al. [21], Sabina et al. [31], Guinovart-Diaz et al. [1], and Rodriguez-Ramos et al. [3]. The estimation of bounds for the overall properties of anisotropic composites has also been considered by many authors. The main overall elastic moduli of fiber composites with transversely isotropic phases are shown to be connected by simple universal relations which are independent of the geometry at given concentration. Upper and lower bounds obtained in terms of phase properties and concentrations were reported by Hill [15]. On the other hand, the plane-strain bulk modulus and the two shear moduli of multiphase transversely isotropic fiber-reinforced materials of arbitrary transverse phase geometry were bounded in terms of phase moduli and phase volume fractions by Hashin [13]. Moreover, expressions and bounds for the five effective elastic moduli of a unidirectional fiber composite, consisting of transversely isotropic fibers and matrix, were derived on the basis of analogies between isotropic and transversely isotropic elasticity equations by Hashin [14]. Bounds of Hashin Shtrikman type and self-consistent estimates for the overall properties of composites, which may be anisotropic, were developed in Willis [36]. In Nemat-Nasser and Hori [27], the Hashin-Shtrikman variational principle is applied to solids with periodic microstructure, and bounds on the overall moduli are obtained by defining energy functionals for the eigenstrain or eigenstress fields in the equivalent homogeneous solid. Finally, bounds for the overall elastic energy of the unit cell, and hence bounds for the overall elasticity parameters of elastic solids with periodic microstructure, were obtained. Moreover, improvable bounds on overall properties of heterogeneous finite solids were also studied by Nemat-Nasser and Hori [26]. In the present paper, the derivationof the upperand lower boundsfor N-phase elastic transversely isotropic composite materials using the procedure of Willis [4] is obtained. The constituents of the two- or three-phase composite studied here are transversely isotropic materials. Unified formulae for the bounds as a function of the aspect ratio of the inclusions are given. Some numerical results for each elastic coefficient are presented. Section 2 starts with the basic equations of Willis and gives the variational bounds explicitly. In Sect. 3 some numerical calculations of the bounds for different geometrical configurations of the inclusions and comparisons with other models are done. Section 4 provides some final remarks. 2 Variational bounds for the anisotropic elastic multiphase composite The problem of estimating bounds for the mean elastic energy density, W (e), in a given sample of material, when it is subjected to loads that generate within it a mean strain field e, has been addressed by Willis [4]. In that chapter, in particular, for a two-phase composite, i.e., a particle-reinforced matrix with anisotropic elastic constituents, he gave a general explicit formulae (21.38) for the upper and lower bounds of the mean energy, which relies on the evaluation of a certain P function depending on the derivatives of the Green s function G. The Willis formulae are given at the outset as follows: where and { N } 1 L U = c r [I (L o L r ) P] 1 r=1 { N } 1 L L = c r [I + (L r L o ) P] 1 r=1 1 2 el Le W (e) 1 2 el Ue, (1) c r [I (L o L r ) P] 1 L r (2) r=1 c r [I + (L r L o ) P] 1 L r (3) are the upper and lower bounds, respectively; the value r = 1 denotes the matrix and r = 2,...,N correspond to the inclusions; c r are the volume fractions of each phase, respectively, such that c 1 + +c N = 1; L r is r=1

3 Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions the anisotropic elastic tensor for r-phase and L o is the anisotropic elastic tensor of the comparison material. P is related to the integral of the Green s function which is defined as a linear operator as follows (see [4]) ( τ) ij (x) = jkl (x, x )τ kl (x )dx, (4) with jkl (x, x ) = 2 G ij x k x l. (5) (ij)(kl) The parenthesis stand for symmetrization, see Willis [4]. Expressions (2) and(3) give an upper and lower bound, respectively, for any L o such that the quadratic forms e(l o L r )e and e(l r L o )e are positive definite for every r. Now, consider that the inclusions are aligned spheroids of aspect ratio δ.letox 3 be the axis of alignment and let the elastic properties of both the inclusions and the matrix be transversely isotropic with its axis parallel to Ox 3. Then both media are characterized by five parameters and the operations indicated in (2) and(3) can be carried out. It is useful to introduce the symbolic notation described in Appendix A. Examples of this can be found in Sabina et al. [32] and Sabina and Willis [33]. In order to work out (2)and(3) it remains to obtain P in the notation of Appendix A, P = (2P 1, P 2, P 3, P 4, 2P 5, 2P 6 ), (6) where P i (i = 1, 2,..., 6) are given in Appendix B. Using the algebraic rules of Appendix A and after some algebraic manipulations, the bounds for the transversely isotropic elastic composite are obtained. They have the following general form: k L k k U, l L l l U, l L l l U, n L n n U, m L m m U, p L p p U, where the particular expressions for the bounds depend on the number of phases. The analytic expressions for the lower and upper bounds are given as follows: ( k B = 1 N a i 4 i ( l B = 1 N d i 4 i ( l B = 1 N e i 4 i ( n B = 1 N d i 4 i (b i k i e i l i ) + (b i k i e i l i ) + e i b i (a i n i 2d i l i 2 ) + N i (b i l i e in i ) + b i a i ) (a i l i 2d i k i ), ) (a i l i 2d i k i ), 2 ) (b i l i 2 e in i ), i ) (a i n i 2d i l i 2 ), i m B = 2 m i f i ( N ) 1, p B = 2 f i p i g i ( N g i ) 1, (7)

4 R. Rodríguez-Ramos et al. where a i = 1 + 4(k i k )P 1 + 2(l i l )P 2, b i = 1 + 2(l i l )P 3 + (n i n )P 4, d i = 2(l i l )P 1 + (n i n )P 2, e i = (l i l )P 4 + 2(k i k )P 3, f i = 1 + 4(m i m )P 5, g i = 1 + 4(p i p )P 6, = 1 2 a ib i d i e i, ( = 1 N a i b i 4 2 d i ) e i. (8) B = L, U denotes the lower and upper bounds, respectively, related to the choice of the comparison body L o. In the calculations that follow the comparison moduli of L o is taken as L o = min(l r ) for the lower bound and L o = max(l r ) for the upper bound, r = 1, 2,...,N. It is not difficult to verify that the bounds in Eq. (7) are identical to those found originally by Hashin and Shtrikman [12] and Hashin [13]. 3 Numerical results In this section, the general expressions for the upper and lower bounds (7) are applied to two- and three-phase composites with different types of inclusions (Fig. 1) for the prediction of their effective elastic properties. 3.1 Two-phase composites (N = 2) It is to be noted that the numerical values of l L, l U, l L,andl U are such that l L = l L and l U = l U, although the given expressions look different. These formulae are provided for completeness and can be used as a further check on the final results. It is worthwhile to note that the upper and lower bounds of m are both equal to the Voigt average when m 1 = m 2. Moreover, the upper and lower role is interchanged when m 2 m 1 changes sign, i.e., for m 2 < m 1, (7) is a lower bound, while for m 1 < m 2,(7) is an upper bound. The same considerations can be applied to p. Fig. 1 Different configuration of composites: a short fibers, b spherical inclusions, c long fibers, and d three-phase fibrous composite

5 Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions In the calculations that follow the comparison moduli of L o is taken as that of the matrix L 2 for the lower bound and as that of the inclusion L 1 for the upper bound. The presentation of the results is divided according to different shapes of inclusions: long cylindrical fibers, spherical inclusions, and short cylindrical fibers (see Fig. 1). The inclusions can be considered in general as ellipsoidal inclusions with different aspect ratios of length to width denoted by δ. Long cylindrical fibers The first example involves ellipsoidal inclusions of the piezoelectric material BaTiO 3 embedded in the piezomagnetic medium CoFe 2 O 4. Only the elastic constants are involved in the calculations of this paper, i.e., the bounds relate only to the elastic properties of both materials. The elastic constants (in the Voigt notation) are given in Table 1 and were taken from Aboudi [1]. Figure 2 displays the five stiffnesses k, l, n, m, and p in GPa plotted versus the inclusion volume fraction c 2 for an aspect ratio δ = 1 and δ =.1. Since the inclusion is softer than the matrix, it shows a characteristic pattern of a weakly reinforced composite. The properties are monotone decreasing functions of c 2. Moreover, the upper and lower bounds calculated using (7) are also shown. Although it is not shown in the figure, it is convenient to comment that, for δ = 1, the effective coefficients k, l, n, m, and p were computed using the asymptotic homogenization method (AHM) [1] and the self-consistent effective field method (SCEF) [21]. All curves agree with the lower and upper bounds (7). These results can be compared with other calculations by Aboudi [1] using the Mori Tanaka method, the generalized method of cells, and a numerical implementation of the asymptotic homogenization method, which also show curves lying on top of each other for each stiffness. The comparison is reasonably good except for a slight maximum near c 2 = 1for m in Aboudi s approach (Fig. 2). This may be due to the numerical implementation of the methods for a rectangular shaped inclusion, which differs from the very long ellipsoidal one used here. The bounds do not overlap for δ =.1, except for m and p. Another combination of materials is now considered: inclusions of PZT-5 embedded in a polymer matrix (see the material data in Table 1). The results (7) are also valid for an isotropic matrix such as the epoxy considered here. Figure 3 shows k, l, m,and n (in GPa) versus c 2, respectively, for a long ellipsoidal inclusion [i.e., δ = 1 in (7)]. This time, the upper and lower bounds are wide but not as wide as the Voigt and Reuss bounds (not shown here). All curves are monotone increasing. The self-consistent effective field (SCEF) solution and finite-element method (FEM) [5 9] lie on top of the lower bound. The asymptotic homogenization method (AHM) also lies on top of the lower bound. Therefore, the lower bound can be considered as an effective value for that composite. In order to validate the obtained formulae (7) an ideal material (stiffer) embedded in a BaTiO 3 matrix (softer) is considered. It was chosen so that m 1 = m 2, i.e., a Hill s solid [15]. All the curves for the effective coefficients k, l, n, and p agree with each other, that is to say, the upper and lower bounds, self-consistent effective field (SCEF), asymptotic homogenization (AHM), and Hill s method curves coincide n UB, δ=1 LB, δ=1 UB, δ=.1 LB, δ= k 1 5 l m p Volume fraction c 2 Fig. 2 Plot of the upper and lower bound (7) of the five effective elastic constants k, l, m, n,and p as a function of the inclusion volume fraction c 2 for an aspect ratio δ = 1 and δ =.1. The materials are inclusion of BaTiO 3 embedded in a matrix of CoFe 2 O 4

6 R. Rodríguez-Ramos et al. Table 1 Elastic material properties in GPa Material C 11 = k + m C 12 = k m C 13 = l C 33 = n C 44 = p BaTiO CoFe 2 O PZT Polymer Isotropinclusion Epoxy matrix I PZT-7A Epoxy matrix II Yttria coating fiber SiC matrix ATH filler PMMA matrix Effective Coefficient k Upper Bound AHM SCEF FEM Lower Bound Effective Coefficient l Effective Coefficient m Effective Coefficient n Volume fraction fiber Volume fraction fiber Fig. 3 Plots of the four effective elastic constants k, l, m,and n versus c 2 for long cylindrical fibers. The materials here are PZT-5 inclusions in an polymer matrix. The upper and lower bound (7) and the results of the self-consistent effective field (SCEF), finite element method (FEM), and asymptotic homogenization methods (AHM) areshown Spherical inclusions As another example, two isotropic materials of Young s moduli E 1 = 3GPa and E 2 = 76 GPa and Poisson s ratios ν 1 =.4 andν 2 =.23 are studied since there are some experimental data available [34] for spherical inclusions (δ = 1). The upper and lower bounds of the axial Young s modulus Ē a (in GPa) and the dimensionless transverse shear modulus m/m 1 computed from (7)usingδ = 1 and finite element method (FEM) [19] are plotted against the inclusion volume fraction c 2 in Fig. 4. The experimental data and FEM results lie between the calculated bounds. One relevant fact is that the FEM results and the experimental data coincide very well. Short cylindrical fibers Figure 5 shows plots the four elastic constants k, m, n, and p (in GPa) versus c 2. The composite is made of isotropinclusions embedded in an epoxy matrix I. The material parameters are in Table 1. The upper and

7 Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions 8 6 UB FEM Exp LB E a 4 m/m Volume fraction c Volume fraction c 2 1 Fig. 4 Plots of the bounds of effective axial Young s modulus and effective property m/m 1 versus c 2 for spherical inclusions. A comparison between the bounds, the FEM results, and the experimental data of Smith for two isotropic materials is shown lower bounds (7) and finite element method (FEM) [2] are shown for aspect ratio δ = 3 (short fibers). In this case, the lower bound can be considered as an effective property for the parameters m and p, respectively. The next example in Fig. 6 considers four elastic constants k, l, n, and p (in GPa) for the materials combination of PZT-7A inclusions in an epoxy matrix II for two fixed values of the inclusion volume fraction c 2 =.3and.8 as a function of the aspect ratio δ. For relatively large values of δ the upper bounds are almost constant for all the elastic constants. The lower bound is likewise almost constant except for n. Near δ =, all curves change very rapidly Upper Bound FEM Lower Bound Effective Coefficient (GPa) k m Effective Coefficient (GPa) n p Volume fraction fiber Volume fraction fiber Fig. 5 Plots of the four effective elastic constants k, m, n,and p versus c 2. The materials here are isotropinclusions in an epoxy matrix. The upper and lower bound (7) and FEM results are shown for aspect ratio δ = 3(short fibers)

8 l R. Rodríguez-Ramos et al. k UB, c =.8 2 PZT 7A/ LB, cepoxy =.8 2 UB, c 2 =.3 LB, c 2 = n 6 p Aspect ratio δ Aspect ratio δ Fig. 6 Plots of the upper and lower bounds of effective moduli k, l, n,and p against the aspect ratio δ for two values of c 2 =.3 and.8. The materials are PZT-7A inclusions in an epoxy matrix 3.2 Three-phase composites (N = 3) In the calculations that follow the comparison moduli of L o are taken as L o = min(l 1, L 2 ) for the lower bound and L o = max(l 1, L 2 ) for the upper bound. The presentation of the results is divided according to the different shapes of the inclusions as before, i.e, long cylindrical fibers, spherical inclusions, and short cylindrical fibers (Fig. 1). Long cylindrical fibers As an example of the results obtained for a three-phase composite, in Fig. 7a, b, we display the normalized shear moduli (transversal and longitudinal) for a coated-fiber composite versus the fiber volume fraction c 3. They have been calculated from the analytical expressions for square and hexagonal symmetries reported by Guinovart-Díaz et al. [11] andcompared with the correspondingupperand lower boundsmoduliobtainedfrom (7), for an aspect ratio of δ = 1 and interphase volume fraction of c 2 =.5. The properties of the constituent materials (SiC matrix/yttria coating/carbon fiber) were taken from Benveniste et al. [3](see Table 1). In Fig.7a we observe that the effective shear modulus m/m 1 is decreasing and stays within the bounds when the volume fraction of fiber increases. For small values of c 3, the effective coefficient m/m 1 lies on the upper bound for both symmetries, while for high concentrations of fiber it lies on the lower bound for square symmetry. As illustrated in Fig. 7b, the other effective properties lie within the bounds and are very close to the corresponding upper bound for all values of the fiber volume fraction. The effective properties and the bound coefficients of the coated-fiber composites have been calculated using the present models for other systems of materials that appear in Benveniste et al. [3] and agree quite well with the effective properties calculated by using the Mori and Tanaka approach in Benveniste et al. [3]. Due to the very small thickness of the coating in the systems considered, the effective properties are almost equal to the effective properties obtained in the lower bound and in the uncoated fiber case. The bound formulae are functions of the aspect ratio δ and they are used for the estimations of the effective properties of particulate reinforced composites. Spherical inclusions Figure 8 presents the effective Young s modulus for spherical inclusions. A comparison between the models reported in Nie and Basaran [28] and Ju and Chen [18] with the lower bound (7) for the three-phase composite

9 Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions Upper Bound AHM Hex AHM Square Lower Bound m / m 1.5 p /p Fiber volume fraction c 3 Fiber volume fraction c 3 Fig. 7 Bounds and AHM results for the normalized effective shear moduli (transversal and longitudinal) for a three-phase composite (SiC matrix/yttria coating/carbon fiber) as a function of volume fraction of filler c 3 is shown. The interphase thickness is c 2 =.1 and the fiber volume fraction is c 3 =.48. The composite is made of Poli-metil metacrilato (PMMA) matrix/alumina trihydrate (ATH) filler. The material parameters are listed in Table 1. The effective Young s modulus and the lower bound can be calculated from the well-known expression E = 9km/(3k + m) for an isotropic composite. Both models are very close to the lower bound, which gives a good estimation for this property. The bounds and the effective properties of the bulk and shear normalized properties k/k 1 and µ/µ 1, respectively, for transversely isotropic elastic composite as a function of interphase region elastic modulus q = E 3 /E 2 and ν 2 = ν 3 were studied as well. The theoretical models 12 1 Nie Basaran (25) Ju and Chen (1994) Lower Bound SCS 8 E (GPa) Interface parameter log(q) Fig. 8 Behavior of the effective Young s modulus for spherical inclusions in a PMMA matrix/ath filler composite as a function of the interphase region elastic modulus q = E 3 /E 2

10 R. Rodríguez-Ramos et al. for the calculations of the above properties were taken from Nie and Basaran [28] and Ju and Chen [18]. The comparison between these two models and the bounds (not shown) illustrates that both models are consistent with the bounds and are near to the lower bound. The bulk modulus in the model of Ju and Chen for the lower values of the parameter q is not consistent with the bounds (it falls outside of the bounds). This indicates that the model presented by Nie Basaran gives a better approximation of the overall bulk modulus, as was mentioned recently in the work reported by Nie and Basaran [28]. However, the shear modulus in both mentioned models falls within the bounds. A comparison for three-phase spherical inclusions with the numerical results by FEM [19,2]isshownin Table 2. The numerical results show good concordance between the FEM results and the variational bounds. In Fig. 9, the behavior of the axial shear modulus p of the composite as a function of the aspect ratio δ and the volume fraction c 3 is illustrated. The composite has the constituents SiC matrix/yttria coating/carbon inclusion (see Table 1). The volume fraction of the coating is c 2 =.5. The aspect ratio δ influences the behavior of the bounds of the composite. The shaded and the mesh surfaces indicated the lower and upper bound, respectively, for the value of the axial shear effective coefficient p with respect to δ (aspect ratio) and c 3 (volume fraction). Notice that for values of δ 2 (short fibers) the bounds are tighter than for values of δ 2. Table 2 Comparison between effective properties using FEM and variational bounds [lower bound (LB) and upper bound (UB)] in GPa Sphere Interphase C 11 C 12 volume fraction volume fraction LB FEM UB LB FEM UB C 13 C C 44 C p (GPa) c δ 8 1 Fig. 9 Surface of the effective axial shear modulus p of the three-phase composite SiC matrix/yttria coating/carbon inclusion as a function of the aspect ratio δ and volume fraction of inner inclusion c 3

11 Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions 4 Conclusions Explicit formulae for the variational upper and lower bounds for transversely isotropic ellipsoidal inclusions embedded in a transversely isotropic matrix are obtained. The inclusions are all aligned in the x 3 direction. The geometrical configuration of the inclusions exerts a significant contribution on the effective properties of the composite and therefore on their bounds. The method used is based on that of Willis [4], who derived a tensor formula for the bounds. Using the symbolic notation and the algebra of Hill [15], it was possible to obtain explicit formulae for the bounds of the effective elastic constants k, l, m, n and p. Several examples for different two- and three-phase composites are shown. Acknowledgments This work was sponsored by CoNACYT project no F and by DFG Graduiertenkolleg 828 Micro Macro Interactions in Structured Media and Particle Systems. The provisions of the Basic Sciences Program Project CITMA no. 9/24 and the Department of Basic Science of the Monterrey Institute of Technology, Campus of México State are also acknowledged. The authors gratefully acknowledge the suggestions and comments given by Prof. John Willis. Thanks are due to Ms. Ana Pérez Arteaga and Fernando Rodríguez for computational support. Appendix A: Transversely isotropic symmetry considerations In the case that the fourth-order tensor L of elastic moduli has transversely isotropic symmetry relative, say, to the Ox 3 axis, then, it can be represented symbolically using the notation employed by Hill [16] in the form L = ( 2k, l, l, n, 2m, 2p ), where the appropriate elastic constants 2k, l, l, n,2m,and 2p appear in the constitutive relation between stress σ and strain e, viz.σ = Le. It is written explicitly as follows 1 2 (σ 11 + σ 22 ) = k(e 11 + e 22 ) + le 33, σ 33 = l (e 11 + e 22 ) + ne 33, σ 11 σ 22 = 2m(e 11 e 22 ), σ 12 = 2me 12, σ 13 = 2pe 13, σ 23 = 2pe 23. It is assumed that l = l since the diagonal symmetry is lost, in general, in a triple product of transversely isotropic tensors. Moreover, let L 1 = (2k, l, l, n, 2m, 2p) and L 2 = ( 2k, l, l, n, 2m, 2p ) be two fourthorder transversely isotropic tensors; then the following algebraic rules hold L 1 + L 2 = ( 2(k 1 + k 2 ), l 1 + l 2, l 1 + l 2, n 1 + n 2, 2(m 1 + m 2 ), 2(p 1 + p 2 ) ), L 1 L 2 = ( 4k 1 k 2 + 2l 1 l 2, 2l 1 k 2 + n 1 l 2, l 1 n 2 + 2k 1 l 2, 2l 1l 2 + n ) 1n 2, 4m 1 m 2, 4p 1 p 2. Then, it follows that the inverse of L, the compliance tensor M, isgivenby ( n M = L 1 = 2, l 2, l 2, k ), 1 2m, 1, 2p where = kn ll. The identity tensor I in this notation is I = (1,,, 1, 1, 1). In particular, if L is isotropic, with Lamé moduli λ, µ, then it can also be represented in this way as L = (2(λ + µ),λ,λ,λ+ 2µ, 2µ, 2µ).

12 R. Rodríguez-Ramos et al. Appendix B This appendix is concerned with the definition of the components of P in (7) and the corresponding Green s function. A suitable representation of the Green s function for the time-reduced elastodynamic equation developed by Willis [37,38] appears slightly modified in Sabina et al. [32]. In the formulae below, two terms are given. The first one represents the static Green function, which is used here; the second term, proportional to ω is the term that includes the contribution to dynamics. Thus the complete formulae are given, so that the contribution of the static term is obtained by setting the circular frequency to zero. G ij (x) = 1 8π 2 3 n=1 ξ =1 ds U i n(ξ)u n j (ξ) ( cn 2 δ(ξ x) + iϖ ) exp (iϖ ξ x c n ). (1) 2c n Here, Ui n and c n are the polarization and wave speed of a plane wave propagating in the direction of the unit vector ξ in the reference medium with properties L o and ρ o. Thus, [ (Lo ) ijkl ξ j ξ l cn 2 (ρ ] o) ik U n k =. (2) Note that here the vectors Ui n are orthonormalized relative to the tensor ρ o. The associated wavenumber k n = w/c n. In the case of transverse isotropy, the three solutions of (2) correspond to the three different wave types. The eigenvectors and eigenvalues of (2)aregiveninSabinaetal.[32]. The integrations of the derivatives of the Green s function to obtain P are carried out over the spheroidal inclusions x1 2 a 2 + x2 2 a 2 + x2 3 c 2 = 1 (3) of aspect ratio δ = c/a. The components of P are given below as: 2P 1 2k = 1 2ρ I 1 P 2 P 3 l = l = (1 u 2 ) [ m 2 1 H 1(ζ(u)) + m 2 3 H 3(ζ(u)) ] du, 1 2(ρ I ρ III ) 1/2 P 4 n = 1 1 u 2 [ m 2 3 ρ H 1(ζ(u)) + m 2 1 H 3(ζ(u)) ] du, III 2P 5 2m = 1 4ρ I 2P 6 2p = (1 u 2 )um 1 m 3 [H 1 (ζ(u)) H 3 (ζ(u))] du, (1 u 2 ) [ m 2 1 H 1(ζ(u)) + H 2 (ζ(u)) + m 2 3 H 3(ζ(u)) ] du, {[ u 2 m2 1 + (1 u 2 ) m2 3 ρ 1 ρ III ] +2(1 u 2 ) 1/2 m 1 m 3 u (ρ I ρ III ) 1/2 H 1 (ζ(u)) + u2 H 2 (ζ(u)) ρ I [ ] + u 2 m2 3 + (1 u 2 ) m2 1 2(1 u 2 ) 1/2 m 1 m 3 u ρ 1 ρ III (ρ I ρ III ) 1/2 } H 3 (ζ(u)) du. (4)

13 Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions Here, ζ(u) = [ 1 + (δ 2 1)u 2] 1/2, m 1 (u) = ρ 1/2 [ I po (1 u 2 ) + n o u 2 ρ III c1 2 ] D 1/2, m 3 (u) = ρ 1/2 III (l o + p o)(1 u 2 ) 1/2 u D 1/2, D = ρ I [ po + (n o p o )u 2 ρ III c 2 1] 2 + ρiii (ló + p o) 2 (1 u 2 )u 2, H n (ζ ) = δ ε(k n aζ) c 2 n ζ 3, 3(1 iz) ε(z) = z 3 e iz (sin z z cos z), (5) where p o, n o, l o, ρ I,andρ III are properties of the reference medium L o, and letting ϖ. The functions H n (ζ ) and ε(z) in the limit of ϖ are H n (ζ ) = δ c 2 n ζ 3, ε(z) = 1. References 1. Aboudi, J.: Micromechanical analysis of fully coupled electro-magneto-thermo-elastic multiphase composites. Smart Mater. Struct. 1, (21) 2. Benveniste, Y.: A new approach to the application of Mori Tanaka s theory in composite materials. Mech. Mater. 6, (1987) 3. Benveniste, Y., Dvorak, G.J., Chen, T.: Stress fields in composites with coated inclusions. Mech. Mater. 7, (1989) 4. Beran, M.J., Molyneux, J.: Use of classical variational principles to determine bounds for the effective bulk modulus in heterogeneous media. Q. Appl. Math. 24, (1966) 5. Berger, H., Gabbert, U., Koeppe, H., Rodriguez-Ramos, R., Bravo-Castillero, J., Guinovart-Diaz, R., Otero, J.A., Maugin, G.A.: Finite element and asymptotic homogenization methods applied to smart composite materials. Comput. Mech. 33, (23) 6. Berger, H., Kari, S., Gabbert, U., Rodriguez-Ramos, R., Guinovart-Diaz, R., Otero, J.A., Bravo-Castillero, J.: An analytical and numerical approach for calculating effective material coefficients of piezoelectric fiber composites. Int. J. Solids Struct. 42, (25a) 7. Berger, H., Kari, S., Gabbert, U., Rodriguez-Ramos, R., Bravo-Castillero, J., Guinovart-Diaz, R.: A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fiber composites. Mater. Sci. Eng. A 412, 53 6 (25b) 8. Berger, H., Kari, S., Gabbert, U., Rodriguez-Ramos, R., Bravo-Castillero, J., Guinovart-Diaz, R.: Calculation of effective coefficients for piezoelectric fiber composites based on a general numerical homogenization technique. Compos. Struct. 71, (25c) 9. Berger, H., Kari, S., Gabbert, U., Rodriguez-Ramos, R., Bravo-Castillero, J., Guinovart-Diaz, R., Sabina, F.J., Maugin, G.A.: Unit cell models of piezoelectric fibre composites for numerical and analytical calculation of effective properties. Smart Mater. Struct. 15, (26) 1. Guinovart-Díaz, R., Bravo-Castillero, J., Rodríguez-Ramos, R., Sabina, F.J.: Closed-form expressions for the effective coefficients of fiber-reinforced composite with transversely isotropic constituents. I: Elastic and hexagonal symmetry. J. Mech. Phys. Solids 49, (21) 11. Guinovart-Díaz, R., Rodríguez-Ramos, R., Bravo-Castillero, J., Sabina, F.J., Maugin, G.A.: Closed-form thermo-elastic moduli of a periodic three-phase fiber-reinforced composite. J. Therm. Stresses 28, (25) 12. Hashin, Z., Shtrikman, S.: A variational approach to the theory of the elastic behaviour of multiphase materials. J. Mech. Phys. Solids 11, (1963) 13. Hashin, Z.: On elastic behaviour of fibre reinforced materials of arbitrary transverse phase geometry. J. Mech. Phys. Solids 13, (1965) 14. Hashin, Z.: Analysis of properties of fiber composites with anisotropic constituents. J. Appl. Mech. 46, (1979) 15. Hill, R.: Theory of mechanical properties of fibre-strengthened materials: I. Elast. Behav. J. Mech. Phys. Solids 12, (1964) 16. Hill, R.: A self-consistent mechanics of composite materials. J. Mech. Phys. Solids 13, (1965a) 17. Hill, R.: Continuum micromechanics of elastoplastic polycrystals. J. Mech. Phys. Solids 13, (1965b) 18. Ju, J.W., Chen, T.: Effective elastic moduli of two phase composites containing randomly dispersed spherical inhomogeneities. Act. Mech. 13, (1994) 19. Kari, S., Berger, H., Rodriguez-Ramos, R., Gabbert, U.: Computational evaluation of effective material properties of composites reinforced by randomly distributed spherical particles. Comp. Struct. 77, (27a) 2. Kari, S., Berger, H., Gabbert, U.: Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites. Comp. Mater. Sci. 39, (27b)

14 R. Rodríguez-Ramos et al. 21. Levin, V.M., Rakovskaja, Krecher, W.S.: The effective thermoelectroelastic properties of microinhomogeneous materials. Int. J. Solids Struct. 36, (1999) 22. McLaughlin, R.: A study of the differential scheme for composite materials. Int. J. Eng. Sci. 15, (1977) 24. Milton, G.W.: Bounds on the elastic and transport properties of two-component composites. J. Mech. Phys. Solids 3, (1982) 24. Milton, G.W., Phan-Thien, N.: New bounds on the effective moduli of two-component materials. Proc. Roy. Soc. Lond. A 38, (1982) 25. Mori, T., Tanaka, K.: Average stress in the matrix and average elastic energy of materials with misfitting inclusions. Acta Metall. Mater. 21, (1973) 26. Nemat-Nasser, S., Hori, M.: Universal bounds for overall properties of linear and nonlinear heterogeneous solids. J. Eng. Mat. Tech. 117, (1995) 27. Nemat-Nasser, S., Hori, M.: Micromechanics: overall properties of heterogeneous materials, second ed. North Holland, Amsterdam (1999) 28. Nie, S., Basaran, C.: A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds. Int. J. Solids Struct. 42, (25) 29. Pobedrya, B.E.: Mechanics of composite materials. Moscow State University Press, Moscow (in Russian) (1984) 3. Rodríguez-Ramos, R., Sabina, F.J., Guinovart-Díaz, R., Bravo-Castillero, J.: Closed-form expressions for the effective coefficients of fibre-reinforced composite with transversely isotropic constituents. I: Elastic and square symmetry. Mech. Mater. 33, (21) 31. Sabina, F.J., Bravo-Castillero, J., Guinovart-Díaz, R., Rodríguez-Ramos, R., Valdiviezo-Mijangos, O.: Overall behavior of two-dimensional periodic composites. Int. J. Solids Struct. 39, (22) 32. Sabina, F.J., Smyshlyaev, V.P., Willis, J.R.: Self-consistent analysis of waves in a matrix-inclusion composite-i. Aligned spheroidal inclusions. J. Mech. Phys. Solids 41(1), (1993) 33. Sabina, F.J., Willis, J.R.: Self-consistent analysis of waves in rocks containing arrays of cracks. In: Fjaer, E., Holt, R.M., Rathore, J.S. (eds.) Seismic Anisotropy Tulsa. Society of Exploration Geophysicists, Oklahoma, pp (1996) 34. Smith, S.C.: Experimental values for the elastic constants of a particulate-filled glassy polymer. J. Res. NBS 8, (1976) 35. Torquato, S.: Random heterogeneous media: microstructure and improved bounds on effective properties. Appl. Mech. Rev. 44, (1991) 36. Willis, J.R.: Bounds and self-consistent estimates for the overall properties of anisotropic composites. J. Mech. Phys. Solids 25, (1977) 37. Willis, J.R.: A polarization approach to the scattering of elastic waves-i. Scattering by a single inclusion. J. Mech. Phys. Solids 28, (198a) 38. Willis, J.R.: A polarization approach to the scattering of elastic waves-ii. Multiple scattering from inclusions. J. Mech. Phys. Solids 28, (198b) 39. Willis, J.R.: Variational and related method for the properties of composites. Adv. Appl. Mech. 21, 1 78 (1981) 4. Willis, J.R.: Bounds for the overall properties of anisotropic composites. In: Mura, T., Koya, T. (eds.) Variational methods in mechanics, chap. 21. Oxford University Press, Oxford (1992)

Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites

Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites Computational Materials Science 39 (2007) 198 204 www.elsevier.com/locate/commatsci Numerical evaluation of effective material properties of randomly distributed short cylindrical fibre composites S. Kari

More information

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method

Module 7: Micromechanics Lecture 34: Self Consistent, Mori -Tanaka and Halpin -Tsai Models. Introduction. The Lecture Contains. Self Consistent Method Introduction In this lecture we will introduce some more micromechanical methods to predict the effective properties of the composite. Here we will introduce expressions for the effective properties without

More information

International Journal of Engineering Science

International Journal of Engineering Science International Journal of Engineering Science 63 (2013) 1 9 Contents lists available at SciVerse ScienceDirect International Journal of Engineering Science journal homepage: www.elsevier.com/locate/ijengsci

More information

A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites

A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites Materials Science and Engineering A 412 (2005) 53 60 A comprehensive numerical homogenisation technique for calculating effective coefficients of uniaxial piezoelectric fibre composites Harald Berger a,,

More information

PLEASE SCROLL DOWN FOR ARTICLE

PLEASE SCROLL DOWN FOR ARTICLE This article was downloaded by:[otto-von-guericke-universitaet] On: 0 September 007 Access Details: [subscription number 77067] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered

More information

Finite Element of Composite Materials with Partial Debonding along Boundary. Kohei HIGUCHI, Tetsuo IWAKUMA

Finite Element of Composite Materials with Partial Debonding along Boundary. Kohei HIGUCHI, Tetsuo IWAKUMA Vol. 8 (2005 8 ) Finite Element of Composite aterials with Partial Debonding along Boundary Kohei HIGUCHI, Tetsuo IWAKUA JFE PhD 980-8579 6-6-06 Damages in composite materials are mostly governed by interfacial

More information

International Journal of Engineering Science

International Journal of Engineering Science International Journal of Engineering Science 49 (2011) 322 332 Contents lists available at ScienceDirect International Journal of Engineering Science ournal homepage: www.elsevier.com/locate/iengsci On

More information

Homogenization methods for multi-phase elastic composites: Comparisons and benchmarks

Homogenization methods for multi-phase elastic composites: Comparisons and benchmarks TECHNISCHE MECHANIK, 3, 4, (), 374 386 submitted: October 6, 9 Homogenization methods for multi-phase elastic composites: Comparisons and benchmarks B. Klusemann, B. Svendsen Usually homogenization methods

More information

Analysis of high loss viscoelastic composites

Analysis of high loss viscoelastic composites Analysis of high loss viscoelastic composites by C. P. Chen, Ph.D. and R. S. Lakes, Ph.D. Department of Engineering Physics Engineering Mechanics Program; Biomedical Engineering Department Materials Science

More information

EVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD

EVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD U.P.B. Sci. Bull., Series D, Vol. 79, Iss. 3, 2017 ISSN 1454-2358 EVALUATION OF EFFECTIVE MECHANICAL PROPERTIES OF COMPLEX MULTIPHASE MATERIALS WITH FINITE ELEMENT METHOD Mohamed said BOUTAANI 1, Salah

More information

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains

Module 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model

More information

RVE Analysis of Light Weight Carbon Nanotubes Embedded Piezoelectric Fibre Composites

RVE Analysis of Light Weight Carbon Nanotubes Embedded Piezoelectric Fibre Composites Nanoscience and Nanotechnology 2016, 6(1): 11-16 DOI: 10.5923/j.nn.20160601.03 RVE Analysis of Light Weight Carbon Nanotubes Embedded Piezoelectric Fibre Composites V. K. Srivastava 1,*, H. Berger 2, U.

More information

Interfacial effects in electromagnetic coupling within piezoelectric phononic crystals

Interfacial effects in electromagnetic coupling within piezoelectric phononic crystals Acta Mech Sin (29) 25:95 99 DOI 1.17/s149-8-21-y RESEARCH PAPER Interfacial effects in electromagnetic coupling within pieoelectric phononic crystals F. J. Sabina A. B. Movchan Received: 14 July 28 / Accepted:

More information

McMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007

McMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007 McMAT 2007 Micromechanics of Materials Austin, Texas, June 3 7, 2007 RANDOM POLYCRYSTALS OF GRAINS WITH CRACKS: MODEL OF ELASTIC BEHAVIOR FOR FRACTURED SYSTEMS James G. Berryman Earth Sciences Division

More information

Variational Asymptotic Method for Unit Cell Homogenization

Variational Asymptotic Method for Unit Cell Homogenization Variational Asymptotic Method for Unit Cell Homogenization Wenbin Yu 1 and Tian Tang 2 1 Utah State University, Logan, Utah 84322-4130 Wenbin.Yu@usu.edu 2 Utah State University, Logan, Utah 84322-4130

More information

A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds

A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds International Journal of Solids and Structures 42 (2005) 4179 4191 www.elsevier.com/locate/ijsolstr A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial

More information

Finite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements

Finite Element Modeling of Residual Thermal Stresses in Fiber-Reinforced Composites Using Different Representative Volume Elements Proceedings o the World Congress on Engineering 21 Vol II WCE 21, June 3 - July 2, 21, London, U.K. Finite Element Modeling o Residual Thermal Stresses in Fiber-Reinorced Composites Using Dierent Representative

More information

Authors: Correspondence: ABSTRACT:

Authors: Correspondence: ABSTRACT: Multi-Scale Modeling of the Impact and Failure of Fiber Reinforced Polymer Structures using DIGIMAT to LS-DYNA Interface Authors: L. Adam, R. Assaker & R. Ramaya e-xstream engineering S.A. 7, Rue du Bosquet.

More information

Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading

Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading Indian Journal of Engineering & Materials Sciences Vol. 15, October 2008, pp. 382-390 Micromechanical analysis of FRP hybrid composite lamina for in-plane transverse loading K Sivaji Babu a *, K Mohana

More information

EFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL

EFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL AMS Subject Classification Index: 74Q20 EFFECTIVE CHARACTERISTICS OF POROUS MEDIA AS A FUNCTION OF POROSITY LEVEL Irini DJERAN-MAIGRE and Sergey V.KUZNETSOV INSA de Lyon, URGC 34 Avenue des Arts, 69621

More information

Fundamentals of Linear Elasticity

Fundamentals of Linear Elasticity Fundamentals of Linear Elasticity Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI bluebox.ippt.pan.pl/ tzielins/ Institute of Fundamental Technological Research of the Polish Academy

More information

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach

Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Mechanical and Thermal Properties of Coir Fiber Reinforced Epoxy Composites Using a Micromechanical Approach Sandhyarani Biswas Department of Mechanical Engineering, N.I.T Rourkela, INDIA Abstract: Now-a-days,

More information

Department of Mechanical Engineering, National Institute of Technology Rourkela, Rourkela , India

Department of Mechanical Engineering, National Institute of Technology Rourkela, Rourkela , India Composites, Article ID 629175, 12 pages http://dx.doi.org/10.1155/2014/629175 Research Article Effect of Fiber Geometry and Representative Volume Element on Elastic and Thermal Properties of Unidirectional

More information

Fig. 1. Circular fiber and interphase between the fiber and the matrix.

Fig. 1. Circular fiber and interphase between the fiber and the matrix. Finite element unit cell model based on ABAQUS for fiber reinforced composites Tian Tang Composites Manufacturing & Simulation Center, Purdue University West Lafayette, IN 47906 1. Problem Statement In

More information

A simple plane-strain solution for functionally graded multilayered isotropic cylinders

A simple plane-strain solution for functionally graded multilayered isotropic cylinders Structural Engineering and Mechanics, Vol. 24, o. 6 (2006) 000-000 1 A simple plane-strain solution for functionally graded multilayered isotropic cylinders E. Pan Department of Civil Engineering, The

More information

Three-dimensional thermoelastic deformations of a functionally graded elliptic plate

Three-dimensional thermoelastic deformations of a functionally graded elliptic plate JCOMB Composites: Part B () 9 6 www.elsevier.com/locate/compositesb Three-dimensional thermoelastic deformations of a functionally graded elliptic plate Z.-Q. Cheng a, R.C. Batra b, * a Department of Modern

More information

Soufiane Belhouideg. Keywords Composite materials, homogenization method, Eshelby s inclusion problem.

Soufiane Belhouideg. Keywords Composite materials, homogenization method, Eshelby s inclusion problem. Prediction of effective mechanical properties of composite materials using homogenization approach: Application to tungsten fiber reinforced bulk metallic glass matrix composite Soufiane Belhouideg Abstract

More information

Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites

Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites International OPEN ACCESS Journal Of Modern Engineering Research (IJMER) Continuum Modeling Techniques to Determine Mechanical Properties of Nanocomposites Sonali Gholap 1, Dr. Dhananjay R. Panchagade

More information

Hysteretic equivalent continuum model of nanocomposites with interfacial stick-slip

Hysteretic equivalent continuum model of nanocomposites with interfacial stick-slip Hysteretic equivalent continuum model of nanocomposites with interfacial stick-slip Giovanni Formica 1, Michela Talò, Walter Lacarbonara 1 Dipartimento di Architettura, Università degli Studi Roma Tre

More information

Statistical continuum mechanics analysis of an elastic two-isotropic-phase composite material

Statistical continuum mechanics analysis of an elastic two-isotropic-phase composite material Composites: Part B 31 (2000) 39 46 www.elsevier.com/locate/compositesb Statistical continuum mechanics analysis of an elastic two-isotropic-phase composite material S. Lin, H. Garmestani* Mechanical Engineering

More information

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING

MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING MICROMECHANICAL ANALYSIS OF FRP COMPOSITES SUBJECTED TO LONGITUDINAL LOADING N. Krishna Vihari 1, P. Phani Prasanthi 1, V. Bala Krishna Murthy 2* and A. Srihari Prasad 3 1 Mech. Engg. Dept., P. V. P. Siddhartha

More information

FINITE ELEMENT ANALYSIS APPLIED TO EVALUATION OF EFFECTIVE MATERIAL COEFFICIENTS FOR PIEZOELECTRIC FIBER COMPOSITES

FINITE ELEMENT ANALYSIS APPLIED TO EVALUATION OF EFFECTIVE MATERIAL COEFFICIENTS FOR PIEZOELECTRIC FIBER COMPOSITES 29 Brazilian Symposium on Aerospace Eng. & Applications Copyright 29 by AAB September 14-16, 29, S. J. Campos, SP, Brazil FINITE ELEMENT ANALYSIS APPLIED TO EALUATION OF EFFECTIE MATERIAL COEFFICIENTS

More information

EFFECTIVE PROPERTIES OF THREE-PHASE ELECTRO-MAGNETO-ELASTIC MULTIFUNCTIONAL COMPOSITE MATERIALS. A Thesis JAE SANG LEE

EFFECTIVE PROPERTIES OF THREE-PHASE ELECTRO-MAGNETO-ELASTIC MULTIFUNCTIONAL COMPOSITE MATERIALS. A Thesis JAE SANG LEE EFFECTIVE PROPERTIES OF THREE-PHASE ELECTRO-MAGNETO-ELASTIC MULTIFUNCTIONAL COMPOSITE MATERIALS A Thesis by JAE SANG LEE Submitted to the Office of Graduate Studies of Texas A&M University in partial fulfillment

More information

Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics

Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Predicting Elastic Properties of Unidirectional SU8/ZnO Nanocomposites using COMSOL Multiphysics Neelam Mishra 1, and Kaushik Das *1 1 School of Minerals Metallurgical and Materials Engineering, Indian

More information

Prediction of Elastic Constants on 3D Four-directional Braided

Prediction of Elastic Constants on 3D Four-directional Braided Prediction of Elastic Constants on 3D Four-directional Braided Composites Prediction of Elastic Constants on 3D Four-directional Braided Composites Liang Dao Zhou 1,2,* and Zhuo Zhuang 1 1 School of Aerospace,

More information

Stiffness Predictions for Unidirectional Short-Fiber Composites: Review and Evaluation

Stiffness Predictions for Unidirectional Short-Fiber Composites: Review and Evaluation Stiffness Predictions for Unidirectional Short-Fiber Composites: Review and Evaluation Charles L. Tucker III Department of Mechanical and Industrial Engineering University of Illinois at Urbana-Champaign

More information

A multiscale approach for composite materials as multifield continua

A multiscale approach for composite materials as multifield continua Materials Science Forum Vols. 539-543 (27) pp. 2551-2556 online at http://www.scientific.net (27) Trans Tech Publications, Switzerland A multiscale approach for composite materials as multifield continua

More information

Supplementary figures

Supplementary figures Supplementary figures Supplementary figure 1: The fang model (longitudinal section) used for the stiffness analysis. The base of the fang is considered fixed (i.e. zero displacement constraint) and the

More information

ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS

ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS ASSESSMENT OF MIXED UNIFORM BOUNDARY CONDITIONS FOR PREDICTING THE MACROSCOPIC MECHANICAL BEHAVIOR OF COMPOSITE MATERIALS Dieter H. Pahr and Helmut J. Böhm Institute of Lightweight Design and Structural

More information

P program, version 2. Theoretical manual & User guide

P program, version 2. Theoretical manual & User guide Documentation for H E L P program, version 2 Theoretical manual & User guide Jan Vorel 1 May 22, 2015 1 Czech Technical University in Prague, Faculty of Civil Engineering, Department of Mechanics Contents

More information

A Survey of Neural Computing of Effective Properties of Random Composites

A Survey of Neural Computing of Effective Properties of Random Composites A Survey of Neural Computing of Effective Properties of Random Composites ABRAHAM I. BELTZER, VICTOR A. GOTLIB Holon Academic Institute of Technology POB 305, Holon 580 ISRAEL TADANOBU SATO Disaster Prevention

More information

CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles

CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles CHEM-C2410: Materials Science from Microstructures to Properties Composites: basic principles Mark Hughes 14 th March 2017 Today s learning outcomes To understand the role of reinforcement, matrix and

More information

Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses

Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses PROC. ITB Eng. Science Vol. 38 B, No. 2, 2006, 159-170 159 Effects of Fracture Parameters in an Anisotropy Model on P-Wave Azimuthal Amplitude Responses Fatkhan Program Studi Teknik Geofisika FIKTM-ITB

More information

Accepted Manuscript. Numerical simulation of the guided Lamb wave propagation in particle reinforced

Accepted Manuscript. Numerical simulation of the guided Lamb wave propagation in particle reinforced Accepted Manuscript Numerical simulation of the guided Lamb wave propagation in particle reinforced composites Ralf Weber, Seyed Mohammad Hossein Hosseini, Ulrich Gabbert PII: S0263-8223(12)00180-8 DOI:

More information

Fatigue Damage Development in a Steel Based MMC

Fatigue Damage Development in a Steel Based MMC Fatigue Damage Development in a Steel Based MMC V. Tvergaard 1,T.O/ rts Pedersen 1 Abstract: The development of fatigue damage in a toolsteel metal matrix discontinuously reinforced with TiC particulates

More information

DESIGN OF LAMINATES FOR IN-PLANE LOADING

DESIGN OF LAMINATES FOR IN-PLANE LOADING DESIGN OF LAMINATES FOR IN-PLANOADING G. VERCHERY ISMANS 44 avenue F.A. Bartholdi, 72000 Le Mans, France Georges.Verchery@m4x.org SUMMARY This work relates to the design of laminated structures primarily

More information

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations

16.21 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive Equations 6.2 Techniques of Structural Analysis and Design Spring 2003 Unit #5 - Constitutive quations Constitutive quations For elastic materials: If the relation is linear: Û σ ij = σ ij (ɛ) = ρ () ɛ ij σ ij =

More information

CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement

CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement CHEM-E2200: Polymer blends and composites Fibre architecture and principles of reinforcement Mark Hughes 19 th September 2016 Outline Fibre architecture Volume fraction and the rule of mixtures Principle

More information

Combined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites

Combined Isotropic-Kinematic Hardening Laws with Anisotropic Back-stress Evolution for Orthotropic Fiber-Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution for Orthotropic Fiber- Reinforced Composites Combined Isotropic-Kinematic Hardening Laws with Antropic Back-stress Evolution

More information

PHASE VELOCITY AND ATTENUATION OF SH WAVES IN A FIBER-

PHASE VELOCITY AND ATTENUATION OF SH WAVES IN A FIBER- PHASE VELOCITY AND ATTENUATION OF SH WAVES IN A FIBER- REINFORCED COMPOSITE Ruey-Bin Yang and Ajit K. Mal Department of Mechanical, Aerospace and Nuclear Engineering University of California, Los Angeles,

More information

International Journal of Solids and Structures

International Journal of Solids and Structures Accepted Manuscript Micromechanical method for effective piezoelectric properties and electromechanical fields in multi-coated long fiber composites George Chatzigeorgiou, Ali Javili, Fodil Meraghni PII:

More information

Fracture Mechanics of Composites with Residual Thermal Stresses

Fracture Mechanics of Composites with Residual Thermal Stresses J. A. Nairn Material Science & Engineering, University of Utah, Salt Lake City, Utah 84 Fracture Mechanics of Composites with Residual Thermal Stresses The problem of calculating the energy release rate

More information

Porous Fibre-Reinforced Materials Under Large Deformations

Porous Fibre-Reinforced Materials Under Large Deformations Porous Fibre-Reinforced Materials Under Large Deformations Salvatore Federico Department of Mechanical and Manufacturing Engineering, The University of Calgary 2500 University Drive NW, Calgary, Alberta,

More information

Development of a code to generate randomly distributed short fiber composites to estimate mechanical properties using FEM

Development of a code to generate randomly distributed short fiber composites to estimate mechanical properties using FEM International Journal of Theoretical and Applied Mechanics. ISSN 0973-6085 Volume 12, Number 4 (2017) pp. 863-872 Research India Publications http://www.ripublication.com Development of a code to generate

More information

STEEL FIBER CURVATURE IN CONCRETE COMPOSITES: MODULUS PREDICTIONS USING EFFECTIVE STEEL FIBER PROPERTIES

STEEL FIBER CURVATURE IN CONCRETE COMPOSITES: MODULUS PREDICTIONS USING EFFECTIVE STEEL FIBER PROPERTIES American Journal of Applied Sciences 11 (1): 145-151, 2014 ISSN: 1546-9239 2014 Science Publication doi:10.3844/ajassp.2014.145.151 Published Online 11 (1) 2014 (http://www.thescipub.com/ajas.toc) STEEL

More information

Bounding of effective thermal conductivities of multiscale materials by essential and natural boundary conditions

Bounding of effective thermal conductivities of multiscale materials by essential and natural boundary conditions PHYSICAL REVIEW B VOLUME 54, NUMBER 1 1 JULY 1996-I Bounding of effective thermal conductivities of multiscale materials by essential and natural boundary conditions M. Ostoja-Starzewski Institute of Paper

More information

A direct evaluation of the Fabric Tensor in anisotropic porous media

A direct evaluation of the Fabric Tensor in anisotropic porous media A direct evaluation of the Fabric Tensor in anisotropic porous media Maria Cristina Pernice 1, Luciano Nunziante 1, Massimiliano Fraldi 1,2 1 Department of Structural Engineering, University of Naples

More information

Evaluation of Effective Elastic Moduli Using Micromechanics

Evaluation of Effective Elastic Moduli Using Micromechanics IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Evaluation of Effective Elastic Moduli Using Micromechanics To cite this article: Misba Mehdi Mrs. et al 2018 IOP Conf. Ser.:

More information

Viscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites

Viscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites Viscoelastic Damping Characteristics of Indium-Tin/SiC Particulate Composites HyungJoo Kim, Colby C. Swan Ctr. for Computer-Aided Design, Univ. of Iowa Roderic Lakes Engineering Physics, Univ. of Wisconsin

More information

Received 5 April 2003; received in revised form 29 May 2003

Received 5 April 2003; received in revised form 29 May 2003 Materials Science and Engineering A360 (2003) 339/344 www.elsevier.com/locate/msea Connection between elastic moduli and thermal conductivities of anisotropic short fiber reinforced thermoplastics: theory

More information

DEVELOPMENT OF A MICROMECHANICAL MODEL IN INTERACTION WITH PARAMETERS RELATED TO THE MICROSTRUCTURE OF CARBON/EPOXY COMPOSITES

DEVELOPMENT OF A MICROMECHANICAL MODEL IN INTERACTION WITH PARAMETERS RELATED TO THE MICROSTRUCTURE OF CARBON/EPOXY COMPOSITES ECCM16-16 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, 22-26 June 2014 DEVELOPMENT OF A MICROMECHANICAL MODEL IN INTERACTION WITH PARAMETERS RELATED TO THE MICROSTRUCTURE OF CARBON/EPOXY

More information

ME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites

ME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-

More information

(This is a sample cover image for this issue. The actual cover is not yet available at this time.)

(This is a sample cover image for this issue. The actual cover is not yet available at this time.) (This is a sample cover image for this issue. The actual cover is not yet available at this time. This article appeared in a journal published by Elsevier. The attached copy is furnished to the author

More information

3D XFEM Modelling of Imperfect Interfaces

3D XFEM Modelling of Imperfect Interfaces 3D XFEM Modelling of Imperfect Interfaces Elena Benvenuti 1, Giulio Ventura 2, Antonio Tralli 1, Nicola Ponara 1 1 Dipartimento di Ingegneria, Universit di Ferrara, Italy E-mail: elena.benvenuti@unife.it,

More information

A Comparative Study of Analytical and Numerical Evaluation of Elastic Properties of Short Fiber Composites

A Comparative Study of Analytical and Numerical Evaluation of Elastic Properties of Short Fiber Composites IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS A Comparative Study of Analytical and Numerical Evaluation of Elastic Properties of Short Fiber Composites To cite this article:

More information

USING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE

USING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE Volume II: Fatigue, Fracture and Ceramic Matrix Composites USING A HOMOGENIZATION PROCEDURE FOR PREDICTION OF MATERIAL PROPERTIES AND THE IMPACT RESPONSE OF UNIDIRECTIONAL COMPOSITE A. D. Resnyansky and

More information

Microstructure design of a two phase composite using two-point correlation functions

Microstructure design of a two phase composite using two-point correlation functions Brigham Young University BYU ScholarsArchive All Faculty Publications 2004-01-01 Microstructure design of a two phase composite using two-point correlation functions Brent L. Adams b_l_adams@byu.edu H.

More information

in materials/structural modellings and designs Liping Liu

in materials/structural modellings and designs Liping Liu Closed-form solutions and their applications in materials/structural modellings and designs Liping Liu University it of Houston Why composites and heterogeneous structures? A material has many different

More information

AN EFFECTIVE SOLUTION OF THE COMPOSITE (FGM S) BEAM STRUCTURES

AN EFFECTIVE SOLUTION OF THE COMPOSITE (FGM S) BEAM STRUCTURES Engineering MECHANICS, Vol. 15, 2008, No. 2, p. 115 132 115 AN EFFECTIVE SOLUTION OF THE COMPOSITE (FGM S) BEAM STRUCTURES Justín Murín, Vladimír Kutiš* The additive mixture rules have been extended for

More information

A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites.

A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites. A coupled field finite element model to predict actuation properties of piezoelectrically actuated bistable composites. P.F.Giddings, C.R.Bowen, H.A.Kim University of Bath, UK Dept. Mech. ng, University

More information

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f

MATERIAL MECHANICS, SE2126 COMPUTER LAB 4 MICRO MECHANICS. E E v E E E E E v E E + + = m f f. f f MATRIAL MCHANICS, S226 COMPUTR LAB 4 MICRO MCHANICS 2 2 2 f m f f m T m f m f f m v v + + = + PART A SPHRICAL PARTICL INCLUSION Consider a solid granular material, a so called particle composite, shown

More information

Aspect Ratio Requirements for Nanotube-Reinforced, Polymer-Matrix Composites

Aspect Ratio Requirements for Nanotube-Reinforced, Polymer-Matrix Composites Aspect Ratio Requirements for Nanotube-Reinforced, Polymer-Matrix Composites J.A.Nairn Wood Science and Engineering, Oregon State University, Corvallis, OR 97330, USA Abstract A fiber s efficiency in a

More information

Moving screw dislocations in piezoelectric bimaterials

Moving screw dislocations in piezoelectric bimaterials phys stat sol (b) 38 No 1 10 16 (003) / DOI 10100/pssb00301805 Moving screw dislocations in piezoelectric bimaterials Xiang-Fa Wu *1 Yuris A Dzenis 1 and Wen-Sheng Zou 1 Department of Engineering Mechanics

More information

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density

The Rotating Inhomogeneous Elastic Cylinders of. Variable-Thickness and Density Applied Mathematics & Information Sciences 23 2008, 237 257 An International Journal c 2008 Dixie W Publishing Corporation, U. S. A. The Rotating Inhomogeneous Elastic Cylinders of Variable-Thickness and

More information

Effects of geometry and properties of fibre and matrix on overall. composite parameters

Effects of geometry and properties of fibre and matrix on overall. composite parameters Int. Journal of Applied Sciences and ngineering Research, ol. 3, Issue 4, 2014 www.ijaser.com 2014 by the authors Licensee IJASR- Under Creative Commons License 3.0 editorial@ijaser.com Research article

More information

Interaction of Two Parallel Short Fibers in the Matrix at Loss of Stability

Interaction of Two Parallel Short Fibers in the Matrix at Loss of Stability Copyright c 2006 Tech Science Press CMES, vol.13, no.3, pp.165-169, 2006 Interaction of Two Parallel Short Fibers in the Matrix at Loss of Stability A.N.Guz and V.A.Dekret 1 Abstract: Stability problem

More information

Quadratic and cubic monocrystalline and polycrystalline materials: their stability and mechanical properties

Quadratic and cubic monocrystalline and polycrystalline materials: their stability and mechanical properties Journal of Physics: Conference Series Quadratic and cubic monocrystalline and polycrystalline materials: their stability and mechanical properties To cite this article: C Jasiukiewicz et al 010 J. Phys.:

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD

DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD ECCM6-6 TH EUROPEAN CONFERENCE ON COMPOSITE MATERIALS, Seville, Spain, -6 June 4 DYNAMIC RESPONSE OF SYNTACTIC FOAM CORE SANDWICH USING A MULTIPLE SCALES BASED ASYMPTOTIC METHOD K. V. Nagendra Gopal a*,

More information

Computational Analysis for Composites

Computational Analysis for Composites Computational Analysis for Composites Professor Johann Sienz and Dr. Tony Murmu Swansea University July, 011 The topics covered include: OUTLINE Overview of composites and their applications Micromechanics

More information

MICROMECHANICS AND HOMOGENIZATION. Materials Containing Coated and Uncoated Spherical Inhomogeneities

MICROMECHANICS AND HOMOGENIZATION. Materials Containing Coated and Uncoated Spherical Inhomogeneities MICROMECHANICS AND HOMOGENIZATION Materials Containing Coated and Uncoated Spherical Inhomogeneities Václav Nežerka, 2012 i Acknowledgement I would like to thank my supervisor, Jan Zeman, who supported

More information

Back Analysis of Measured Displacements of Tunnels

Back Analysis of Measured Displacements of Tunnels Rock Mechanics and Rock Engineering 16, 173--180 (1983) Rock Mechanics and Rock Engineering 9 by Springer-Verlag 1983 Back Analysis of Measured Displacements of Tunnels By S. Sakurai and K. Takeuchi Kobe

More information

Effect of fibre shape on transverse thermal conductivity of unidirectional composites

Effect of fibre shape on transverse thermal conductivity of unidirectional composites Sādhanā Vol. 4, Part 2, April 25, pp. 53 53. c Indian Academy of Sciences Effect of fibre shape on transverse thermal conductivity of unidirectional composites B RAGHAVA RAO,, V RAMACHANDRA RAJU 2 and

More information

MODELLING THE THERMOELASTIC PROPERTIES OF SHORT FIBRE COMPOSITES WITH ANISOTROPIC PHASES

MODELLING THE THERMOELASTIC PROPERTIES OF SHORT FIBRE COMPOSITES WITH ANISOTROPIC PHASES MODELLING THE THERMOELASTIC PROPERTIES OF SHORT FIBRE COMPOSITES WITH ANISOTROPIC PHASES P.J.Hine*, C.D.Price*, B.Whiteside $, A.M.Cunha # and I.M.Ward* * - IRC in Polymer Science and Technology, University

More information

Interfacial partial debonding and its in uence on the elasticity of a two-phase composite

Interfacial partial debonding and its in uence on the elasticity of a two-phase composite Mechanics of Materials 32 (2000) 695±709 www.elsevier.com/locate/mechmat Interfacial partial debonding and its in uence on the elasticity of a two-phase composite S.F. Zheng, M. Denda, G.J. Weng * Department

More information

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1

3 2 6 Solve the initial value problem u ( t) 3. a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 Math Problem a- If A has eigenvalues λ =, λ = 1 and corresponding eigenvectors 1 3 6 Solve the initial value problem u ( t) = Au( t) with u (0) =. 3 1 u 1 =, u 1 3 = b- True or false and why 1. if A is

More information

Macroscopic theory Rock as 'elastic continuum'

Macroscopic theory Rock as 'elastic continuum' Elasticity and Seismic Waves Macroscopic theory Rock as 'elastic continuum' Elastic body is deformed in response to stress Two types of deformation: Change in volume and shape Equations of motion Wave

More information

Rayleigh waves in magneto-electro-elastic half planes

Rayleigh waves in magneto-electro-elastic half planes Acta Mech 202, 127 134 (2009) DOI 10.1007/s00707-008-0024-8 W. J. Feng E. Pan X. Wang J. Jin Rayleigh waves in magneto-electro-elastic half planes Received: 23 August 2007 / Revised: 17 February 2008 /

More information

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE

A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE A FINITE ELEMENT MODEL TO PREDICT MULTI- AXIAL STRESS-STRAIN RESPONSE OF CERAMIC MATRIX COMPOSITES WITH STRAIN INDUCED DAMAGE Daxu Zhang and D. R. Hayhurst School of Mechanical, Aerospace and Civil Engineering,

More information

Archetype-Blending Multiscale Continuum Method

Archetype-Blending Multiscale Continuum Method Archetype-Blending Multiscale Continuum Method John A. Moore Professor Wing Kam Liu Northwestern University Mechanical Engineering 3/27/2014 1 1 Outline Background and Motivation Archetype-Blending Continuum

More information

Exploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description. Igor Dobovšek

Exploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description. Igor Dobovšek Exploring Piezoelectric Properties of Wood and Related Issues in Mathematical Description Igor Dobovšek University of Ljubljana Faculty of Mathematics and Physics Institute of Mathematics Physics and Mechanics

More information

CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES

CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES CHAPTER 4 MODELING OF MECHANICAL PROPERTIES OF POLYMER COMPOSITES 4. Introduction Fillers added to polymer matrices as additives are generally intended for decreasing the cost (by increase in bulk) of

More information

Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate

Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate Buckling Behavior of 3D Randomly Oriented CNT Reinforced Nanocomposite Plate Outline Introduction Representative Volume Element (RVE) Periodic Boundary Conditions on RVE Homogenization Method Analytical

More information

Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996).

Cellular solid structures with unbounded thermal expansion. Roderic Lakes. Journal of Materials Science Letters, 15, (1996). 1 Cellular solid structures with unbounded thermal expansion Roderic Lakes Journal of Materials Science Letters, 15, 475-477 (1996). Abstract Material microstructures are presented which can exhibit coefficients

More information

3D and Planar Constitutive Relations

3D and Planar Constitutive Relations 3D and Planar Constitutive Relations A School on Mechanics of Fibre Reinforced Polymer Composites Knowledge Incubation for TEQIP Indian Institute of Technology Kanpur PM Mohite Department of Aerospace

More information

LA-UR-01-3491 Approved for public release; distribution is unlimited. Title: APPLYING MICRO-MECHANICS TO FINITE ELEMENT SIMULATIONS OF SPLIT HOPKINSON PRESSURE BAR EXPERIMENTS ON HIGH EXPLOSIVES Author(s):

More information

P ticles of one material in a polymer matrix of another

P ticles of one material in a polymer matrix of another A Gene ra I i ze d Se If - C o n s i s t e n t M or i -Ta na ka Sc h e m e for Prediction of the Effective Elastic Moduli of Hybrid Multiphase Particulate Composites L. H. DAI, 2. P. HUANG, andr. WANG

More information

Exercise: concepts from chapter 8

Exercise: concepts from chapter 8 Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic

More information

Methodology for the evaluation of yield strength and hardening behavior of metallic materials by indentation with spherical tip

Methodology for the evaluation of yield strength and hardening behavior of metallic materials by indentation with spherical tip JOURNAL OF APPLIED PHYSICS VOLUME 94, NUMBER 1 1 JULY 2003 Methodology for the evaluation of yield strength and hardening behavior of metallic materials by indentation with spherical tip Dejun Ma Department

More information

ANALYTICAL FRACTURE MECHANICS ANALYSIS OF THE PULL-OUT TEST INCLUDING THE EFFECTS OF FRICTION AND THERMAL STRESSES

ANALYTICAL FRACTURE MECHANICS ANALYSIS OF THE PULL-OUT TEST INCLUDING THE EFFECTS OF FRICTION AND THERMAL STRESSES Advanced Composite Letters, Vol. 9, No. 6, 373-383 000 Full Article ANALYTICAL FRACTURE MECHANICS ANALYSIS OF THE PULL-OUT TEST INCLUDING THE EFFECTS OF FRICTION AND THERMAL STRESSES John A. Nairn Material

More information