Search results

The purpose of this paper is to propose the generalized integral transform technique (GITT) to the solution of convection-diffusion problems with nonlinear boundary conditions by employing the corresponding nonlinear eigenvalue problem in the construction of the expansion basis.

The original nonlinear boundary condition coefficients in the problem formulation are all incorporated into the adopted eigenvalue problem, which may be itself integral transformed through a representative linear auxiliary problem, yielding a nonlinear algebraic eigenvalue problem for the associated eigenvalues and eigenvectors, to be solved along with the transformed ordinary differential system. The nonlinear eigenvalues computation may also be accomplished by rewriting the corresponding transcendental equation as an ordinary differential system for the eigenvalues, which is then simultaneously solved with the transformed potentials.

An application on one-dimensional transient diffusion with nonlinear boundary condition coefficients is selected for illustrating some important computational aspects and the convergence behavior of the proposed eigenfunction expansions. For comparison purposes, an alternative solution with a linear eigenvalue problem basis is also presented and implemented.

This novel approach can be further extended to various classes of nonlinear convection-diffusion problems, either already solved by the GITT with a linear coefficients basis, or new challenging applications with more involved nonlinearities.

This paper aims to develop an efficient numerical method for nonlinear transient heat conduction problems with local radiation boundary conditions and nonlinear heat sources.

Based on the physical characteristic of the transient heat conduction and the distribution characteristic of the Green’s function, a quasi-superposition principle is presented for the transient heat conduction problems with local nonlinearities. Then, an efficient method is developed, which indicates that the solution of the original nonlinear problem can be derived by solving some nonlinear problems with small structures and a linear problem with the original structure. These problems are independent of each other and can be solved simultaneously by the parallel computing technique.

Within a small time step, the nonlinear thermal loads can only induce significant temperature responses of the regions near the positions of the nonlinear thermal loads, whereas the temperature responses of the remaining regions are very close to zero. According to the above physical characteristic, the original nonlinear problem can be transformed into some nonlinear problems with small structures and a linear problem with the original structure.

An efficient and accurate numerical method is presented for transient heat conduction problems with local nonlinearities, and some numerical examples demonstrate the high efficiency and accuracy of the proposed method.

The purpose of this paper is to describe how a generalized single-system-single-solve (GS4-1) computational framework, previously developed for linear first order transient systems, can be properly extended for use in nonlinear counterparts, with particular applications to time dependent Burgers’ equation, which is well-known to serve as a simplified model of fluid dynamics, for illustrations of the essential concepts.

The framework permits, for a very general family of time integrators where traditional methods are a subset, much needed desirable features including second order time accuracy, robustness and unconditional stability, zero-order overshoot behavior, and additionally, a selective control of high frequency damping for both the primary variable and its time derivative. The latter, which is a new, key desirable feature not available in past/existing methods to-date, allows for different amounts of high frequency damping for both the primary variable and its time derivative to ensure physically accurate solutions of these variables. This is in contrast to having only limited control of these numerical dampings, often indiscriminately, as in some past developments which can lead to numerical instabilities in the time derivative variable. The extension of the framework to nonlinear problems, as described in this paper, is achieved via the use of a normalized time weighted residual approach, which naturally allows the time discretization of the transient nonlinear systems as being the natural extensions of the linear systems.

The primary aim is also to demonstrate the advantage of the selective control feature inherit in the present numerical methodologies for these nonlinear first order transient systems as in the linear counterparts.

The authors wish to tackle the challenges to further enable extensions to nonlinear first order transient systems that frequently arise in fluid dynamics problems; this is the focus of this paper. The primary wish is to demonstrate the ability of the GS4-1 framework for nonlinear first order transient systems as seen in the linear transient counterparts; while on one hand the authors show that an equal amount of high frequency damping (i.e. ρ _∞ = ρ _∞ s) leads to non-physical instability in the time derivative variable for a minimal damping required to obtain acceptable solution of the primary variable, on the other hand, the authors particularly demonstrate how this instability can be easily tuned off via the selective control feature (i.e. ρ _∞ ≠ ρ _∞ s) offered by the developed framework; thereby, demonstrating its robustness and superiority.

