Deformation and failure of a film/substrate system subjected to spherical indentation: Part 1. Experimental validation of stresses and strains derived using Hankel transform technique in an elastic film/substrate system

Our concern here is to rationalize experimental observations of failure modes brought about by indentation of hard thin ceramic films deposited on metallic substrates. By undertaking this exercise, we would like to evolve an analytical framework that can be used for designs of coatings. In Part I of the paper we develop an algorithm and test it for a model system. Using this analytical framework we address the issue of failure of columnar TiN films in Part II [J. Mater. Res. 21, 783 (2006)] of the paper. In this part, we used a previously derived Hankel transform procedure to derive stress and strain in a birefringent polymer film glued to a strong substrate and subjected to spherical indentation. We measure surface radial strains using strain gauges and bulk film stresses using photo elastic technique (stress freezing). For a boundary condition based on Hertzian traction with no film interface constraint and assuming the substrate constraint to be a function of the imposed strain, the theory describes the stress distributions well. The variation in peak stresses also demonstrates the usefulness of depositing even a soft film to protect an underlying substrate.

Artificial Bee Colony (ABC) for multi-objective design optimization of composite structures

In this paper, we present a generic method/model for multi-objective design optimization of laminated composite components, based on Vector Evaluated Artificial Bee Colony (VEABC) algorithm. VEABC is a parallel vector evaluated type, swarm intelligence multi-objective variant of the Artificial Bee Colony algorithm (ABC). In the current work a modified version of VEABC algorithm for discrete variables has been developed and implemented successfully for the multi-objective design optimization of composites. The problem is formulated with multiple objectives of minimizing weight and the total cost of the composite component to achieve a specified strength. The primary optimization variables are the number of layers, its stacking sequence (the orientation of the layers) and thickness of each layer. The classical lamination theory is utilized to determine the stresses in the component and the design is evaluated based on three failure criteria: failure mechanism based failure criteria, maximum stress failure criteria and the tsai-wu failure criteria. The optimization method is validated for a number of different loading configurations-uniaxial, biaxial and bending loads. The design optimization has been carried for both variable stacking sequences, as well fixed standard stacking schemes and a comparative study of the different design configurations evolved has been presented. Finally the performance is evaluated in comparison with other nature inspired techniques which includes Particle Swarm Optimization (PSO), Artificial Immune System (AIS) and Genetic Algorithm (GA). The performance of ABC is at par with that of PSO, AIS and GA for all the loading configurations. (C) 2009 Elsevier B.V. All rights reserved.

Modelling of symmetric laminates under extension

A higher-order theory of laminated composites under in-plane loads is developed. The displacement field is expanded in terms of the thickness co-ordinate, satisfying the zero shear stress condition at the surfaces of the laminate. Using the principle of virtual displacement, the governing equations and boundary conditions are established. Numerical results for interlaminar stresses arising in the case of symmetric laminates under uniform extension have been obtained and are compared with similar results available in the literature.

Analysis of edge delaminations in laminates through combined use of quasi-three-dimensional, eight-noded, two-noded and transition elements

The use of appropriate finite elements in different regions of a stressed solid can be expected to be economical in computing its stress response. This concept is exploited here in studying stresses near free edges in laminated coupons. The well known free edge problem of [0/90], symmetric laminate is considered to illustrate the application of the concept. The laminate is modelled as a combination of three distinct regions. Quasi-three-dimensional eight-noded quadrilateral isoparametric elements (Q3D8) are used at and near the free edge of the laminate and two-noded line elements (Q3D2) are used in the region away from the free edge. A transition element (Q3DT) provides a smooth inter-phase zone between the two regions. Significant reduction in the problem size and hence in the computational time and cost have been achieved at almost no loss of accuracy.

Curved composite beams—optimum lay-up for buckling by ranking

Instability of laminated curved composite beams made of repeated sublaminate construction is studied using finite element method. In repeated sublaminate construction, a full laminate is obtained by repeating a basic sublaminate which has a smaller number of plies. This paper deals with the determination of optimum lay-up for buckling by ranking of such composite curved beams (which may be solid or sandwich). For this purpose, use is made of a two-noded, 16 degress of freedom curved composite beam finite element. The displacements u, v, w of the element reference axis are expressed in terms of one-dimensional first-order Hermite interpolation polynomials, and line member assumptions are invoked in formulation of the elastic stiffness matrix and geometric stiffness matrix. The nonlinear expressions for the strains, occurring in beams subjected to axial, flexural and torsional loads, are incorporated in a general instability analysis. The computer program developed has been used, after extensive checking for correctness, to obtain optimum orientation scheme of the plies in the sublaminate so as to achieve maximum buckling load for typical curved solid/sandwich composite beams.

Optimum Design of Composite Laminates for Strength by Ranking

