Dual finite element methods in homogenization for elastic–plastic fibrous composite material

The problem of homogenization for a periodic, elastic–perfectly plastic, fiber reinforced, composite material is considered. The overall mechanical behavior of the material is described using the anisotropic model of elastic–plastic body with kinematic hardening. The appropriate initial–boundary value problem, set for one repeatable cell of the composite, is solved in order to find the effective constitutive relations. The cell problem is solved using the finite element method formulated in two dual forms: in displacements and in stresses. Stress functions are used in the latter formulation.     

Micromechanics-based elastic-damage analysis of laminated composite structures

Based on the micromechanics-based constitutive model, derived in our preceding work [Lee, H.K., Pyo, S.H., 2009. A 3D-damage model for fiber-reinforced brittle composites with microcracks and imperfect interfaces. Journal of Engineering Mechanics-ASCE, in press, doi:10.1061/(ASCE)EM.1943.7889.0000039.], incorporating a multi-level damage model and a continuum damage model, the overall elastic behavior and damage evolution of laminated composite structures are studied in detail. The constitutive model is implemented into the finite element program ABAQUS using a user-subroutine UMAT to solve boundary value problems of the composite structures. The validity of the implemented constitutive model is assured by comparing the predicted stress–strain curves with experimental data available in literature under uniaxial tension with various fiber orientations. The results show that the proposed micromechanics-based constitutive model accurately predict the overall elastic-damage behavior of laminated composite structures having different material compositions, thicknesses and boundary conditions.     

Analytical solution to a mixed boundary value elastic problem of a roller-guided panel of laminated composite

This study presents an analytical solution to elastic field in a roller-guided panel of symmetric cross-ply laminated composite material. The mixed boundary value two-dimensional plane stress elasticity problem is formulated in terms of a single displacement potential function. This reduces the problem to the solution of a single fourth order partial differential equation of equilibrium as the other equilibrium equation is satisfied automatically. The solution is obtained in terms of an infinite Fourier series. To present some numerical results, a panel of glass/epoxy laminated composite is considered and different components of stress and displacement at different sections of the panel are presented graphically. To justify the present analytical solution, it is compared with the finite element solution obtained by using the commercial software ANSYS. It is found that the two solutions agree well with each other. This ensures that the formulation developed in this study based on the displacement potential approach can be used to obtain analytical solution of an elastic field in structural elements of laminated composite under any mode of boundary conditions prescribed in terms of either stress, displacement or any combination of these.     

Axisymmetric Vibrations and Dissipative Heating of a Laminated Inelastic Disk

The coupled thermomechanical dynamic behavior of an inhomogeneous body is investigated for a partial case where a laminated inelastic disc is subject to forced axisymmetric vibrations and dissipative heating. The problem is solved in complete and approximate formulations. In the former case, the behavior of the material is described using generalized flow theory. In the latter case, the behavior of the material is characterized by complex moduli. The spatial distributions of the field quantities and the temperature– and amplitude–frequency characteristics of the disc are analyzed. The results are compared.     

Function theoretic solutions to problems of orthotropic elasticity

The plane strain problem for a two dimensional orthotropic elastic body is investigated. In particular analytic representations for the solution of the displacement boundary value problem and the stress boundary value problems are found. To this end, the Navier equations are reduced by means of composite transformations to normal form. These are the so-called equations for bianalytic function of the type (k). The generalized Cauchy integral formula for this function theory is used to obtain representation formulae. A simplified method to solve these problems by bianalytic function theory is given for certain situations of plane strain for an orthotropic elastic body. AMS (MOS): 35A20, 35CO5, 35G15, 35J55.Applied Mathematics Institute Technical Report No. 140A, July 1983.The work of this author was supported in part by grant no. DE-AC01-81ER-10967 from the Department of Energy.     

A continuum theory for viscoelastic composite beams

