Logarithmic matrix norms in motion stability problems

The problem of the stability of the motions of mechanical systems, described by non-linear non-autonomous systems of ordinary differential equations, is considered. Using the logarithmic matrix norm method, and constructing a reference system, the sufficient conditions for the asymptotic and exponential stability of unperturbed motion and for the stabilization of progammed motions of such systems are obtained. The problem of the asymptotic stability of a non-conservative system with two degrees of freedom is solved, taking for parametric disturbances into account. Examples of the solution of the problem of stabilizing programmed motions – for an inverted double pendulum and for a two-link manipulator on a stationary base – are considered.     

Linear systems with a quadratic integral

It is shown that a linear system of n differential equations with constant coefficients, at least one of whose integrals is a non-degenerate quadratic form, may be reduced to a canonical system of Hamiltonian equations. In particular, n is even and the phase flow preserves the standard measure; if the index of the quadratic integral is odd, the trivial solution is unstable, and so on. For the case n = 4 the stability conditions are given a geometrical form. The general results are used to investigate small oscillations of non-holonomic systems, and also the problem of the stability of invariant manifolds of non-linear systems that have Morse functions as integrals.     

The onset of resonance processes in the Couette–Taylor problem close to the point of codimension-2 bifurcation

The bifurcations on passing around the point of intersection of two neutral curves (points of codimension-2 bifurcation) are considered in the Couette–Taylor problem of the fluid motion between rotating cylinders. The secondary modes in a small neighbourhood of a point of codimension-2 bifurcation are studied using a system of non-linear amplitude equations in a central manifold. The steady-state solutions of the amplitude systems, to which secondary periodic modes of the travelling-wave type, non-linear mixtures of travelling waves and unsteady two-, three- and four-frequency quasiperiodic solutions of the system of Navier–Stokes equations correspond, are analysed. A numerical analysis of the conditions for the existence and stability of irrotationally symmetric steady-state fluid flows between unidirectionally rotating cylinders is carried out.     

Periodic and almost periodic solutions of strongly non-linear impulsive systems

The problem of the existence and stability of periodic and almost periodic solutions of strongly non-linear impulsive systems is investigated. The Poincaré method [1] is justified for the case of an isolated generating solution. A dynamical system consisting of a bead on a vibrating surface is considered as an example.

The small parameter method for investigating systems with discontinuous solutions was previously applied [2, 3] to the case when the periodic solution is non-isolated.

A method is used below for reducing the investigation of a system of equations with impulsive actions on surfaces to equations with fixed moments of inpulsive action.     

Problems of the partial stability and detectability of dynamical systems

The conditions under which uniform stability (uniform asymptotic stability) with respect to a part of the variables of the zero equilibrium position of a non-linear non-stationary system of ordinary differential equations signifies uniform stability (uniform asymptotic stability) of this equilibrium position with respect the other, larger part of the variables, which include an additional group of coordinates of the phase vector, are established. These conditions include the condition for uniform asymptotic stability of the zero equilibrium position of the “reduced” subsystem of the original system with respect to the additional group of variables. Since within the conditions obtained the stability with respect to the remaining unmeasured coordinates of the phase vector remains undetermined or is investigated additionally, partial zero-detectability of the original system occurs in this case, and the conditions obtained supplement the series of known results from partial stability theory. The application of the results obtained to problems of the partial stabilization of non-linear controlled systems, particularly to the problem of stabilizing an asymmetric rigid body relative to an assigned direction in an inertial space, is considered. The partial detectability of linear systems with constant coefficients is also investigated.     

The stability and stabilization of the motion of non-conservative mechanical systems

Two stability problems are solved. In the first, the stability of mechanical systems, on which dissipative, gyroscopic, potential and positional non-conservative forces (systems of general form) act, is investigated. The stability is considered in the case when the potential energy has a maximum at equilibrium. The condition for asymptotic stability is obtained by constructing Lyapunov's function. In the second problem, the possibility of stabilizing a gyroscopic system with two degrees of freedom up to asymptotic stability using non-linear dissipative and positional non-conservative forces is investigated. Stability of the gyroscopic system is achieved by gyroscopic stabilization. The stability conditions are obtained in terms of the system parameters. Cases when the gyroscopic stabilization is disrupted by these non-linear forces are indicated.     

Numerical solutions for convection–diffusion equation using El-Gendi method

The method of El-Gendi [El-Gendi SE. Chebyshev solution of differential integral and integro-differential equations. J Comput 1969;12:282–7; Mihaila B, Mihaila I. Numerical approximation using Chebyshev polynomial expansions: El-gendi’s method revisited. J Phys A Math Gen 2002;35:731–46] is presented with interface points to deal with linear and non-linear convection–diffusion equations.The linear problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using three-level time scheme.The non-linear problem is reduced to three systems of ordinary differential. Each one of these systems is, then, solved using three-level time scheme. Numerical results for Burgers’ equation and modified Burgers’ equation are shown and compared with other methods. The numerical results are found to be in good agreement with the exact solutions.     

