Transient response of piezoelectric composites caused by the sudden formation of localized defects

AbstractA continuum model is presented which is capable of generating the transient electroelastic field in piezoelectric composites of periodic microstructure, caused by the sudden appearance of localized defects. These defects are simulated by associating to every one of the ten piezoelectric parameters of the constituents a distinct damage variable. This procedure enables the modeling of localized cracks, soft and stiff inclusions and cavities. As a result, the constitutive equations of the piezoelectric phases appear in a specific form that includes eigen-electromechanical field variables which represent these defects. The method of solution is based on the combination of two distinct approach. In the first one, the representative cell method is employed according to which the periodic composite, which is discretized into several cells, is reduced to a problem of a single cell in the discrete Fourier transform domain. The resulting coupled elastodynamic and electric equations, initial, boundary and interfacial conditions in the transform domain are solved by employing a wave propagation in piezoelectric composite analysis which forms the second approach. The method of solution is verified by comparison with an analytical solution for the transient response of a piezoelectric material with a semi-infinite mode III-crack. Several applications are presented for the sudden formation of cracks in homogeneous and layered piezoelectric materials which are subjected to various types of electromechanical loading, and for the sudden appearance of a cavity. The effect of electromechanical coupling on the dynamic response is discussed

Micromechanics of metal matrix composites using the Generalized Method of Cells model (GMC) user's guide

A user's guide for the program gmc.f is presented. The program is based on the generalized method of cells model (GMC) which is capable via a micromechanical analysis, of predicting the overall, inelastic behavior of unidirectional, multi-phase composites from the knowledge of the properties of the viscoplastic constituents. In particular, the program is sufficiently general to predict the response of unidirectional composites having variable fiber shapes and arrays

Micromechanical Prediction of the Effective Coefficients of Thermo-Piezoelectric Multiphase Composites

The micromechanical generalized method of cells model is employed for the prediction of the effective elastic, piezoelectric, dielectric, pyroelectric and thermal-expansion constants of multiphase composites with embedded piezoelectric materials. The predicted effective constants are compared with other micromechanical methods available in the literature and good agreements are obtained

Micromechanical analysis of thermo-inelastic multiphase short-fiber composites

A micromechanical formulation is presented for the prediction of the overall thermo-inelastic behavior of multiphase composites which consist of short fibers. The analysis is an extension of the generalized method of cells that was previously derived for inelastic composites with continuous fibers, and the reliability of which was critically examined in several situations. The resulting three dimensional formulation is extremely general, wherein the analysis of thermo-inelastic composites with continuous fibers as well as particulate and porous inelastic materials are merely special cases

A new and general formulation of the parametric HFGMC micromechanical method for two and three-dimensional multi-phase composites

AbstractThe recent two-dimensional (2D) parametric formulation of the high fidelity generalized method of cells (HFGMC) reported by the authors is generalized for the micromechanical analysis of three-dimensional (3D) multiphase composites with periodic microstructure. Arbitrary hexahedral subcell geometry is developed to discretize a triply periodic repeating unit-cell (RUC). Linear parametric-geometric mapping is employed to transform the arbitrary hexahedral subcell shapes from the physical space to an auxiliary orthogonal shape, where a complete quadratic displacement expansion is performed. Previously in the 2D case, additional three equations are needed in the form of average moments of equilibrium as a result of the inclusion of the bilinear terms. However, the present 3D parametric HFGMC formulation eliminates the need for such additional equations. This is achieved by expressing the coefficients of the full quadratic polynomial expansion of the subcell in terms of the side or face average-displacement vectors. The 2D parametric and orthogonal HFGMC are special cases of the present 3D formulation. The continuity of displacements and tractions, as well as the equilibrium equations, are imposed in the average (integral) sense as in the original HFGMC formulation. Each of the six sides (faces) of a subcell has an independent average displacement micro-variable vector which forms an energy-conjugate pair with the transformed average-traction vector. This allows generating symmetric stiffness matrices along with internal resisting vectors for the subcells which enhances the computational efficiency. The established new parametric 3D HFGMC equations are formulated and solution implementations are addressed. Several applications for triply periodic 3D composites are presented to demonstrate the general capability and varsity of the present parametric HFGMC method for refined micromechanical analysis by generating the spatial distributions of local stress fields. These applications include triply periodic composites with inclusions in the form of a cavity, spherical inclusion, ellipsoidal inclusion, and discontinuous aligned short fiber. A 3D repeating unit-cell for foam material composite is simulated

Nonlinear micromechanical formulation of the high fidelity generalized method of cells

AbstractThe recent High Fidelity Generalized Method of Cells (HFGMC) micromechnical modeling framework of multiphase composites is formulated in a new form which facilitates its computational efficiency that allows an effective multiscale material–structural analysis. Towards this goal, incremental and total formulations of the governing equations are derived. A new stress update computational method is established to solve for the nonlinear material constituents along with the micromechanical equations. The method is well-suited for multiaxial finite increments of applied average stress or strain fields. Explicit matrix form of the HFGMC model is presented which allows an immediate and convenient computer implementation of the offered method. In particular, the offered derivations provide for the residual field vector (error) in its incremental and total forms along with an explicit expression for the Jacobian matrix. This enables the efficient iterative computational implementation of the HFGMC as a stand alone. Furthermore, the new formulation of the HFGMC is used to generate a nested local-global nonlinear finite element analysis of composite materials and structures. Applications are presented to demonstrate the efficiency of the proposed approach. These include the behavior of multiphase composites with nonlinearly elastic, elastoplastic and viscoplastic constituents