In this paper we describe a method for the optimum design of fiber rein forced composite laminates for strength by ranking. The software developed based on this method is capable of designing laminates for strength; which are subjected to inplane and/or bending loads and optionally hygrothermal loads. Symmetric laminates only are considered which are assumed to be made of repeated sublaminate construction. Various layup schemes are evaluated based on the laminated plate theory and quadratic failure cri terion for the given mechanical and hygrothermal loads. The optimum layup sequence in the sublaminate and the number of such sublaminates required are obtained. Further, a ply-drop round-off scheme is adopted to arrive at an optimum laminate thickness. As an example, a family of 0/90/45/ -45 bi-directional lamination schemes are examined for dif ferent types of loads and the gains in optimising the ply orientations in a sublaminate are demonstrated.

An improved iterative finite element solution for pin joints

Accurate, reliable and economical methods of determining stress distributions are important for fastener joints. In the past the contact stress problems in these mechanically fastened joints using interference or push or clearance fit pins were solved using both inverse and iterative techniques. Inverse techniques were found to be most efficient, but at times inadequate in the presence of asymmetries. Iterative techniques based on the finite element method of analysis have wider applications, but they have the major drawbacks of being expensive and time-consuming. In this paper an improved finite element technique for iteration is presented to overcome these drawbacks. The improved iterative technique employs a frontal solver for elimination of variables not requiring iteration, by creation of a dummy element. This automatically results in a large reduction in computer time and in the size of the problem to be handled during iteration. Numerical results are compared with those available in the literature. The method is used to study an eccentrically located pin in a quasi-isotropic laminated plate under uniform tension.

Atmospheric boundary-layer across Hadley and Ferrel cells over the Indian Ocean

A Southern Ocean Pilot cruise covering the latitudes from 10 degrees N to 56 degrees S in the open Indian Ocean was carried out during January February 2004. Surface and upper air data collected during this cruise are reported here. It is shown that the broad features of the atmosphere, in particular that of temperature, follow the tropical and mid-latitude weather expected during January February in this region. However, the atmospheric boundary-layer shows large variations, both in its height and structure between tropics and high latitudes. Strong influence of the surface heat flux on boundary layer structure is clearly seen. Humidity field reveals several local maxima and minima, suggesting a laminated atmosphere with air from different sources moving almost unmixed in adjacent layers.

Finite element estimation of elastic interlaminar stresses in laminates

Low interlaminar strength and the consequent possibility of interlaminar failures in composite laminates demand an examination of interlaminar stresses and/or strains to ensure their satisfactory performance. As a first approximation, these stresses can be obtained from thickness-wise integration of ply equilibrium equations using in-plane stresses from the classical laminated plate theory. Implementation of this approach in the finite element form requires evaluation of third and fourth order derivatives of the displacement functions in an element. Hence, a high precision element developed by Jayachandrabose and Kirkhope (1985) is used here and the required derivatives are obtained in two ways. (i) from direct differentiation of element shape functions; and (ii) by adapting a finite difference technique applied to the nodal strains and curvatures obtained from the finite element analysis. Numerical results obtained for a three-layered symmetric and a two-layered asymmetric laminate show that the second scheme is quite effective compared to the first scheme particularly for the case of asymmetric laminates.

Recovery of tropical rainforest avifauna in relation to vegetation succession following shifting cultivation in Mizoram, north-east India

1. Recovery of rainforest bird community structure and composition, in relation to forest succession after slash-and-burn shifting cultivation or jhum was studied in Mizoram, north-east India. Replicate fallow sites abandoned after shifting cultivation 1, 5, 10, 25 and approximate to 100 years ago, were compared with primary evergreen and semi-evergreen forest using transect and quadrat sampling. 2. Vegetation variables such as woody plant species richness, tree density and vertical stratification increased with fallow age in a rapid. nun-linear, asymptotic manner. Principal components analysis of vegetation variables summarized 92.8% of the variation into two axes: PC1 reflecting forest development and woody plant succession (variables such as tree density, woody plant species richness), and PC2 depicting bamboo density, which increased from 1 to 25 years and declined thereafter. 3. Bird species richness, abundance and diversity, increased rapidly and asymptotically during succession paralleling vegetation recovery as shown by positive correlations with fallow age and PC1 scores of sites. Bamboo density reflected by PC2 had a negative effect on bird species richness and abundance. 4. The bird community similarity (Morisita index) of sites with primary forest also increased asymptotically with fallow age indicating sequential species turnover during succession. Bird community similarity of sites with primary forest (or between sites) was positively correlated with both physiognomic and floristic similarities with primary forest (or between sites). 5. The number of bird species in guilds associated with forest development and woody plants (canopy insectivores, frugivores: bark feeders) was correlated with PCI scores of the sites. Species in other guilds (e. g. granivores, understorey insectivores) appeared to dominate during early and mid-succession. 6. The non-linear relationships imply that fallow periods less than a threshold of 25 years for birds, and about 50-75 years for woody plants, are likely to cause substantial community alteration. 7. As 5-10-year rotation periods or jhum cycles prevail in many parts of north-east India. there is a need to protect and conserve tracts of late-successional and primary forest.

