The stability analysis for nonlinear differentialalgebraic systems is addressed using tools from classical control theory. Sufficient stability conditions relying on matrix inequalities are established via Lyapunov Di...The stability analysis for nonlinear differentialalgebraic systems is addressed using tools from classical control theory. Sufficient stability conditions relying on matrix inequalities are established via Lyapunov Direct Method. In addition, a novel interpretation of differential-algebraic systems as feedback interconnection of a purely differential system and an algebraic system allows reducing the stability analysis to a smallgain-like condition. The study of stability properties for constrained mechanical systems, for a class of Lipschitz differential-algebraic systems and for an academic example is used to illustrate the theory.展开更多
Numerical methods for Differential-Algebraic systems with discontinuous right-hand sides is discussed. A class of continuous Rosenbrock methods are constructed, and numerical experiments show that the continuous Rosen...Numerical methods for Differential-Algebraic systems with discontinuous right-hand sides is discussed. A class of continuous Rosenbrock methods are constructed, and numerical experiments show that the continuous Rosenbrock methods are effective. Applying the methods, a fast and high-precision numerical algorithm is given to deal with typical discontinuous parts, which occur frequently in differential-algebraic systems(DAS).展开更多
In this paper, the stable problem for differential-algebraic systems is investigated by a convex op-timization approach. Based on the Lyapunov functional method and the delay partitioning approach, some delay and its ...In this paper, the stable problem for differential-algebraic systems is investigated by a convex op-timization approach. Based on the Lyapunov functional method and the delay partitioning approach, some delay and its time-derivative dependent stable criteria are obtained and formulated in the form of simple linear matrix inequalities (LMIs). The obtained criteria are dependent on the sizes of delay and its time-derivative and are less conservative than those produced by previous approaches.展开更多
A class of parallel Runge-Kutta Methods for differential-algebraic equations of index 2are constructed for multiprocessor system. This paper gives the order conditions and investigatesthe convergence theory for such m...A class of parallel Runge-Kutta Methods for differential-algebraic equations of index 2are constructed for multiprocessor system. This paper gives the order conditions and investigatesthe convergence theory for such methods.展开更多
Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some eq...Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some equivalent system,,; are given. Using them the paper shows how to state properly initial and boundary conditions for these DAEs. The existence and uniqueness theory of the solution of the initial and boundary value problems for higher index DAEs are proposed.展开更多
An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction...An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations (DAEs) where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula (BDF) using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.展开更多
In this article, linear regular index 2 DAEs A(t)[D(t)x(t)]' + B(t)x(t) = q(t) are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular inde...In this article, linear regular index 2 DAEs A(t)[D(t)x(t)]' + B(t)x(t) = q(t) are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.展开更多
Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As...Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As far as we known, the nonnegativity of solutions of the nonlinear fractional differential-algebraic equations is still not studied. In this article, we investigate the nonnegativity of solutions of the equations. Firstly, we discuss the existence of nonnegative solutions of the equations, and then we show that the nonnegative solution can be approached by a monotone waveform relaxation sequence provided the initial iteration is chosen properly. The choice of initial iteration is critical and we give a method of finding it. Finally, we present an example to illustrate the efficiency of our method.展开更多
The objective of dynamical system learning tasks is to forecast the future behavior of a system by leveraging observed data.However,such systems can sometimes exhibit rigidity due to significant variations in componen...The objective of dynamical system learning tasks is to forecast the future behavior of a system by leveraging observed data.However,such systems can sometimes exhibit rigidity due to significant variations in component parameters or the presence of slow and fast variables,leading to challenges in learning.To overcome this limitation,we propose a multiscale differential-algebraic neural network(MDANN)method that utilizes Lagrangian mechanics and incorporates multiscale information for dynamical system learning.The MDANN method consists of two main components:the Lagrangian mechanics module and the multiscale module.The Lagrangian mechanics module embeds the system in Cartesian coordinates,adopts a differential-algebraic equation format,and uses Lagrange multipliers to impose constraints explicitly,simplifying the learning problem.The multiscale module converts high-frequency components into low-frequency components using radial scaling to learn subprocesses with large differences in velocity.Experimental results demonstrate that the proposed MDANN method effectively improves the learning of dynamical systems under rigid conditions.展开更多
Two cubical 3D electric circuits with single and double capacitors and twelve ohmic resistors are considered. The resistors are the sides of the cube. The circuit is fed with a single internal emf. The charge on the c...Two cubical 3D electric circuits with single and double capacitors and twelve ohmic resistors are considered. The resistors are the sides of the cube. The circuit is fed with a single internal emf. The charge on the capacitor(s) and the current distributions of all twelve sides of the circuit(s) vs. time are evaluated. The analysis requires solving twelve differential-algebraic intertwined symbolic equations. This is accomplished by applying a Computer Algebra System (CAS), specifically Mathematica. The needed codes are included. For a set of values assigned to the elements, the numeric results are depicted.展开更多
In this article,we investigate a fractional-order singular Leslie-Gower prey-predator bioeconomic model,which describes the interaction bet ween populations of prey and predator,and takes into account the economic int...In this article,we investigate a fractional-order singular Leslie-Gower prey-predator bioeconomic model,which describes the interaction bet ween populations of prey and predator,and takes into account the economic interest.We firstly obtain the solvability condition and the st ability of the model sys tem,and discuss the singularity induced bifurcation phenomenon.Next,we introduce a st ate feedback controller to elimina te the singularity induced bifurcation phenomenon,and discuss the optimal control problems.Finally,numerical solutions and their simulations are considered in order to illustrate the theoretical results and reveal the more complex dynamical behavior.展开更多
In this paper, Lyapunov-like exponential stability and unstability of differentialalgebraic equation are considered from the viewpoint of stability of system motion, and the criteria of exponential stability and unsta...In this paper, Lyapunov-like exponential stability and unstability of differentialalgebraic equation are considered from the viewpoint of stability of system motion, and the criteria of exponential stability and unstability of nonlinear nonautonomous differential-algebraic equation are given by using Lyapunov-like function similar to ordinary differential equation.展开更多
A nonlinear mathematical model for the analysis of large deformation of frame structures with discontinuity conditions and initial displacements, subject to dynamic loads is formulated with arc-coordinates. The differ...A nonlinear mathematical model for the analysis of large deformation of frame structures with discontinuity conditions and initial displacements, subject to dynamic loads is formulated with arc-coordinates. The differential quadrature element method (DQEM) is then applied to discretize the nonlinear mathematical model in the spatial domain, An effective method is presented to deal with discontinuity conditions of multivariables in the application of DQEM. A set of DQEM discretization equations are obtained, which are a set of nonlinear differential-algebraic equations with singularity in the time domain. This paper also presents a method to solve nonlinear differential-algebra equations. As application, static and dynamical analyses of large deformation of frames and combined frame structures, subjected to concentrated and distributed forces, are presented. The obtained results are compared with those in the literatures. Numerical results show that the proposed method is general, and effective in dealing with disconti- nuity conditions of multi-variables and solving differential-algebraic equations. It requires only a small number of nodes and has low computation complexity with high precision and a good convergence property.展开更多
Projected Runge-Kutta (R-K) methods for constrained Hamiltonian systems are proposed. Dynamic equations of the systems, which are index-3 differential-algebraic equations (DAEs) in the Heisenberg form, are establi...Projected Runge-Kutta (R-K) methods for constrained Hamiltonian systems are proposed. Dynamic equations of the systems, which are index-3 differential-algebraic equations (DAEs) in the Heisenberg form, are established under the framework of Lagrangian multipliers. R-K methods combined with the technique of projections are then used to solve the DAEs. The basic idea of projections is to eliminate the constraint violations at the position, velocity, and acceleration levels, and to preserve the total energy of constrained Hamiltonian systems by correcting variables of the position, velocity, acceleration, and energy. Numerical results confirm the validity and show the high precision of the proposed method in preserving three levels of constraints and total energy compared with results reported in the literature.展开更多
Linear multistep methods and one-leg methods are applied to a class of index-2 nonlinear differential-algebraic equations with a variable delay.The corresponding convergence results are obtained and successfully confi...Linear multistep methods and one-leg methods are applied to a class of index-2 nonlinear differential-algebraic equations with a variable delay.The corresponding convergence results are obtained and successfully confirmed by some numerical examples.The results obtained in this work extend the corresponding ones in literature.展开更多
This paper develops a class of general one-step discretization methods for solving the index-1 stochastic delay differential-algebraic equations. The existence and uniqueness theorem of strong solutions of index-1 equ...This paper develops a class of general one-step discretization methods for solving the index-1 stochastic delay differential-algebraic equations. The existence and uniqueness theorem of strong solutions of index-1 equations is given. A strong convergence criterion of the methods is derived, which is applicable to a series of one-step stochastic numerical methods. Some specific numerical methods, such as the Euler-Maruyama method, stochastic ^-methods, split-step ^-methods are proposed, and their strong convergence results are given. Numerical experiments further illustrate the theoretical results.展开更多
In this paper,we propose a novel Legendre neural network combined with the extreme learning machine algorithm to solve variable coefficients linear delay differential-algebraic equations with weak discontinuities.Firs...In this paper,we propose a novel Legendre neural network combined with the extreme learning machine algorithm to solve variable coefficients linear delay differential-algebraic equations with weak discontinuities.First,the solution interval is divided into multiple subintervals by weak discontinuity points.Then,Legendre neural network is used to eliminate the hidden layer by expanding the input pattern using Legendre polynomials on each subinterval.Finally,the parameters of the neural network are obtained by training with the extreme learning machine.The numerical examples show that the proposed method can effectively deal with the difficulty of numerical simulation caused by the discontinuities.展开更多
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 ...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.展开更多
Based on constructing programmed constraint and constraint perturbation equation, a kinematics and dynamics numerical simulation model is established for virtual mechanism, in which the difference scheme guarantee pre...Based on constructing programmed constraint and constraint perturbation equation, a kinematics and dynamics numerical simulation model is established for virtual mechanism, in which the difference scheme guarantee precision in simulation procedure and its mtmerical solutions satisfy programmed manifold stability. A crank-piston mechanism in a car engine, a steering mechanism and a suspension mechanism are simulated in a virtual environment, then comparing the simulation results with those obtained in ADAMS under the same circumstances proved the solver valid.展开更多
文摘The stability analysis for nonlinear differentialalgebraic systems is addressed using tools from classical control theory. Sufficient stability conditions relying on matrix inequalities are established via Lyapunov Direct Method. In addition, a novel interpretation of differential-algebraic systems as feedback interconnection of a purely differential system and an algebraic system allows reducing the stability analysis to a smallgain-like condition. The study of stability properties for constrained mechanical systems, for a class of Lipschitz differential-algebraic systems and for an academic example is used to illustrate the theory.
文摘Numerical methods for Differential-Algebraic systems with discontinuous right-hand sides is discussed. A class of continuous Rosenbrock methods are constructed, and numerical experiments show that the continuous Rosenbrock methods are effective. Applying the methods, a fast and high-precision numerical algorithm is given to deal with typical discontinuous parts, which occur frequently in differential-algebraic systems(DAS).
文摘In this paper, the stable problem for differential-algebraic systems is investigated by a convex op-timization approach. Based on the Lyapunov functional method and the delay partitioning approach, some delay and its time-derivative dependent stable criteria are obtained and formulated in the form of simple linear matrix inequalities (LMIs). The obtained criteria are dependent on the sizes of delay and its time-derivative and are less conservative than those produced by previous approaches.
文摘A class of parallel Runge-Kutta Methods for differential-algebraic equations of index 2are constructed for multiprocessor system. This paper gives the order conditions and investigatesthe convergence theory for such methods.
基金Project supported by the National Natural Science Foundation of China by Jiangsu Provincial Natural Science Foundation
文摘Linear differential-algebraic equations (DAEs) with time-varying coefficients A(t)x(1)(t) + B(t)x(t) = q(t), which are tractable with a higher index. are discussed. Their essential properties are investigated. Some equivalent system,,; are given. Using them the paper shows how to state properly initial and boundary conditions for these DAEs. The existence and uniqueness theory of the solution of the initial and boundary value problems for higher index DAEs are proposed.
文摘An efficient numerical approach for the general thermomechanical problems was developed and it was tested for a two-dimensional thermoelasticity problem. The main idea of our numerical method is based on the reduction procedure of the original system of PDEs describing coupled thermomechanical behavior to a system of Differential Algebraic Equations (DAEs) where the stress-strain relationships are treated as algebraic equations. The resulting system of DAEs was then solved with a Backward Differentiation Formula (BDF) using a fully implicit algorithm. The described procedure was explained in detail, and its effectiveness was demonstrated on the solution of a transient uncoupled thermoelastic problem, for which an analytical solution is known, as well as on a fully coupled problem in the two-dimensional case.
基金Project supported by the Foundation for the Authors of the National Excellent Doctoral Thesis Award of China (200720)
文摘In this article, linear regular index 2 DAEs A(t)[D(t)x(t)]' + B(t)x(t) = q(t) are considered. Using a decoupling technique, initial condition and boundary condition are properly formulated. Regular index 1 DAEs are obtained by a regularization method. We study the behavior of the solution of the regularization system via asymptotic expansions. The error analysis between the solutions of the DAEs and its regularization system is given.
基金supported by the Natural Science Foundation of China(NSFC)under grant 11501436Young Talent fund of University Association for Science and Technology in Shaanxi,China(20170701)
文摘Nonlinear fractional differential-algebraic equations often arise in simulating integrated circuits with superconductors. How to obtain the nonnegative solutions of the equations is an important scientific problem. As far as we known, the nonnegativity of solutions of the nonlinear fractional differential-algebraic equations is still not studied. In this article, we investigate the nonnegativity of solutions of the equations. Firstly, we discuss the existence of nonnegative solutions of the equations, and then we show that the nonnegative solution can be approached by a monotone waveform relaxation sequence provided the initial iteration is chosen properly. The choice of initial iteration is critical and we give a method of finding it. Finally, we present an example to illustrate the efficiency of our method.
基金supported by the National Natural Science Foundations of China(Nos.12172186 and 11772166).
文摘The objective of dynamical system learning tasks is to forecast the future behavior of a system by leveraging observed data.However,such systems can sometimes exhibit rigidity due to significant variations in component parameters or the presence of slow and fast variables,leading to challenges in learning.To overcome this limitation,we propose a multiscale differential-algebraic neural network(MDANN)method that utilizes Lagrangian mechanics and incorporates multiscale information for dynamical system learning.The MDANN method consists of two main components:the Lagrangian mechanics module and the multiscale module.The Lagrangian mechanics module embeds the system in Cartesian coordinates,adopts a differential-algebraic equation format,and uses Lagrange multipliers to impose constraints explicitly,simplifying the learning problem.The multiscale module converts high-frequency components into low-frequency components using radial scaling to learn subprocesses with large differences in velocity.Experimental results demonstrate that the proposed MDANN method effectively improves the learning of dynamical systems under rigid conditions.
文摘Two cubical 3D electric circuits with single and double capacitors and twelve ohmic resistors are considered. The resistors are the sides of the cube. The circuit is fed with a single internal emf. The charge on the capacitor(s) and the current distributions of all twelve sides of the circuit(s) vs. time are evaluated. The analysis requires solving twelve differential-algebraic intertwined symbolic equations. This is accomplished by applying a Computer Algebra System (CAS), specifically Mathematica. The needed codes are included. For a set of values assigned to the elements, the numeric results are depicted.
文摘In this article,we investigate a fractional-order singular Leslie-Gower prey-predator bioeconomic model,which describes the interaction bet ween populations of prey and predator,and takes into account the economic interest.We firstly obtain the solvability condition and the st ability of the model sys tem,and discuss the singularity induced bifurcation phenomenon.Next,we introduce a st ate feedback controller to elimina te the singularity induced bifurcation phenomenon,and discuss the optimal control problems.Finally,numerical solutions and their simulations are considered in order to illustrate the theoretical results and reveal the more complex dynamical behavior.
文摘In this paper, Lyapunov-like exponential stability and unstability of differentialalgebraic equation are considered from the viewpoint of stability of system motion, and the criteria of exponential stability and unstability of nonlinear nonautonomous differential-algebraic equation are given by using Lyapunov-like function similar to ordinary differential equation.
基金Project supported by Shanghai Pujiang Program(No.07pj14073)and Shanghai Leading Academic Discipline Project(No.Y0103)
文摘A nonlinear mathematical model for the analysis of large deformation of frame structures with discontinuity conditions and initial displacements, subject to dynamic loads is formulated with arc-coordinates. The differential quadrature element method (DQEM) is then applied to discretize the nonlinear mathematical model in the spatial domain, An effective method is presented to deal with discontinuity conditions of multivariables in the application of DQEM. A set of DQEM discretization equations are obtained, which are a set of nonlinear differential-algebraic equations with singularity in the time domain. This paper also presents a method to solve nonlinear differential-algebra equations. As application, static and dynamical analyses of large deformation of frames and combined frame structures, subjected to concentrated and distributed forces, are presented. The obtained results are compared with those in the literatures. Numerical results show that the proposed method is general, and effective in dealing with disconti- nuity conditions of multi-variables and solving differential-algebraic equations. It requires only a small number of nodes and has low computation complexity with high precision and a good convergence property.
基金Project supported by the National Natural Science Foundation of China(No.11432010)the Doctoral Program Foundation of Education Ministry of China(No.20126102110023)+2 种基金the 111Project of China(No.B07050)the Fundamental Research Funds for the Central Universities(No.310201401JCQ01001)the Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University(No.CX201517)
文摘Projected Runge-Kutta (R-K) methods for constrained Hamiltonian systems are proposed. Dynamic equations of the systems, which are index-3 differential-algebraic equations (DAEs) in the Heisenberg form, are established under the framework of Lagrangian multipliers. R-K methods combined with the technique of projections are then used to solve the DAEs. The basic idea of projections is to eliminate the constraint violations at the position, velocity, and acceleration levels, and to preserve the total energy of constrained Hamiltonian systems by correcting variables of the position, velocity, acceleration, and energy. Numerical results confirm the validity and show the high precision of the proposed method in preserving three levels of constraints and total energy compared with results reported in the literature.
基金This work is supported by NSF of China(10971175)Specialized Research Fund for the Doctoral Program of Higher Education of China(20094301110001)+2 种基金Program for Changjiang Scholars and Innovative Research Team in University of China(IRT1179)NSF of Hunan Province(10JJ7001)the Aid Program for Science and Technology Innovative Research Team in Higher Educational Institutions of Hunan Province,and Fund Project of Hunan Province Education Office(11C1220).
文摘Linear multistep methods and one-leg methods are applied to a class of index-2 nonlinear differential-algebraic equations with a variable delay.The corresponding convergence results are obtained and successfully confirmed by some numerical examples.The results obtained in this work extend the corresponding ones in literature.
文摘This paper develops a class of general one-step discretization methods for solving the index-1 stochastic delay differential-algebraic equations. The existence and uniqueness theorem of strong solutions of index-1 equations is given. A strong convergence criterion of the methods is derived, which is applicable to a series of one-step stochastic numerical methods. Some specific numerical methods, such as the Euler-Maruyama method, stochastic ^-methods, split-step ^-methods are proposed, and their strong convergence results are given. Numerical experiments further illustrate the theoretical results.
基金supported by the National Natural Science Foundation of China(No.11971412)the Natural Science Foundation of Hunan Province of China(No.2018JJ2378)Scientific Research Fund of Hunan Provincial Science and Technology Department(No.2018WK4006).
文摘In this paper,we propose a novel Legendre neural network combined with the extreme learning machine algorithm to solve variable coefficients linear delay differential-algebraic equations with weak discontinuities.First,the solution interval is divided into multiple subintervals by weak discontinuity points.Then,Legendre neural network is used to eliminate the hidden layer by expanding the input pattern using Legendre polynomials on each subinterval.Finally,the parameters of the neural network are obtained by training with the extreme learning machine.The numerical examples show that the proposed method can effectively deal with the difficulty of numerical simulation caused by the discontinuities.
基金supported by the National Natural Science Foundation of China(10772136)Shanghai Leading Academic Discipline Project(B302)The authors wish to thank Dr.Guyue Jiao for the literary suggestions on the manuscript
文摘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.
基金This project is supported by National Natural Science Foundation of China (No.50305033, No.60375020)National Basic Research Program of China(973 Program, No.2004CB719400, No.2002CB312106)Provincial Natural Science Foundation of Zhejiang, China(No.Y105430).
文摘Based on constructing programmed constraint and constraint perturbation equation, a kinematics and dynamics numerical simulation model is established for virtual mechanism, in which the difference scheme guarantee precision in simulation procedure and its mtmerical solutions satisfy programmed manifold stability. A crank-piston mechanism in a car engine, a steering mechanism and a suspension mechanism are simulated in a virtual environment, then comparing the simulation results with those obtained in ADAMS under the same circumstances proved the solver valid.