A PARAMETRIC STUDY ON THE AXISYMMETRIC MODES OF VIBRATION OF MULTI-LAYERED SPHERICAL SHELLS WITH LIQUID CORES OF RELEVANCE TO HEAD IMPACT MODELLING

The free vibration of spheres composed of inviscid compressible liquid cores surrounded by spherical layers of linear elastic, homogeneous and isotropic materials are studied using three-dimensional elasticity equations. The exact three-dimensional equations are first derived for an N -layered sphere with a liquid core and an extensive parametric study is then presented for the first few natural frequencies of the spheroidal modes of vibration. Non-dimensional frequency parameters are compared with values obtained using lower order membrane and shell theories. It is shown that for a remarkably wide range of geometric and material parameters, which encompasses values typical for the human head, the first ovalling mode of a fluid-filled shell behaves like a membrane filled with incompressible fluid and a simple closed-form expression is derived which closely approximates the natural frequencies obtained using the exact three-dimensional equations.     

Exact spectral elements for follower tension buckling by power series

The spectral dynamic stiffness method using exact solutions of the governing equations as shape functions has been popular for vibration and dynamic stability analyses of framed structures consisting of uniform members. Since non-uniform members do not generally have closed form solutions, special cases only have been considered. However, exact solutions are still possible for generally non-uniform members using power series. The paper studies the exact dynamic stability of columns with distributed axial force by power series. Both uniform and distributed, compression and tension, and conservative and non-conservative axial forces are considered. Interaction diagrams of various kinds of axial loads on the natural frequencies including different intensities of the distributed loads and degree of tangency are given. Follower tension buckling is reported for the first time. It is found that the power series outperforms the dynamic stiffness method in terms of versatility in applications and numerical stability at the very low and high ends of the frequency spectrum.     

New analytical solutions for nonlinear physical models of the coupled Higgs equation and the Maccari system via rational exp <Emphasis Type="BoldItalic">(−φ(η))</Emphasis>-expansion method

In this article, a variety of solitary wave solutions are found for some nonlinear equations. In mathematical physics, we studied two complex systems, the Maccari system and the coupled Higgs field equation. We construct sufficient exact solutions for nonlinear evolution equations. To study travelling wave solutions, we used a fractional complex transform to convert the particular partial differential equation of fractional order into the corresponding partial differential equation and the rational exp (?φ(η))-expansion method is implemented to find exact solutions of nonlinear equation. We find hyperbolic, trigonometric, rational and exponential function solutions using the above equation. The results of various studies show that the suggested method is very effective and can be used as an alternative for finding exact solutions of nonlinear equations in mathematical physics. A comparative study with the other methods gives validity to the technique and shows that the method provides additional solutions. Graphical representations along with the numerical data reinforce the efficacy of the procedure used. The specified idea is very effective, pragmatic for partial differential equations of fractional order and could be protracted to other physical phenomena.     

Radial anomalous diffusion in an annulus

In this study, time fractional radial diffusion has been modeled in cylindrical coordinates in order to analyze the anomalous diffusion in an annulus. By using an integral transform technique, the analytical solution of the concentration distribution formula is obtained. The establishing of the concentration distribution is found to be controlled by the fractional derivative α, and the influences of α on the concentration field, the total amount diffused and the quantity of mass passing through the inner wall are presented graphically and studied in detail. Asymptotic expressions for the exact solutions are developed in order to explain the numerical results at small and large time, respectively, and the physical mechanism explanation for the paradoxical behavior shown in the numerical results is given.     

Higher accuracy analytical approximations to a nonlinear oscillator with discontinuity by He's homotopy perturbation method

He's homotopy perturbation method is used to calculate higher-order approximate periodic solutions of a nonlinear oscillator with discontinuity for which the elastic force term is proportional to sgn(x). We find He's homotopy perturbation method works very well for the whole range of initial amplitudes, and the excellent agreement of the approximate frequencies and periodic solutions with the exact ones has been demonstrated and discussed. Only one iteration leads to high accuracy of the solutions with a maximal relative error for the approximate period of less than 1.56% for all values of oscillation amplitude, while this relative error is 0.30% for the second iteration and as low as 0.057% when the third-order approximation is considered. Comparison of the result obtained using this method with those obtained by different harmonic balance methods reveals that He's homotopy perturbation method is very effective and convenient.     

