T.S. Ramamurthy

Generalized modified crack closure integral (GMCCI) and its application to interface crack problems

It is well known that at an interface crack tip, mode mixity exists even under single mode loading situations. The modified crack closure integral (MCCI) technique together with finite element analysis has been extensively used in the past to estimate the mixed mode strain energy release rate (SERR) components at these crack tips. One of the aspects which attracted considerable attention in these problems is the convergence of SERR components as the virtual crack extension size is progressively decreased to zero. Such studies have proved to be computationally expensive as the problem has to be solved each time for various sizes of virtual crack extension. Also it was not possible to estimate the effect of preventing the crack face interpenetration on SERR estimates for large sizes of virtual crack extension greater than the contact zone size. To overcome these in this paper, a generalized modified crack closure integral (GMCCI) algorithm is proposed for four and eight-noded isoparametric quadrilateral elements, which can estimate the SERR components for several sizes of virtual crack extension through a single finite element analysis. This is achieved through a numerical integration of the crack closure integral (CCI) and for various sizes of virtual crack extension covering a large number of equal or unequal size elements ahead and behind the crack tip. The procedure developed is demonstrated with application to typical interface crack problems under remote tension and shear.

Elasto-plastic contact stress analysis of joints subjected to cyclic loading

Rigorous elastic-plastic finite element analysis of joints subjected to cyclic loading is carried out. An incremental-iterative algorithm is developed in a modular form combining elasto-plastic material behaviour and contact stress analysis. For the case of the interference fit, the analysis sequentially carries out insertion of the pin and application of the load on the joint, covering possible initiation of separation (and/or yielding) and progressively the receding/advancing contact at the pin-plate interface. Deformations of both the plate and the pin are considered in the analysis. Numerical examples are presented for the case of an interference fit pin in a large plate under remote cyclic tension, and for an interference fit pin lug joint subjected to cyclic loading. A detailed study is carried out for the latter problem considering the effect of change in contact/separation at the pin-plate interface on local stresses, strains and redistribution of these stresses with the spread of a plastic zone. The results of the study are a useful input for the estimation of the fatigue life of joints. Copyright (C) 1996 Elsevier Science Ltd

Delamination tolerance studies in laminated composite panels

Determination of levels of tolerance in delaminated composite panels is an important issue in composite structures technology. The primary intention is to analyse delaminated composite panels and estimate Strain Energy Release Rate (SERR) parameters at the delamination front to feed into acceptability criteria. Large deformation analysis is necessary to cater for excessive rotational deformations in the delaminated sublaminate. Modified Virtual Crack Closure Integral (MVCCI) is used to estimate all the three SERR components at the delamination front from the finite element output containing displacements, strains and stresses. The applied loading conditions are particularly critical and compressive loading on the panel could lead to buckling of the delaminated sublaminate and consequent growth of delamination. Numerical results are presented for circular delamination of various sizes and delamination at various interfaces (varying depth-wise location) between the base- and the sub-laminates. Numerical data are also presented on the effect of bi-axial loading and in particular on compressive loading in both directions. The results can be used to estimate delamination tolerance at various depths (or at various interfaces) in the laminate.

Analysis of a cracked lug loaded by an interference fit pin

In the present study, a lug joint fitted with an interference fit (oversized) pin is considered with radial through cracks situated at diametrically opposite points perpendicular to the loading direction. A finite element contact stress algorithm is developed with linear elastic assumptions to deal with varying partial contact/separation at the pin-plate interface using a marching solution. Stress Intensity Factor (SIF) at the crack tips is evaluated using the Modified Crack Closure Integral (MCCI) method. The effect of change in crack length and edge distance on the load-contact relation, SIFs and stress distributions are studied. A rigorous plane stress elasticity solution of the pin-plate interface at the crack mouth confirmed the existence of the stress concentration leading to a local peak in the radial stress at the crack mouth and provided a method of estimating it quantitatively. Copyright (C) 1996 Elsevier Science Ltd.

Elastic analysis of pin joints in plates under some combined pin and plate loads

Joints are primary sources of weakness in structures. Pin joints are very common and are used where periodic disassembly of components is needed. A circular pin in a circular hole in an infinitely large plate is an abstraction of such a pin joint. A two-dimensional plane-stress analysis of such a configuration is carried out, here, subjected to pin-bearing and/or biaxial-plate loading. The pin is assumed to be rigid compared to the plate material. For pin load the reactive stresses at the edges of the infinite plate tend to zero though their integral over the external boundary equals to the pin load. The pin-hole interface is unbonded and so beyond some load levels the plate separates from the pin and the extent of separation is a non-linear function of load level. The problem is solved by inverse technique where the extent of contact is specified and the causative loads are evaluated directly. In the situations where combined load is acting the separation-contact zone specification generally needs two parameters (angles) to be specified. The present report deals with analysing such a situation in metallic (or isotropic) plates. Numerical results are provided for parametric representation and the methodology is demonstrated.