Gives a bibliographical review of the finite element methods (FEMs) applied for the linear and nonlinear, static and dynamic analyses of basic structural elements from the theoretical as well as practical points of view. The range of applications of FEMs in this area is wide and cannot be presented in a single paper; therefore aims to give the reader an encyclopaedic view on the subject. The bibliography at the end of the paper contains 2,025 references to papers, conference proceedings and theses/dissertations dealing with the analysis of beams, columns, rods, bars, cables, discs, blades, shafts, membranes, plates and shells that were published in 1992‐1995.

This paper aims to provide a well-behaved nonlinear scheme and accelerating iteration for the nonlinear convection diffusion equation with fundamental properties illustrated.

A nonlinear finite difference scheme is studied with fully implicit (FI) discretization used to acquire accurate simulation. A Picard–Newton (PN) iteration with a quadratic convergent ratio is designed to realize fast solution. Theoretical analysis is performed using the discrete function analysis technique. By adopting a novel induction hypothesis reasoning technique, the L^∞ (H¹) convergence of the scheme is proved despite the difficulty because of the combination of conservative diffusion and convection operator. Other properties are established consequently. Furthermore, the algorithm is extended from first-order temporal accuracy to second-order temporal accuracy.

Theoretical analysis shows that each of the two FI schemes is stable, its solution exists uniquely and has second-order spatial and first/second-order temporal accuracy. The corresponding PN iteration has the same order of accuracy and quadratic convergent speed. Numerical tests verify the conclusions and demonstrate the high accuracy and efficiency of the algorithms. Remarkable acceleration is gained.

The numerical method provides theoretical and technical support to accelerate resolving convection diffusion, non-equilibrium radiation diffusion and radiation transport problems.

The FI schemes and iterations for the convection diffusion problem are proposed with their properties rigorously analyzed. The induction hypothesis reasoning method here differs with those for linearization schemes and is applicable to other nonlinear problems.

To provide high torques needed to move a robot’s links, electric actuators are followed by a transmission system with a high transmission rate. For instance, gear ratios of 100:1 are often used in the joints of a robotic manipulator. This results into an actuator with large mechanical impedance (also known as nonback-drivable actuator). This in turn generates high contact forces when collision of the robotic mechanism occur and can cause humans’ injury. Another disadvantage of electric actuators is that they can exhibit overheating when constant torques have to be provided. Comparing to electric actuators, pneumatic actuators have promising properties for robotic applications, due to their low weight, simple mechanical design, low cost and good power-to-weight ratio. Electropneumatically actuated robots usually have better friction properties. Moreover, because of low mechanical impedance, pneumatic robots can provide moderate interaction forces which is important for robotic surgery and rehabilitation tasks. Pneumatic actuators are also well suited for exoskeleton robots. Actuation in exoskeletons should have a fast and accurate response. While electric motors come against high mechanical impedance and the risk of causing injuries, pneumatic actuators exhibit forces and torques which stay within moderate variation ranges. Besides, unlike direct current electric motors, pneumatic actuators have an improved weight-to-power ratio and avoid overheating problems.

The aim of this paper is to analyze a nonlinear optimal control method for electropneumatically actuated robots. A two-link robotic exoskeleton with electropneumatic actuators is considered as a case study. The associated nonlinear and multivariable state-space model is formulated and its differential flatness properties are proven. The dynamic model of the electropneumatic robot is linearized at each sampling instance with the use of first-order Taylor series expansion and through the computation of the associated Jacobian matrices. Within each sampling period, the time-varying linearization point is defined by the present value of the robot’s state vector and by the last sampled value of the control inputs vector. An H-infinity controller is designed for the linearized model of the robot aiming at solving the related optimal control problem under model uncertainties and external perturbations. An algebraic Riccati equation is solved at each time-step of the control method to obtain the stabilizing feedback gains of the H-infinity controller. Through Lyapunov stability analysis, it is proven that the robot’s control scheme satisfies the H-infinity tracking performance conditions which indicate the robustness properties of the control method. Moreover, global asymptotic stability is proven for the control loop. The method achieves fast convergence of the robot’s state variables to the associated reference trajectories, and despite strong nonlinearities in the robot’s dynamics, it keeps moderate the variations of the control inputs.

