ON DOUBLE PEAK PROBABILITY DENSITY FUNCTIONS OF DUFFING OSCILLATOR TO COMBINED DETERMINISTIC AND RANDOM EXCITATIONS

The principal resonance of Duffing oscillator to combined deterministic and random external excitation was investigated. The random excitation was taken to be white noise or harmonic with separable random amplitude and phase. The method of multiple scales was used to determine the equations of modulation of amplitude and phase. The one peak probability density function of each of the two stable stationary solutions was calculated by the linearization method. These two one-peak-density functions were combined using the probability of realization of the two stable stationary solutions to obtain the double peak probability density function. The theoretical analysis are verified by numerical results.     

Similarity solutions to viscous flow and heat transfer of nanofluid over nonlinearly stretching sheet

The boundary-layer flow and heat transfer in a viscous fluid containing metallic nanoparticles over a nonlinear stretching sheet are analyzed. The stretching velocity is assumed to vary as a power function of the distance from the origin. The governing partial differential equation and auxiliary conditions are reduced to coupled nonlinear ordinary differential equations with the appropriate corresponding auxiliary conditions. The resulting nonlinear ordinary differential equations (ODEs) are solved numerically. The effects of various relevant parameters, namely, the Eckert number Ec, the solid volume fraction of the nanoparticles φ, and the nonlinear stretching parameter n are discussed. The comparison with published results is also presented. Different types of nanoparticles are studied. It is shown that the behavior of the fluid flow changes with the change of the nanoparticles type.     

Differential quadrature time element method for structural dynamics

An accurate and efficient differential quadraturetime element method(DQTEM) is proposed for solving ordinary differential equations(ODEs),the numerical dissipationand dispersion of DQTEM is much smaller than that of thedirect integration method of single/multi steps.Two methodsof imposing initial conditions are given,which avoids thetediousness when derivative initial conditions are imposed,and the numerical comparisons indicate that the first method,in which the analog equations of initial displacements andvelocities are used to directly replace the differential quadrature(DQ) analog equations of ODEs at the first and the lastsampling points,respectively,is much more accurate thanthe second method,in which the DQ analog equations ofinitial conditions are used to directly replace the DQ analogequations of ODEs at the first two sampling points.On thecontrary to the conventional step-by-step direct integrationschemes,the solutions at all sampling points can be obtainedsimultaneously by DQTEM,and generally,one differentialquadrature time element may be enough for the whole timedomain.Extensive numerical comparisons validate the efficiency and accuracy of the proposed method.     

Differential transformation method for studying flow and heat transfer due to stretching sheet embedded in porous medium with variable thickness,variable thermal conductivity,and thermal radiation

This article presents a numerical solution for the flow of a Newtonian fluid over an impermeable stretching sheet embedded in a porous medium with the power law surface velocity and variable thickness in the presence of thermal radiation. The flow is caused by non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The governing partial differential equations （PDEs） are transformed into a system of coupled non-linear ordinary differential equations （ODEs） with appropriate boundary conditions for various physical parameters. The remaining system of ODEs is solved numerically using a differential transformation method （DTM）. The effects of the porous parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter, and the Prandtl number on the flow and temperature profiles are presented. Moreover, the local skin-friction and the Nusselt numbers are presented. Comparison of the obtained numerical results is made with previously published results in some special cases, with good agreement. The results obtained in this paper confirm the idea that DTM is a powerful mathematical tool and can be applied to a large class of linear and non-linear problems in different fields of science and engineering.     

SEISMIC RANDOM VIBRATION ANALYSIS OF LOCALLY NONLINEAR STRUCTURES

A nonlinear seismic analysis method for complex frame structures subjected to stationary random ground excitations is proposed. The nonlinear elasto-plastic behaviors may take place only on a small part of the structure. The Bouc-Wen differential equation model is used to model the hysteretic characteristics of the nonlinear components. The Pseudo Excitation Method (PEM) is used in solving the linearized random differential equations to replace the solution of the less efficient Lyapunov equation. Numerical results of a real bridge show that .the method proposed is effective for practical engineering analysis.     

NONLINEAR VIBRATION OF THE COUPLED STRUCTURE OF SUSPENDED-CABLE-STAYED BEAM-1:2 INTERNAL RESONANCE