The role of the interface in foam-bearing glass-epoxy composites subjected to impact loading

foam, either stacked together as three layers (MC) or inserted at three different positions (3L) while arranging the stacking sequence during the fabrication of glass fiber-epoxy composites, form the subject of investigation. This stacking variation resulted in a different interfacial area between these foam materials and the glass-epoxy regions in the laminates. This area in designed to be maximum for the 3L variety. The energy of impact being high enough to cause development of the crack in the samples, how the change in interfacial area affects the traverse of the crack front and the failure feature of the laminated composite are reported in the form of photomacrographs in this work. The results point to significant changes for the impact data, like for instance the peak load attained by the different samples, through thickness crack propagation and tensile fracture features on the non-impacted end for the plain variety, separation about the mid-zone for the MC laminates and two or more layer separations for the 3L variety. The separation for the foam-bearing systems occur invariably at the interface and here again one of the (two identical) interfaces only is chosen for the separation.

Micro and macroscopic features of impact failed defect induced carbon-epoxy composites

Laminated composite structures are susceptible to damage under impacts with attendant properly degradation. While studies on damage tolerance behaviour are emphasised and the findings reported, the citations correlating impacts with the fracture features are limited. In the present study, therefore, attempts have been made to depict how the transition of the fracture features take place depending on the type and extent of defect introduced onto the carbon-epoxy system. The test specimens were subjected to differing levels of low energy pendulum impacts with a view to have specimens with varying levels of intial impacts history. Into such specimens, additional defect in the form of slits of varying depths were introduced by a mechanical process. The test coupons were then allowed to fail by impact. The fracture surface was studied under scanning electron microscope. The fractographic features that appear, based on the induced/inserted defects, are presented in this paper. It was noticed that the energy absorbed for final fracture could be associated with the defect introduced into the system. It was also observed that the size of the mechanically inserted defect had a significant influence on the features of the fracture surface.

FE analysis of piezo-laminate composites under thermal loads

A shear flexible 4-noded finite element formulation, having five mechanical degrees of freedom per node, is presented for modeling the dynamic as well as the static thermal response of laminated composites containing distributed piezoelectric layers. This element has been developed to have one electrical degree of freedom per piezoelectric layer. The mass, stiffness and thermo-electro-mechanical coupling effects on the actuator and sensor layers have been considered. Numerical studies have been conducted to investigate both the sensory and active responses on piezoelectric composite beam and plate structures. It is. concluded that both the thermal and pyroelectric effects are important and need to be considered in the precision distributed control of intelligent structures.

On the utility of an instrumented impact set-up to characterise the response of discontinuous buffer strips in carbon–epoxy systems in dry and wet conditions

The use of an instrumented impact test set-up to evaluate the influence of water ingress on the impact response of a carbon–epoxy (C–E) laminated composite system containing discontinuous buffer strips (BS) has been examined. The data on the BS-free C–E sample in dry conditions are used as reference to compare with the data derived from those immersed in water. The work demonstrated the utility of an instrumented impact test set-up in characterising the response, first owing to the architectural difference due to introduction of buffer strips and then due to the presence of an additional phase in the form of water ingressed into the sample. The presence of water was found to enhance the energy absorption characteristics of the C–E system with BS insertions. It was also noticed that with an increasing number of BS layer insertions, the load–time plots displayed characteristic changes. The ductility indices (DI) were found to display a lower value for the water immersed samples compared to the dry ones.

Fully Comprehensive Geometrically Non-Linear Dynamic Analysis of Multi-Body Beam Systems with Elastic Couplings

This paper is concerned with the dynamic analysis of flexible,non-linear multi-body beam systems. The focus is on problems where the strains within each elastic body (beam) remain small. Based on geometrically non-linear elasticity theory, the non-linear 3-D beam problem splits into either a linear or non-linear 2-D analysis of the beam cross-section and a non-linear 1-D analysis along the beam reference line. The splitting of the three-dimensional beam problem into two- and one-dimensional parts, called dimensional reduction,results in a tremendous savings of computational effort relative to the cost of three-dimensional finite element analysis,the only alternative for realistic beams. The analysis of beam-like structures made of laminated composite materials requires a much more complicated methodology. Hence, the analysis procedure based on Variational Asymptotic Method (VAM), a tool to carry out the dimensional reduction, is used here.The analysis methodology can be viewed as a 3-step procedure. First, the sectional properties of beams made of composite materials are determined either based on an asymptotic procedure that involves a 2-D finite element nonlinear analysis of the beam cross-section to capture trapeze effect or using strip-like beam analysis, starting from Classical Laminated Shell Theory (CLST). Second, the dynamic response of non-linear, flexible multi-body beam systems is simulated within the framework of energy-preserving and energy-decaying time integration schemes that provide unconditional stability for non-linear beam systems. Finally,local 3-D responses in the beams are recovered, based on the 1-D responses predicted in the second step. Numerical examples are presented and results from this analysis are compared with those available in the literature.