Averaging method for the solution of non-linear differential equations with periodic non-harmonic solutions

While Krylov and Bogolyubov used harmonic functions in their averaging method for the approximate solution of weakly non-linear differential equations with oscillatory solution, we apply a similar averaging technique using Jacobi elliptic functions. These functions are also periodic and are exact solutions of strongly non-linear differential equations. The method is used to solve non-linear differential equations with linear and non-linear small dissipative terms and/or with time dependent parameters. It is also shown that quite general dissipative terms can be transformed into time-dependent parameters. As a special example, the Langevin (collisional) equation of motion of electrons in a neutralizing ion background under the influence of a time and space-dependent electric field is presented. The method may also be used for non-linear control theory, dynamic and parametric stabilization of non-linear oscillations in plasma physics, etc.     

Asymptotic solutions to the non-linear cantilever elastica

In this work it is shown that by a series of admissible functional transformations the constructed higher-order strongly non-linear differential equation (ODE), describing the elastica of a cantilever due to a terminal generalized concentrated, as well as to a lateral uniformly distributed loading, is reduced to a first-order non-linear integrodifferential equation consisting of the first intermediate integral of the original equation. The absence of exact analytic solutions in terms of known (tabulated) functions of the above reduced equation leads to the conclusion that there are no exact analytic solutions of this complicated elastica problem. In the limits of small values of the slope parameter of the deflected elastica, we expand asymptotically the above integrodifferential equation to non-linear ODEs of the generalized Emden–Fowler types, exact analytic solutions of which are constructed in parametric form.     

Exact solution of the multi-cracked Euler–Bernoulli column

The use of distributions (generalized functions) is a powerful tool to treat singularities in structural mechanics and, besides providing a mathematical modelling, their capability of leading to closed form exact solutions is shown in this paper. In particular, the problem of stability of the uniform Euler–Bernoulli column in presence of multiple concentrated cracks, subjected to an axial compression load, under general boundary conditions is tackled. Concentrated cracks are modelled by means of Dirac’s delta distributions. An integration procedure of the fourth order differential governing equation, which is not allowed by the classical distribution theory, is proposed. The exact buckling mode solution of the column, as functions of four integration constants, and the corresponding exact buckling load equation for any number, position and intensity of the cracks are presented. As an example a parametric study of the multi-cracked simply supported and clamped–clamped Euler–Bernoulli columns is presented.     

A class of coupled nonlinear Schr(o)dinger equations: Painlevé property,exact solutions,and application to atmospheric gravity waves

The Painlev(e) integrability and exact solutions to a coupled nonlinear Schrodinger (CNLS) equation applied in atmospheric dynamics are discussed. Some parametric restrictions of the CNLS equation are given to pass the Painleve test. Twenty periodic cnoidal wave solutions are obtained by applying the rational expansions of fundamental Jacobi elliptic functions. The exact solutions to the CNLS equation are used to explain the generation and propagation of atmospheric gravity waves.     

A class of coupled nonlinear Schrödinger equations: Painlevé property,exact solutions,and application to atmospheric gravity waves

The Painleve integrability and exact solutions to a coupled nonlinear Schrodinger （CNLS） equation applied in atmospheric dynamics are discussed. Some parametric restrictions of the CNLS equation are given to pass the Painleve test. Twenty periodic cnoidal wave solutions are obtained by applying the rational expansions of fundamental Jacobi elliptic functions. The exact solutions to the CNLS equation are used to explain the generation and propagation of atmospheric gravity waves.     

Note on the capillary curve

The plane capillary curve of a stationary system liquid-gas enclosed between two horizontal plates is investigated. Starting from the Laplace-Young equation of the theory of capillarity it is shown that the governing equation of the problem can be reduced to the equation of motion of a simple pendulum oscillating with finite amplitude. Integration of this non-linear differential equation is performed by means of the Jacobian elliptic functions.     

Constructing transient response probability density of non-linear system through complex fractional moments

The probability density function for transient response of non-linear stochastic system is investigated through the stochastic averaging and Mellin transform. The stochastic averaging based on the generalized harmonic functions is adopted to reduce the system dimension and derive the one-dimensional Itô stochastic differential equation with respect to amplitude response. To solve the Fokker–Plank–Kolmogorov equation governing the amplitude response probability density, the Mellin transform is first implemented to obtain the differential relation of complex fractional moments. Combining the expansion form of transient probability density with respect to complex fractional moments and the differential relations at different transform parameters yields a set of closed-form first-order ordinary differential equations. The complex fractional moments which are determined by the solution of the above equations can be used to directly construct the probability density function of system response. Numerical results for a van der Pol oscillator subject to stochastically external and parametric excitations are given to illustrate the application, the convergence and the precision of the proposed procedure.     

Prediction of the response of non-linear oscillators under stochastic parametric and external excitations

