Effect of matrix cracking and material uncertainty on composite plates

A laminated composite plate model based on first order shear deformation theory is implemented using the finite element method.Matrix cracks are introduced into the finite element model by considering changes in the A, B and D matrices of composites. The effects of different boundary conditions, laminate types and ply angles on the behavior of composite plates with matrix cracks are studied.Finally, the effect of material property uncertainty, which is important for composite material on the composite plate, is investigated using Monte Carlo simulations. Probabilistic estimates of damage detection reliability in composite plates are made for static and dynamic measurements. It is found that the effect of uncertainty must be considered for accurate damage detection in composite structures. The estimates of variance obtained for observable system properties due to uncertainty can be used for developing more robust damage detection algorithms. (C) 2010 Elsevier Ltd. All rights reserved.

Wave propagation analysis in laminated composite plates with transverse cracks using the wavelet spectral finite element method

This paper presents a newly developed wavelet spectral finite element (WFSE) model to analyze wave propagation in anisotropic composite laminate with a transverse surface crack penetrating part-through the thickness. The WSFE formulation of the composite laminate, which is based on the first-order shear deformation theory, produces accurate and computationally efficient results for high frequency wave motion. Transverse crack is modeled in wavenumber-frequency domain by introducing bending flexibility of the plate along crack edge. Results for tone burst and impulse excitations show excellent agreement with conventional finite element analysis in Abaqus (R). Problems with multiple cracks are modeled by assembling a number of spectral elements with cracks in frequency-wavenumber domain. Results show partial reflection of the excited wave due to crack at time instances consistent with crack locations. (C) 2014 Elsevier B.V. All rights reserved.

The effect of shear strength on the ballistic response of laminated composite plates

The ballistic performance of clamped circular carbon fibre reinforced polymer (CFRP) and Ultra High Molecular Weight Polyethylene (UHMWPE) fibre composite plates of equal areal mass and 0/90 lay-up were measured and compared with that of monolithic 304 stainless steel plates. The effect of matrix shear strength upon the dynamic response was explored by testing: (i) CFRP plates with both a cured and uncured matrix and (ii) UHMWPE laminates with identical fibres but with two matrices of different shear strength. The response of these plates when subjected to mid-span, normal impact by a steel ball was measured via a dynamic high speed shadow moiré technique. Travelling hinges emanate from the impact location and travel towards the supports. The anisotropic nature of the composite plate results in the hinges travelling fastest along the fibre directions and this results in square-shaped moiré fringes in the 0/90 plates. Projectile penetration of the UHMWPE and the uncured CFRP plates occurs in a progressive manner, such that the number of failed plies increases with increasing velocity. The cured CFRP plate, of high matrix shear strength, fails by cone-crack formation at low velocities, and at higher velocities by a combination of cone-crack formation and communition of plies beneath the projectile. On an equal areal mass basis, the low shear strength UHMWPE plate has the highest ballistic limit followed by the high matrix shear strength UHMWPE plate, the uncured CFRP, the steel plate and finally the cured CFRP plate. We demonstrate that the high shear strength UHMWPE plate exhibits Cunniff-type ballistic limit scaling. However, the observed Cunniff velocity is significantly lower than that estimated from the laminate properties. The data presented here reveals that the Cunniff velocity is limited in its ability to characterise the ballistic performance of fibre composite plates as this velocity is independent of the shear properties of the composites: the ballistic limit of fibre composite plates increases with decreasing matrix shear strength for both CFRP and UHMWPE plates. © 2013 Elsevier Masson SAS. All rights reserved.

Fasteners in Composite Plates - Effects of Interfacial Friction

The effects of tangential friction at pin—hole interfaces are appropriately modelled for the analysis of fasteners in large composite (orthotropic) plate loaded along its edges. The pin—hole contact could be of interference, clearance or neat fit. When the plate load is monotonically increased, interference fits give rise to receding contact, whereas clearance fits result in advancing contact. In either case, the changing contact situations lead to non-linear moving boundary value problems. The neat fit comes out as a special case in which the contact and separation regions are invariant with the applied load level and so the problem remains linear. The description of boundary conditions in the presence of tangential friction, will depend on whether the problem is one of advancing or receding contact, advancing contact presenting a special problem. A model is developed for the limiting case of a rigid pin and an ideally rough interface (infinitely large friction coefficient). The non-linearity resulting from the continuously varying proportions of contact and separation at the interface, is handled by an “Inverse Formulation” which was successfully applied earlier by the authors for smooth (zero friction) interfacial conditions. The additional difficulty introduced by advancing contact is handled by adopting a “Marching Solution”. The modelling and the procedure are illustrated in respect of symmetric plate load cases. Numerical results are presented bringing out the effects of interfacial friction and plate orthotropy on load-contact relations and plate stresses.