Through the Galerkin method the nonlinear ordinary differential equations （ODEs） in time are obtained from the nonlinear partial differential equations （PDEs） to describe the mo- tion of the coupled structure of a suspended-cable-stayed beam. In the PDEs, the curvature of main cables and the deformation of cable stays are taken into account. The dynamics of the struc- ture is investigated based on the ODEs when the structure is subjected to a harmonic excitation in the presence of both high-frequency principle resonance and 1：2 internal resonance. It is found that there are typical jumps and saturation phenomena of the vibration amplitude in the struc- ture. And the structure may present quasi-periodic vibration or chaos, if the stiffness of the cable stays membrane and frequency of external excitation are disturbed.     

Flow of micropolar fluid between two orthogonally moving porous disks

The unsteady,laminar,incompressible,and two-dimensional flow of a micropolar fluid between two orthogonally moving porous coaxial disks is considered.The extension of von Karman’s similarity transformations is used to reduce the governing partial differential equations(PDEs) to a set of non-linear coupled ordinary differential equations(ODEs) in the dimensionless form.The analytical solutions are obtained by employing the homotopy analysis method(HAM).The effects of various physical parameters such as the expansion ratio and the permeability Reynolds number on the velocity fields are discussed in detail.     

GURTIN VARIATIONAL PRINCIPLE AND FINITE ELEMENT SIMULATION FOR DYNAMICAL PROBLEMS OF FLUID-SATURATED POROUS MEDIA

Based on the theory of porous media,a general Gurtin variational principle for theinitial boundary value problem of dynamical response of fluid-saturated elastic porous media isdeveloped by assuming infinitesimal deformation and incompressible constituents of the solid andfluid phase.The finite element formulation based on this variational principle is also derived.Asthe functional of the variational principle is a spatial integral of the convolution formulation,thegeneral finite element discretization in space results in symmetrical differential-integral equationsin the time domain.In some situations,the differential-integral equations can be reduced to sym-metrical differential equations and,as a numerical example,it is employed to analyze the reflectionof one-dimensional longitudinal wave in a fluid-saturated porous solid.The numerical results canprovide further understanding of the wave propagation in porous media.     

On the ratio of expectation crossings of random-excited dielectric elastomer balloon

The ratio of expectation crossings of dielectric elastomer balloon excited by random pressure is analytically evaluated in this letter.The Mooney–Rivlin model is adopted to describe the constitutive relation while the random pressure is described by Gaussian white noise.Through a specific transformation,the stochastic differential equations for the total energy and phase are derived.With the application of the stochastic averaging,the system total energy is then approximated by a one-dimensional diffusion process.Solving the associated Fokker–Planck–Kolmogorov(FPK)equation yields the stationary probability density of the system total energy.The ratio of expectation crossings is then derived based on the joint stationary probability density of stretch ratio and its ratio of change.The efficacy and accuracy of the proposed procedure are verified by comparing with the results from Monte Carlo simulation(MCS).     

ELASTIC MEDIA WITH RANDOMLY DISTRIBUTED DEFFCTS

In this paper,the elastic field in a solid with randomly distributed defects is derivedThese defects are composed of cavities and microcracks,whose locations,ortentation andsize are random variables.The Random Point Field Model is proposed proposed to describe therandom defects,and the basic equations for elastic field in a random defect medium aredeveloped.Two examples are studied in detail.One is a solid with random microcracks andthe other is a solid with ellipsoidal cavities.     

Random wave propagation in a viscoelastic layered half space

An extremely efficient and accurate solution method is presented for the propagation of stationary random waves in a viscoelastic, transversely isotropic and stratified half space. The efficiency and accuracy are obtained by using the pseudo excitation method (PEM) with the precise integration method (PIM). The solid is multi-layered and located above a semi-infinite space. The excitation sources form a random field which is stationary in the time domain. PEM is used to transform the random wave equation into deterministic equations. In the frequency-wavenumber domain, these equations are ordinary differential equations which can be solved precisely by using PIM. The power spectral densities (PSDs) and the variances of the ground responses can then be computed. The paper presents the full theory and gives results for instructive examples. The comparison between the analytical solutions and the numerical results confirms that the algorithm presented in this paper has exceptionally high precision. In addition, the numerical results presented show that: surface waves are very important for the wave propagation problem discussed; the ground displacement PSDs and variances are significant over bigger regions in the spatial domain when surface waves exist; and as the depth of the source increases the ground displacement PSDs decrease and the regions over which they have significant effect become progressively more restricted to low frequencies while becoming more widely distributed in the spatial domain.     

