On an analytical-numerical procedure for the analysis of cylindrical shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

A Spectral Iteration Technique For The Analysis Of Yagi-Uda Arrays

The paper present a spectral iteration technique for the analysis of linear arrays of unequally spaced dipoles of unequal lengths. As an example, the Yagi-Uda array is considered for illustration. Analysis is carried out in both the spatial as well as the spectral domains, the two being linked by the Fourier transform. The fast Fourier transform algorithm is employed to obtain an iterative solution to the electric field integral equation and the need for matrix inversion is circumvented. This technique also provides a convenient means for testing the satisfaction of the boundary conditions on the array elements. Numerical comparison of the input impedance and radiation pattern have been made with results deduced elsewhere by other methods. The computational efficency of this technique has been found to be significant for large arrays.

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.

Transmission loss analysis of rectangular expansion chamber with arbitrary location of inlet/outlet by means of Green's functions

Transmission loss of a rectangular expansion chamber, the inlet and outlet of which are situated at arbitrary locations of the chamber, i.e., the side wall or the face of the chamber, are analyzed here based on the Green's function of a rectangular cavity with homogeneous boundary conditions. The rectangular chamber Green's function is expressed in terms of a finite number of rigid rectangular cavity mode shapes. The inlet and outlet ports are modeled as uniform velocity pistons. If the size of the piston is small compared to wavelength, then the plane wave excitation is a valid assumption. The velocity potential inside the chamber is expressed by superimposing the velocity potentials of two different configurations. The first configuration is a piston source at the inlet port and a rigid termination at the outlet, and the second one is a piston at the outlet with a rigid termination at the inlet. Pressure inside the chamber is derived from velocity potentials using linear momentum equation. The average pressure acting on the pistons at the inlet and outlet locations is estimated by integrating the acoustic pressure over the piston area in the two constituent configurations. The transfer matrix is derived from the average pressure values and thence the transmission loss is calculated. The results are verified against those in the literature where use has been made of modal expansions and also numerical models (FEM fluid). The transfer matrix formulation for yielding wall rectangular chambers has been derived incorporating the structural–acoustic coupling. Parametric studies are conducted for different inlet and outlet configurations, and the various phenomena occurring in the TL curves that cannot be explained by the classical plane wave theory, are discussed.

Surface acoustic waves of finite amplitude excited by a monochromatic line source

A technique for obtaining a uniformly valid solution to the problem of nonlinear propagation of surface acoustic waves excited by a monochromatic line source is presented. The method of solution is an extension of the method of strained coordinates wherein both the dependent and independent variables are expanded in perturbation series. A special transformation is proposed for the independent variables so as to make the expansions uniformly valid and also to satisfy all the boundary conditions. This perturbation procedure, carried out to the second order, yields a solution containing a second harmonic surface wave whose amplitude and phase exhibit an oscillatory variation along the direction of propagation. In addition, the solution also contains a second harmonic bulk wave of constant amplitude but varying phase propagating into the medium.

Indentation strength of a piezoelectric ceramic: Experiments and simulations

The spherical indentation strength of a lead zirconate titanate (PZT) piezoelectric ceramic was investigated under poled and unpoled conditions and with different electrical boundary conditions (arising through the use of insulating or conducting indenters). Experimental results show that the indentation strength of the poled PZT is higher than that of the unpoled PZT. The strength of a poled PZT under a conducting indenter is higher than that under an insulating indenter. Poling direction (with respect to the direction of indentation loading) did not significantly affect the strength of material. Complementary finite element analysis (FEA) of spherical indentation of an elastic, linearly coupled piezoelectric half-space is conducted for rationalizing the experimental observations. Simulations show marked dependency of the contact stress on the boundary conditions. In particular, contact stress redistribution in the Coupled problem leads to a change in the fracture initiation, from Hertzian cracking in the unpoled material to Subsurface damage initiation in poled PZT. These observations help explain the experimental ranking of strength the PZT in different material conditions or under different boundary conditions.

Heat transfer with solid liquid phase change initiating at a point in a semi infinite mold

Short-time analytical solutions of temperature and moving boundary in two-dimensional two-phase freezing due to a cold spot are presented in this paper. The melt occupies a semi-infinite region. Although the method of solution is valid for various other types of boundary conditions, the results in this paper are given only for the prescribed flux boundary conditions which could be space and time dependent. The freezing front propagations along the interior of the melt region exhibit well known behaviours but the propagations along the surface are of new type. The freezing front always depends on material parameters. Several interesting results can be obtained as particular cases of the general results.

Load Transfer from a Smooth Elastic Pin to a Large Sheet

The plane problem of load transfer from an elastic interference or clearance fit pin to a large elastic sheet with a perfectly smooth interface is solved. As the load on the pin is monotonically increased, the pin-hole interface is in partial contact above certain critical load in interference fit and throughout the loading range in clearance fit.Such situations result in mixed boundary-value problems with moving boundaries and the arc of contact varies nonlinearly with applied load. These problems are analyzed by an inverse technique in which the arcs of contact/separation are prescribed and the causative loads are evaluated. A direct method of analysis is adopted using biharmonic polar trigonometric stress functions and a simple collocation method for satisfying the boundary conditions. A unified analytical formulation is achieved for interference and clearance fits. The solutions for the linear problem of push fits are inherent in the unified analysis. Numerical results highlighting the effects of pin and sheet elasticity parameters are presented.