Summary A dynamical continuum theory is developed for laminated composite beams. Starting with an assumed displacement- and temperature field, the one-dimensional approximate theory is consistently constructed within the frame of the three-dimensional theory of linear, nonisothermal, anisotropic, coupled viscoelasticity. Each constituent of the beam may possess different constant thickness and mechanical properties. All dynamic interactions between the adjacent constituents are included. Further, the effects of transverse shear and normal strains and rotatory inertia as well as those of cross-sectional distortion are all taken into account. The resulting equations consist of the macroscopic beam equations of motion and heat conduction, the kinematical relations, the initial and boundary conditions and the constitutive equations, and they govern the extensional, flexural and torsional motions of laminated composite beams. The special cases of constituents which made of either isotropic thermoviscoelastic or anisotropic thermoelastic materials are discussed briefly.Supported by the Office of Naval Research.With 1 figure     

Homogenization of nonlinear problems in the mechanics of composites

A further development of the homogenization method is proposed to solve the physically nonlinear equilibrium problems for the laminated plates or the plates made of functionally graded materials. In the linear case, according to this method, the corresponding solution is a superposition of the solution to the global problem in the entire domain and the solution to the local problem in a representative domain, e.g., in a periodicity cell. In the nonlinear case, such a superposition is not valid, which complicates the application of the homogenization method. In order to eliminate this difficulty, it is possible to combine the homogenization method and the linearization method when solving a boundary value problem or a variational problem. In the mechanics of deformable solids, the constitutive relations can be considered as equations with respect to velocities or the stress and strain differentials in time or in the loading parameter. When these equations are linear with respect to velocities, it is possible to use the homogenization method. In this paper such an approach is illustrated by the example of a symmetric laminated plate bent under a uniformly distributed time-dependent load.     

Couple stress theory for solids

By relying on the definition of admissible boundary conditions, the principle of virtual work and some kinematical considerations, we establish the skew-symmetric character of the couple-stress tensor in size-dependent continuum representations of matter. This fundamental result, which is independent of the material behavior, resolves all difficulties in developing a consistent couple stress theory. We then develop the corresponding size-dependent theory of small deformations in elastic bodies, including the energy and constitutive relations, displacement formulations, the uniqueness theorem for the corresponding boundary value problem and the reciprocal theorem for linear elasticity theory. Next, we consider the more restrictive case of isotropic materials and present general solutions for two-dimensional problems based on stress functions and for problems of anti-plane deformation. Finally, we examine several boundary value problems within this consistent size-dependent theory of elasticity.     

A note on incremental equations for a new class of constitutive relations for elastic bodies

Some new classes of constitutive relations for elastic bodies have been proposed in the literature, wherein the stresses and strains are obtained from implicit constitutive relations. A special case of the above relations corresponds to a class of constitutive equations where the linearized strain tensor is given as a nonlinear function of the stresses. For such constitutive equations we consider the problem of decomposing the stresses into two parts: one corresponds to a time-independent solution of the boundary value problem, plus a small (in comparison with the above) time-dependent stress tensor. The effect of this initial time-independent stress in the propagation of a small wave motion is studied for an infinite medium.     

Stress-resultant plasticity theories for composite laminated plates

Constitutive laws are presented for the inelastic analysis of laminated composite plates. The implications of using an elastoplastic theory, applied in a stress-resultant formulation, are discussed and investigated. Two different stress-resultant plasticity theories are proposed, both of which overlook the matrix and fiber inelastic behavior and describe the inelastic response of the laminate as a function of overall laminate properties. Results from numerical experiments with the proposed models are compared with results obtained using a micromechanical elastoplastic composite constitutive model.     

Accurate equations for laminated composite deep thick shells

This paper derives accurate equations of elastic deformation for laminated composite deep, thick shells. The equations include shells with a pre-twist and accurate force and moment resultants which are considerably different than those used for plates. This is due to the fact that the stresses over the thickness of the shell have to be integrated on a trapezoidal-like cross-section of a shell element to obtain the stress resultants. Numerical results are obtained and showed that accurate stress resultants are needed for laminated composite deep thick shells, especially if the curvature is not spherical. A consistent set of equations of motion, energy functionals and boundary conditions are also derived. These may be used in obtaining exact solutions or approximate ones like the Ritz or finite element methods.     

