A new approximate solution technique for randomly excited non-linear oscillators

A new technique is proposed to obtain an approximate probability density for the response of a non-linear oscillator under Gaussian white noise excitations. The random excitations may be either multiplicative (also known as parametric) or additive (also known as external), or both. In this new technique, the original non-linear oscillator is replaced by another oscillator belonging to the class of generalized stationary potential for which the exact solution is obtainable. The replacement oscillator is selected on the basis that the average energy dissipation remains unchanged. Examples are given to illustrate the application of the new procedure. In one of the examples, the new procedure leads to a better approximation than that obtained by stochastic averaging.     

A new approximate solution technique for randomly excited non-linear oscillators—II

The method of weighted residuals is applied to the reduced Fokker-Planck equation associated with a non-linear oscillator, which is subjected to both additive and multiplicative Gaussian white noise excitations. A set of constraints are deduced for obtaining an approximate stationary probability density for the system response. One of the constraints coincides with the previously proposed criterion of dissipation energy balancing, and the others are useful for calculating the equivalent conservative force. It is shown that these constraints imply certain relationships among certain statistical moments; their imposition guarantees that such moments computed from the approximate probability density satisfy the corresponding exact equations derived from the original equation of motion. Moreover, the well-known procedure of stochastic linearization and its improved version of partial linearization are shown to be special cases of this scheme, and they are less accurate since the approximations are not chosen from the entire set of the solution pool of generalized stationary potential. Applications of the scheme are illustrated by examples, and its accuracy is substantiated by Monte Carlo simulation results.     

EPC procedure for PDF solution of non-linear oscillators excited by Poisson white noise

The stationary probability density function (PDF) solution of the responses of non-linear stochastic oscillators subjected to Poisson pulses is analyzed. The PDF solutions are obtained by the exponential-polynomial closure (EPC) method. To assess the effectiveness of the solution procedure numerically, non-linear oscillators are analyzed with different impulse arrival rates, degree of oscillator non-linearity and excitation intensity. Numerical results show that the PDFs obtained with the EPC method yield good agreement with those obtained from Monte Carlo simulation when the polynomial order is 4 or 6. It is also observed that the EPC procedure is the same as the equivalent linearization procedure under Gaussian white noise in the case of the polynomial order being 2.     

A probabilistic linearization method for non-linear systems subjected to additive and multiplicative excitations

A general method to obtain approximate solutions for the random response of non-linear systems subjected to both additive and multiplicative Gaussian white noises is presented. Starting from the concept of linearization, the proposed method of “Probabilistic Linearization” (PL) is based on the replacement of the Fokker–Planck equation of the original non-linear system with an equivalent one relative to a linear system subjected to additive excitation only. By means of the general scheme of the weighted residuals, the unknown coefficients of the equivalent system are determined. Assuming a Gaussian probability density function of the response process and by choosing the weighting functions in a suitable way, the equivalence of the proposed method, called “Gaussian Probabilistic Linearization” (GPL), with the “Gaussian Stochastic Linearization” (GSL) applied to the coefficients of the Itô differential rule is evidenced. In addition, the generalization of the proposed method, called “Generalized Gaussian Probabilistic Linearization” (GGPL), is presented. Numerical applications show as, varying the choice of the weighting functions, it is possible to obtain different linearizations, with a variable degree of accuracy. For the two examples considered, different suitable combinations of the weighting functions lead to different equivalent linear systems, all characterized by the exact solution in terms of variance.     

A new method for the non-linear deflection analysis of an infinite beam resting on a non-linear elastic foundation

The aim of this paper is to develop a new method of analyzing the non-linear deflection behavior of an infinite beam on a non-linear elastic foundation. Non-linear beam problems have traditionally been dealt with by semi-analytical approaches that involve small perturbations or by numerical methods, such as the non-linear finite element method. In this paper, in contrast, a transformed non-linear integral equation that governs non-linear beam deflection behavior is formulated to develop a new method for non-linear solutions. The proposed method requires an iteration to solve non-linear problems, but is fairly simple and straightforward to apply. It also converges quickly, whereas traditional non-linear solution procedures are generally quite complex in application. Mathematical analysis of the proposed method is performed. In addition, illustrative examples are presented to demonstrate the validity of the method developed in the present study.     