Spectral methods for some singularly perturbed third order ordinary differential equations

Spectral methods with interface point are presented to deal with some singularly perturbed third order boundary value problems of reaction-diffusion and convection-diffusion types. First, linear equations are considered and then non-linear equations. To solve non-linear equations, Newton’s method of quasi-linearization is applied. The problem is reduced to two systems of ordinary differential equations. And, then, each system is solved using spectral collocation methods. Our numerical experiments show that the proposed methods are produce highly accurate solutions in little computer time when compared with the other methods available in the literature.      

1:2 and 1:3 internal resonance active absorber for non-linear vibrating system

Vibration and dynamic chaos should be controlled in either structures or machines. An active vibration absorber for suppressing the vibration of the non-linear plant when subjected to external and parametric excitations is studied in the presence of one-to-two and one-to-three internal resonance. The main attention is focused on the study of the active control and stability of two systems, which can be used to reduce vibrations due to rotor blade flapping motion. The method of multiple scale perturbation technique is applied to determine four first-order non-linear ordinary differential equations that govern the modulation of the amplitudes and phases in the presence of internal resonance of the two systems with quadratic and cubic order of control. These equations are used to determine the steady state solutions and their stability. The stability study of non-linear periodic solution for two cases (1:2 and 1:3 internal resonance) and the stability of the obtained numerical solution are investigated using frequency, force-response curves and phase-plane method. Also, effects of some parameters on the steady state solution of the vibrating system are investigated and reported in this paper. Variation of some parameters leads to the bending of the frequency, force-response curves and hence to the jump phenomenon occurrence. The reported results are compared to the available published work.     

Hyperbolic system of balance laws modeling dissipative ionic crystals

Based on the classical continuum theory of electroelasticity which includes polarization gradients as independent variables, we propose a constitutive model for ionic crystals accounting for both ionic and electronic contributions to polarization. Dissipation is modeled via internal variables which satisfy suitable evolution equations and the consequences of the second law of thermodynamics are exploited to cast the non-linear problem in the form of a symmetric hyperbolic system of balance laws. The stability of perturbations with respect to unstrained, unpolarized states is discussed. A set of linear equations is also derived for the fully electromagnetic problem which generalizes previous results.     

The non-linear Saint-Venant problem of the torsion,stretching and bending of a naturally twisted rod

Problems on large stretching, torsional and bending deformations of a naturally twisted rod, loaded with end forces and moments, are considered from the point of view of the non-linear three-dimensional theory of elasticity. Particular solutions of the equations of elastostatics are found, which are two-parameter families of finite deformations and which possess the property that, for these deformations, the initial system of three-dimensional non-linear equations reduces to a system of equations with two independent variables. The use of these equations enables one to reduce certain Saint-Venant problems for a naturally twisted rod to two-dimensional non-linear boundary-value problems for a planar domain in the form of the cross-section of a rod. Different formulations of the two-dimensional boundary-value problem for the cross-section are proposed, which differ in the choice of the unknown functions. A non-linear problem of the torsion and stretching of a circular cylinder with helical anisotropy, which is reduced to ordinary differential equations, is considered as a special case.     

Detection of maximum loadability limits and weak buses using Chaotic PSO considering security constraints

In the current research chaotic search is used with the optimization technique for solving non-linear complicated power system problems because Chaos can overcome the local optima problem of optimization technique. Power system problem, more specifically voltage stability, is one of the practical examples of non-linear, complex, convex problems. Smart grid, restructured energy system and socio-economic development fetch various uncertain events in power systems and the level of uncertainty increases to a great extent day by day. In this context, analysis of voltage stability is essential. The efficient method to assess the voltage stability is maximum loadability limit (MLL). MLL problem is formulated as a maximization problem considering practical security constraints (SCs). Detection of weak buses is also important for the analysis of power system stability. Both MLL and weak buses are identified by PSO methods and FACTS devices can be applied to the detected weak buses for the improvement of stability. Three particle swarm optimization (PSO) techniques namely General PSO (GPSO), Adaptive PSO (APSO) and Chaotic PSO (CPSO) are presented for the comparative study with obtaining MLL and weak buses under different SCs. In APSO method, PSO-parameters are made adaptive with the problem and chaos is incorporated in CPSO method to obtain reliable convergence and better performances. All three methods are applied on standard IEEE 14 bus, 30 bus, 57 bus and 118 bus test systems to show their comparative computing effectiveness and optimization efficiencies.     

Unique solvability of a non-linear non-local boundary-value problem for systems of non-linear functional differential equations

General conditions for the unique solvability of a non-linear nonlocal boundary-value problem for systems of non-linear functional differential equations are obtained.     

Criteria for the functional controllability and invertibility of non-linear systems

In the development of investigations on inverse problems [1, 2], criteria for the functional controllability and invertibility of non-linear systems of equations with an output are obtained. The solution is based on the construction of an inverse system for which the input action of the initial system is the output. An identification problem is considered which corresponds to the problem of invertibility with an unknown initial state. The properties of λ-invertibility and λ-identifiability, which arise in cases when the output signal is known in a set of trajectories, are investigated.     

