Construction of approximate analytical solutions to strongly nonlinear damped oscillators

An analytical approximate method for strongly nonlinear damped oscillators is proposed. By introducing phase and amplitude of oscillation as well as a bookkeeping parameter, we rewrite the governing equation into a partial differential equation with solution being a periodic function of the phase. Based on combination of the Newton’s method with the harmonic balance method, the partial differential equation is transformed into a set of linear ordinary differential equations in terms of harmonic coefficients, which can further be converted into systems of linear algebraic equations by using the bookkeeping parameter expansion. Only a few iterations can provide very accurate approximate analytical solutions even if the nonlinearity and damping are significant. The method can be applied to general oscillators with odd nonlinearities as well as even ones even without linear restoring force. Three examples are presented to illustrate the usefulness and effectiveness of the proposed method.     

Accurate Higher-Order Analytical Approximate Solutions to Large-Amplitude Oscillating Systems with a General Non-Rational Restoring Force

A new approximate analytical approach for accurate higher-order nonlinear solutions of oscillations with large amplitude is presented in this paper. The oscillatory system is subjected to a non-rational restoring force. This approach is built upon linearization of the governing dynamic equation associated with the method of harmonic balance. Unlike the classical harmonic balance method, simple linear algebraic equations instead of nonlinear algebraic equations are obtained upon linearization prior to harmonic balancing. This approach also explores large parameter regions beyond the classical perturbation methods which in principle are confined to problems with small parameters. It has significant contribution as there exist many nonlinear problems without small parameters. Through some examples in this paper, we establish the general approximate analytical formulas for the exact period and periodic solution which are valid for small as well as large amplitudes of oscillation.     

Analytical approximation to large-amplitude oscillation of a non-linear conservative system

This paper deals with non-linear oscillation of a conservative system having inertia and static non-linearities. By combining the linearization of the governing equation with the method of harmonic balance, we establish analytical approximate solutions for the non-linear oscillations of the system. Unlike the classical harmonic balance method, linearization is performed prior to proceeding with harmonic balancing, thus resulting in a set of linear algebraic equations instead of one of non-linear algebraic equations. Hence, we are able to establish analytical approximate formulas for the exact frequency and periodic solution. These analytical approximate formulas show excellent agreement with the exact solutions, and are valid for small as well as large amplitudes of oscillation.     

Nonlinear vibration of a two-mass system with nonlinear stiffnesses

An analytical approach is developed for nonlinear free vibration of a conservative, two-degree-of-freedom mass–spring system having linear and nonlinear stiffnesses. The main contribution of the proposed approach is twofold. First, it introduces the transformation of two nonlinear differential equations of a two-mass system using suitable intermediate variables into a single nonlinear differential equation and, more significantly, the treatment a nonlinear differential system by linearization coupled with Newton’s method and harmonic balance method. New and accurate higher-order analytical approximate solutions for the nonlinear system are established. After solving the nonlinear differential equation, the displacement of two-mass system can be obtained directly from the governing linear second-order differential equation. Unlike the common perturbation method, this higher-order Newton–harmonic balance (NHB) method is valid for weak as well as strong nonlinear oscillation systems. On the other hand, the new approach yields simple approximate analytical expressions valid for small as well as large amplitudes of oscillation unlike the classical harmonic balance method which results in complicated algebraic equations requiring further numerical analysis. In short, this new approach yields extended scope of applicability, simplicity, flexibility in application, and avoidance of complicated numerical integration as compared to the previous approaches such as the perturbation and the classical harmonic balance methods. Two examples of nonlinear two-degree-of-freedom mass–spring system are analyzed and verified with published result, exact solutions and numerical integration data.     

On the harmonic balance with linearization for asymmetric single degree of freedom non-linear oscillators