Stochastic response determination of U-OWC energy harvesters: a statistical linearization solution treatment accounting for intermittent wave excitation

Nonlinear Dynamics - A novel statistical linearization technique is developed for determining approximately the response statistics and the power output of U-Oscillating Water Column (U-OWC) energy...     

A one-dimensional moving grid solution for the coupled non-linear equations governing multiphase flow in porous media. 1: Model development

A straightforward moving grid finite element method is developed to solve the one-dimensional coupled system of non-linear partial differential equations (PDEs) governing two- and three-phase flow in porous media. The method combines features from a number of self-adaptive grid techniques. These techniques are the equidistribution, the moving grid finite element and the local grid refinement/coarsening methods. Two equidistribution criteria, based on solution gradient and curvature, are employed and nodal distributions are computed iterativcly. Using the developed approach, an intermingle-free nodal distribution is guaranteed. The method involves examination of a single representative gradient to facilitate the application of moving grid algorithms to solve a non-linear coupled set of PDEs and includes a feature to limit mass balance error during nodal redistribution. The finite element part of the developed algorithm is verified against an existing finite difference model. A numerical simulation example involving a single-front two-phase flow problem is presented to illustrate model performance. Additional simulation examples are given in Part 2 of this paper. These examples include single and double moving fronts in two- and three-phase flow systems incorporating source/sink terms. Simulation sensitivity to the moving grid parameters is also explored in Part 2.     

A deformational and configurational framework for geometrically non-linear continuum thermomechanics coupled to diffusion

A general framework encompassing both the (conventional) deformational and configurational settings of continuum mechanics is presented. A systematic application of balance principles over a migrating control volume in the undeformed configuration of the continuum body yields the system of governing equations in the bulk, on the surface and on a coherent interface within the continuum. The equations governing the response of the bulk agree with those of the conventional deformational approach. The localised balance equations are expressed in the configurational setting using a pull-back operator and reformulated in terms of the Eshelby stress. The configurational expression of the dissipation elucidates the energy loss associated with configurational changes. The general framework is introduced by considering the problem of coupled deformation, heat conduction and species diffusion within a geometrically non-linear continuum body intersected by a coherent interface. The nature of the coupling is emphasised throughout the presentation and via an example.     

A general solution for one dimensional chemo-mechanical coupled hydrogel rod

Smart hydrogels are environmentally sensitive hydrogels, which can produce a sensitive response to external stimuli, and often exhibit the characteristics of multi filed coupling. In this paper, a hydrogel rod under chemo-mechanical coupling was analytically studied based on a poroelastical model. The already known constitutive and governing equations were simplified into the one dimensional case, then two different boundary conditions were considered. The expressions of concentration, displacement, chemical potential and stress related to time were obtained in a series form. Examples illustrate the interaction mechanism of chemical and mechanical effect. It was found that there was a balance state in the diffusion of concentration and the diffusion process could lead to the expansion or the stress change of the hydrogel rod.     

A one-dimensional moving grid solution for the coupled non-linear equations governing multiphase flow in porous media. 2: Example simulations and sensitivity analysis

This paper presents numerical examples for the moving grid finite element algorithm derived in Part Ito solve the non-linear coupled set of PDEs governing immiscible multiphase flow in porous media in one dimension. Examples include single- and double-front simulations for two- and three-phase flow regimes and incorporating a mass sink. The modelling approach is shown to achieve significant savings in computation time and memory allocation when compared with fixed grid solutions of equivalent accuracy. This work includes sensitivity analyses for the parameters which are incorporated in the grid adaptation method, including the curvature weights, artificial viscosity and artificial repulsive force. It is found that the curvature weights are exponential functions of the negative ratio of the square root of the domain length to the number of discrete nodes. These weighting parameters are also shown to depend upon the shape of the front. On the basis of the examined simulations, it is recommended that artificial viscosity be neglected in the solution of the coupled non-linear set of PDEs governing multiphase flow in porous media. Similarly, use of a repulsive force is found to be unnecessary in simulations involving the migration of two liquid phases. For multiphase flows incorporating a gas phase it is recommended to use a non-zero value for the repulslive force to avoid development of an ill-conditioned nodal distribution matrix. An equation to evaluate the repulsive force under these circumstances is suggested.     