ARC-length method for differential equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｒｄｉｎａｒｙａｎｄｐａｒｔｉａｌｄｉｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｏｆｃｏｎｔｉｎｕｕｍｐｒｏｂｌｅｍａｒｅｏｆｔｅｎｗｉｔｈｃｅｒｔａｉｎｔｙｐｅｓｏｆｓｉｎｇｕｌａｒｉｔｙａｓｓｔｉｆｐｒｏｐｅｒｔｙ,ｏｒ...     

Analytical solution for bending problem of moderately thick composite annular sector plates with general boundary conditions and loadings using multi-term extended Kantorovich method

An analytical solution for bending of composite sector plates is presented using multi-term extended Kantorovich method (MTEKM). The governing equations are derived using the displacement field of the first-order shear deformation theory and converted into two sets of coupled ordinary differential equations (ODEs) utilizing MTEKM. Next, an analytical iterative procedure is presented for solving the derived sets of ODEs based on state-space method. Numerous examples are studied by the present method, and as special cases, solid sector and rectangular plates are also investigated. Next, the results obtained by the present method are compared to those of finite element method and other results available in the literature. It is found that the present method has a high convergence rate as well as good accuracy in all cases.     

Stationary Solutions for Stochastic Differential Equations Driven by Lévy Processes

In the paper, stationary solutions of stochastic differential equations driven by Lévy processes are considered. And the existence of these stationary solutions follows from the theory of random dynamical systems and their attractors. Moreover, under a one-sided Lipschitz continuity condition and a temperedness condition, Itô and Marcus stochastic differential equations driven by Lévy processes are proved to have stationary solutions. Besides, continuous dependence of stationary solutions on drift coefficients of these equations is presented.     

Neural network as a function approximator and its application in solving differential equations

A neural network(NN) is a powerful tool for approximating bounded continuous functions in machine learning. The NN provides a framework for numerically solving ordinary differential equations(ODEs) and partial differential equations(PDEs)combined with the automatic differentiation(AD) technique. In this work, we explore the use of NN for the function approximation and propose a universal solver for ODEs and PDEs. The solver is tested for initial value problems and boundary value problems of ODEs, and the results exhibit high accuracy for not only the unknown functions but also their derivatives. The same strategy can be used to construct a PDE solver based on collocation points instead of a mesh, which is tested with the Burgers equation and the heat equation(i.e., the Laplace equation).     

Solution of the incompressible Navier–Stokes equations by the method of lines

The numerical method of lines (NUMOL) is a numerical technique used to solve efficiently partial differential equations. In this paper, the NUMOL is applied to the solution of the two‐dimensional unsteady Navier–Stokes equations for incompressible laminar flows in Cartesian coordinates. The Navier–Stokes equations are first discretized (in space) on a staggered grid as in the Marker and Cell scheme. The discretized Navier–Stokes equations form an index 2 system of differential algebraic equations, which are afterwards reduced to a system of ordinary differential equations (ODEs), using the discretized form of the continuity equation. The pressure field is computed solving a discrete pressure Poisson equation. Finally, the resulting ODEs are solved using the backward differentiation formulas. The proposed method is illustrated with Dirichlet boundary conditions through applications to the driven cavity flow and to the backward facing step flow. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulation of electrically driven jet using Chebyshev collocation method

The model of electrically driven jet is governed by a series of quasi 1D dimensionless partial differential equations (PDEs). Following the method of lines, the Chebyshev collocation method is employed to discretize the PDEs and obtain a system of differential-algebraic equations (DAEs). By differentiating constrains in DAEs twice, the system is transformed into a set of ordinary differential equations (ODEs) with invariants. Then the implicit differential equations solver “ddaskr” is used to solve the ODEs and post-stabilization is executed at the end of each step. Results show the distributions of radius, linear charge density, stretching ratio and also the horizontal velocity at a time point. Meanwhile, the spiral and expanding projections to X-Y plane of the jet centerline suggest the occurring of bending instability.     

Lyapunov, Bohl and Sacker-Sell Spectral Intervals for Differential-Algebraic Equations

Lyapunov and exponential dichotomy spectral theory is extended from ordinary differential equations (ODEs) to nonautonomous differential-algebraic equations (DAEs). By using orthogonal changes of variables, the original DAE system is transformed into appropriate condensed forms, for which concepts such as Lyapunov exponents, Bohl exponents, exponential dichotomy and spectral intervals of various kinds can be analyzed via the resulting underlying ODE. Some essential differences between the spectral theory for ODEs and that for DAEs are pointed out. It is also discussed how numerical methods for computing the spectral intervals associated with Lyapunov and Sacker-Sell (exponential dichotomy) can be extended from those methods proposed for ODEs. Some numerical examples are presented to illustrate the theoretical results.