The method of equivalent external excitation is derived to predict the stationary variances of the states of non-linear oscillators subjected to both stochastic parametric and external excitations. The oscillator is interpreted as one which is excited solely by an external zero-mean stochastic process. The Fokker-Planck-Kolmogorov equation is then applied to solve for the density functions and match the stationary variances of the states. Four examples which include polynomial, non-polynomial, and Duffing type non-linear oscillators are used to illustrate this approach. The validity of the present approach is compared with some exact solutions and with Monte Carlo simulations.     

Solutions and conservation laws of a (3+1)-dimensional Zakharov–Kuznetsov equation

In this paper, we study a nonlinear evolution partial differential equation, namely the (3+1)-dimensional Zakharov–Kuznetsov equation. Kudryashov method together with Jacobi elliptic function method is used to obtain the exact solutions of the (3+1)-dimensional Zakharov–Kuznetsov equation. Furthermore, the conservation laws of the (3+1)-dimensional Zakharov–Kuznetsov equation are obtained by using the multiplier method.     

Free vibration of non-uniform column using DQM

Most of engineering problems are governed by a set of partial differential equations with proper boundary conditions. The present work is concerned with free vibration analysis of non-uniform column resting on elastic foundation and subjected to follower force. The used method of solution is the differential quadrature method (DQM). Formulation of the problem is introduced. The results obtained and compared with the exact solution and traditional numerical techniques such as finite element method. The parametric study is used to investigate the effect of column geometry on the natural frequencies, the mode shapes and the critical load.     

Internal-external resonance of beams on non-linear viscoelastic foundation traversed by moving load

Vibration of a finite Euler–Bernoulli beam, supported by non-linear viscoelastic foundation traversed by a moving load, is studied and the Galerkin method is used to discretize the non-linear partial differential equation of motion. Subsequently, the solution is obtained for different harmonics using the Multiple Scales Method (MSM) as one of the perturbation techniques. Free vibration of a beam on non-linear foundation is investigated and the effects of damping and non-linear stiffness of the foundation on the responses are examined. Internal-external resonance condition is then stated and the frequency responses of different harmonics are obtained by MSM. Different conditions of the external resonance are studied and a parametric study is carried out for each case. The effects of damping and non-linear stiffness of the foundation as well as the magnitude of the moving load on the frequency responses are investigated. Finally, a thorough local stability analysis is performed on the system.     

Stochastic finite element analysis of non-linear plane trusses

—This study considers the responses of geometrically and materially non-linear plane trusses under random excitations. The stress-strain law in the inelastic range is based on an explicit differential equation model. After a total Lagrangian finite element discretization, the nodal displacements satisfy a system of stochastic non-linear ordinary differential equations with right-hand-sides given by random functions of time. The exact solution of the above stochastic differential equation is generally difficult to obtain. To seek an approximate solution with good accuracy and reasonable computational effort, the stochastic linearization method is used to find the first and second statistical moments (i.e. the mean vector and the one-time covariance matrix) of the nodal displacements. Results of simple structures under Gaussian white-noise excitation indicate that the proposed method has good accuracy (generally underestimates the r.m.s. stationary response by 5–14%) and requires only a small fraction of the computation time of the time-history Monte-Carlo method.     

On axisymmetric flow of Sisko fluid over a radially stretching sheet

This study presents an analysis of the axisymmetric flow of a non-Newtonian fluid over a radially stretching sheet. The momentum equations for two-dimensional flow are first modeled for Sisko fluid constitutive model, which is a combination of power-law and Newtonian fluids. The general momentum equations are then simplified by invoking the boundary layer analysis. Then a non-linear ordinary differential equation governing the axisymmetric boundary layer flow of Sisko fluid over a radially stretching sheet is obtained by introducing new suitable similarity transformations. The resulting non-linear ordinary differential equation is solved analytically via the homotopy analysis method (HAM). Closed form exact solution is then also obtained for the cases n=0 and 1. Analytical results are presented for the velocity profiles for some values of governing parameters such as power-law index, material parameter and stretching parameter. In addition, the local skin friction coefficient for several sets of the values of physical parameter is tabulated and analyzed. It is shown that the results presented in this study for the axisymmetric flow over a radially non-linear stretching sheet of Sisko fluid are quite general so that the corresponding results for the Newtonian fluid and the power-law fluid can be obtained as two limiting cases.     

Nonlinear flexural waves and chaos behavior in finite-deflection Timoshenko beam

Based on the Timoshenko beam theory, the finite-deflection and the axial inertia are taken into account, and the nonlinear partial differential equations for flexural waves in a beam are derived. Using the traveling wave method and integration skills, the nonlinear partial differential equations can be converted into an ordinary differential equation. The qualitative analysis indicates that the corresponding dynamic system has a heteroclinic orbit under a certain condition. An exact periodic solution of the nonlinear wave equation is obtained using the Jacobi elliptic function expansion. When the modulus of the Jacobi elliptic function tends to one in the degenerate case, a shock wave solution is given. The small perturbations are further introduced, arising from the damping and the external load to an original Hamilton system, and the threshold condition of the existence of the transverse heteroclinic point is obtained using Melnikov＇s method. It is shown that the perturbed system has a chaotic property under the Smale horseshoe transform.     