Symbolic computation of limit cycles associated with Hilbert’s 16th problem

This paper is concerned with the practical complexity of the symbolic computation of limit cycles associated with Hilbert’s 16th problem. In particular, in determining the number of small-amplitude limit cycles of a non-linear dynamical system, one often faces computing the focus values of Hopf-type critical points and solving lengthy coupled polynomial equations. These computations must be carried out through symbolic computation with the aid of a computer algebra system such as Maple or Mathematica, and thus usually gives rise to very large algebraic expressions. In this paper, efficient computations for the focus values and polynomial equations are discussed, showing how to deal with the complexity in the computation of non-linear dynamical systems.     

Solución de la cinemática directa de la plataforma Gough-Stewart general usando un algoritmo híbrido de optimización

The direct kinematics problem for parallel robots can be stated as follows: given values of the joint variables, the corresponding Cartesian variable values, the pose of the end-effector, must be found. Most of the times the direct kinematics problem involves the solution of a system of non-linear equations. The most efficient methods to solve such kind of equations assume convexity in a cost function which minimum is the solution of the non-linear system. In consequence, the capacity of such methods depends on the knowledge about an starting point which neighboring region is convex, hence the method can find the global minimum. This article propose a method based on probabilistic learning about an adequate starting point for the Dogleg method which assumes local convexity of the function. The proposed method efficiently avoids the local minima, without need of human intervention or apriori knowledge, thus it shows a more robust performance than the simple Dogleg method or other gradient based methods. To demonstrate the performance of the proposed hybrid method, numerical experiments and the respective discussion are presented. The proposal can be extended to other structures of closed-kinematics chains, to the general solution of systems of non-linear equations, and to the minimization of non-linear functions.     

Non-linear deformation properties of materials with stochastically distributed anisotropic inclusions

In the present contribution, the problem of non-linear deformation of materials with stochastically distributed anisotropic inclusions is considered on the basis of the methods of mechanics of stochastically non-homogeneous media. The homogenization model of materials of stochastic structure with physically non-linear components is developed for the case of a matrix which is strengthened by unidirectional ellipsoidal inclusions. It is assumed that the matrix is isotropic, deforms non-linearly; inclusions are linear-elastic and have transversally-isotropic symmetry of physical and mechanical properties. Stochastic differential equations of physically non-linear elasticity theory form the underlying equations. Transformation of these equations into integral equations by using the Green's function and application of the method of conditional moments allow us to reduce the problem to a system of non-linear algebraic equations. This system of non-linear algebraic equations is solved by the Newton-Raphson method. On the analytical as well as the numerical basis, the algorithm for determination of the non-linear effective characteristics of such a material is introduced. The non-linear behavior of such a material is caused by the non-linear matrix deformations. On the basis of the numerical solution, the dependences of homogenized Poisson's coefficients on macro-strains and the non-linear stress-strain diagrams for a material with randomly distributed unidirectional ellipsoidal pores are predicted and discussed for different volume fractions of pores. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Tropical differential equations

Tropical differential equations are introduced and an algorithm is designed which tests solvability of a system of tropical linear differential equations within the complexity polynomial in the size of the system and in the absolute values of its coefficients. Moreover, we show that there exists a minimal solution, and the algorithm constructs it (in case of solvability). This extends a similar complexity bound established for tropical linear systems. In case of tropical linear differential systems in one variable a polynomial complexity algorithm for testing its solvability is designed.We prove also that the problem of solvability of a system of tropical non-linear differential equations in one variable is NP-hard, and this problem for arbitrary number of variables belongs to NP. Similar to tropical algebraic equations, a tropical differential equation expresses the (necessary) condition on the dominant term in the issue of solvability of a differential equation in power series.     

On solvability of non-linear semi-periodic boundary-value problem for system of hyperbolic equations

We study a non-linear semi-periodic boundary-value problem for a system of hyperbolic equations with mixed derivative. At that, the semi-periodic boundary-value problem for a system of hyperbolic equations is reduced to an equivalent problem, consisting of a family of periodic boundary-value problems for ordinary differential equations and functional relation. When solving a family of periodic boundary-value problems of ordinary differential equations we use the method of parameterization. This approach allowed to establish sufficient conditions for the existence of an isolated solution of non-linear semi-periodic boundary-value problem for a system of hyperbolic equations.     

Application to differential transformation method for solving systems of differential equations

In this paper, we present an analytical solution for different systems of differential equations by using the differential transformation method. The convergence of this method has been discussed with some examples which are presented to show the ability of the method for linear and non-linear systems of differential equations. We begin by showing how the differential transformation method applies to a non-linear system of differential equations and give two examples to illustrate the sufficiency of the method for linear and non-linear stiff systems of differential equations. The results obtained are in good agreement with the exact solution and Runge–Kutta method. These results show that the technique introduced here is accurate and easy to apply.