Three-dimensional microscopic interlaminar analysis of cross-ply laminates based on a homogenization theory

In this study, a method for three-dimensional microscopic interlaminar analysis of cross-ply laminates is developed based on a homogenization theory to analyze microscopic interactions between unidirectional long fiber-reinforced laminae. For this, a unit cell of a cross-ply laminate, which includes interlaminar areas, is defined under the assumption that each lamina in the laminate has a transversely square fiber array. Then, showing that the laminate has a point-symmetric internal structure, the symmetry is utilized to introduce half of the unit cell as the domain of analysis. Moreover, the domain of analysis is divided into substructures using a substructure method combined with the homogenization theory, significantly reducing the computational costs. The present method is then applied to the analysis of interlaminar stress distributions in a carbon fiber/epoxy cross-ply laminate subjected to in-plane uniaxial tension. It is shown that microscopic shear stress noticeably occurs at the interface between the 0°- and 90°-plies. It is also shown that the microscopic interaction between the two plies is observed only in the vicinity of the interface.     

Thermodynamic driving force in ferroelectric crystals with a rank-2 laminated domain pattern, and a study of enhanced electrostriction

The thermodynamic driving force for domain growth in a rank-2 laminated ferroelectric crystal is derived in this article, and we used it, together with a homogenization theory, to study the issue of enhanced electrostrictive actuation recently reported by Burcsu et al. [2004. Large electrostrictive actuation of barium titanate single crystals. J. Mech. Phys. Solids 52, 823-846]. We derived this force from the reduction of Gibbs free energy with respect to the increase of domain concentration. It is shown that both the free energy and the thermodynamic force consist of three parts: the first arises from the difference in M₀ and M₁, the linear electromechanical compliances of the parent and product domains, respectively, at a given level of applied stress and electric field, the second stems from the electromechanical work associated with the change of spontaneous strain and spontaneous polarization during domain switch, and the third from the internal energy due to the distribution of polarizations strain and electric polarization inside the crystal. We prove that the first term is substantially lower than the second one, and the third one is identically zero with compatible domain pattern. The second one is, however, not exactly equal to the commonly written sum of the products of stress with strain, and electric field with polarization during switch, unless both domains have identical moduli in the common global axes. We also show that, with compatible domain patterns and when M₁=M₀, this driving force is identical to Eshelby's driving force acting on a flat interface due to the jump of energy-momentum tensor. Applications of the theory to a BaTiO₃ crystal subjected to a fixed axial compression and decreasing electric field from the [0 0 1] state reveal that the crystal undergoes a three-stage switching process: (i) the 0→90° switch to form a rank-1 laminate, (ii) the 0→180° switch inside the 0° domain to form a laminate I with a concurrent 90°→−90° switch inside the 90° domain to form laminate II, creating a rank-2-laminated domain pattern, and (iii) finally the 90→180° switch. It is the exchange of stability between the 0, 90°, and 180° domains under compression and electric field that is the origin of the enhanced actuation. We illustrate these intrinsic features by showing the evolution of these domains, and demonstrate how the reported large actuation strain can be attained with a rank-2 laminate.     

The influence of a compressive stress on the nonlinear response of ferroelectric crystals

The origin of nonlinearity in a ferroelectric crystal is domain reorientation, and such a process can be affected by the presence of a compressive stress. In this article we examine how a superimposed compression affects the evolution of new domain and how it changes the shape of the hysteresis loop. We start out by considering the thermodynamic driving force for domain reorientation, and then use a dual-phase homogenization theory to calculate the overall response. To uncover the influence of a compressive stress, the theory is used to calculate the hysteresis loop between the electric displacement D and the electric field E of a BaTiO₃ crystal, first without and then with a compression, using a two-consecutive 90° switch model (i.e. 0° → 90° → 180°). It is found that, from the initial 0° position, the compressive stress will increase the thermodynamic driving force and promote an earlier onset of the 90° domain, but its presence will cause a significant delay for the reorientation process to pass through the intermediate 90° state in route to its final 180° configuration. The D vs. E loop then exhibits a more round shape and a lesser steep slope near the coercive field. The delayed passage and more rounded shape are found to be consistent with a recent experimental observation [Burcsu et al., 2004. J. Mech. Phys. Solids 52, 823–846].     