Non-linear beam-type vibrations of long cylindrical shells

This paper presents a closed-form solution of the problem of free beam-type vibration of a long cylindrical shell subject to uniform axial tension, uniform internal pressure, and elastic axial restraint. The shell is assumed to be clamped at the ends. The bending moment-curvature relationship employed in this paper is non-linear due to the effect of ovalization (flattening) of initially circular cross sections. Another non-linear effect taken into account is the stretching of the shell axis which is caused by the axial restraint.The analysis results in a cubic differential equation; the frequency of the solution of this equation is found exactly using elliptic integrals.     

Accurate analytical perturbation approach for large amplitude vibration of functionally graded beams

The present work derives the accurate analytical solutions for large amplitude vibration of thin functionally graded beams. In accordance with the Euler–Bernoulli beam theory and the von Kármán type geometric non-linearity, the second-order ordinary differential equation having odd and even non-linearities can be formulated through Hamilton's principle and Galerkin's procedure. This ordinary differential equation governs the non-linear vibration of functionally graded beams with different boundary constraints. Building on the original non-linear equation, two new non-linear equations with odd non-linearity are to be constructed. Employing a generalised Senator–Bapat perturbation technique as an ingenious tool, two newly formulated non-linear equations can be solved analytically. By selecting the appropriate piecewise approximate solutions from such two new non-linear equations, the analytical approximate solutions of the original non-linear problem are established. The present solutions are directly compared to the exact solutions and the available results in the open literature. Besides, some examples are selected to confirm the accuracy and correctness of the current approach. The effects of boundary conditions and vibration amplitudes on the non-linear frequencies are also discussed.     

Exact parametric analytic solutions of the elastica ODEs for bars including effects of the transverse deformation

We succeed in constructing exact parametric analytic solutions for the non-linear ordinary differential equations governing the elastica response of a cantilever due to a generalized end loading by taking into account the effects of transverse deformation. Application to the case of the eccentric buckling of a cantilever by taking into account the above influences is developed.     

Effects of radiation and heat source on MHD flow of a viscoelastic liquid and heat transfer over a stretching sheet

We study the MHD flow and also heat transfer in a viscoelastic liquid over a stretching sheet in the presence of radiation. The stretching of the sheet is assumed to be proportional to the distance from the slit. Two different temperature conditions are studied, namely (i) the sheet with prescribed surface temperature (PST) and (ii) the sheet with prescribed wall heat flux (PHF). The basic boundary layer equations for momentum and heat transfer, which are non-linear partial differential equations, are converted into non-linear ordinary differential equations by means of similarity transformation. The resulting non-linear momentum differential equation is solved exactly. The energy equation in the presence of viscous dissipation (or frictional heating), internal heat generation or absorption, and radiation is a differential equation with variable coefficients, which is transformed to a confluent hypergeometric differential equation using a new variable and using the Rosseland approximation for the radiation. The governing differential equations are solved analytically and the effects of various parameters on velocity profiles, skin friction coefficient, temperature profile and wall heat transfer are presented graphically. The results have possible technological applications in liquid-based systems involving stretchable materials.     

Exact solutions for buckling of non-uniform columns under axial concentrated and distributed loading

In this paper, the buckling problem of non-uniform columns subjected to axial concentrated and distributed loading is studied. The expression for describing the distribution of flexural stiffness of a non-uniform column is arbitrary, and the distribution of axial forces acting on the column is expressed as a functional relation with the distribution of flexural stiffness and vice versa. The governing equation for buckling of a non-uniform column with arbitrary distribution of flexural stiffness or axial forces is reduced to a second-order differential equation without the first-order derivative by means of functional transformations. Then, this kind of differential equation is reduced to Bessel equations and other solvable equations for 12 cases, several of which are important in engineering practice. The exact solutions that represent a class of exact functional solutions for the buckling problem of non-uniform columns subjected to axial concentrated and distributed loading are obtained. In order to illustrate the proposed method, a numerical example is given in the last part of this paper.     

Flow distribution in parallel and reverse flow manifolds

Numerical predictions of the flow distribution in parallel and reverse flow manifold systems have been obtained by a novel finite-difference procedure. A one-dimensional elliptic formulation is proposed. An iterative numerical scheme solves the differential equations for momentum in the longitudinal direction and continuity, while in the lateral direction an integral equation for momentum is solved. The predicted pressure distributions compare well with the available experimental results. A parametric computation has been performed to demonstrate the effect of the governing non-dimensional parameters on the flow distribution.