Vibrations of a Beam in Variable Contact With a Flat Surface

We study small vibrations of cantilever beams contacting a rigid surface. We study two cases: the first is a beam that sags onto the ground due to gravity, and the second is a beam that sticks to the ground through reversible adhesion. In both cases, the noncontacting length varies dynamically. We first obtain the governing equations and boundary conditions, including a transversality condition involving an end moment, using Hamilton's principle. Rescaling the variable length to a constant value, we obtain partial differential equations with time varying coefficients, which, upon linearization, give the natural frequencies of vibration. The natural frequencies for the first case (gravity without adhesion) match that of a clamped-clamped beam of the same nominal length; frequencies for the second case, however, show no such match. We develop simple, if atypical, single degree of freedom approximations for the first modes of these two systems, which provide insights into the role of the static deflection profile, as well as the end moment condition, in determining the first natural frequencies of these systems. Finally, we consider small transverse sinusoidal forcing of the first case and find that the governing equation contains both parametric and external forcing terms. For forcing at resonance, w find that either the internal or the external forcing may dominate.

Combined Use of Solid of Revolution, Thin Shell, and Interphase Elements for Analysis of Cylindrical Shells

Near the boundaries of shells, thin shell theories cannot always provide a satisfactory description of the kinematic situation. This imposes severe limitations on simulating the boundary conditions in theoretical shell models. Here an attempt is made to overcome the above limitation. Three-dimensional theory of elasticity is used near boundaries, while thin shell theory covers the major part of the shell away from the boundaries. Both regions are connected by means of an “interphase element.” This method is used to study typical static stress and natural vibration problems

On an analytical-numerical procedure for the analysis of cylindrical-shells with arbitrarily oriented cracks

An analytical-numerical procedure for obtaining stress intensity factor solutions for an arbitrarily oriented crack in a long, thin circular cylindrical shell is presented. The method of analysis involves obtaining a series solution to the governing shell equation in terms of Mathieu and modified Mathieu functions by the method of separation of variables and satisfying the crack surface boundary conditions numerically using collocation. The solution is then transformed from elliptic coordinates to polar coordinates with crack tip as the origin through a Taylor series expansion and membrane and bending stress intensity factors are computed. Numerical results are presented and discussed for the pressure loading case.

Some studies on free vibration of composite laminates

Free vibration analysis is carried out to study the vibration characteristics of composite laminates using the modified shear deformation, layered, composite plate theory and employing the Rayleigh-Ritz energy approach. The analysis is presented in a unified form so as to incorporate all different combinations of laminate boundary conditions and with full coverage with regard to the various design parameters of a laminated plate. A parametric study is made using a beam characteristic function as the admissible function for the numerical calculations. The numerical results presented here are for an example case of fully clamped boundary conditions and are compared with previously published results. The effect of parameters, such as the aspect ratio of plates, ply-angle, number of layers and also the thickness ratios of plies in laminates on the frequencies of the laminate, is systematically studied. It is found that for anti-symmetric angle-ply or cross-ply laminates unique numerical values of the thickness ratios exist which improve the vibration characteristics of such laminates. Numerical values of the non-dimensional frequencies and nodal patterns, using the thickness ratio distribution of the plies, are then obtained for clamped laminates, fabricated out of various commonly used composite materials, and are presented in the form of the design curves.

Transverse vibrations of hybrid laminated plates

In this paper, an attempt is made to obtain the free vibration response of hybrid, laminated rectangular and skew plates. The Galerkin technique is employed to obtain an approximate solution of the governing differential equations. It is found that this technique is well suited for the study of such problems. Results are presented in a graphical form for plates with one pair of opposite edges simply supported and the other two edges clamped. The method is quite general and can be applied to any other boundary conditions.

Application of partition method to vibration problems of plates

The “partition method” or “sub-domain method” consists of expressing the solution of a governing differential equation, partial or ordinary, in terms of functions which satisfy the boundary conditions and setting to zero the error in the differential equation integrated over each of the sub-domains into which the given domain is partitioned. In this paper, the use of this method in eigenvalue problems with particular reference to vibration of plates is investigated. The deflection of the plate is expressed in terms of polynomials satisfying the boundary conditions completely. Setting the integrated error in each of the subdomains to zero results in a set of simultaneous, linear, homogeneous, algebraic equations in the undetermined coefficients of the deflection series. The algebraic eigenvalue problem is then solved for eigenvalues and eigenvectors. Convergence is examined in a few typical cases and is found to be satisfactory. The results obtained are compared with existing results based on other methods and are found to be in very good agreement.

A refined analysis of composite laminates

The purpose of this paper is to develop a sufficiently accurate analysis, which is much simpler than exact three-dimensional analysis, for statics and dynamics of composite laminates. The governing differential equations and boundary conditions are derived by following a variational approach. The displacements are assumed piecewise linear across the thickness and the effects of transverse shear deformations and rotary inertia are included. A procedure for obtaining the general solution of the above governing differential equations in the form of hyperbolic-trigonometric series is given. The accuracy of the present theory is assessed by obtaining results for free vibrations and flexure of simply supported rectangular laminates and comparing them with results from exact three-dimensional analysis.