Multi-level modeling of woven glass/epoxy composite for multilayer printed circuit board applications

The objective of this paper is to develop a hybrid homogenization method to predict the elastic properties of a common woven glass/epoxy composite substrate for multilayer circuit board applications. Comprehensive high resolution 3D finite element (FE) models of a quarter of the repeated unit cell (RUC) for the woven glass/epoxy composite were developed based on different micromechanical schemes. . Specifically, four different micromechanics schemes were investigated: self-consistent, Mori–Tanaka, three-phase approach and composite cylinder assemblage (CCA). The element based strain concentration matrices were determined and used to obtain the homogenized woven glass/epoxy composite properties via a specially developed MATLAB code. Attention was further devoted to the predictions of the homogenized elastic moduli of the multilayer printed circuit board (PCB). The results from our simulations, based on Mori–Tanaka and CCA, are in good agreement with existing experimental results, indicating that the newly proposed homogenization scheme can be used as a design tool to predict the overall properties of woven composite materials typically used in multilayer PCB applications.     

Stress resultants and moments around holes in unsymmetrical composite laminates subjected to remote uniform loading

Stress resultants and moments around the holes of infinite unsymmetrical composite laminates under remote uniform loading are analyzed using the complex potential approach developed by Chen and Shen [Mech. Res. Commun. 28 (4) (2001) 423; 28 (5) (2001) 513]. In the analysis, the authors’ previous work has been modified to determine the complex constants B_k appeared in the complex potential functions of Chen and Shen [Mech. Res. Commun. 28 (5) (2001) 513]. Herein B_k are the unknown quantities related to the remote loading conditions. The effect of bending extension coupling has been discussed for several laminates such as [0/90]_T, [45/−45]_T and [0/90/45/−45]_nT. Results indicates that the coupling between bending and extension significantly affects the stresses around the hole of an unsymmetrical laminate containing one sub-laminate only. Such the effect rapidly decreases as the number of the sub-laminates increases. The results confirm the validity of the developed complex potential approach.     

Cyclic behavior of unidirectional and cross-ply titanium matrix composites

Relatively simple and efficient micromechanical models are used to obtain the uniaxial response of SCS-6/Timetal 21S with [0]₄ and [0/90]_s laminates when subjected to isothermal and thermomechanical fatigue (TMF) loadings. Features of the modeling that are required to obtain the accurate deformation behavior for this class of materials under these loadings are highlighted. To this end, a comparison is made between the concentric cylinder model and the uniaxial stress model for representing the [0] laminate. The axial stresses from the two models are very similar under mechanical loading. The greatest differences appear under thermal loading alone. The differences on the composite response between a time-independent elastic-plastic and a viscoplastic matrix constitutive model are also examined. The latter is based on the Bodner-Partom unified constitutive model. The [0/90] laminate is treated by adding a parallel element with smeared [90] ply properties to the [0] model and invoking axial strain compatibility as well as stress equilibrium. The proposed constitutive law for the [90] ply includes both matrix viscoplasticity and fiber/matrix separation damage and is based on damage mechanics concepts. The effect of cyclic frequency on TMF behavior is examined. The in-phase TMF life is shown to be very sensitive to frequency due to the relaxation of matrix stress and the attendant increase in fiber stress.     

On the use of crack-tip-opening displacement to predict the fracture strength of notched graphite/epoxy laminates