An exact solution for the natural frequencies and mode shapes of an immersed elastically restrained wedge beam carrying an eccentric tip mass with mass moment of inertia

In general, the exact solutions for natural frequencies and mode shapes of non-uniform beams are obtainable only for a few types such as wedge beams. However, the exact solution for the natural frequencies and mode shapes of an immersed wedge beam is not obtained yet. This is because, due to the “added mass” of water, the mass density of the immersed part of the beam is different from its emerged part. The objective of this paper is to present some information for this problem. First, the displacement functions for the immersed part and emerged part of the wedge beam are derived. Next, the force (and moment) equilibrium conditions and the deflection compatibility conditions for the two parts are imposed to establish a set of simultaneous equations with eight integration constants as the unknowns. Equating to zero the coefficient determinant one obtains the frequency equation, and solving the last equation one obtains the natural frequencies of the immersed wedge beam. From the last natural frequencies and the above-mentioned simultaneous equations, one may determine all the eight integration constants and, in turn, the corresponding mode shapes. All the analytical solutions are compared with the numerical ones obtained from the finite element method and good agreement is achieved. The formulation of this paper is available for the fully or partially immersed doubly tapered beams with square, rectangular or circular cross-sections. The taper ratio for width and that for depth may also be equal or unequal.     

Convergence studies on static and dynamic analyses of plates by using the U-transformation and the finite difference method

This paper presents the exact static and dynamic analyses of simply supported rectangular plates. The analytical solutions for displacements, buckling loads, natural frequencies and dynamic responses are obtained by using the double U-transformation method and the finite difference method. After an equivalent system with cyclic periodicity in two directions is established, the difference governing equation for such an equivalent system can be uncoupled by applying the double U-transformation. Then the exact analytical finite difference solutions, the exact error expressions and the exact convergence rates are derived. These results cannot be obtained if other methods are used instead.     

Solution for fractional Zakharov–Kuznetsov equations by using two reliable techniques

In this paper, the approximated analytical solutions for nonlinear dispersive fractional Zakharov-Kuznetsov (FZK(a, b, c)) equations are obtained with the help of two novel techniques, called fractional natural decomposition method (FNDM) and q-homotopy analysis transform method (q-HATM). In order to validate and illustrate the efficiency of the proposed techniques, we consider two special cases FZK(2, 2, 2) and FZK(3, 3, 3) in FZK(a, b, c). The numerical simulations have been conducted in order to verify the proposed techniques are reliable and accurate. The outcomes are revealed through the plots and tables. The comparison between the obtained solutions with the exact solutions exhibits that, both the featured schemes are efficient and effective in solving nonlinear complex problems.     

Ground-State Entropy of ±J Ising Lattices by Monte Carlo Simulations

An accurate numerical calculation of the ground-state entropy associated to two-dimensional ±J Ising lattices is presented. The method is based on the use of the thermodynamic integration method. Total energy is calculated by means of the Monte Carlo method. Then the entropy (or degeneracy) of a state of interest is obtained by using thermodynamic integration starting at a known reference state. Results for small sizes are compared to exact values obtained by exhaustive scanning of the entire ground-state manifold, which serves as a test for the reliability of the simulation model developed here. The close agreement between simulated and exact results for energy and remnant entropy supports the validity of the technique used for describing the properties of ±J Ising lattices at the fundamental level. Finally, the results are extrapolated in order to estimate tendencies for larger systems.     

Stability and bifurcation of an axially moving beam tuned to three-to-one internal resonances

This study analyzed the nonlinear vibration of an axially moving beam subject to periodic lateral force excitations. Attention is paid to the fundamental and subharmonic resonances, since the excitation frequency is close to the first two natural frequencies of the system. The incremental harmonic balance (IHB) method was used to evaluate the nonlinear dynamic behaviour of the axially moving beam. The stability and bifurcations of the periodic solutions for given parameters were determined by the multivariable Floquet theory using Hsu’s method. The solutions obtained from the IHB method agreed very well with those obtained from numerical integration. Furthermore, numerical examples are given to illustrate the effects of the three-to-one internal resonance on the response of the system.     

Exact solutions of a nonpolynomially nonlinear Schrödinger equation

