International Journal of Solids and Structures
|
|
- Damon Owen
- 5 years ago
- Views:
Transcription
1 Accepted Manuscript Micromechanical method for effective piezoelectric properties and electromechanical fields in multi-coated long fiber composites George Chatzigeorgiou, Ali Javili, Fodil Meraghni PII: S-7683(8) DOI: Reference: SAS To appear in: International Journal of Solids and Structures Received date: 7 April 8 Revised date: 3 September 8 Accepted date: 4 September 8 Please cite this article as: George Chatzigeorgiou, Ali Javili, Fodil Meraghni, Micromechanical method for effective piezoelectric properties and electromechanical fields in multi-coated long fiber composites, International Journal of Solids and Structures (8), doi: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
2 Micromechanical method for effective piezoelectric properties and electromechanical fields in multi-coated long fiber composites George Chatzigeorgiou a, Ali Javili b, Fodil Meraghni a, Abstract a Arts et Métiers ParisTech, LEM3-UMR 739 CNRS, 4 Rue Augustin Fresnel 5778 Metz, France b Department of Mechanical Engineering, Bilkent University, 68 Ankara, Turkey This paper proposes a micromechanical framework for identifying the macroscopic behavior of multi-coated long fiber composites, as well as the average electromechanical microscopic fields of all phases (matrix, fibers, coating layers), generated upon known macroscopic conditions. The work aims at developing a unified micromechanical approach that provides an analytical solution standing for non-coated and multi-coated long fiber composites with transversely isotropic piezoelectric behavior. The proposed method solves specific boundary value problems and utilizes the Mori-Tanaka homogenization scheme, in which the dilute strain and electric field concentration tensors are obtained analytically with the help of an extended composite cylinders method that accounts for coupled electromechanical fields. The capabilities of this homogenization strategy are illustrated with the help of numerical examples, and comparisons with known solutions from the literature for non-coated and coated fiber piezoelectric composites are provided. Keywords: piezoelectricity; fiber composites; multi-coated fibers; composite cylinders method.. Introduction Piezoelectric materials are very attractive in applications involving the design of sensors, actuators, transducers, etc. due to their uniue capability to convert electrical into mechanical energy. Using piezoelectric ceramics, like PZT, in bulk form is not always convenient, mainly due to Corresponding author. addresses: georges.chatzigeorgiou@ensam.eu (George Chatzigeorgiou), ajavili@bilkent.edu.tr (Ali Javili), fodil.meraghni@ensam.eu (Fodil Meraghni) Preprint submitted to International Journal of Solids and Structures September 5, 8
3 their increased weight. To avoid such issues, an efficient solution is to combine these materials with non-piezoelectric polymers in the form of composites. These advanced composite materials opened new horizons in the development of new transducers and sensors with high strength, low thermal expansion coefficients, increased thermal conductivities and decreased dielectric constants. During the last 3 years several models have been proposed in the literature to study the piezoelectric and the combined thermo-magneto-electro-elastic response of composites. Dunn and Taya (993) have introduced an Eshelby-type approach by identifying appropriate Eshelby and concentration tensors for the combined electromechanical response. This techniue later was extended to account for other phenomena, linear (Li and Dunn, 998) and nonlinear (Hossain et al., 5), while Zou et al. () identified Eshelby tensors for arbitrary shaped piezoelectric inclusions. Benveniste (994), based on the initial framework of Benveniste and Dvorak (99), studied the macroscopic response of piezoelectric composites using partially the composite cylinders method. Aboudi () developed a computational method for coupled electro-magneto-thermo-elastic composites and Lee et al. (5) have proposed numerical and Eshelby-based analytical strategies for three phase electro-magneto-elastic composites. Piezoelectric composites (Berger et al., 5, 6; Maruccio et al., 5), electro-magneto-thermo-elastic composites (Bravo-Castillero et al., 9) and piezoelastic plates (Kalamkarov and Kolpakov, ) have been studied by using the periodic homogenization theory. A homogenization approach for studying piezoelectric composites with periodic and random microstructure was proposed by Spinelli and Lopez-Pamies (4). Sharma et al. (7) identified a theoretical framework that describes the conditions under which non-piezoelectric materials can be used for the design of piezoelectric nanocomposites. Ray and Batra (9) developed a micromechanical scheme for studying piezoelectric composites with suare cross section fibers, and similar techniue was developed later for magneto-electro-elastic with suare cross section fibers (Pakam and Arockiarajan, 4). Koutsawa et al. (), using Hill s interfacial operators, proposed a self consistent scheme for studying thermo-electro-elastic properties of composites with multi-coated ellipsoidal particles. Wang et al. (4) developed a micromechanical method for piezoelectric composites with imperfect interfaces between the ellipsoidal particles and the matrix phase, using the concept of euivalent particle. Piezoelectric
4 composites with imperfect interfaces have been also studied by Gu et al. (4, 5). The goal of this work is to develop a unified micromechanical approach aimed at providing an analytical solution that stands for non-coated and multi-coated long fiber composites with transversely isotropic piezoelectric behavior. Studying the coating between a fiber and a matrix in composites is of vital importance. A lot of nonlinear deformation mechanisms, like plasticity and/or martensite transformation, occur freuently at a small region close to the fibers and lead to an interaction with the local damage of the fiber/matrix interface (Payandeh et al.,, ). Developing appropriate computational tools that identify the mechanical and electric fields at the proximity of the fibers can assist in the design of more accurate damage and failure criteria for piezoelectric composites. The developed approach is based on solving specific boundary value problems, extending the composite cylinders model of Hashin and Rosen (964). This effort can be considered as a generalization of the Dvorak and Benveniste (99); Benveniste (994) methodology, providing analytical expressions of the dilute strain-electric field coupled concentration tensors, which can be utilized in classical micromechanical techniues, like Mori-Tanaka or self consistent. The advantage of such information is that it permits to identify not only the overall response of the composite, but also the various average electromechanical fields generated at the matrix, the fiber and the coating layers for known macroscopic electromechanical conditions. For non-coated fibers the obtained dilute tensors are euivalent with those of Dunn and Taya (993). To the best of the authors knowledge, the only available framework in the literature that computes concentration tensors for piezoelectric composites with coated fibers is the one of Koutsawa et al. (), but it is based on certain approximations. The approach discussed herein does not reuire such approximations, since the solution proposed in this work for the Eshelby s inhomogeneity problem is exact. The organization of the manuscript is as follows: After the introduction, section starts with a small recall on the piezoelectricity concepts, and then it describes the Mori-Tanaka type micromechanical framework and the general procedure for obtaining the dilute concentration tensors for composites with transversely isotropic piezoelectric material constituents (matrix, fiber, coating layers). Section 3 presents the case of non-coated fiber composites and discusses the consistency of the approach with published results from the open literature. Section 4 discusses the case of 3
5 coated fibers with one coating layer. Numerical examples and comparisons with existing finite element calculations and other micromechanical approaches are demonstrated in section 5. The section also includes a study in which, for given macroscopic strain and electric field and known volume fractions of the phases, all the electromechanical fields at every phase and at the overall composite are computed (Table ). The paper finishes with a section giving some concluding remarks and future developments. For the purpose of the paper s completeness, two appendices at the end of the article provide the piezoelectric Eshelby tensor of Dunn and Taya (993) and explain briefly the framework of Koutsawa et al. ().. Micromechanical approach for a coated fiber/matrix piezoelectric composite.. General concepts from piezoelectricity and notations In a linear piezoelectric material, the constitutive law that describes the relation between the stress tensor σ, the electric displacement vector D, the strain tensor ε and the electric field vector E is written in the following indicial form (double indices denote summation): σ i j = C i jkl ε kl e i jm E m, D i = e imn ε mn + κ i j E j. In tensorial notation, the above relations are expressed as σ = C:ε e E, D = e T :ε + κ E. () In the above relations C is the fourth order elasticity tensor, κ is the second order permittivity moduli tensor and e is the third order piezoelectric moduli tensor. The strain tensor ε and electric field vector E are related with the displacement vector u and the electric potential φ respectively through the relations ε = gradu + gradu T, E = gradφ. () Moreover, the euilibrium and electrostatic euations read divσ =, divd =. (3) 4
6 Unless stated otherwise, the Voigt notation will be adopted in the seuel to represent second, third and fourth order tensors. The convention used in strains and stresses to replace the indices i j to a single index k is the following: i j = k =, i j = k =, i j = 33 k = 3, i j = k = 4, i j = 3 k = 5, i j = 3 k = 6... Application of the Mori-Tanaka method N Figure : Unidirectional multi-coated fiber composite. The coating has N distinct layers. Let s consider a N + -phase composite, consisting of a matrix, denoted as, and infinitely long, multi-coated, cylindrical fibers. The fibers are denoted with the index, while the N coating layers are denoted with the indices,3,...,n (Figure ). All phases are assumed to be made by linear piezoelectric materials. The phases are characterized by their piezoelectric moduli C, e and κ, =,,, 3,..., N. According to the usual micromechanics arguments, a Representative Volume Element (RVE) with total volume V is sufficient to describe the overall response of the composite. Each phase 5
7 Ω Ω N r N r Ω Ω Figure : Typical RVE of a unidirectional multi-coated fiber composite. in the RVE occupies the space Ω (Figure ), has its own volume V and volume fraction c for =,,, 3,..., N. The classical volume summation rule states that N V = V, = r N c = = N = V V =. The macroscopic strain ε, stress σ, electric field E and electric displacement D corresponding to the RVE are computed from the volume averages of their microscopic counterparts at all phases, i.e. N N N N ε = c ε, σ = c σ, E = c E, D = c D. (4) = = In these expressions ε, σ, E and D denote the average strain, stress, electric field and electric displacement in the th phase respectively: ε = ε(x) dv, σ = σ(x) dv, V Ω V Ω The constitutive law of each phase states that = = E = E(x) dv, D = D(x) dv. (5) V Ω V Ω σ = C :ε e E, D = e T :ε + κ E, =,,, 3,..., N. (6) 6
8 The goal of homogenization is to identify a similar type of constitutive law for the macroscopic uantities, i.e. σ = C:ε e E, D = e T :ε + κ E, (7) where C, κ and e are the macroscopic elasticity, permittivity and piezoelectric tensors respectively. The law (7) represents the overall behavior of the composite. According to the Mori-Tanaka approach, the average fields in the fiber or the coatings and the matrix phase ( = ) are connected to each other through the dilute concentration tensors: T mm ε = T mm :ε + T me E, E = T em :ε + T ee E, =,, 3,..., N. (8) are fourth order tensors, written as 6 6 matrices, T me are third order tensors, written as 6 3 matrices, T em are third order tensors, written as 3 6 matrices and T ee are second order tensors, written as 3 3 matrices. Combining (8) with (4),3 and after some algebra, the following relations are obtained: where and ε = A mm :ε + A me E, E = A em :ε + A ee E, =,,, 3,..., N, (9) A mm = Π mm Π me Π ee Π em, A me = A mm :Π me Π ee, A em = Π ee Π em : A mm, A ee = Πee Π ee :Π em : A me, () A mm = T mm : A mm + T me A em A em = T em : A mm + T ee A em for =,, 3,..., N. In the above relations,, Ame, Aee = T mm : A me + T me A ee Π mm = c I + N c T mm, Π me = = N c T me, Π em = =, = T em : A me + T ee A ee, () N c T em, Π ee = c I + = N = c T ee, () Contrarily to the electromechanical fields, the macroscopic moduli are generally not eual to the volume averages of their microscopic counterparts. The Voigt bound, which considers such relation, often provides a poor estimate of the real response. 7
9 while I and I, with I i jkl = δ ikδ jl + δ il δ jk, I i j = δ i j, δ i j : Kronecker delta, denote the symmetric fourth order and second order identity tensors respectively. Substituting (9) in (6) gives σ = C : A mm ε A em ε e A ee C : A me E, D = e T : A mm + κ A em ε + e T : A me + κ A ee E, (3) for =,,, 3,..., N. Implementing these results in (4),4 and comparing with (7) eventually yields C = e T = N N c C : A mm ε A em, e = c e A ee C : A me, = = N N c e T : A mm + κ A em, κ = c e T : A me + κ A ee. (4) = Of course, the forms of the dilute concentration tensors T mm, T me, T em and T ee should respect the compatibility between the euations (4) and (4) 3. It is worth mentioning that it is possible to construct one single dilute concentration tensor that combines the four tensors of the expressions (8) (see for instance Dunn and Taya, 993). The advantage of using separate dilute tensors for each field (mechanical and electrical) and the couplings arising from them is that the scheme can be extended more naturally to account for additional mechanisms. Indeed, one can obtain dilute concentration tensors for thermomechanical response (Benveniste et al., 99; Chatzigeorgiou et al., 8), or for inelastic response, like in the case of the Transformation Field Analysis approach (Dvorak and Benveniste, 99; Dvorak, 99). =.3. Dilute piezoelectric concentration tensors For the identification of the dilute concentration tensors, a single linear piezoelectric coated fiber is assumed to be embedded inside the linear piezoelectric matrix, as shown in Figure 3. 8
10 u = ε x, φ = E x Ω r N Ω N Ω N Ω N Ω Ω r r Figure 3: Coated fiber inside a matrix, subjected to linear displacement and electric potential at far distance. All phases are piezoelectric materials. The fiber occupies the space Ω with volume V, its coatings the spaces Ω with volumes V, =, 3,..., N and the matrix occupies the space Ω, which is extended to far distance from the fiber. The matrix is subjected to linear displacement u = ε x and linear electric potential φ = E x at far distance from the fiber (r ). The interface between each phase is considered perfect. The interface between the phase and the phase + for =,..., N is denoted as Ω, while the interface between the last coating layer and the matrix is denotes as Ω N. The various electromechanical fields generated at every phase depend on the spatial position, i.e. u () (x), ε () (x), σ () (x), φ () (x), E () (x), D () (x), x Ω. (5) In infinitely long fiber composites, it is more convenient to describe all the necessary euations in cylindrical coordinates, by describing the position vector in terms of the radius r, the angle θ and the longitudinal position z (Figure 4). For transversely isotropic piezoelectric matrix, fiber and Ω Ω 9
11 z r θ y x Figure 4: Cylindrical coordinate system. coatings, the constitutive law (6) for each phase is written in Voigt notation as σ () rr σ () θθ σ () zz σ () rθ σ () rz σ () θz = K tr + µ tr K tr µ tr K tr µ tr l K tr + µ tr l l l n µ tr µ ax µ ax ε () rr ε () θθ ε () zz ε () rθ ε () rz ε () θz e 3 e 3 e 33 e 5 e 5 E r () E () θ E z (), (6)
12 ε () rr D () r D () θ D () z = e 5 e 5 e 3 e 3 e 33 ε () θθ ε () zz ε () rθ ε () rz ε () θz + κ κ κ 33 E r () E () θ E z () The strain tensor and the electric field vector at each phase are given by the expressions ε () rr = u() r r, ε () rz = u() z + u() r r E r () z, = φ() r, ε() θθ ε() θz = r = r E() θ u () θ + u() r θ r, u () z + u() θ θ = r φ () θ, z, ε() zz = u() z z, ε() rθ = u() θ E() z while the euilibrium and electrostatic euations are written as σ () rr r + r σ () rθ r σ () rz r D () r r σ () rθ + σ() rr θ + r + r + r σ () r θθ + σ() rz z σ () θθ + σ() rθ + σ() θz θ r z σ () θz + σ() rz + σ() zz θ r z D () θ + D() r + D() z θ r z r + r u () r u() θ θ r,. (7) = φ() z, (8) =, =, =, =. (9) The fiber is considered to have radius r = r and every coating has external radius r (Figure 3). The interface conditions between the phase and the phase + for =,, 3,..., N are
13 expressed as u () r (r, θ, z) = u (+) r (r, θ, z), u () θ (r, θ, z) = u (+) θ (r, θ, z), u () z (r, θ, z) = u (+) z (r, θ, z), σ () rr (r, θ, z) = σ (+) rr (r, θ, z), σ () rθ (r, θ, z) = σ (+) rθ (r, θ, z), σ () rz (r, θ, z) = σ (+) rz (r, θ, z), φ () (r, θ, z) = φ (+) (r, θ, z), D () r (r, θ, z) = D (+) r (r, θ, z). () Additionally, the interface conditions between the N th coating and the matrix are written as σ (N) u (N) r (r N, θ, z) = u () r (r N, θ, z), u (N) θ (r N, θ, z) = u () θ (r N, θ, z), u (N) z (r N, θ, z) = u () z (r N, θ, z), rr (r N, θ, z) = σ () rr (r N, θ, z), σ (N) rθ (r N, θ, z) = σ () rθ (r N, θ, z), φ (N) (r N, θ, z) = φ () (r N, θ, z), σ (N) rz (r N, θ, z) = σ () rz (r N, θ, z), D (N) r (r N, θ, z) = D () r (r N, θ, z). () Following the Eshelby s methodology, the dilute concentration tensors provide the relationship between the average strain ε and electric field E inside the phase ( =,, 3,..., N), and the strain ε and electric field E at the far field: ε = T mm :ε + T me E, E = T em :ε + T ee E. () The Mori-Tanaka approach assumes that the strain ε and the electric field E applied in the far field of this Eshelby-type problem correspond to the average strain ε and the average electric field E of the matrix in the RVE of Figure, which are used in euations (8). Using the divergence theorem, one obtains for the fiber the average fields ε = V ε () dv = Ω V E = V Ω E () dv = V Ω u() n + n u () ds, Ω φ () nds, (3) where n is the unit vector of the interface Ω. Due to the homothetic topology of the coated fiber, the unit vector is the same at all interfaces. For all the coating layers, the average strain and electric field are defined as ε = ε () dv = V Ω V Ω u() n + n u () ds V Ω u( ) n + n u ( ) ds, E = E () dv = φ () nds φ ( ) nds. (4) V Ω V Ω Ω
14 Obtaining the dilute piezoelectric tensors is not an easy task. For certain, ellipsoidal-type, forms of inclusions without coating, Dunn and Taya (993) have computed the Eshelby tensors for the combined strain/electric field system (see Appendix A for infinitely long fibers). In the present work, the dilute concentration tensors are computed directly for long coated fiber piezoelectric composites with transversely isotropic behavior at every phase. To achieve such goal, analytical solutions are utilized of similar boundary value problems with those of the composite cylinders method proposed by Hashin (99), taking into account the effects caused by the combined presence of mechanical and electric fields. In the pure mechanical problem, similar techniues have been utilized in the literature to obtain dilute (Benveniste et al., 989) and semi-dilute (Chatzigeorgiou et al., ) stress concentration tensors, as well as dilute strain concentration tensors (Wang et al., 6) for coated fiber composites. In cylindrical coordinates, the surface element in a surface of constant radius r (a vertical cylinder) is ds r = rdθdz and the surface element in a surface of constant z is ds z = rdrdθ. For an arbitrary tensor Q(r, θ, z) and a cylinder of radius r and length L, the sum of surface integrals with the general form F = Lπr L π L Q(r, θ, z)dθdz + Lπr π r Q(r, θ, L) Q(r, θ, L) rdrdθ, (5) is reuired for the computations of the average uantities (3) and (4). The three normal vectors in cylindrical coordinates are expressed as cos θ n = sin θ, n = sin θ cos θ, n 3 = The displacements of the phases are represented in matrix form as u () = u () n + u () r θ. (6) n + u () z n 3. (7) Another important point to be mentioned is that in long fiber composites with, at most, transversely isotropic phases (axis of symmetry: the direction of fibers), the dilute strain concentration tensors For infinitely long cylinder, L. To avoid infinite values, the division by volume takes care of L. 3
15 present transverse isotropy. In Voigt notation, they take the form T mm T mm T mm 44 T mm 3 T me 3 T mm T mm 44 T mm T mm 3 T me 3 T mm = T em = T mm 44 T mm 55 T mm T em 5 T em 5 55, T me, T ee = = T me 5 T me T ee 5 T ee To obtain the unknown terms of these tensors, four types of boundary value problems should be examined. In the seuel two cases are studied: a fiber/matrix composite without coating and a fiber/matrix composite with only one coating layer. 3. First case: non-coated fibers The first studied case considers non-coated piezoelectric fibers embedded in a piezoelectric matrix. This is a well examined problem in the literature and the scope of this section is to illustrate that the proposed methodology produces reliable results. 3.. Axial shear strain and in-plane electric field conditions For this case the following displacement vector and electric potential are applied: u (x,y,z) = βx., (8), φ (x,y,z) = γx. (9) 4
16 These correspond to the following strain tensor (in classical tensorial form) and electric field vector β/ γ ε = β/, E = In cylindrical coordinates, the displacement and electric potential are transformed to u (r,θ,z) = βr cos θ., φ (r,θ,z) = γr cos θ. For these boundary conditions, the displacement vector and electric potential at the matrix ( = ) and the fiber ( = ) are given by the general expressions u () z (r, θ) = βru z () (r) cos θ, u () r = u () θ =, U z () (r) = Ξ () + Ξ () r/r, φ () (r, θ) = γrφ () (r) cos θ, Φ () (r) = Ξ () 3 + Ξ () 4 r/r. (3) where Ξ () i, i =,, 3, 4, are unknown constants. These general expressions lead to stresses and electric displacements that satisfy the euilibrium and electrostatic euations (9). The important stresses and electric displacements for identifying the unknown constants are σ () rz (r, θ) = β Σ m() rz (r) + γ Σ e() rz (r) cos θ, D () r (r, θ) = β D m() r (r) + γ D e() r (r) cos θ, Σ m() rz (r) = µ ax Ξ () Ξ (), Σ e() r/r rz (r) = e 5 Ξ () 3 Ξ () 4, r/r D m() r (r) = e 5 Ξ (), D e() r/r r (r) = κ Ξ () 3 Ξ () 4. r/r The boundary conditions that should be satisfied in this boundary value problem are u () Ξ () z and φ () finite at r = Ξ () = Ξ () 4 =, u () z (r, θ) = βr cos θ and φ () (r, θ) = γr cos θ Ξ () = Ξ () 3 =. (3) Considering these results, the interface conditions u () z (r, θ) = u () z (r, θ), σ () rz (r, θ) = σ () rz (r, θ), φ () (r, θ) = φ () (r, θ), D () r (r, θ) = D () r (r, θ), 5
17 construct the linear system K β I m + γ I e Ξ m + Ξ e = β F β + γ F γ, with K = Ξ m = F β = µ ax µ ax e 5 e 5, I m = e 5 e 5 κ κ T T Ξ () Ξ (), Ξ e = Ξ () 3 Ξ () 4, T T µ ax e 5, Fγ = e 5 κ., I e = The solution of this linear system is written in the form Ξ m = I m K F β + γβ K F γ, Ξ e = I e β γ K F β + K F γ. After some algebra, the following results are obtained: Ξ () = B + γ β B, Ξ () 3 = β γ B 3 + B 4, e 5 e5 + e 5 + µ ax κ + κ B = e5 + e 5 + κ + κ µ ax e 5 µ ax B 3 = e 5 µ ax e5 + e 5 + κ + κ µ ax + µax + µax e 5 κ e 5 κ, B = e5 + e 5 + κ + κ µ ax +, µax e 5 e5 + e 5 + κ µ ax, B 4 = + µax e5 + e 5 + κ + κ µ ax +. µax, 6
18 Implementing the transformations (5), (6) and (7) in (3) yields the average strain and electric field inside the fiber, β/ γ/ ε = U z () (r ) ε ε = B + B β/ γ/ β γ E = Φ () (r ) E E = B 3 + B 4. Comparing these results with the expressions () and (8) it becomes clear that T mm 55 = ε xz e 5 e5 + e 5 + µ ax κ + κ = B β/ = γ= e5 + e 5 + κ + κ µ ax +, µax T me 5 = ε xz e 5 κ e 5 κ = B γ/ = β= e5 + e 5 + κ + κ µ ax +, µax T em 5 = E x e 5 µ ax = B β 3 = e 5 µ ax γ= e5 + e 5 + κ + κ µ ax +, µax T ee = E x e 5 e5 + e 5 + κ µ ax = B γ 4 = + µax β= e5 + e 5 +. (3) κ + κ µ ax 3.. Transverse shear strain conditions For this case the following displacement vector is applied: u (x,y,z) = which corresponds to the strain tensor (in classical tensorial form) βy βx β ε = β + µax,, (33). 7
19 In cylindrical coordinates, the displacement is transformed to the vector u (r,θ,z) = βr sin θ βr cos θ For these boundary conditions, the displacements at each phase are given by the general expressions u () r (r, θ) = β r U r () (r) sin θ, u () θ U () U () (r, θ) = β r U() θ (r) cos θ, r (r) = Ktr µ tr r/r K tr + µ tr Ξ () + Ξ () r/r 4 Ξ() 3 + Ktr + µ tr µ tr r/r Ξ() 4, θ (r) = r/r Ξ () + Ξ () + r/r 4 Ξ() 3 + r/r Ξ() 4, (34) where the unknowns that need to be defined are the Ξ () i, i =,, 3, 4. These expressions lead to stresses and electric displacements that satisfy the euilibrium and electrostatic euations (9). The important stresses for the calculations are given by σ () rr (r, θ) = β Σ () rr (r) sin θ, σ () rθ (r, θ) = β Σ() rθ (r) cos θ, rr (r) = µ tr Ξ () + 6µ tr r/r 4 Ξ() 3 4K tr r/r Ξ() Σ () Σ () µ tr K tr + µ tr rθ (r) = 6K tr r/r Ξ () + µ tr Ξ (). 4, 6µ tr r/r 4 Ξ() 3 + K tr r/r Ξ() The boundary conditions that should be satisfied in this boundary value problem are u () u () r, u () θ finite at r = Ξ () 3 = Ξ () 4 =, r (r, θ) = βr sin θ and u () θ (r, θ) = βr cos θ Ξ() =, Ξ () =. (35) 4. Considering these results, the interface conditions u () r (r, θ) = u () r (r, θ), u () θ (r, θ) = u () θ (r, θ), σ () rr (r, θ) = σ () rr (r, θ), σ () rθ (r, θ) = σ () rθ (r, θ), 8
20 construct the linear system K Ξ = F, with K = K tr µtr K tr + µtr Ktr µ tr 6K tr µtr K tr + µtr + µtr µ tr 6µ tr 4K tr µ tr 6µ tr K tr The solution of this linear system gives Ξ () =, Ξ () = µ tr, Ξ = K tr µ tr µ tr µtr + Ktr + µtr Ξ () Ξ () Ξ () 3 Ξ () 4 + µtr., F = Implementing the transformations (5), (6) and (7) in (3) yields the average strain inside the fiber, ε = U () r (r ) + U () θ (r ) ε ε = Ξ () ε. Comparing these results with the expressions () and (8) it becomes clear that T mm 44 = ε xy β = µ tr Ξ() K tr = + µtr. (36) µ tr 3.3. Plane strain and axial electric field conditions µ tr µtr + Ktr + µtr µ tr µ tr For this case the following displacement vector and electric potential are applied: u (x,y,z) = βx βy, φ (x,y,z) = γz. (37) 9.
21 These correspond to the following strain tensor (in classical tensorial form) and electric field vector β ε = β, E = In cylindrical coordinates, the displacement and electric potential are transformed to βr u (r,θ,z) = γ., φ (r,θ,z) = γz. For these boundary conditions, the displacement vector and electric potential at the matrix ( = ) and the fiber ( = ) are given by the general expressions u () r (r) = βru r () (r), u () θ = u () z =, U r () (r) = Ξ () + Ξ () r/r, φ () (z) = γz, where Ξ () i, i =,, are unknown constants. These general expressions lead to stresses and electric displacements that satisfy the euilibrium and electrostatic euations (9). The important stresses for identifying the unknown constants are σ () rr (r) = β Σ m() rr (r) + γ Σ e() rr, Σ m() rr (r) = K tr Ξ () µ tr Ξ (), Σ e() r/r rr = e () The boundary conditions that should be satisfied in this boundary value problem are u () r finite at r = D () =, u () Considering these results, the interface conditions 3. (38) r (r ) = βr D () =. (39) u () r (r ) = u () r (r ), σ () rr (r ) = σ () rr (r ), construct the linear system K Ξ = F β + γ β F γ,
22 with K = K tr µ tr, Ξ = The solution of this linear system gives Ξ () Ξ () Ξ (), F β = = B + γ β B, B = Ktr + µtr K tr + µtr K tr, F γ =, B = e 3 e 3 K tr + µtr e 3 e 3.. (4) Implementing the transformations (5), (6) and (7) in (3) yields the average strain inside the fiber, β γ ε = U r () (r ) ε ε = B β + B γ The general form of the dilute concentration tensors () and (8) permits to write ε xx = T mm T mm 44 β + T me 3 γ. From (4) and (4), it becomes clear that T mm = K tr + µtr K tr + T mm 3.4. Hydrostatic strain conditions + µtr 44 For this case the following displacement vector is applied: u (x,y,z) = K tr + µtr which corresponds to the strain tensor (in classical tensorial form). (4), T me 3 = e 3 e 3. (4) βx βy βz β ε = β β, (43).
23 In cylindrical coordinates, the displacement is transformed to the vector u (r,θ,z) = For these boundary conditions, the displacement vector and electric potential at the matrix ( = ) and the fiber ( = ) are given by the general expressions u () r βr βz. (r) = βru r () (r), u () θ =, u () z (z) = βz, U r () (r) = Ξ () + Ξ () r/r, (44) where Ξ () i, i =,, are unknown constants. These general expressions lead to stresses and electric displacements that satisfy the euilibrium and electrostatic euations (9). The important stresses for identifying the unknown constants are σ () rr (r) = β Σ () rr (r), Σ () rr (r) = K tr Ξ () µ tr Ξ () + l r/r. The boundary conditions that should be satisfied in this boundary value problem are u () r finite at r = D () =, u () Considering these results, the interface conditions construct the linear system with K = u () K tr µ tr r (r ) = βr D () =. (45) r (r ) = u () r (r ), σ () rr (r ) = σ () rr (r ),, Ξ = K Ξ = F, Ξ () Ξ (), F = K tr + l l.
24 The solution of this linear system gives Ξ () = Ktr + µtr K tr + µtr + l l K tr +. (46) µtr Implementing the transformations (5), (6) and (7) in (3) yields the average strain inside the fiber, ε = U r () (r ) ε ε = Ξ () ε. (47) The general form of the dilute concentration tensors () and (8) permits to write From (46), (47) and (4), it becomes clear that 3.5. Macroscopic properties ε xx = T mm T mm 44 + T mm 3 β. l l T mm 3 = K tr + µtr. (48) The dilute concentration tensors (8), whose terms are given by (3), (36), (4) and (48), are euivalent to the global piezoelectric concentration tensor provided by Dunn and Taya (993) (see Appendix A). Assuming that the volume fraction of fibers in the composite is c and applying the Mori-Tanaka scheme, as described in subsection., the following macroscopic properties are obtained: K tr = cktr K tr µtr + Ktr + µtr cµtr ck tr cktr + Ktr + µtr l = c K tr l + c K tr, µ tr = + µtr + µtr l ck tr cktr + Ktr + µtr µ tr µ tr cµ tr, n = µ tr µtr + Ktr cµtr + µtr e 3 = ce 3 K tr + µtr + ce3 K tr + µtr ck tr cktr + Ktr +, e 33 = µtr µ tr cµtr + µtr + cµtr + K tr µ tr Z ck tr cktr + Ktr + µtr Z + cµtr + µtr cµtr, µ ax = W, W ck tr cktr + Ktr +, e 5 = W 3, µtr W, κ = W 4 W, κ 33 = Z 3 ck tr cktr + Ktr +, µtr 3
25 where Z = c l l K tr Ktr n n + n cn + cn K tr + µtr c l l + n K tr Ktr, Z = c K tr Ktr e33 e 33 +c l l e 3 e 3 + e33 K tr Z 3 = c e3 e 3 + K tr + l l e 3 e 3 + e33 ce 33 + ce 33 K tr + µtr Ktr +c e3 e 3 + κ 33 K tr µtr κ33 κ 33 Ktr, W = + c e 5 + c e 5 e 5 + c e 5 + κ + cκ + κ cκ µ ax, + κ33 cκ 33 + cκ 33 K tr + cµax + µax cµax W = c e 5 e 5 µ ax + c e 5 + e 5 µ ax + 4ce 5 µ ax + κ + cκ + κ cκ µ ax µ ax cµax + µax + cµax W 3 = c e e 5 e c e 5 + 4cκ µ ax + c κ µ ax + cµax + µax cµax + κ µ ax + µtr,, e5 cµax + µax + cµax W 4 = c e 5 e 5 κ + c e 5 + e 5 κ + 4ce 5 κ + κ cκ + κ + cκ µ ax κ. + cµax + µax cµax e5, The obtained solution respects the compatibility condition between euations (4) and (4) 3. The above properties agree with the findings of Benveniste (994). The transverse bulk modulus K tr is given by the same expression, while for the transverse shear modulus µ tr, Benveniste also concluded that it depends exclusively on the mechanical properties. For calculating µ tr he proposed the use of the Christensen and Lo (979) approach, which is a generalized self consistent method and provides uite similar results with Mori-Tanaka. The rest of the properties satisfy the universal relations (Benveniste and Dvorak, 99) K tr Ktr = Ktr Ktr = cl + cl l l l l l cn () + cn () n = ce + ce 3 e 3 3 ce () 33 + ce() 33 e, 33 K tr Ktr = Ktr Ktr cl + cl l = e 3 e 3 e 3 e 3 ce () 33 + ce() 33 e = e 3 ce 3 ce 3, 33 cκ 33 + cκ 33 κ 33 4
26 and µ ax e 5 κ Det µ ax e 5 κ µ ax e 5 κ There are two main advantages of the proposed methodology compared to that of Benveniste (994): a) The computation of the strain and electric field concentration tensors allow the identification of the microscopic electromechanical fields for known macroscopic fields. Such information can be very valuable for detecting possible inelastic mechanisms due to high stresses in the matrix or the fiber. b) The dilute concentration tensors can be used into the Mori-Tanaka scheme for obtaining the macroscopic response of composites with more complex microstructure, for instance with two different types of fibers, or the same fibers with random orientation (Chatzigeorgiou et al., ; Seidel et al., 4). 4. Second case: fibers with one coating layer In this second studied case, the composite consists of coated fibers and matrix. Each fiber has radius r = r and is coated with a layer that has external radius r = r. In the seuel, the ratio ρ = r /r will be utilized. Considering as c the volume fraction of the fibers, the coating layers and the matrix volume fraction can be computed through the ratio ρ. Indeed, the volume fractions of all the phases (fiber, coating, matrix) are =. c = c, c = c ρ c, c = c c = c ρ. (49) The dilute concentration tensors are computed with similar boundary value problems as those discussed in the previous section. Here, only the essential different points are going to be presented. For the sake of the procedure s simplicity, only one coating layer is considered for this section. Nevertheless, the approach is general and the addition of more coating layers is straightforward. 5
27 4.. Axial shear strain and in-plane electric field conditions Similarly to the non-coated fibers, the infinite matrix is subjected to the far field displacement vector and electric potential of (9). For all phases (matrix, fiber, coating) the general expressions (3) for these fields introduce unknowns (4 constants per phase). After eliminating four from the boundary conditions (3), the rest of the unknowns are identified from the interface conditions that hold between the fiber and the coating and between the coating and the matrix. The produced system of linear euations take the form with K = K β I m + γ I e Ξ m + Ξ e = β F β + γ F γ, ρ ρ µ ax µ ax µ ax e 5 e 5 e 5 µ ax ρ µ ax ρ µ ax e 5 ρ e 5 ρ e 5 ρ ρ e 5 e 5 e 5 κ κ κ e 5 ρ e 5 ρ e 5 κ ρ κ ρ κ, 6
28 I m = Ξ m = F β =, I e = T Ξ () Ξ () Ξ () Ξ (), Ξ e = Ξ () 3 Ξ () 3 Ξ () T T µ ax e 5, Fγ = e 5 κ. The solution of this linear system is expressed as Ξ m = I m K F β + γβ K F γ, Ξ e = I e β γ K F β + K F γ., 4 Ξ () 4 Considering the dilute tensors, the important terms Ξ (), Ξ() 3, Ξ() and Ξ () 4 can be written in the general form Ξ () = B + γ β B, Ξ () 3 = β γ B 3 + B 4, Ξ () = B 5 + γ β B 6, Ξ () 4 = β γ B 7 + B 8. Implementing the transformations (5), (6) and (7) in (3) yields the average strain and electric field inside the fiber, β/ γ/ ε = U z () (r ) ε ε = B + B β/ γ/ β γ E = Φ () (r ) E E = B 3 + B 4., T, On the other hand, implementing the transformations (5), (6) and (7) in (4) and taking into account the displacement and electric potential continuity conditions yields the average strain and 7
29 electric field inside the coating, ε = ρ U() z (r ) ρ ρ U() z (r ) β/ ε γ/ ε = + ρ B 5 ρ B ρ + ρ B 6 ρ B ρ β/ γ/ E = ρ Φ() (r ) ρ ρ Φ() (r ) E β γ E = ρ B 7 ρ B 3 ρ + + ρ B 8 ρ B 4 ρ. Comparing these results with the expressions () and (8) it becomes clear that T mm T mm 55 = B, T me 5 = B, T em 5 = B 3, T ee = B 4, 55 = + ρ B 5 ρ B, T me ρ 5 = ρ B 6 ρ B, ρ T em 4.. Transverse shear strain conditions 5 = ρ B 7 ρ B 3, T ee ρ = + ρ B 8 ρ B 4. (5) ρ For the transverse shear term of the dilute concentration tensor, the infinite matrix is subjected to the far field displacement vector of (33). For all phases (matrix, fiber, coating) the general expressions (34) introduce unknowns (4 constants per phase). After eliminating four from the boundary conditions (35), the rest of the unknowns are identified from the interface conditions that hold between the fiber and the coating and between the coating and the matrix. The produced system of linear euations take the form K Ξ = F,, 8
30 with K tr µtr K tr + µtr Ktr µtr K tr + µtr Ktr + µtr µ tr Ξ = K = K tr µtr ρ K tr + µtr ρ 4 ρ Ktr + µtr µ tr ρ 4 ρ Ktr + µtr µ tr /ρ ρ 4 ρ ρ 4 ρ µ tr µ tr 6K trµtr K tr + µtr µ tr 6Ktr µtr K tr + µtr µ tr 6µ tr 4K tr 6µ tr K tr µ tr 6ρ 4 µ tr 4ρ K tr 6ρ 4 µ tr 4ρ K tr 6K trµtr ρ K tr + µtr µ tr Ξ () Ξ () Ξ () Ξ () Ξ () 3 Ξ () 4 Ξ () 3 Ξ () 4 6ρ 4 µ tr ρ K tr 6ρ 4 µ tr ρ K tr T, F = µ tr µ tr The solution of this linear system is expressed as Ξ = K F. Implementing the transformations (5), (6) and (7) in (3) and (4) and taking into account the displacement continuity conditions yields the average strain inside the fiber and the coating, ε = U () r (r ) + U () θ (r ) 3K tr ε ε = 4K tr + Ξ () µtr + Ξ () ε, ε = U () ρ r (r ) + U () θ (r ) ρ U () ρ r (r ) + U () θ (r ) ε ε = + ρ Ktr + µtr ρ µ tr Ξ () 3K tr 4 ρ 4K tr Ξ () ρ Ξ () ε. + µtr Comparing these results with the expressions () and (8) it becomes clear that, T. T mm 44 = T mm 44 = 3K tr 4K tr + µtr ρ Ξ () + Ξ (), + ρ Ktr + µtr µ tr Ξ () 9 3K tr 4 ρ 4K tr + µtr Ξ () ρ Ξ (). (5)
31 4.3. Plane strain and axial electric field conditions Similarly to the corresponding boundary value problem of the previous section, the infinite matrix is subjected to the far field displacement vector and electric potential of (37). For all phases (matrix, fiber, coating) the general expressions (38) introduce 6 unknowns ( constants per phase). After eliminating two from the boundary conditions (39), the rest of the unknowns are identified from the interface conditions that hold between the fiber and the coating and between the coating and the matrix. The produced system of linear euations take the form with K = Ξ = F γ = K Ξ = F β + γ β F γ, ρ ρ K tr K tr µ tr K tr ρ µ tr ρ µ tr T T, Fβ = K tr, Ξ () Ξ () Ξ () Ξ () e 3 e 3 e 3 e 3 T. The solution of this linear system is expressed as Ξ = K F β + γ β K F γ. Considering the dilute tensors, the important terms Ξ () and Ξ () can be written in the general form Ξ () = B + γ β B, Ξ () = B 3 + γ β B 4., 3
32 Implementing the transformations (5), (6) and (7) in (3) and (4) and taking into account the displacement continuity conditions yields the average strain inside the fiber and the coating, β γ ε = U r () (r ) ε = B β + B γ, ε = ρ U() r (r ) ρ ρ U() r (r ) ε β γ = + ρ B 3 ρ B ρ β + ρ B 4 ρ B ρ γ. The expressions of the dilute concentration tensors () and (8) permit to write ε xx = T mm T mm 44 β + T me 3 γ, ε xx = T mm T mm 44 β + T me 3 γ. From these results it becomes clear that = B + T mm 44, T me 3 = B, T mm = + ρ B 3 ρ B + ρ T mm 44, T me 3 = ρ B 4 ρ B. (5) ρ T mm 4.4. Hydrostatic strain conditions For the remaining terms of the dilute mechanical concentration tensors, the infinite matrix is subjected to the far field displacement vector (43). For all phases (matrix, fiber, coating) the general expressions (44) introduce 6 unknowns ( constants per phase). After eliminating two from the boundary conditions (45), the rest of the unknowns are identified from the interface conditions that hold between the fiber and the coating and between the coating and the matrix. The produced system of linear euations take the form K Ξ = F, 3
33 with K = Ξ = ρ ρ, K tr K tr µ tr K tr ρ µ tr ρ µ tr T T, F = l l K tr + l l. Ξ () Ξ () Ξ () Ξ () The solution of this linear system is expressed as Ξ = K F. Implementing the transformations (5), (6) and (7) in (3) and (4) and taking into account the displacement continuity conditions yields the average strain inside the fiber and the coating, ε = U r () (r ) ε = Ξ () ε, ε = ρ U() r (r ) ρ ρ U() r (r ) ε = + ρ Ξ () ρ Ξ () ε ρ. The general form of the dilute concentration tensors () and (8) permits to write ε xx = T mm T mm 44 + T mm 3 β, ε xx = T mm T mm 44 + T mm 3 β. From these results it becomes clear that 3 = Ξ () T mm + T mm 44, T mm 3 = + ρ Ξ () ρ Ξ () T mm ρ + T mm 44. (53) T mm The dilute concentration tensors (8), whose terms are given by (5), (5), (5) and (53), can be directly utilized in the Mori-Tanaka scheme, described in subsection.. 5. Numerical examples To illustrate the capabilities of the proposed methodology, several numerical studies are performed. In the following examples the matrix phase is considered a non-piezoelectric epoxy, the 3
34 fibers are PZT-7A and the coating layers (if they exist) are made of PZT-5. The piezoelectric properties for these materials are taken from Berger et al. (5, 6) and they are summarized in Table. In the first example, no coating phase is considered. Several macroscopic piezoelectric properties of the fiber composite as a function of the fibers volume fraction are illustrated in Figure 5. In the same plots the numerical results of Berger et al. (6), which solve periodic unit cell problems with hexagonal arrangement of the fibers through finite element calculations, are utilized as reference values. As shown in the Figure, except from a slight under-prediction in the transverse shear modulus at high fiber content (above 5% volume fraction), the analytical method provides excellent accuracy in identifying the macroscopic response of the composite. As expected, the absolute values of the electromechanical properties show significant increase, almost exponential, at above 5% fiber volume fraction. material parameters material n l K tr µ tr µ ax e 3 e 33 e 5 κ κ 33 epoxy PZT-7A PZT units GPa GPa GPa GPa GPa Cm Cm Cm ncv m ncv m Table : Material parameters for the matrix (epoxy), the fiber (PZT-7A) and the coating layer (PZT-5). The values for the first two materials are taken from Berger et al. (6), while for the third material are taken from Berger et al. (5). In the second parametric study, the same macroscopic properties are compared with those provided by coated fiber composites. The comparisons consider the same fiber content c, while the coating and the matrix have volume fractions that depend on the radii ratio ρ = r /r between the fiber and the coating/matrix interface, as given by the expressions (49). Since the coating is stronger than the matrix phase, increase in the absolute values of the electromechanical properties is expected for the coated fiber composite. In Figure 6 two coating thicknesses are considered, one with ρ =.95 and the second with ρ =.9. In the same Figure, the predictions of the 33
35 3 K tr (GPa) µ ax (GPa) e5 (Cm ) FE hexagonal proposed method fiber volume fraction (%) 8 4 (a) FE hexagonal proposed method fiber volume fraction (%) (c) FE hexagonal proposed method µ tr (GPa) e3 (Cm ) FE hexagonal proposed method fiber volume fraction (%) κ (ncv m ) fiber volume fraction (%) (b) FE hexagonal proposed method fiber volume fraction (%) (d) FE hexagonal proposed method (e) (f) fiber volume fraction (%) Figure 5: Macroscopic piezoelectric properties of fiber composites without coating: (a) K tr, (b) µ tr, (c) µ ax, (d) e 3, (e) e 5 and (f) κ. Comparison between the proposed method and the numerical results of Berger et al. (6) using finite element computations in periodic hexagonal unit cells. 34
36 proposed method are compared with those obtained by the Mori-Tanaka scheme, when the dilute concentration tensors are obtained through the Koutsawa et al. () framework. Details about this approach are given in Appendix B. The obtained results show that the coupled piezoelectric parameters are the properties with the highest increase. At 6% fibers volume fraction and ρ =.9, both e 3 and e 5 become more than 4% higher (in absolute value) compared to the corresponding parameters of the non-coated fiber composite. For the same ρ and fibers volume fraction, the gain in mechanical properties does not exceed 55% and in electrical permittivity does not exceed 66%. The proposed methodology and the one based on the Mori-Tanaka with the dilute concentration tensors of the Koutsawa et al. () theory (denoted as KBBNC in Figure 6) render exactly the same predictions for all the macroscopic properties, with the exception of the transverse shear modulus. For the µ tr there is an insignificant difference (less than.%), which is due to the assumption made by Koutsawa et al. () in the usage of the iterfacial operators (see euation (B.7) and the discussion in Appendix B). It is worth mentioning that, in the case of multiple coatings, such assumption affects all the macroscopic properties and it can cumulate a non-negligible error in the macroscopic response. On the other hand, the methodology proposed in this work does not suffer from this issue, since the Eshelby s inhomogeneity problem is solved analytically and the dilute concentration tensors are computed through exact averaging. As already mentioned in the previous sections, the main advantage of identifying dilute concentration tensors is that they provide information about the average microscopic fields of the various phases, when the macroscopic strain and electric field are known. To demonstrate this capability, it is considered that at a coated fiber composite with fibers volume fraction 45% and 35
37 K tr (GPa) µ ax (GPa) e5 (Cm ) no coating ρ=.95 ρ=.9 ρ=.95, KBBNC ρ=.9, KBBNC fiber volume fraction (%) (a) no coating ρ=.95 ρ=.9 ρ=.95, KBBNC ρ=.9, KBBNC fiber volume fraction (%) (c) no coating ρ=.95 ρ=.9 ρ=.95, KBBNC ρ=.9, KBBNC fiber volume fraction (%) µ tr (GPa) e3 (Cm ) (e) κ (ncv m ) no coating ρ=.95 ρ=.9 ρ=.95, KBBNC ρ=.9, KBBNC fiber volume fraction (%) (b).3 no coating ρ=.95.4 ρ=.9 ρ=.95, KBBNC.5 ρ=.9, KBBNC fiber volume fraction (%) (d) no coating ρ=.95 ρ=.9 ρ=.95, KBBNC ρ=.9, KBBNC fiber volume fraction (%) Figure 6: Macroscopic piezoelectric properties of fiber composites: (a) K tr, (b) µ tr, (c) µ ax, (d) e 3, (e) e 5 and (f) κ. Comparison between non-coated fibers and coated fibers with ρ =.95 and ρ =.9. The dashed lines correspond to the prediction of the current methodology and the points to the predictions of the Mori-Tanaka scheme with the dilute concentration tensors according to the Koutsawa et al. () theory. (f) 36
38 ρ =.9, the following uniform macroscopic strain tensor and electric field vector are applied: ε = ε xx ε yy ε zz ε xy ε xz ε yz = 4, E = 4 6 E x E y E z = 4 GV For these macroscopic conditions, the average microscopic fields at all phases (matrix, fiber, coating), as well as the macroscopic stresses and electric displacements are summarized in Table. 6. Conclusions The article presented a micromechanical approach for identifying dilute concentration tensors for piezoelectric multi-coated, long fiber composites, which can be utilized in classical homogenization schemes, such as Mori-Tanaka and self consistent. The procedure shown here for noncoated and coated fibers with one layer is applicable to fibers with multiple coating layers, by using the boundary value problems and the analytical formalism described in sections 3 and 4. The extensive comparisons with known analytical and computational results from the literature illustrated the method s reliability and capabilities. The main advantage of the proposed methodology is that it solves analytically the Eshelby s inhomogeneity problem, which in turn permits the exact computation of the dilute concentration piezoelectric tensors. The described procedure is also easily extendable to account for other types of mechanisms, like thermomechanical behavior or inelastic response through a transformation field analysis type approach. In addition, the proposed approach provides new capabilities in studying piezoelectric composite materials and permits, with proper modifications, to account for imperfect mechanical (Chatzigeorgiou et al., 7) and/or dielectric (Chatzigeorgiou et al., 5) interfaces between the material constituents. m. 37
39 component macroscopic matrix fiber coating units ε xx ε yy ε zz ε xy ε xz ε yz σ xx MPa σ yy MPa σ zz MPa σ xy MPa σ xz MPa σ yz MPa E x GV/m E y GV/m E z GV/m D x C/m D y C/m D z C/m Table : Macroscopic and average microscopic fields in a coated fiber composite with fibers volume fraction 45% and ρ =.9. 38
40 Appendix A. Eshelby tensor for infinitely long cylinder in piezoelectric material In a piezoelectric material, the constitutive law euations () can be expressed in the global matrix form Σ = L Z, (A.) where Σ and Z are the stress/electric displacement and strain/electric field tensors respectively, written as 9 vectors Σ = σ D = σ xx σ yy σ zz σ xy σ xz σ yz D x D y T D z, Z = ε E = ε xx ε yy ε zz ε xy ε xz ε yz E x E y T E z, (A.) and L is the symmetric global piezoelectric tensor, written as 9 9 matrix L = C e T e κ. (A.3) Dunn and Taya (993) have obtained analytical expressions for the Eshelby tensor, S, considering a transversely isotropic piezoelectric material, written in the grobal form (A.), and infinitely 39
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 informationPredicting 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 informationInternational 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 informationModule 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 informationMathematics and Mechanics of Solids. A systematic study on homogenization and the utility of circular simplified RVE
athematics and echanics of Solids A systematic study on homogenization and the utility of circular simplified RVE Journal: athematics and echanics of Solids anuscript ID S--0 anuscript Type: Original anuscript
More informationLecture 8. Stress Strain in Multi-dimension
Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]
More informationAuthors: 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 information16.20 HANDOUT #2 Fall, 2002 Review of General Elasticity
6.20 HANDOUT #2 Fall, 2002 Review of General Elasticity NOTATION REVIEW (e.g., for strain) Engineering Contracted Engineering Tensor Tensor ε x = ε = ε xx = ε ε y = ε 2 = ε yy = ε 22 ε z = ε 3 = ε zz =
More informationPEAT 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 informationME 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 informationA 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 informationA Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core
Commun. Theor. Phys. 56 774 778 Vol. 56, No. 4, October 5, A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core JIANG
More informationNumerical analyses of cement-based piezoelectric smart composites
Numerical analyses of cement-based piezoelectric smart composites *Jan Sladek 1, Pavol Novak 2, Peter L. Bishay 3, and Vladimir Sladek 1 1 Institute of Construction and Architecture, Slovak Academy of
More informationDYNAMIC 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 informationFig. 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 informationP 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 informationFinite 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 informationStiffness 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 information3 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 informationMicromechanical 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 informationBasic concepts to start Mechanics of Materials
Basic concepts to start Mechanics of Materials Georges Cailletaud Centre des Matériaux Ecole des Mines de Paris/CNRS Notations Notations (maths) (1/2) A vector v (element of a vectorial space) can be seen
More informationMacroscopic 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 informationELECTROMECHANICAL RESPONSE OF PIEZOELECTRIC FOAMS
18 TH INTRNATIONAL CONFRNC ON COMPOSIT MATRIALS LCTROMCHANICAL RSPONS OF PIZOLCTRIC FOAMS K.S. Challagulla 1 *, T.A. Venkatesh 1 School of ngineering, Laurentian University, Sudbury, Canada, Department
More informationMICROMECHANICS 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* I / 2c (2 DI )/(2 + D! AXIAL SHEAR MODULUS OF A FIBER-REINFORCED COMPOSITE WITH RANDOM FIBER CROSS-SECTIONS. I and 2. I ntn. J. Math. & Math. Si.
I ntn. J. Math. & Math. Si. Vol. 5 No. (198) 393-401 33 AXIAL SHEAR MODULUS OF A FIBER-REINFORCED COMPOSITE WITH RANDOM FIBER CROSS-SECTIONS $.K. BOSE Department of Mathematics Regional Engineering College
More informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationDEVELOPMENT 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 informationContinuum Mechanics. Continuum Mechanics and Constitutive Equations
Continuum Mechanics Continuum Mechanics and Constitutive Equations Continuum mechanics pertains to the description of mechanical behavior of materials under the assumption that the material is a uniform
More informationMechanics PhD Preliminary Spring 2017
Mechanics PhD Preliminary Spring 2017 1. (10 points) Consider a body Ω that is assembled by gluing together two separate bodies along a flat interface. The normal vector to the interface is given by n
More informationCellular 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 informationConstitutive Equations
Constitutive quations David Roylance Department of Materials Science and ngineering Massachusetts Institute of Technology Cambridge, MA 0239 October 4, 2000 Introduction The modules on kinematics (Module
More informationEFFECTIVE 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 informationDepartment 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 informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationOn the uniformity of stresses inside an inhomogeneity of arbitrary shape
IMA Journal of Applied Mathematics 23) 68, 299 311 On the uniformity of stresses inside an inhomogeneity of arbitrary shape Y. A. ANTIPOV Department of Mathematics, Louisiana State University, Baton ouge,
More informationInternational 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 informationModelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases
Modelling and numerical simulation of the wrinkling evolution for thermo-mechanical loading cases Georg Haasemann Conrad Kloß 1 AIMCAL Conference 2016 MOTIVATION Wrinkles in web handling system Loss of
More informationFINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS
FINITE ELEMENT ANALYSIS OF COMPOSITE MATERIALS Ever J. Barbero Department of Mechanical and Aerospace Engineering West Virginia University USA CRC Press Taylor &.Francis Group Boca Raton London New York
More informationPLEASE 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 informationSoufiane 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 informationAnalytical Solution for Electro-mechanical Behavior of Piezoelectric Rotating Shaft Reinforced by BNNTs Under Nonaxisymmetric
Journal of Solid Mechanics Vol. No. () pp. 339-35 Analytical Solution for Electro-mechanical Behavior of Piezoelectric Rotating Shaft Reinforced by BNNTs Under Nonaxisymmetric Internal Pressure A. Ghorbanpour
More informationA 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 informationThree-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 informationContinuum 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 information3D Elasticity Theory
3D lasticity Theory Many structural analysis problems are analysed using the theory of elasticity in which Hooke s law is used to enforce proportionality between stress and strain at any deformation level.
More information16.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 informationPROPERTY CALCULATION SYSTEM FOR INJECTION MOLDING AND COMPRESSION MOLDING OF FIBER-FILLED POLYMER COMPOSITES
PROPERTY CALCULATION SYSTEM FOR INJECTION MOLDING AND COMPRESSION MOLDING OF FIBER-FILLED POLYMER COMPOSITES X. Jin *, J. Wang, S. Han Autodesk Inc., 353 N. Triphammer Rd., Ithaca, NY, 485, U.S.A. Corning
More information3D 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 informationMICROMECHANICAL 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 informationHomogenization 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 informationVIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS
Journal of Engineering Science and Technology Vol. 12, No. 12 (217) 3398-3411 School of Engineering, Taylor s University VIBRATION CONTROL OF RECTANGULAR CROSS-PLY FRP PLATES USING PZT MATERIALS DILEEP
More informationMATERIAL 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 informationDistributed: Wednesday, March 17, 2004
MASSACHUSETTS INSTITUTE OF TECHNOLOGY DEPARTMENT OF MECHANICAL ENGINEERING CAMBRIDGE, MASSACHUSETTS 019.00 MECHANICS AND MATERIALS II QUIZ I SOLUTIONS Distributed: Wednesday, March 17, 004 This quiz consists
More informationAspect 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 informationEnhancement of magnetoelectric coupling in multiferroic composites via FEM simulation
Enhancement of magnetoelectric coupling in multiferroic composites via FEM simulation *Artjom Avakian 1), Andreas Ricoeur 2) 1), 2) Institute of Mechanics, University of Kassel, Kassel 34125, Germany 1)
More informationGeneralized Plane Piezoelectric Problem: Application to Heterostructure Nanowires
Generalized Plane Piezoelectric Problem: Application to Heterostructure Nanowires H. T. Mengistu and A. Garcia-Cristobal Institute of Material Science, University of Valencia, P.O.Box 22085, E-46071, Valencia,
More informationwhere C f = A ρ g fluid capacitor But when squeezed, h (and hence P) may vary with time even though V does not. Seems to imply C f = C f (t)
ENERGY-STORING COUPLING BETWEEN DOMAINS MULTI-PORT ENERGY STORAGE ELEMENTS Context: examine limitations of some basic model elements. EXAMPLE: open fluid container with deformable walls P = ρ g h h = A
More informationTABLE OF CONTENTS. Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA
Mechanics of Composite Materials, Second Edition Autar K Kaw University of South Florida, Tampa, USA TABLE OF CONTENTS 1. INTRODUCTION TO COMPOSITE MATERIALS 1.1 Introduction... 1.2 Classification... 1.2.1
More informationModeling of Axisymmetric Waves in a Piezoelectric Circular Fiber Coated with Thin Film
Mechanics and Mechanical Engineering Vol. 21, No. 3 2017) 637 648 c Lodz University of Technology Modeling of Axisymmetric Waves in a Piezoelectric Circular Fiber Coated with Thin Film R. Selvamani Department
More information3D Compression Molding
Autodesk Simulation Moldflow Insight 2014 3D Compression Molding Executive summary In this work, the simulation results from a program developed for the three-dimensional analysis of compression molding
More informationElastic parameters prediction under dynamic loading based on the. unit cell of composites considering end constraint effect
Elastic parameters prediction under dynamic loading based on the unit cell of composites considering end constraint effect Wang Meng 1,, Fei Qingguo 1,, Zhang Peiwei 1, (1. Institute of Aerospace Machinery
More informationA short review of continuum mechanics
A short review of continuum mechanics Professor Anette M. Karlsson, Department of Mechanical ngineering, UD September, 006 This is a short and arbitrary review of continuum mechanics. Most of this material
More informationFundamentals 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 informationMoving 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 informationNumerical 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 informationPiezoelectric Bimorph Response with Imperfect Bonding Conditions
Copyright c 28 ICCES ICCES, vol.6, no.3, pp.5-56 Piezoelectric Bimorph Response with Imperfect Bonding Conditions Milazzo A., Alaimo A. and Benedetti I. Summary The effect of the finite stiffness bonding
More informationEVALUATION 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 informationCHAPTER OBJECTIVES CHAPTER OUTLINE. 4. Axial Load
CHAPTER OBJECTIVES Determine deformation of axially loaded members Develop a method to find support reactions when it cannot be determined from euilibrium euations Analyze the effects of thermal stress
More informationCOUPLING MOISTURE DIFFUSION AND INTERNAL MECHANICAL STATES IN POLYMERS AND COMPOSITES A THERMODYNAMICAL APPROACH
18 TH INTERNATIONAL CONFERENCE ON COMPOSITE MATERIALS COUPLING MOISTURE DIFFUSION AND INTERNAL MECHANICAL STATES IN POLYMERS AND COMPOSITES A B.E. Sar 1, S. Fréour 1, P. Davies 2 et F. Jacquemin 1 * 1
More informationOriginal Article Effective properties of a piezocomposite containing shape memory alloy and inert inclusions
Continuum Mech. Thermodyn. () : 87 Digital Object Identifier (DOI).7/s675 Original Article Effective properties of a piezocomposite containing shape memory alloy and inert inclusions B. Jiang, R. C. Batra
More informationMODELLING 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 informationStress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina
Chapter 2 Macromechanical Analysis of a Lamina Stress, Strain Stress strain relationships for different types of materials Stress strain relationships for a unidirectional/bidirectional lamina Islamic
More informationThe New Boundary Condition Effect on The Free Vibration Analysis of Micro-beams Based on The Modified Couple Stress Theory
International Research Journal of Applied and Basic Sciences 2015 Available online at www.irjabs.com ISSN 2251-838X / Vol, 9 (3): 274-279 Science Explorer Publications The New Boundary Condition Effect
More informationMultiscale analyses of the behaviour and damage of composite materials
Multiscale analyses of the behaviour and damage of composite materials Presented by Didier BAPTISTE ENSAM, LIM, UMR CNRS 8006 151 boulevard de l hôpital l 75013 PARIS, France Research works from: K.Derrien,
More informationCHAPTER 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 informationExercise: 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 informationMechanics of Materials
Mechanics of Materials 3 (1 47 6 Contents lists available at SciVerse ScienceDirect Mechanics of Materials journal homepage: www.elsevier.com/locate/mechmat Effective properties of a novel composite reinforced
More informationNUMERICAL MODELING ON ELECTRIC RESPONSE OF FIBRE- ORIENTATION OF COMPOSITES WITH PIEZOELECTRICITY
www.arpapress.com/volumes/vol15issue2/ijrras_15_2_08.pdf NUMRICAL MODLIN ON LCTRIC RSPONS OF FIBR- ORINTATION OF COMPOSITS WITH PIZOLCTRICITY Cheuk-Yu Lee, Qing H. Qin and uy Walpole Research School of
More informationIn this section, thermoelasticity is considered. By definition, the constitutive relations for Gradθ. This general case
Section.. Thermoelasticity In this section, thermoelasticity is considered. By definition, the constitutive relations for F, θ, Gradθ. This general case such a material depend only on the set of field
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationIntelligent Materials and their applications. Autors: Jana Pintea, Ionut Balan INCDIE ICPE-CA Bucharest, ROMANIA
Intelligent Materials and their applications Autors: Jana Pintea, Ionut Balan INCDIE ICPE-CA Bucharest, ROMANIA What are intelligent materials? The intelligence of a material is a characteristic difficult
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationFINITE ELEMENT AND EXPERIMENTAL STUDY OF NOVEL CONCEPT OF 3D FIBRE CELL STRUCTURE
FINITE ELEMENT AND EXPERIMENTAL STUDY OF NOVEL CONCEPT OF 3D FIBRE CELL STRUCTURE M. Růžička, V. Kulíšek 2, J. Had, O. Prejzek Department of Mechanics, Biomechanics and Mechatronics, Faculty of Mechanical
More informationMechanical 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 informationOptimal Location of an Active Segment of Magnetorheological Fluid Layer in a Sandwich Plate
Acta Montanistica Slovaca Ročník 16 (2011), číslo 1, 95-100 Optimal Location of an Active Segment of Magnetorheological Fluid Layer in a Sandwich Plate Jacek Snamina 1 Abstract: In the present study a
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Standard Solids and Fracture Fluids: Mechanical, Chemical Effects Effective Stress Dilatancy Hardening and Stability Mead, 1925
More informationLAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS
XXII. LAMINATION THEORY FOR THE STRENGTH OF FIBER COMPOSITE MATERIALS Introduction The lamination theory for the elastic stiffness of fiber composite materials is the backbone of the entire field, it holds
More information3.2 Hooke s law anisotropic elasticity Robert Hooke ( ) Most general relationship
3.2 Hooke s law anisotropic elasticity Robert Hooke (1635-1703) Most general relationship σ = C ε + C ε + C ε + C γ + C γ + C γ 11 12 yy 13 zz 14 xy 15 xz 16 yz σ = C ε + C ε + C ε + C γ + C γ + C γ yy
More informationComputational 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 informationMultipole moments. November 9, 2015
Multipole moments November 9, 5 The far field expansion Suppose we have a localized charge distribution, confined to a region near the origin with r < R. Then for values of r > R, the electric field must
More informationMECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS. Vinu Unnithan and J. N. Reddy
MECHANICS OF CARBON NANOTUBE BASED COMPOSITES WITH MOLECULAR DYNAMICS AND MORI TANAKA METHODS Vinu Unnithan and J. N. Reddy US-South American Workshop: Mechanics and Advanced Materials Research and Education
More informationComposite materials: mechanical properties
Composite materials: mechanical properties A computational lattice model describing scale effects @ nano-scale Luciano Colombo In collaboration with: Pier Luca Palla and Stefano Giordano Department of
More informationOn the Numerical Modelling of Orthotropic Large Strain Elastoplasticity
63 Advances in 63 On the Numerical Modelling of Orthotropic Large Strain Elastoplasticity I. Karsaj, C. Sansour and J. Soric Summary A constitutive model for orthotropic yield function at large strain
More informationMechanics of Biomaterials
Mechanics of Biomaterials Lecture 7 Presented by Andrian Sue AMME498/998 Semester, 206 The University of Sydney Slide Mechanics Models The University of Sydney Slide 2 Last Week Using motion to find forces
More informationNUMERICAL EVALUATION OF A TEFLON BASED PIEZOELECTRIC SENSOR EFFECTIVITY FOR THE MONITORING OF EARLY AGE COCRETE STRENGTHING
NUMERICAL EVALUATION OF A TEFLON BASED PIEZOELECTRIC SENSOR EFFECTIVITY FOR THE MONITORING OF EARLY AGE COCRETE STRENGTHING Evangelos V. Liarakos Postdoctoral researcher School of Architecture, Technical
More informationEffect of Thermal Stresses on the Failure Criteria of Fiber Composites
Effect of Thermal Stresses on the Failure Criteria of Fiber Composites Martin Leong * Institute of Mechanical Engineering Aalborg University, Aalborg, Denmark Bhavani V. Sankar Department of Mechanical
More informationChapter 2. Rubber Elasticity:
Chapter. Rubber Elasticity: The mechanical behavior of a rubber band, at first glance, might appear to be Hookean in that strain is close to 100% recoverable. However, the stress strain curve for a rubber
More informationComposites Design and Analysis. Stress Strain Relationship
Composites Design and Analysis Stress Strain Relationship Composite design and analysis Laminate Theory Manufacturing Methods Materials Composite Materials Design / Analysis Engineer Design Guidelines
More informationThe 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 informationMicromechanics Analysis of Fiber Reinforced Composite
Micromechanics Analysis of Fiber Reinforced Composite Paresh K Katariya Student, ME-CAD/CAM, Department of mechanical engineering, Shree S ad Vidya Mandal Institute of Technology, Bharuch, Gujarat, India
More information