The fracture strength and crack-opening displacement of notched graphite/epoxy laminates were measured experimentally using the center-cracked tension-specimen geometry. Four replicate tests were conducted for a variety of laminate stacking sequences, thicknesses, and notch lengths. Most laminates exhibited extensive notch-tip damage prior to fracture. Values of crack-tip-opening displacement (CTOD) at fracture were estimated from values of crack-opening displacement measured at the crack center line. CTOD was independent of specimen crack length for the [0/±45/90] _s , [0/±45/90]_15s , [0/±45] _s , [0/±45/]_15s , and [0/90]_24s laminates. In addition, notched laminate strength was accurately predicted using a Dugdale-type model along with the estimated CTOD.Paper was presented at V International Congress on Experimental Mechanics held in Montreal, Quebec, Canada on June 10–15, 1984.     

Singular full-field stresses in composite laminates with open holes

A method of the superposition of a hybrid and displacement approximation was developed to provide the accurate stress fields in a multilayered composite laminate, including the singular neighborhood of the ply interface and the hole edge. Asymptotic analysis was used to derive the hybrid stress functions. The displacement approximation is based on the polynomial B-spline functions. The method provides the determination of the coefficient of the singular term along with convergent stress components including the singular regions. Reissner’s variational principle was employed. Simple [45/−45]_s and quasi-isotropic IM7/5250 [45/90/−45/0]_s laminates were analyzed. Uniaxial loading and residual stress calculation (quasi-isotropic laminate) were considered. A convergence study showed that accurate values of the coefficient of the singular term of the asymptotic stress expansion could be obtained with coarse out-of-plane and in-plane subdivisions. The interaction between the singular terms on the neighboring interfaces was found to be important for the convergence with coarse subdivisions. Converged transverse interlaminar stress components as a function of the distance from the hole edge, were shown for all examples.     

BEM SOLUTION FOR THE PROBLEM OF FLUX OF A MULTICOMPONENT MIXTURE OF GASES OUT OF A MULTILAYER LANDFILL

A two-dimensional numerical model for convection–diffusion flow of a multigas mixture through a multilayer porous medium was developed with the aim to be used for evaluation of emissions of gases from landfills. The proposed model is based on the boundary element–dual reciprocity method. Time-independent one-dimensional analytical solutions for a multilayer domain were found for the cases of a single gas and a two-gas mixture and used to verify the accuracy of the model. Although the proposed technique is a simple one, consisting only of boundary integrals, it was found that the technique can be applied with satisfactory accuracy to the problem at which it was initially aimed.     

Motion of an incompressible conducting liquid in a strong self-consistent magnetic field

One of the common regimes of operation of many laboratory and industrial magnetohydrodynamic (MHD) devices using liquid metals as working medium is the regime for which the Alfvén number A, the ratio of the magnetic and kinetic energy densities, appreciably exceeds unity. For example, for a typical MHD device [1] with characteristic length 0.1 m of the working region, velocity 1 m/sec of the medium, and magnetic induction 1 T (the medium is molten sodium at temperature 330°C) the Alfvén number is A - 900. To simplify the investigation of the processes in such devices, one can use the approximation of a strong magnetic field proposed by Somov and Syrovatskii [2] to describe certain types of hydrodynamic flows of a dissipationless plasma in a magnetic field. In the present paper, the approach to the analysis of the self-consistent magnetohydrodynamic problem in this asymptotic approximation is extended to the case of an incompressible liquid with finite conductivity. A study is made of the closed reduced system of MHD equations obtained from the complete model in the zeroth order in the small parameter A^–1, in which the magnetic field is a force-free field. An investigation is made of the free diffusion of force-free magnetic field with constant coefficient a of proportionality between the current density and the magnetic induction in a spatially unbounded liquid, and the kinematic properties of a velocity field of the liquid in which the force-free nature of the magnetic field is maintained during the damping process are determined. It is shown that the complete class of such velocity fields is represented by the group of rigid-body motions of the liquid.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 3–9, January–February, 1991.     

Numerical simulation of large amplitude oscillatory shear of a high-density polyethylene melt using the MSF model