Flexure of composite plates

A new higher order shear deformation theory of laminated composite plates is developed. The basic displacement variables in this theory are two partial normal displacements and two in-plane displacement parameters. The governing equations are presented in the form of four simultaneous partial differential equations. The shear deformation theories of Bhimareddy and Stevens, and of Reddy are special cases of this formulation. In their models, transverse shear strains will become zero at points in the plate where displacements are constrained to be zero such as those on fixed edges. This limitation has been overcome in the present formulation.

Deflection analysis of hybrid laminated composite plates

In this paper an attempt is made to obtain deflections of hybrid, laminated, rectangular and skew composite plates. Analysis is performed by employing the Galerkin technique. Numerical results have been obtained for two types of layups employing Kevlar/epoxy and Boron/epoxy laminae. It is observed that for a given aspect ratio the rigidity of the skew plate increases with an increase in the skew angle. Further, for a specified deflection, the hybrid laminates turn out to be lighter.

Recent studies on the behaviour of interference fit pins in composite plates

Composites are finding increasing application in many advanced engineering fields like aerospace, marine engineering, hightech sports equipment, etc., due to their high specific strength and/or specific stiffness values. The use of composite components in complex situations like airplane wing root or locations of concentrated load transfer is limited due to the lack of complete understanding of their behaviour in the region of joints. Joints are unavoidable in the design and manufacture of complex structures. Pin joints are one of the most commonly used methods of connection. In regions of high stresses like airplane wing root joints interference fit pins are used to increase its fatigue life and thereby increase the reliability of the whole structure. The present contribution is a study on the behaviour of the interference fit pin in a composite plate subjected to both pull and push type of loads. The interference fit pin exhibits partial contact/separation under the loads and the contact region is a non-linear function of the load magnitude. This non-linear behaviour is studied by adopting the inverse technique and some new results are presented in this paper.

A geodesic-based triangulation technique for damage location in metallic and composite plates

Lamb-wave-based damage detection methods using the triangulation technique are not suitable for handling structures with complex shapes and discontinuities as the parametric/analytical representation of these structures is very difficult. The geodesic concept is used along with the triangulation technique to overcome the above problem. The present work is based on the fundamental fact that a wave takes the minimum energy path to travel between two points on any multiply connected surface and this reduces to the shortest distance path or geodesic. The geodesics are computed on the meshed surface of the structure using the fast marching method. The wave response matrix of the given sensor configuration for the healthy and the damaged structure is obtained experimentally. The healthy and damage response matrices are compared and their difference gives the time information about the reflection of waves from the damage. A wavelet transform is used to extract the arrival time information of the wave scattered by the damage from the acquired Lamb wave signals. The computed geodesics and time information are used in the ellipse algorithm of triangulation formulation to locate the loci of possible damage location points for each actuator-sensor pair. The results obtained for all actuator-sensor pairs are combined and the intersection of multiple loci gives the damage location result. Experiments were conducted in aluminum and composite plate specimens to validate this method.

Stability of laminated composite plates with cut-outs

Critical buckling loads of laminated fibre-reinforced plastic square panels have been obtained using the finite element method. Various boundary conditions, lay-up details, fibre orientations, cut-out sizes are considered. A 36 degrees of freedom triangular element, based on the classical lamination theory (CLT) has been used for the analysis. The performance of this element is validated by comparing results with some of those available in literature. New results have been given for several cases of boundary conditions for [0°/ ± 45°/90°]s laminates. The effect of fibre-orientation in the ply on the buckling loads has been investigated by considering [±?]6s laminates.

Large deformation finite element analysis of laminated circular composite plates