A New and General Formulation of the Parametric HFGMC Micromechanical Method for Three-Dimensional Multi-Phase Composites

The recent two-dimensional (2-D) parametric formulation of the high fidelity generalized method of cells (HFGMC) reported by the authors is generalized for the micromechanical analysis of three-dimensional (3-D) multiphase composites with periodic microstructure. Arbitrary hexahedral subcell geometry is developed to discretize a triply periodic repeating unit-cell (RUC). Linear parametric-geometric mapping is employed to transform the arbitrary hexahedral subcell shapes from the physical space to an auxiliary orthogonal shape, where a complete quadratic displacement expansion is performed. Previously in the 2-D case, additional three equations are needed in the form of average moments of equilibrium as a result of the inclusion of the bilinear terms. However, the present 3-D parametric HFGMC formulation eliminates the need for such additional equations. This is achieved by expressing the coefficients of the full quadratic polynomial expansion of the subcell in terms of the side or face average-displacement vectors. The 2-D parametric and orthogonal HFGMC are special cases of the present 3-D formulation. The continuity of displacements and tractions, as well as the equilibrium equations, are imposed in the average (integral) sense as in the original HFGMC formulation. Each of the six sides (faces) of a subcell has an independent average displacement micro-variable vector which forms an energy-conjugate pair with the transformed average-traction vector. This allows generating symmetric stiffness matrices along with internal resisting vectors for the subcells which enhances the computational efficiency. The established new parametric 3-D HFGMC equations are formulated and solution implementations are addressed. Several applications for triply periodic 3-D composites are presented to demonstrate the general capability and varsity of the present parametric HFGMC method for refined micromechanical analysis by generating the spatial distributions of local stress fields. These applications include triply periodic composites with inclusions in the form of a cavity, spherical inclusion, ellipsoidal inclusion, discontinuous aligned short fiber. A 3-D repeating unit-cell for foam material composite is simulated

Determining Shear Stress Distribution in a Laminate

A "simplified shear solution" method approximates the through-thickness shear stress distribution within a composite laminate based on an extension of laminated beam theory. The method does not consider the solution of a particular boundary value problem; rather, it requires only knowledge of the global shear loading, geometry, and material properties of the laminate or panel. It is thus analogous to lamination theory in that ply-level stresses can be efficiently determined from global load resultants at a given location in a structure and used to evaluate the margin of safety on a ply-by-ply basis. The simplified shear solution stress distribution is zero at free surfaces, continuous at ply boundaries, and integrates to the applied shear load. The method has been incorporated within the HyperSizer commercial structural sizing software to improve its predictive capability for designing composite structures. The HyperSizer structural sizing software is used extensively by NASA to design composite structures. In the case of through-thickness shear loading on panels, HyperSizer previously included a basic, industry-standard, method for approximating the resulting shear stress distribution in sandwich panels. However, no such method was employed for solid laminate panels. The purpose of the innovation is to provide an approximation of the through-thickness shear stresses in a solid laminate given the through-thickness shear loads (Qx and Qy) on the panel. The method was needed for implementation within the HyperSizer structural sizing software so that the approximated ply-level shear stresses could be utilized in a failure theory to assess the adequacy of a panel design. The simplified shear solution method was developed based on extending and generalizing bi-material beam theory to plate-like structures. It is assumed that the through-thickness shear stresses arise due to local bending of the laminate induced by the through-thickness shear load, and by imposing equilibrium both vertically and horizontally, the through-thickness shear stress distribution can be calculated. The resulting shear stresses integrate to the applied shear load, are continuous at the ply interfaces, and are zero at the laminate-free surfaces. If both Qx and Qy shear loads are present, it is assumed that they act independently and that their effects can be superposed. The calculated shear stresses can be rotated within each ply to the principal material coordinates for use in a ply-level failure criterion. The novelty of the simplified shear solution method is its simplicity and the fact that it does not require solution of a particular boundary value problem. The advantages of the innovation are that an approximation of the though-thickness shear stress distribution can be quickly determined for any solid laminate or solid laminate region within a stiffened panel

Analysis of Fiber Clustering in Composite Materials Using High-Fidelity Multiscale Micromechanics

A new multiscale micromechanical approach is developed for the prediction of the behavior of fiber reinforced composites in presence of fiber clustering. The developed method is based on a coupled two-scale implementation of the High-Fidelity Generalized Method of Cells theory, wherein both the local and global scales are represented using this micromechanical method. Concentration tensors and effective constitutive equations are established on both scales and linked to establish the required coupling, thus providing the local fields throughout the composite as well as the global properties and effective nonlinear response. Two nondimensional parameters, in conjunction with actual composite micrographs, are used to characterize the clustering of fibers in the composite. Based on the predicted local fields, initial yield and damage envelopes are generated for various clustering parameters for a polymer matrix composite with both carbon and glass fibers. Nonlinear epoxy matrix behavior is also considered, with results in the form of effective nonlinear response curves, with varying fiber clustering and for two sets of nonlinear matrix parameters