Rotating beams and non-rotating beams with shared eigenpair for pinned-free boundary condition

In this paper we look for a rotating beam, with pinned-free boundary conditions, whose eigenpair (frequency and mode-shape) is same as that of a uniform non-rotating beam for a particular mode. It is seen that for any given mode, there exists a flexural stiffness function (FSF) for which the ith mode eigenpair of a rotating beam with uniform mass distribution, is identical to that of a corresponding non-rotating beam with same length and mass distribution. Inserting these derived FSF's in a finite element code for a rotating pinned-free beam, the frequencies and mode shapes of a non-rotating pinned-free beam are obtained. For the first mode, a physically realistic equivalent rotating beam is possible, but for higher modes, the FSF has internal singularities. Strategies for addressing these singularities in the FSF for finite element analysis are provided. The proposed functions can be used as test functions for rotating beam codes and also for targeted destiffening of rotating beams.

Boundary perturbations and the Helmholtz equation in three dimensions

We propose an analytic perturbative scheme in the spirit of Lord Rayleigh's work for determining the eigenvalues of the Helmholtz equation in three dimensions inside an arbitrary boundary where the eigenfunction satisfies either the Dirichlet boundary condition or the Neumann boundary condition. Although numerous works are available in the literature for arbitrary boundaries in two dimensions, to the best of our knowledge the formulation in three dimensions is proposed for the first time. In this novel prescription, we have expanded the arbitrary boundary in terms of spherical harmonics about an equivalent sphere and obtained perturbative closed-form solutions at each order for the problem in terms of corrections to the equivalent spherical boundary for both the boundary conditions. This formulation is in parallel with the standard time-independent Rayleigh-Schrodinger perturbation theory. The efficacy of the method is tested by comparing the perturbative values against the numerically calculated eigenvalues for spheroidal, superegg and superquadric shaped boundaries. It is shown that this perturbation works quite well even for wide departure from spherical shape and for higher excited states too. We believe this formulation would find applications in the field of quantum dots and acoustical cavities.

Spin populations in one-dimensional alternant spin systems: A theoretical approach

Spin-density maps, deduced from polarized neutron diffraction experiments, for both the pair and chain compounds of the system Mn2+Cu2+ have been reported recently. These results have motivated us to investigate theoretically the spin populations in such alternant mixed-spin systems. In this paper, we report our studies on the one-dimensional ferrimagnetic systems (S-A,S-B)(N) where hi is the number of AB pairs. We have considered all cases in which the spin Sri takes on allowed values in the range I to 7/2 while the spin S-B is held fixed at 1/2. The theoretical studies have been carried out on the isotropic Heisenberg model, using the density matrix renormalization group method. The effect of the magnitude of the larger spin SA On the quantum fluctuations in both A and B sublattices has been studied as a function of the system size N. We have investigated systems with both periodic and open boundary conditions, the latter with a view to understanding end-of-chain effects. The spin populations have been followed as a function of temperature as well as an applied magnetic field. High-magnetic fields are found to lead to interesting re-entrant behavior. The ratio of spin populations P-A-P-B is not sensitive to temperature at low temperatures.

Static analysis of infinite laminated cylindrical shell panel

This paper presents an approximate three-dimensional elasticity solution for an infinitely long, cross-ply laminated circular cylindrical shell panel with simply supported boundary conditions, subjected to an arbitrary discontinuous transverse loading. The solution is based on the principal assumption that the ratio of the thickness of the lamina to its middle surface radius is negligible compared to unity. The validity of this assumption and the range of application of this approximate solution have been established through a comparison with an exact solution. Results of classical and first-order shear deformation shell theories have been compared with the results of the present solution to bring out the accuracy of these theories. It is also shown that for very shallow shell panels the definition of a thin shell should be based on the ratio of thickness to chord width rather than the ratio of thickness to mean radius.

Transonic Properties of Accretion Disk Around Compact Objects

An accretion flow is necessarily transonic around a black hole.However, around a neutron star it may or may not be transonic, depending on the inner disk boundary conditions influenced by the neutron star. I will discuss various transonic behavior of the disk fluid in general relativistic (or pseudo general relativistic) framework. I will address that there are four types of sonic/critical point. possible to form in an accretion disk. It will be shown that how the fluid properties including location of sonic point's vary with angular momentum of the compact object which controls the overall disk dynamics and outflows.

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.

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.

Flow through a plateau border of cellular foam

Experimentally measured average velocities through plateau borders of stationary cellular foam, when compared with those calculated with the assumption of rigid Plateau Border walls, show that the assumption of rigid walls severely underestimates the velocities. An analysis of the situation wherein plateau border walls have velocities, as decided by the surface viscosity of the system, is presented here. The plateau border is idealized as a pipe of equilateral triangular cross-section with vertices of the triangle having zero velocity. The pertinent form of Navier-Stoke's equations with inhomogeneous boundary conditions and its solution through a procedure of successive approximations is presented in dimensionless form. The solution reduces to the known solution of slow steady flow through a triangular pipe, when surface viscosity is infinite. Results indicate that the assumption of rigid plateau border walls is valid only when value of the inverse of dimensionless surface viscosity is less than 0.044. Beyond that the assumption severely underestimates the flow and the effect of nonrigidity of the wall must be considered.

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.