This paper deals with analytical approximation of non-linear oscillations of conservative asymmetric single degree of freedom systems, using the method of harmonic balance with linearization. This technique which consists of linearizing the governing equations prior to harmonic balance permits us to avoid solving complicated non-linear algebraic equations. But it could be applied only to symmetric oscillations for which it proves to be very simple and effective. This restriction is due to the fact that the method requires an appropriate initial approximate solution as input. Such a solution could not be readily identified for nonsymmetric oscillations, contrary the symmetric case where the fundamental harmonic works well. For these nonsymmetric oscillations, we propose in this paper to consider an initial approximation which consists of a small bias plus the fundamental harmonic. By expanding the corresponding harmonic balance equations respectively to first and second order in the bias, we are able to easily determine the bias and thus the required initial approximate solution that yields consistent solution at higher order. We use three examples to illustrate the proposed approach and reveal its simplicity and its very good convergence.     

Thermomechanically coupled non-linear vibration of plates

An analytical method is presented for the analysis of large amplitude thermomechanically coupled vibrations of rectangular elastic thin plates with various boundary conditions. The field of temperature and deformation are assumed to be coupled, and the transverse and longitudinal deformations are mutually dependent. The fundamental equations of non-linear flexural vibration of a plate stemming from Berger's analysis are coupled with the energy equation. Based on one-term approximate solution technique, the system of non-linear equations is solved by employing the methods of Galerkin and successive approximations. The analytical solutions are compared with those for the linear case and from the numerical analysis to investigate the influence of thermomechanical coupling and large amplitude on the period of plate lateral vibration.     

A New Method for Approximate Analytical Solutions to Nonlinear Oscillations of Nonnatural Systems

This paper deals with nonlinear oscillations of a conservative,nonnatural, single-degree-of-freedom system with odd nonlinearity. Bycombining the linearization of the governing equation with the method ofharmonic balance, we establish approximate analytical solutions for thenonlinear oscillations of the system. Unlike the classical harmonicbalance method, the linearization is performed prior to proceeding withharmonic balancing thus resulting in linear algebraic equations insteadof nonlinear algebraic equations. Hence, we are able to establish theapproximate analytical formulas for the exact period and periodicsolution. These approximate solutions are valid for small as well aslarge amplitudes of oscillation. Two examples are presented toillustrate that the proposed formulas can give excellent approximateresults.     

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.     

电动力学中衰减机械系统的非线性振荡方法

衰减机械系统的非线性振荡可用来研究长约瑟夫逊结的电动力学方程式，而这方程式等同于弱衰减机械系统的非线性振荡。本文应用的方法是将控制方程线性化及结合谐波平衡法（线性谐波平衡法）而产生色散关系，再把平均法应用在弱非线性的耗散系统中得到非常准确的瞬变反应。在此提出的方法不仅考虑能量耗散，而且利用简单的线性代数等式关系来代替冗长及复杂的分析近似解。     

A comparative study of pressure correction and block-implicit finite volume algorithms on parallel computers

A comparison of a new parallel block-implicit method and the parallel pressure correction procedure for the solution of the incompressible Navier–Stokes equations is presented. The block-implicit algorithm is based on a pressure equation. The system of non-linear equation s is solved by Newton's method. For the solution of the linear algebraic systems the Bi-CGSTAB algorithm with incomplete lower–upper (ILU) decomposition of the matrix is applied. Domain decomposition serves as a strategy for the parallelization of the algorithms. Different algorithms for the parallel solution of the linear system of algebraic equations in conjunction with the pressure correction procedure are proposed. Three different flows are predicted with the parallel algorithms. Results and efficiency data of the block-implicit method are compared with the parallel version of the pressure correction algorithm. The block-implicit method is characterized by stable convergence behaviour, high numerical efficiency, insensitivity to relaxation parameters and high spatial accuracy. © 1997 John Wiley & Sons, Ltd.     

Approximate analytical solution in slow-fast system based on modified multi-scale method