In this paper, a novel solution has been proposed for the nonlinear optimal control problem of robotic exoskeletons with electropneumatic actuators. As a case study, the dynamic model of a two-link lower-limb robotic exoskeleton with electropneumatic actuators has been considered. The dynamic model of this robotic system undergoes first approximate linearization at each iteration of the control algorithm around a temporary operating point. Within each sampling period, this linearization point is defined by the present value of the robot’s state vector and by the last sampled value of the control inputs vector. The linearization process relies on first-order Taylor series expansion and on the computation of the associated Jacobian matrices. The modeling error which is due to the truncation of higher-order terms from the Taylor series is considered to be a perturbation which is asymptotically compensated by the robustness of the control algorithm. To stabilize the dynamics of the electropneumatically actuated robot and to achieve precise tracking of reference setpoints, an H-infinity (optimal) feedback controller is designed. Actually, the proposed H-infinity controller for the model of the two-link electropneumatically actuated exoskeleton achieves the solution of the associated optimal control problem under model uncertainty and external disturbances. This controller implements a min-max differential game taking place between: (i) the control inputs which try to minimize a cost function which comprises a quadratic term of the state vector’s tracking error and (ii) the model uncertainty and perturbation inputs which try to maximize this cost function. To select the stabilizing feedback gains of this H-infinity controller, an algebraic Riccati equation is being repetitively solved at each time-step of the control method. The global stability properties of the H-infinity control scheme are proven through Lyapunov analysis.

Pneumatic actuators are characterized by high nonlinearities which are due to air compressibility, thermodynamics and valves behavior and thus pneumatic robots require elaborated nonlinear control schemes to ensure their fast and precise positioning. Among the control methods which have been applied to pneumatic robots, one can distinguish differential geometric approaches (Lie algebra-based control, differential flatness theory-based control, nonlinear model predictive control [NMPC], sliding-mode control, backstepping control and multiple models-based fuzzy control). Treating nonlinearities and fault tolerance issues in the control problem of robotic manipulators with electropneumatic actuators has been a nontrivial task.

The novelty of the proposed control method is outlined as follows: preceding results on the use of H-infinity control to nonlinear dynamical systems were limited to the case of affine-in-the-input systems with drift-only dynamics. These results considered that the control inputs gain matrix is not dependent on the values of the system’s state vector. Moreover, in these approaches the linearization was performed around points of the desirable trajectory, whereas in the present paper’s control method the linearization points are related with the value of the state vector at each sampling instance as well as with the last sampled value of the control inputs vector. The Riccati equation which has been proposed for computing the feedback gains of the controller is novel, so is the presented global stability proof through Lyapunov analysis. This paper’s scientific contribution is summarized as follows: (i) the presented nonlinear optimal control method has improved or equally satisfactory performance when compared against other nonlinear control schemes that one can consider for the dynamic model of robots with electropneumatic actuators (such as Lie algebra-based control, differential flatness theory-based control, nonlinear model-based predictive control, sliding-mode control and backstepping control), (ii) it achieves fast and accurate tracking of all reference setpoints, (iii) despite strong nonlinearities in the dynamic model of the robot, it keeps moderate the variations of the control inputs and (iv) unlike the aforementioned alternative control approaches, this paper’s method is the only one that achieves solution of the optimal control problem for electropneumatic robots.

The use of electropneumatic actuation in robots exhibits certain advantages. These can be the improved weight-to-power ratio, the lower mechanical impedance and the avoidance of overheating. At the same time, precise positioning and accurate execution of tasks by electropneumatic robots requires the application of elaborated nonlinear control methods. In this paper, a new nonlinear optimal control method has been developed for electropneumatically actuated robots and has been specifically applied to the dynamic model of a two-link robotic exoskeleton. The benefit from using this paper’s results in industrial and biomedical applications is apparent.