Secondary torsional moment deformation effect in inelastic nonuniform torsion of bars of doubly symmetric cross section by BEM

In this paper a boundary element method is developed for the inelastic nonuniform torsional problem of simply or multiply connected prismatic bars of arbitrarily shaped doubly symmetric cross section, taking into account the secondary torsional moment deformation effect. The bar is subjected to arbitrarily distributed or concentrated torsional loading along its length, while its edges are subjected to the most general torsional boundary conditions. A displacement based formulation is developed and inelastic redistribution is modeled through a distributed plasticity model exploiting three dimensional material constitutive laws and numerical integration over the cross sections. An incremental–iterative solution strategy is adopted to resolve the elastic and plastic part of stress resultants along with an efficient iterative process to integrate the inelastic rate equations. The one dimensional primary angle of twist per unit length, a two dimensional secondary warping function and a scalar torsional shear correction factor are employed to account for the secondary torsional moment deformation effect. The latter is computed employing an energy approach under elastic conditions. Three boundary value problems with respect to (i) the primary warping function, (ii) the secondary warping one and (iii) the total angle of twist coupled with its primary part per unit length are formulated and numerically solved employing the boundary element method. Domain discretization is required only for the third problem, while shear locking is avoided through the developed numerical technique. Numerical results are worked out to illustrate the method, demonstrate its efficiency and wherever possible its accuracy.     

Damage analysis and dynamic response of elasto-plastic laminated composite shallow spherical shell under low velocity impact

Based on the elasto-plastic mechanics, the damage analysis and dynamic response of an elasto-plastic laminated composite shallow spherical shell under low velocity impact are carried out in this paper. Firstly, a yielding criterion related to spherical tensor of stress is proposed to model the mixed hardening orthotropic material, and accordingly an incremental elasto-plastic damage constitutive relation for the laminated shallow spherical shell is founded when a strain-based Hashin failure criterion is applied to assess the damage initiation and propagation. Secondly, using the presented constitutive relations and the classical nonlinear shell theory, a series of incremental nonlinear motion equations of orthotropic moderately thick laminated shallow spherical shell are obtained. The questions are solved by using the orthogonal collocation point method, Newmark method and iterative method synthetically. Finally, a modified elasto-plastic contact law is developed to determine the normal contact force and the effect of damage, geometrical parameters, elasto-plastic contact and boundary conditions on the contact force and the dynamic response of the structure under low velocity impact are investigated.     

Analytical solution for the cylindrical bending vibration of piezoelectric composite plates

An analytical solution for the cylindrical bending vibrations of linear piezoelectric laminated plates is obtained by extending the Stroh formalism to the generalized plane strain vibrations of piezoelectric materials. The laminated plate consists of homogeneous elastic or piezoelectric laminae of arbitrary thickness and width. Fourier basis functions for the mechanical displacements and electric potential that identically satisfy the equations of motion and the charge equation of electrostatics are used to solve boundary value problems via the superposition principle. The coefficients in the infinite series solution are determined from the boundary conditions at the edges and continuity conditions at the interfaces between laminae, which are satisfied in the sense of Fourier series. The formulation admits different boundary conditions at the edges of the laminated piezoelectric composite plate. Results for laminated elastic plates with either distributed or segmented piezoelectric actuators are presented for different sets of boundary conditions at the edges.     

Parameters of the Edge Effect as a Function of the Mechanical Characteristics of a Composite with an Interfacial Crack

The paper studies the dependence of edge effects in a laminated composite of regular structure with a periodic system of symmetric interfacial cracks on the mechanical characteristics of the composite components. Loading the composite induces constant strain in the reinforcement direction. The problem is solved in exact formulation using the linear elastic equations for piecewise-homogeneous media and a criterion for quantitative evaluation of edge effects. An approximate solution of the problem is found by the mesh approach. A difference scheme for the mixed problem under consideration is derived. The edge-effect zones for normal stresses are constructed. The maximum values of the normal stress concentration factors are determined     