A nonlinear generalisation of Schrödinger's equation had previously been obtained using information-theoretic arguments. The nonlinearities in that equation were of a nonpolynomial form, equivalent to the occurrence of higher-derivative nonlinear terms at all orders. Here we construct some exact solutions to that equation in 1+1$1+1$ dimensions. On the half-line, the solutions resemble (exponentially damped) Bloch waves even though no external periodic potential is included. The solutions are nonperturbative as they do not reduce to solutions of the linear theory in the limit that the nonlinearity parameter vanishes. An intriguing feature of the solutions is their infinite degeneracy: for a given energy, there exists a very large arbitrariness in the normalisable wavefunctions. We also consider solutions to a q-deformed version of the nonlinear equation and discuss a natural discretisation implied by the nonpolynomiality. Finally, we contrast the properties of our solutions with other solutions of nonlinear Schrödinger equations in the literature and suggest some possible applications of our results in the domains of low-energy and high-energy physics.     

Solutions of the full- and quasi-vector wave equations based on H-field for optical waveguides by using mapped Galerkin method

The mapped Galerkin method in solving the full-vector and quasi-vector wave equations in terms of transverse magnetic fields (H-formulation) for optical waveguides with step-index profiles is described. By transforming the whole x-y space onto a unit square and using two-dimensional Fourier series expansion, the modal distributions and propagation constants for optical waveguides are obtained in the absence of boundary truncation. Results for step-index circular fiber, buried rectangular waveguide, and optical rib waveguide are presented. Solutions are good agreed with exact solutions and numerical results by using vector nonlinear iterative method, Fourier operator transform method, and vector beam propagation method.     

Dynamic analysis of a multi-span uniform beam carrying a number of various concentrated elements

This paper employs the numerical assembly method (NAM) to determine the “exact” frequency–response amplitudes of a multiple-span beam carrying a number of various concentrated elements and subjected to a harmonic force, and the exact natural frequencies and mode shapes of the beam for the case of zero harmonic force. First, the coefficient matrices for the intermediate concentrated elements, pinned support, applied force, left-end support and right-end support of a beam are derived. Next, the overall coefficient matrix for the whole vibrating system is obtained using the numerical assembly technique of the conventional finite element method (FEM). Finally, the exact dynamic response amplitude of the forced vibrating system corresponding to each specified exciting frequency of the harmonic force is determined by solving the simultaneous equations associated with the last overall coefficient matrix. The graph of dynamic response amplitudes versus various exciting frequencies gives the frequency–response curve for any point of a multiple-span beam carrying a number of various concentrated elements. For the case of zero harmonic force, the above-mentioned simultaneous equations reduce to an eigenvalue problem so that natural frequencies and mode shapes of the beam can also be obtained.     

Exact and Asymptotic Solutions for the Field of a Point Source in a Transversely Isotropic Medium

Exact and asymptotic solutions are obtained for the acoustic field generated by an isotropic pulsed point source in an infinite transversely isotropic elastic medium. The exact solution for the displacement field is obtained in the form of a double integral over the horizontal slowness and the frequency by using the method of integral transforms. The calculation of the integral over the horizontal slowness by the method of stationary phase reduces the exact solution to an asymptotic solution that is convenient for numerical calculations. Formulas are given for calculating the spreading factors and the wave fronts of quasi-longitudinal qP-waves and quasi-transverse qSV-waves. With the formulas obtained, the displacement field of a point source is investigated for a particular transversely isotropic medium.     

Generalized asymptotic expansions for coupled wavenumbers in fluid-filled cylindrical shells

Analytical expressions are found for the coupled wavenumbers in an infinite fluid-filled cylindrical shell using the asymptotic methods. These expressions are valid for any general circumferential order (n). The shallow shell theory (which is more accurate at higher frequencies) is used to model the cylinder. Initially, the in vacuo shell is dealt with and asymptotic expressions are derived for the shell wavenumbers in the high- and the low-frequency regimes. Next, the fluid-filled shell is considered. Defining a relevant fluid-loading parameter μ, we find solutions for the limiting cases of small and large μ. Wherever relevant, a frequency scaling parameter along with some ingenuity is used to arrive at an elegant asymptotic expression. In all cases, Poisson's ratio ν is used as an expansion variable. The asymptotic results are compared with numerical solutions of the dispersion equation and the dispersion relation obtained by using the more general Donnell-Mushtari shell theory (in vacuo and fluid-filled). A good match is obtained. Hence, the contribution of this work lies in the extension of the existing literature to include arbitrary circumferential orders (n).     