A solution for the lateral stability of beams using an electrical analog

This paper shows how an electrical-resistance analog based on the “conjugate-beam” concept and developed to simulate the elastic behavior of beams in the plane of loading may be used to examine the lateral buckling behavior of beams.     

A quasilinearization approach to the solution of elastic beams on nonlinear foundations

In this paper, the solution of a beam on nonlinear elastic foundation whose deflection satisfies the nonlinear boundary value problem (1, 2), is studied by means of the theory of quasilinearization. The problem is formulated in Section 2 where conditions for the existence and uniqueness of the solution are stated. In Section 3, the idea of quasilinearization is introduced and the positivity of an associated linear differential operator is investigated. In Section 4 the usual version of quasilinearization, i.e. The Newton-Raphson-Kantorovich sequence, is presented and conditions under which this sequence is monotonically convergent, are established. In Section 5, an alternative successive approximation scheme whose derivation relies on ideas of quasilinearization, is presented. Finally, an example is solved by numerical procedures based in the methods discussed in previous sections.     

A new method for predicting the maximum vibration amplitude of periodic solution of non-linear system

An original method based on the proposed framework for calculating the maximum vibration amplitude of periodic solution of non-linear system is presented. The problem of determining the worst maximum vibration is transformed into a non-linear optimization problem. The harmonic balance method and the Hill method are selected to construct the general non-linear equality and inequality constraints. The resulting constrained maximization problem is then solved by using the MultiStart algorithm. Finally, the effectiveness of the proposed approach is illustrated through two numerical examples. Numerical examples show that the proposed method can, at much lower cost, give results with higher accuracy as compared with numerical results obtained by a parameter continuation method.     

A novel block cryptosystem based on the coupled chaotic map lattice

In this paper, the authors propose a novel block cryptographic scheme based on the coupled chaotic map lattice. A pseudorandom number generator is constructed with the coupled spatiotemporal chaotic map lattice, which depends on the plaintext. The plaintext block is then encrypted by the classical chaotic masking technique. The pseudorandom sequences produced by the generator have excellent statistical properties; therefore, they are suitable for encrypting messages. Moreover, it can be easily implemented in parallel by hardware. Both theoretical analysis and numerical simulation prove the block cryptographic scheme secure and practical.     

A new efficient method to evaluate exact stiffness and mass matrices of non-uniform beams resting on an elastic foundation

In this paper, a new efficient method to evaluate the exact stiffness and mass matrices of a non-uniform Bernoulli–Euler beam resting on an elastic Winkler foundation is presented. The non-uniformity may result from variable cross-section and/or from inhomogeneous linearly elastic material. It is assumed that there is no abrupt variation in the cross-section of the beam so that the Euler–Bernoulli theory is valid. The method is based on the integration of the exact shape functions which are derived from the solution of the axial deformation problem of a non-uniform bar and the bending problem of a non-uniform beam which are both formulated in terms of the two displacement components. The governing differential equations are uncoupled with variable coefficients and are solved within the framework of the analog equation concept. According to this, the two differential equations with variable coefficients are replaced by two linear ones pertaining to the axial and transverse deformation of a substitute beam with unit axial and bending stiffness, respectively, under ideal load distributions. The key point of the method is the evaluation of the two ideal loads which in this work is achieved by approximating them by two polynomials. More specifically, the axial ideal load is approximated by a linear polynomial while the transverse one by a cubic polynomial. The numerical implementation of the method is simple, and the results are compared favorably to those obtained by exact solutions available in literature.     

A method for finding the invariants and exact solutions of coupled non-linear differential equations with applications to dynamical systems

The polynomial invariants of (a set) non-linear differential equations are found by using a direct approach. The integrability of these invariants deserves the integrability of the given set of coupled differential equations. As applications, the Lorenz and Rikitake sets, among others, are studied. New invariants are obtained.     

A novel series of solution methods for solving nonlinear stiff dynamic problems

Nonlinear Dynamics - A novel type of linear multi-step formulas is proposed for solving initial value problems, such as the problems of multi-body systems and vibration systems and a variety of...