Elastic waves in an electro-microelastic solid

The present study is concerned with the wave propagation in an electro-microelastic solid. The reflection phenomenon of plane elastic waves from a stress free plane boundary of an electro-microelastic solid half-space is studied. The condition and the range of frequency for the existence of elastic waves in an infinite electro-microelastic body are investigated. The constitutive relations and the field equations for an electro-microelastic solid are stemmed from the Eringen’s theory of microstretch elasticity with electromagnetic interactions. Amplitude ratios and energy ratios of various reflected waves are presented when an elastic wave is made incident obliquely at the stress free plane boundary of an electro-microelastic solid half-space. It has been verified that there is no dissipation of energy at the boundary surface during reflection. Numerical computations are performed for a specific model to calculate the phase speeds, amplitude ratios and energy ratios, and the results obtained are depicted graphically. The effect of elastic parameter corresponding to micro-stretch is noticed on reflection coefficients, in particular. Results of Parfitt and Eringen [Parfitt, V.R., Eringen, A.C., 1969. Reflection of plane waves from a flat boundary of a micropolar elastic half-space. J. Acoust. Soc. Am. 45, 1258–1272] have also been reduced as a special case from the present formulation.     

A model of layered prismatic shells

The present paper is devoted to a model for elastic layered prismatic shells which is constructed by means of a suggested in the paper approach which essentially differs from the known approaches for constructing models of laminated structures. Using Vekua’s dimension reduction method after appropriate modifications, hierarchical models for elastic layered prismatic shells are constructed. We get coupled governing systems for the whole structure in the projection of the structure. The advantage of this model consists in the fact that we solve boundary value problems separately for each ply. In addition, beginning with the second ply, we use a solution of a boundary value problem of the preceding ply. We indicate ways of investigating boundary value problems for the governing systems. For the sake of simplicity, we consider the case of two plies, in the zeroth approximation. However, we also make remarks concerning the cases when either the number of plies is more than two or higher-order approximations (hierarchical models) should be applied. As an example, we consider a special case of deformation and solve the corresponding boundary value problem in the explicit form.     

Stresses in rotating heterogeneous viscoelastic composite cylinders with variable thickness

An analytical solution is presented for the rotation problem of a two-layer composite elastic cylinder under a plane strain assumption. The external cylinder has variable-thickness formulation, and is made of a heterogeneous orthotropic material. It contains a fiber-reinforced viscoelastic homogeneous isotropic solid core of uniform thickness. The thickness and elastic properties of the external cylinder are taken as power functions of the radial direction. By the boundary and continuity conditions, the radial displacement and stresses for the rotating composite cylinder are determined. The effective moduli and Illyushin’s approximation methods are used to obtain the viscoelastic solution to the problem. The effects of heterogeneity, thickness variation, constitutive, time parameters on the radial displacement, and stresses are investigated.     

Micromechanics-based constitutive modeling for unidirectional laminated composites

A micromechanics-based constitutive model is developed to predict the effective mechanical behavior of unidirectional laminated composites. A newly developed Eshelby’s tensor for an infinite circular cylindrical inclusion [Cheng, Z.Q., Batra, R.C., 1999. Exact Eshelby tensor for a dynamic circular cylindrical inclusion. J. Appl. Mech. 66, 563–565] is adopted to model the unidirectional fibers and is incorporated into the micromechanical framework. The progressive loss of strength resulting from the partial fiber debonding and the nucleation of microcracks is incorporated into the constitutive model. To validate the proposed model, the predicted effective stiffness of transversely isotropic composites under far field loading conditions is compared with analytical solutions. The constitutive model incorporating the damage models is then implemented into a finite element code to numerically characterize the elastic behavior of laminated composites. Finally, the present predictions on the stress–strain behavior of laminated composite plate containing an open hole is compared with experimental data to verify the predictive capability of the model.