Exact solution of the Dirac–Weyl equation in graphene under electric and magnetic fields

In this paper, we have obtained exact analytical solutions for the bound states of a graphene Dirac electron in magnetic fields with various q-parameters under an electrostatic potential. In order to solve the time-independent Dirac–Weyl equation, the Nikoforov–Uvarov (NU) and Frobenius methods have been used. We have also investigated the thermodynamic properties by using the Hurwitz zeta function method for one of the states. Finally, some of the numerical results are also shown.     

Comparison between the generalized tanh–coth and the (<Emphasis Type="Italic">G</Emphasis>′/<Emphasis Type="Italic">G</Emphasis>)-expansion methods for solving NPDEs and NODEs

In this paper, we find exact solutions of some nonlinear evolution equations by using generalized tanh–coth method. Three nonlinear models of physical significance, i.e. the Cahn–Hilliard equation, the Allen–Cahn equation and the steady-state equation with a cubic nonlinearity are considered and their exact solutions are obtained. From the general solutions, other well-known results are also derived. Also in this paper, we shall compare the generalized tanh–coth method and generalized (G ^′/G )-expansion method to solve partial differential equations (PDEs) and ordinary differential equations (ODEs). Abundant exact travelling wave solutions including solitons, kink, periodic and rational solutions have been found. These solutions might play important roles in engineering fields. The generalized tanh–coth method was used to construct periodic wave and solitary wave solutions of nonlinear evolution equations. This method is developed for searching exact travelling wave solutions of nonlinear partial differential equations. It is shown that the generalized tanh–coth method, with the help of symbolic computation, provides a straightforward and powerful mathematical tool for solving nonlinear problems.     

An exact frequency analysis of annular plates with small core having elastically restrained outer edge and sliding inner edge

The present study deals with an exact analysis of free transverse vibrations of annular plates having small core and sliding inner edge and the outer edge being elastically restrained based on classical plate theory. This study focuses mainly on the influence of variations in the elastic restraint parameters on the fundamental frequencies of plate vibration. The natural frequencies for the first six modes of annular plate vibrations are computed for different materials and varying values of the radius parameter and these natural frequencies may correspond to either axisymmetric and/or non-axisymmetric modes of plate vibration. The extensive data of values of fundamental frequency parameter presented in this paper is believed to be of use in the design of acoustic underwater transducers, ocean and naval structures, compressor and pump elements, offshore platforms. These results may serve as bench mark values for researchers to validate their results obtained using approximate numerical methods.     

IDENTIFICATION OF DAMPING: PART 4, ERROR ANALYSIS

In previous papers (S. ADHIKARI and J. WOODHOUSE 2001 Journal of Sound and Vibration243, 43-61; 63-88; S. ADHIKARI and J. WOODHOUSE 2002 Journal of Sound and Vibration251, 477-490) methods were proposed to obtain the coefficient matrix for a viscous damping model or a non-viscous damping model with an exponential relaxation function, from measured complex natural frequencies and modes. In all these works, it has been assumed that exact complex natural frequencies and complex modes are known. In reality, this will not be the case. The purpose of this paper is to analyze the sensitivity of the identified damping matrices to measurement errors. By using numerical and analytical studies it is shown that the proposed methods can indeed be expected to give useful results from moderately noisy data provided a correct damping model is selected for fitting. Indications are also given of what level of noise in the measured modal properties is needed to mask the true physical behaviour.     

Free vibration analysis of a cracked beam by finite element method

In this paper, the natural frequencies and mode shapes of a cracked beam are obtained using the finite element method. An ‘overall additional flexibility matrix’, instead of the ‘local additional flexibility matrix’, is added to the flexibility matrix of the corresponding intact beam element to obtain the total flexibility matrix, and therefore the stiffness matrix. Compared with analytical results, the new stiffness matrix obtained using the overall additional flexibility matrix can give more accurate natural frequencies than those resulted from using the local additional flexibility matrix. All the elements in the overall additional flexibility matrix are computed by 128-point (1D) or (128×128)-point (2D) Gauss quadrature, and then further best fitted using the least-squares method. The explicit form best-fitted formulas agree very well with the numerical integration results, and are very convenient for use and valuable for further reference. In addition, the authors constructed a shape function that can perfectly satisfy the local flexibility conditions at the crack locations, which can give more accurate vibration modes.