We study the flow response in large amplitude oscillatory shear of the molecular stress function (MSF) model that has recently been proposed by Wagner et al. [M.H. Wagner, P. Rubio, H. Bastian, The molecular stress function model for polydisperse polymer melts with dissipative convective constraint release, J. Rheol. 45 (2001) 1387–1412]. The MSF model is derived from molecular theory and has only two parameters to describe the non-linear material response. The model predictions are analysed in both the frequency and time domain. It shows good agreement with experimental data for a linear high-density polyethylene melt. At low and medium strains, MSF model predictions are in excellent agreement with experimental data and predictions of a six-mode Giesekus model which has six parameters to describe the non-linear material response. At medium strains, the basic Doi–Edwards model, which has no non-linear parameters, already underpredicts the data. At high strains, the MSF model predictions agree slightly better with the experimental data than the Giesekus model. Surprisingly, however, it is the Doi–Edwards model that shows excellent agreement with experimental data at high strains. For the linear melt we consider, it outperforms the models that have non-linear parameters, both in the time and frequency domain.     

Forced oscillations of pressure in tubes of active material

A study is made in the quasione-dimensional inertialess approximation of the axisymmetric flow of a Newtonian fluid in a tube of finite length made of a nonlinear active material with the capability of reducing deformations in response to an increase in tensile stresses [1, 2]. A study is made of the influence of the frequency and amplitude of forced oscillations of pressure at the entrance of the tube on its flow rate characteristics and on the behavior of the tube, depending on its length and certain rheological parameters. The first attempts at a study within the framework of this model of flow for unsteady conditions at the ends of the tube and in the ambient medium are described in [3, 4]. A general solution of this problem for external periodic disturbances of low amplitude is constructed in [5]. The present study gives an analysis of certain results of the numerical solution of an analogous problem for a wide range of variations in the frequency and amplitude of the pressure oscillations at the entrance to the tube.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 88–90, March–April, 1985.     

Parametric excitation of waves at an interface

Experiments on the parametric excitation of waves at a fluid interface show a strong disagreement with theoretical results [1–3], since the latter do not take into account the influence of the second medium. This proves to be especially important at low frequencies. Thus, for a water-air interface with an excitation frequency = 60 sec^–1 the contribution amounts to 10%,and with = 30 sec^–1, even 20%. In this paper the stability of the interface of two viscous, incompressible fluids of finite depth in a variable gravity field is considered. The problem is put in the linear form by making an expansion with respect to the small viscosity and is solved by taking the Laplace transform with respect to time. A second-order integrodifferential equation with periodic coefficients is obtained for the deviation of the interface from the equilibrium position; its solution is sought by the method of averaging [4]. It is shown that the presence of the second fluid significantly raises the threshold of instability.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 167–170, March–April, 1977.     

Fracture of composite laminates containing cracks due to shear loading

The responses of graphite/epoxy [0/90/±45] _s , [±45]_2s , [0/90]_2s and [0/±45]_2s composite laminates with and without center cracks were studied under shear loading using the three-rail shear test. The shear stress/strain relationship, the failure mechanisms and the notched strength were analyzed. Substantial amounts of local buckling were observed in some of the laminates. The present paper shows that shear modulus can be determined accurately using the three-rail shear test with proper interpretation of data. Using the minimum strength model, only one characteristic length was needed to predict accurately the notched strength of a composite laminate under shear and tensile loadings.Paper was presented at the 1987 SEM Spring Conference on Experimental Mechanics held in Houston, TX on June 14–19.     

一种新的宏、细观变形全场测量方法

基于光学干涉和衍射理论,提出了一种新的用于宏、细观面内变形测量的光学方法。利用三束准直相干光之间的干涉,分别得到面内位移U V和U-V的两幅耦合条纹图。借助于FFT相移技术,直接计算出全场的U和V分布。本实验系统中包含有长距离显微镜,三光束系统,数字图像采集和处理系统等。它具有光路简单,避开条纹级次的确定,直接计算得到U、V场等优点。将它用于热塑性复合材料 56层板AS4/PEEK[0/±45/90]7S三点弯曲梁位移场的宏、细观测量。得到了层间的变形规律。本文还利用三维各向异性弹性有限元法对AS4/PEEK[0/±45/90]7S三点弯曲梁进行了数值计算,得到了自由表面的位移场和应变场的数值解。定量分析了层间的变形传递规律,有助于进一步分析AS4/PEEK[0/±45/90]7S层合板宏、细观变形。     