A comparison of the proposed nonlinear optimal (H-infinity) control method against other linear and nonlinear control schemes for electropneumatically actuated robots shows the following: (1) Unlike global linearization-based control approaches, such as Lie algebra-based control and differential flatness theory-based control, the optimal control approach does not rely on complicated transformations (diffeomorphisms) of the system’s state variables. Besides, the computed control inputs are applied directly on the initial nonlinear model of the electropneumatic robot and not on its linearized equivalent. The inverse transformations which are met in global linearization-based control are avoided and consequently one does not come against the related singularity problems. (2) Unlike model predictive control (MPC) and NMPC, the proposed control method is of proven global stability. It is known that MPC is a linear control approach that if applied to the nonlinear dynamics of the electropneumatic robot, the stability of the control loop will be lost. Besides, in NMPC the convergence of its iterative search for an optimum depends on initialization and parameter values selection and consequently the global stability of this control method cannot be always assured. (3) Unlike sliding-mode control and backstepping control, the proposed optimal control method does not require the state-space description of the system to be found in a specific form. About sliding-mode control, it is known that when the controlled system is not found in the input-output linearized form the definition of the sliding surface can be an intuitive procedure. About backstepping control, it is known that it cannot be directly applied to a dynamical system if the related state-space model is not found in the triangular (backstepping integral) form. (4) Unlike PID control, the proposed nonlinear optimal control method is of proven global stability, the selection of the controller’s parameters does not rely on a heuristic tuning procedure, and the stability of the control loop is assured in the case of changes of operating points. (5) Unlike multiple local models-based control, the nonlinear optimal control method uses only one linearization point and needs the solution of only one Riccati equation so as to compute the stabilizing feedback gains of the controller. Consequently, in terms of computation load the proposed control method for the electropneumatic actuator’s dynamics is much more efficient.

Permanent magnet synchronous spherical motors can have wide use in robotics and industrial automation. They enable three-DOF omnidirectional motion of their rotor. They are suitable for several applications, such as actuation in robotics, traction in electric vehicles and use in several automation systems. Unlike conventional synchronous motors, permanent magnet synchronous spherical motors consist of a fixed inner shell, which is the stator, and a rotating outer shell, which is the rotor. Their dynamic model is multivariable and strongly nonlinear. The treatment of the associated control problem is important.

In this paper, the multivariable dynamic model of permanent magnet synchronous spherical motors is analysed, and a nonlinear optimal (H-infinity) control method is developed for it. Differential flatness properties are proven for the spherical motors’ state-space model. Next, the motors’ state-space description undergoes approximate linearization with the use of first-order Taylor series expansion and through the computation of the associated Jacobian matrices. The linearization process takes place at each sampling instance around a time-varying operating point, which is defined by the present value of the motors’ state vector and by the last sampled value of the control input vector. For the approximately linearized model of the permanent magnet synchronous spherical motors, a stabilizing H-infinity feedback controller is designed. To compute the controller’s gains, an algebraic Riccati equation has to be repetitively solved at each time-step of the control algorithm. The global stability properties of the control scheme are proven through Lyapunov analysis. Finally, the performance of the nonlinear optimal control method is compared against a flatness-based control approach implemented in successive loops.

Due to the nonlinear and multivariable structure of the state-space model of spherical motors, the solution of the associated nonlinear control problem is a nontrivial task. In this paper, a novel nonlinear optimal (H-infinity) control approach is proposed for the dynamic model of permanent magnet synchronous spherical motors. The method is based on approximate linearization of the motor’s state-space model with the use of first-order Taylor series expansion and the computation of the associated Jacobian matrices. Furthermore, the paper has introduced a different solution to the nonlinear control problem of the permanent magnet synchronous spherical motor, which is based on flatness-based control implemented in successive loops.