A 48 d.o.f., four-noded quadrilateral laminated composite shell finite element is particularised to a sector finite element and is used for the large deformation analysis of circular composite laminated plates. The strain-displacement relationships for the sector element are obtained by reducing those of the quadrilateral shell finite element by substituting proper values for the geometric parameters. Subsequently, the linear and tangent stiffness matrices are formulated using conventional methods. The Newton-Raphson method is employed as the nonlinear solution technique. The computer code developed is validated by solving an isotropic case for which results are available in the literature. The method is then applied to solve problems of cylindrically orthotropic circular plates. Some of the results of cylindrically orthotropic case are compared with those available in the literature. Subsequently, application is made to the case of laminated composite circular plates having different lay-up schemes. The computer code can handle symmetric/unsymmetric lay-up schemes. The large displacement analysis is useful in estimating the damage in composite plates caused by low-velocity impact.

Large deformation dynamic finite element analysis of delaminated composite plates using contact-impact conditions

A new C-0 composite plate finite element based on Reddy's third order theory is used for large deformation dynamic analysis of delaminated composite plates. The inter-laminar contact is modeled with an augmented Lagrangian approach. Numerical results show that the widely used ``unconditionally stable'' beta-Newmark method presents instability problems in the transient simulation of delaminated composite plate structures with large deformation. To overcome this instability issue, an energy and momentum conserving composite implicit time integration scheme presented by Bathe and Baig is used. It is found that a proper selection of the penalty parameter is very crucial in the contact simulation. (C) 2014 Elsevier Ltd. All rights reserved.

Bending solution of high-order refined shear deformation-theory for rectangular composite plates

A new high-order refined shear deformation theory based on Reissner's mixed variational principle in conjunction with the state- space concept is used to determine the deflections and stresses for rectangular cross-ply composite plates. A zig-zag shaped function and Legendre polynomials are introduced to approximate the in-plane displacement distributions across the plate thickness. Numerical results are presented with different edge conditions, aspect ratios, lamination schemes and loadings. A comparison with the exact solutions obtained by Pagano and the results by Khdeir indicates that the present theory accurately estimates the in-plane responses.

Theoretical analysis on microwave heating of oil–water emulsions supported on ceramic, metallic or composite plates

A detailed theoretical analysis has been carried out to study efficient heating due to microwaves for one-dimensional (1D) oil–water emulsion samples placed on various ceramic, metallic (reflective) and ceramic–metallic composite supports. Two typical emulsion systems are considered such as oil-in-water (o/w) and water-in-oil (w/o). A preliminary study has been carried out via average power vs emulsion thickness diagram to estimate microwave power absorption within emulsion samples for various cases. The maxima in average power, also termed as ‘resonances’, are observed for specific emulsion thicknesses and the two consecutive resonances of significant magnitudes are termed as R1 and R2 modes. For both o/w and w/o emulsions, it is observed that microwave power absorption is enhanced in presence of metallic and composite supports during both R1 and R2 modes. The efficient heating strategies characterized by ‘large heating rates’ with ‘minimal thermal runaway’ i.e. uniform temperature distributions within the sample have been assessed for each type of emulsion. Based on the detailed spatial distributions of power and temperature for various cases, SiC-metallic composite support may be recommended as an optimal heating strategy for o/w samples with higher oil fractions (0.45) whereas metallic and Alumina-metallic composite supports may be favored for samples with smaller oil fractions (=0.3) during R1 mode. For w/o samples, SiC-metallic composite support may be suitable heating strategy for all ranges of water fractions during R1 mode. During R2 mode, metallic and Alumina-metallic composite supports are favored for both o/w and w/o emulsion samples. Current study recommends the efficient way to use microwaves in a single mode waveguide and the heating strategy can be suitably extended for heating of any other emulsions for which dielectric properties are easily measurable or available in the literature.

Multiple cutout optimization in composite plates using evolutionary structural optimization

The optimization of cutouts in composite plates was investigated by implementing a procedure known as Evolutionary Structural Optimization. Perforations were introduced into a finite element mesh of the plate from which one or more cutouts of a predetermined size were evolved. In the examples presented, plates were rejected from around each evolving cutout based on a predefined rejection criterion. The Limiting ply within each plate element around the cutout was determined based on the Tsai-Hill failure criterion. Finite element plates with values below the product of the average Tsai-Hill number and a rejection criterion were subsequently removed. This process was iterated until a steady state was reached and the rejection criterion was then incremented by an evolutionary rate and the above steps repeated until the desired cutout area was achieved. Various plates with differing lay-up and loading parameters were investigated to demonstrate the generality and robustness of this optimization procedure.