Dispersion of a filtration flow

In this paper we shall consider the transport of a dynamically neutral impurity in a porous medium containing random inhomogeneities. The original versions of the equations for the mean impurity concentration [1, 2] were based on the hyphothesis that the random motions obeyed the Markov principle, use being made of the diffusion equations of A. N. Kolmogorov. Later [3, 4] the method of perturbations was used to study the complete system of equations for the impurity concentration and random filtration velocity in the case of a constant, nonrandom porosity; after an averaging process this yields a generalized equation for the average concentration. In the limiting cases of small- and large-scale inhomogeneities in the permeability of the medium, the basic integrodifferential equation may be, respectively, reduced to parabolic and hyperbolic equations of the second order. In the present analysis we shall use the perturbation method to study the transport of an impurity by a flow when the filtration velocity of the latter fluctuates around inhomogeneities in the permeability field, the porosity of the medium in which the flow is taking place also constituting a random field, correlating with the field of permeability. We shall derive equations for the average concentration and should formulate the corresponding boundary-value problems for these equations; we shall also calculate the components of the dispersion tensor and shall consider the equilibrium sorption of an impurity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 65–69, July–August, 1976.The author is grateful to A. I. Shnirel'man for useful discussions.     

均匀化方法在粘弹性多层复合材料中的应用

主要研究了由各向同性线弹性加强体和各向同性线粘弹性基体组成的多层复合材料的问题,在已有的线弹性多层材料的均匀化方法的基础上,应用弹性一粘弹性对应原理,在Carson域中求解粘弹性多层材料的问题。通过Burgers模型表示线粘弹性基体材料,反演得到了多层材料的有效松弛模量和有效泊松比在时间域中的表达式,并且与实验结果和其他结果进行了比较。     

Experimental and analytical comparisons of failure in thermoplastic composite laminates

Failure characteristics of Gr/PEEK were studied, using an experimental investigation and a fully nonlinear ply-by-ply finite-element technique. The stacking sequence of the laminates (with centrally located holes) investigated were: 0, 90, ±45 deg, (0/45/90/−45 deg)_2s and (0±45/90 deg)_2s. The [0 deg] laminate failure was characterized by splitting at the extremities of the hole and along the fibers. The [90 deg] laminates failed in the transverse direction, whereas the [±45 deg] laminates exhibited considerable elongation to failure. In the case of the quasi-isotropic laminates, the failure progression appeared to be due dominantly to matrix cracking followed by fiber failure. Analytical predictions of the failure process showed reasonably good correlation with the experimentally determined data.     

Reflection from a semi-infinite stack of layers using homogenization

A canonical scattering problem is that of a plane wave incident upon a periodic layered medium. Our aim here is to replace the periodic medium by a homogenized counterpart and then to investigate whether this captures the reflection and transmission behaviour accurately at potentially high frequencies.We develop a model based upon high frequency homogenization and compare the reflection coefficients and full fields with the exact solution. For some material properties it is shown that the asymptotic behaviour of the dispersion curves are locally linear near critical frequencies and that low frequency behaviour is replicated at these critical, high, frequencies. The homogenization approach accurately replaces the periodic medium and the precise manner in which this is achieved then opens the way to future numerical implementation of this technique to scattering problems.     

Homogenization of the Navier-Stokes equations in open sets perforated with tiny holes I. Abstract framework,a volume distribution of holes

This paper treats the homogenization of the Stokes or Navier-Stokes equations with a Dirichlet boundary condition in a domain containing many tiny solid obstacles, periodically distributed in each direction of the axes. (For example, in the three-dimensional case, the obstacles have a size of ³ and are located at the nodes of a regular mesh of size .) A suitable extension of the pressure is used to prove the convergence of the homogenization process to a Brinkman-type law (in which a linear zero-order term for the velocity is added to a Stokes or Navier-Stokes equation).