The presented control approaches do not exhibit any limitations, but on the contrary, they have specific advantages. In comparison to global linearization-based control schemes (such as Lie-algebra-based control), they do not make use of complicated changes of state variables (diffeomorphisms) and transformations of the system's state-space description. The computed control inputs are applied directly to the initial nonlinear state-space model of the permanent magnet spherical motor without the intervention of inverse transformations and thus without coming against the risk of singularities.

The motion control problem of spherical motors is nontrivial because of the complicated nonlinear and multivariable dynamics of these electric machines. So far, there have been several attempts to apply nonlinear feedback control to permanent magnet-synchronous spherical motors. However, due to the model’s complexity, few results exist about the associated nonlinear optimal control problem. The proposed nonlinear control methods for permanent magnet synchronous spherical motors make more efficient, precise and reliable the use of such motors in robotics, electric traction and several automation systems.

The treated research topic is central for robotic and industrial automation. Permanent magnet synchronous spherical motors are suitable for several applications, such as actuation in robotics, traction in electric vehicles and use in several automation systems. The solution of the control problem for the nonlinear dynamic model of permanent magnet synchronous spherical motors has many industrial applications and therefore contributes to economic growth and development.

The proposed nonlinear optimal control method is novel compared to past attempts to solve the optimal control problem for nonlinear dynamical systems. Unlike past approaches, in the new nonlinear optimal control method, linearization is performed around a temporary operating point, which is defined by the present value of the system's state vector and by the last sampled value of the control inputs vector and not at points that belong to the desirable trajectory (setpoints). Besides, the Riccati equation which is used for computing the feedback gains of the controller is new, and so is the global stability proof for this control method. Compared to nonlinear model predictive control, which is a popular approach for treating the optimal control problem in industry, the new nonlinear optimal (H-infinity) control scheme is of proven global stability, and the convergence of its iterative search for the optimum does not depend on initial conditions and trials with multiple sets of controller parameters. It is also noteworthy that the nonlinear optimal control method is applicable to a wider class of dynamical systems than approaches based on the solution of state dependent Riccati equations (SDRE). The SDRE approaches can be applied only to dynamical systems which can be transformed into the linear parameter varying form. Besides, the nonlinear optimal control method performs better than nonlinear optimal control schemes, which use approximation of the solution of the Hamilton–Jacobi–Bellman equation by Galerkin series expansions. Furthermore, the second control method proposed in this paper, which is flatness-based control in successive loops, is also novel and demonstrates substantial contribution to nonlinear control for robotics and industrial automation.

The numerical solution of electromagnetic field nonlinear problems requires successive building and solving of linear systems of equations. This is the most time–consuming part, especially for large problems. Both fast linear solvers and efficient nonlinear iterative algorithms, are critical for the overall efficiency of the nonlinear electromagnetic field solver. This paper presents an analysis of a variety of techniques that can be efficiently used to reduce the solution time of nonlinear magnetic field equations in large finite element method (FEM) problems.

The purpose of this paper is to present a new method based on the homotopy analysis method (HAM) with the aim of fast searching and calculating multiple solutions of nonlinear boundary value problems (NBVPs).

A major problem with the previously modified HAM, namely, predictor homotopy analysis method, which is used to predict multiplicity of solutions of NBVPs, is a time-consuming computation of high-order HAM-approximate solutions due to a symbolic variable namely “prescribed parameter”. The proposed new technique which is based on traditional shooting method, and the HAM cuts the dependency on the prescribed parameter.

To demonstrate the computational efficiency, the mentioned method is implemented on three important nonlinear exactly solvable differential equations, namely, the nonlinear MHD Jeffery–Hamel flow problem, the nonlinear boundary value problem arising in heat transfer and the strongly nonlinear Bratu problem.

The more high-order approximate solutions are computable, multiple solutions are easily searched and discovered and the more accurate solutions can be obtained depending on how nonhomogeneous boundary conditions are transcribed to the homogeneous boundary conditions.