Analysis of a finite composite plate with smooth rigid pin

A continuum method of analysis is presented in this paper for the problem of a smooth rigid pin in a finite composite plate subjected to uniaxial loading. The pin could be of interference, push or clearance fit. The plate is idealized to an orthotropic sheet. As the load on the plate is progressively increased, the contact along the pin-hole interface is partial above certain load levels in all three types of fit. In misfit pins (interference or clearance), such situations result in mixed boundary value problems with moving boundaries and in all of them the arc of contact and the stress and displacement fields vary nonlinearly with the applied load. In infinite domains similar problems were analysed earlier by ‘inverse formulation’ and, now, the same approach is selected for finite plates. Finite outer domains introduce analytical complexities in the satisfaction of boundary conditions. These problems are circumvented by adopting a method in which the successive integrals of boundary error functions are equated to zero. Numerical results are presented which bring out the effects of the rectangular geometry and the orthotropic property of the plate. The present solutions are the first step towards the development of special finite elements for fastener joints.

Fracture Analysis of Delaminated Composite Panels in Compression Using Numerically Integrated MVCCI

Delamination tolerance is one of the key issues in laminated composites used in high technology structures. Delamination in these structures could grow due to repeated in-plane compressive loading during operation and the delamination tolerance levels are obtained by comparing mode-I SERR component with the corresponding fatigue threshold values. The Modified Virtual Crack Closure Integral (MVCCI) is used to estimate individual mode SERR components in a mixed-mode situation, which occurs at an embedded delamination. The present work proposed a numerical integration for the first time to evaluate the MVCCI for Hexa20 elements used in three-dimensional finite element analysis. Extensive numerical results are presented including a parametric study on a typical laminate.

Development of special fastener elements

A special finite element (FASNEL) is developed for the analysis of a neat or misfit fastener in a two-dimensional metallic/composite (orthotropic) plate subjected to biaxial loading. The misfit fasteners could be of interference or clearance type. These fasteners, which are common in engineering structures, cause stress concentrations and are potential sources of failure. Such cases of stress concentration present considerable numerical problems for analysis with conventional finite elements. In FASNEL the shape functions for displacements are derived from series stress function solutions satisfying the governing difffferential equation of the plate and some of the boundary conditions on the hole boundary. The region of the plate outside FASNEL is filled with CST or quadrilateral elements. When a plate with a fastener is gradually loaded the fastener-plate interface exhibits a state of partial contact/separation above a certain load level. In misfit fastener, the extent of contact/separation changes with applied load, leading to a nonlinear moving boundary problem and this is handled by FASNEL using an inverse formulation. The analysis is developed at present for a filled hole in a finite elastic plate providing two axes of symmetry. Numerical studies are conducted on a smooth rigid fastener in a finite elastic plate subjected to uniaxial loading to demonstrate the capability of FASNEL.

Analytical solutions for load transfer through smooth elastic pin in an infinite orthotropic plate

Abstract Analytical solutions are developed for stresses in an infinite orthotropic plate with smooth elastic pin under pin and plate loading. The pin could be of interference, push or clearance fit. A two-dimensional plane stress elastic analysis using a complex variable formulation in conjunction with equidistant boundary collocation scheme for satisfying boundary conditions is carried out. An inverse technique is used to obtain loads, stresses and displacements for a specified contact/separation configuration. Numerical results are presented which bring out the effect of pin modulus of elasticity and the degree of orthotropy of plate on bearing stress, hoop stress and load contact behaviour. The results indicate that decrease in the pin modulus of elasticity relative to that of the plate ( E p E 1 ) decreases the bearing stresses at the pin-plate interface.

Thermomechanical generalization of pin joints in finite plates

The understanding of thermoelastic behaviour of joints is significant in order to ensure the integrity of large and complex structures exposed to a thermal environment, particularly in fields such as aerospace and nuclear engineering. Thermomechanical generalization of partial contact behaviour of a pin joint under combined in-plane mechanical loading and on-axis unidirectional heat flow has already been established by the authors for the analytically simpler domains of large plates. This paper successfully extends the on-going investigation to a single pin in a finite rectangular isotropic plate as a two-dimensional abstraction from a practical situation of a multipin fastener joint. The finite element method is used to analyse the joint problem under on-axis thermomechanical loading and unified load-contact relationships are established for a class of loading conditions.