A simple, yet accurate modified multi-scale method (MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method (MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simplified system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking (SP), spiking-quiescence (SP-QS), and quiescence (QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method (HBM). The results obtained by the present method are considerably better than those obtained by traditional methods, quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.     

Approximate analytical solution in slow-fast system based on modi?ed multi-scale method

A simple, yet accurate modi?ed multi-scale method(MMSM) for an approximately analytical solution in nonlinear oscillators with two time scales under forced harmonic excitation is proposed. This method depends on the classical multi-scale method(MSM) and the method of variation of parameters. Assuming that the forced excitation is a constant, one could easily obtain the approximate analytical solution of the simpli?ed system based on the traditional MSM. Then, this solution for the oscillator under forced harmonic excitation could be established after replacing the harmonic excitation by the constant excitation. To certify the correctness and precision of the proposed analytical method, the van der Pol system with two scales subject to slowly periodic excitation is investigated; this system presents rich dynamical phenomena such as spiking(SP),spiking-quiescence(SP-QS), and quiescence(QS) responses. The approximate analytical expressions of the three types of responses are given by the MMSM, and it can be found that the precision of the new analytical method is higher than that of the classical MSM and better than that of the harmonic balance method(HBM). The results obtained by the present method are considerably better than those obtained by traditional methods,quantitatively and qualitatively, particularly when the excitation frequency is far less than the natural frequency of the system.     

多尺度法二义性的一种解释

多尺度法是为解决含小参数系统发展起来的应用最广泛的摄动法之一. 在求解高阶近 似方程时,多尺度法一般只求特解. 用多尺度法求解van der Pol 方程的三阶解时 将出现矛盾. 以van der Pol方程为例,证明了忽略一阶修正量中的一阶谐波 项使得混合偏导数不能交换顺序,从而导致了多尺度法的二义性和另一个数学矛盾. 在求解一阶修正量时采用含有一阶谐波项的全解,消除了二义性和该矛盾. 该 方法所求得的近似解与数值解进行了比较,结果非常吻合,验证了其合理性.     

Momentum‐based approximation of incompressible multiphase fluid flows

We introduce a time stepping technique using the momentum as dependent variable to solve incompressible multiphase problems. The main advantage of this approach is that the mass matrix is time‐independent making this technique suitable for spectral methods. A level set method is applied to reconstruct the fluid properties such as density. We also introduce a stabilization method using an entropy‐viscosity technique and a compression technique to limit the flattening of the level set function. We extend our algorithm to immiscible conducting fluids by coupling the incompressible Navier‐Stokes and the Maxwell equations. We validate the proposed algorithm against analytical and manufactured solutions. Results on test cases such as Newton's bucket problem and a variation thereof are provided. Surface tension effects are tested on benchmark problems involving bubbles. A numerical simulation of a phenomenon related to the industrial production of aluminium is presented at the end of the paper.     

Nonlinear analysis of chatter vibration in a cylindrical transverse grinding process with two time delays using?a?nonlinear time transformation method

A nonlinear dynamic system of cylindrical transverse grinding process is studied in this paper. The system consists of a grinding wheel and a workpiece, which are connected to the base by spring-damper elements, interacting with nonlinear normal forces. This two DOF model includes two time delays originated from the regenerative effects of the workpiece and the grinding wheel. Bifurcation points are located using a numerical algorithm by which we can find all the eigenvalues in a given rectangular region on the complex plane for the delayed differential equations. Supercritical bifurcation has been found for some sets of system parameter values. The amplitudes of the limit cycles are predicted using a nonlinear time transformation method, which is similar to the harmonic balance approach in that a periodic solution is approximated by a Fourier series. However, the main difference is that a nonlinear time ? is introduced in the Fourier series rather than the physical time t. The analytical solutions of stable limit cycles up to the third harmonics are compared with numerical simulations for the retarded system. It is shown that the proposed method gives accurate approximate solutions.     

Localized Basis Function Method for Computing Limit Cycle Oscillations

An alternate approach to the standard harmonic balance method (based on Fourier transforms) is proposed. The proposed method begins with an idea similar to the harmonic balance method, i.e. to transform the initial set of differential equations of the dynamics to a set of discrete algebraic equations. However, as distinct from previous harmonic balance techniques, the proposed method uses a set of basis functions which are localized in time and are not necessarily sinusoidal. Also as distinct from previous harmonic balance methods, the algebraic equations obtained after the transformation of the differential equations of the dynamics are solved in the time domain rather than the frequency domain. Numerical examples are provided to demonstrate the performance of the method for autonomous and forced dynamics of a Van der Pol oscillator.     

Approximate stationary solution and stochastic stability for a class of differential equations with parametric colored noise

This paper aims to study a class of differential equations with parametric Gaussian colored noise. We present the general framework to get the solvability conditions of the approximate stationary probability density function, which is determined by the Fokker-Planck-Kolmogorov (FPK) equations. These equations are derived using the stochastic averaging method and the operator theory with the perturbation technique. An illustrative example is proposed to demonstrate the procedure of our proposed method. The analytical expression of approximate stationary probability density function is obtained. Numerical simulation is carried out to verify the analytical results and excellent agreement can be easily found. The FPK equation for the probability density function of order ε ⁰ is used to examine the almost-sure stability for the amplitude process. Finally, the stability in probability of the amplitude process is investigated by Lin and Cai’s method.     

Approximate Methods for the Linear and Nonlinear Analysis of Plates and Shallow Shells

Abstract

Berger's equations for the large amplitude deformation of membranes are used to produce a simple approximate expression for the large amplitude deflection of plates. The deformation of shallow shells is also considered and two approximate methods are outlined. Several important problems are discussed, the obtained solution being in good agreement with both experimental data and other approximate results. The main advantage of this technique is its ease of application, as it requires comparatively little computational work. A simple approximate formula for computing the fundamental frequency of a vibrating shallow shell is also presented and is shown to yield very accurate values in the case of a shallow dome and a rectangular panel.     

Stability of nonlinear periodic vibrations of 3D beams

The geometrically nonlinear periodic vibrations of beams with rectangular cross section under harmonic forces are investigated using a p-version finite element method. The beams vibrate in space; hence they experience longitudinal, torsional, and nonplanar bending deformations. The model is based on Timoshenko’s theory for bending and assumes that, under torsion, the cross section rotates as a rigid body and is free to warp in the longitudinal direction, as in Saint-Venant’s theory. The theory employed is valid for moderate rotations and displacements, and physical phenomena like internal resonances and change of the stability of the solutions can be investigated. Green’s nonlinear strain tensor and Hooke’s law are considered and isotropic and elastic beams are investigated. The equation of motion is derived by the principle of virtual work. The differential equations of motion are converted into a nonlinear algebraic form employing the harmonic balance method, and then solved by the arc-length continuation method. The variation of the amplitude of vibration in space with the excitation frequency of vibration is determined and presented in the form of response curves. The stability of the solution is investigated by Floquet’s theory.     

Complex variable function method for the scattering of plane waves by an arbitrary hole in a porous medium

Based on the complex variable function method, a new approach for solving the scattering of plane elastic waves by a hole with an arbitrary configuration embedded in an infinite poroelastic medium is developed in the paper. The poroelastic medium is described by Biot's theory. By introducing three potentials, the governing equations for Biot's theory are reduced to three Helmholtz equations for the three potentials. The series solutions of the Helmholtz equations are obtained by the wave function expansion method. Through the conformal mapping method, the arbitrary hole in the physical plane is mapped into a unit circle in the image plane. Integration of the boundary conditions along the unit circle in the image plane yields the algebraic equations for the coefficients of the series solutions. Numerical solution of the resulting algebraic equations yields the displacements, the stresses and the pore pressure for the porous medium. In order to demonstrate the proposed approach, some numerical results are given in the paper.