The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple ba...The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple back-and-forth motion develops into water drinking motion. The solution to a drinking bird equation of motion manifests itself the transition from thermodynamic equilibrium to nonequilibrium irreversible states. The independent solution signifying a nonequilibrium thermal state seems to be constructed as if two independent bifurcation solutions are synthesized, and so, the solution is tentatively termed as the bifurcation-integration solution. The bifurcation-integration solution expresses the transition from mechanical and thermodynamic equilibrium to a nonequilibrium irreversible state, which is explicitly shown by the nonlinear differential equation with time-dependent coefficients (NDE-TC). The analysis established a new theoretical approach to nonequilibrium irreversible states, thermomechanical dynamics (TMD). The TMD method enables one to obtain thermodynamically consistent and time-dependent progresses of thermodynamic quantities, by employing the bifurcation-integration solutions of NDE-TC. We hope that the basic properties of bifurcation-integration solutions will be studied and investigated further in mathematics, physics, chemistry and nonlinear sciences in general.展开更多
In this paper, we study the existence of the transcendental meromorphic solution of the delay differential equations , where a(z) is a rational function, and are polynomials in w(z) with rational c...In this paper, we study the existence of the transcendental meromorphic solution of the delay differential equations , where a(z) is a rational function, and are polynomials in w(z) with rational coefficients, k is a positive integer. Under the assumption when above equations own transcendental meromorphic solutions with minimal hyper-type, we derive the concrete conditions on the degree of the right side of them. Specially, when w(z)=0 is a root of , its multiplicity is at most k. Some examples are given here to illustrate that our results are accurate.展开更多
In this article we derive a general differential equation that describes long-term economic growth in terms of cyclical and trend components. Equation is based on the model of non-linear accelerator of induced investm...In this article we derive a general differential equation that describes long-term economic growth in terms of cyclical and trend components. Equation is based on the model of non-linear accelerator of induced investment. A scheme is proposed for obtaining approximate solutions of nonlinear differential equation by splitting solution into the rapidly oscillating business cycles and slowly varying trend using Krylov-Bogoliubov-Mitropolsky averaging. Simplest modes of the economic system are described. Characteristics of the bifurcation point are found and bifurcation phenomenon is interpreted as loss of stability making the economic system available to structural change and accepting innovations. System being in a nonequilibrium state has a dynamics with self-sustained undamped oscillations. The model is verified with economic development of the US during the fifth Kondratieff cycle (1982-2010). Model adequately describes real process of economic growth in both quantitative and qualitative aspects. It is one of major results that the model gives a rough estimation of critical points of system stability loss and falling into a crisis recession. The model is used to forecast the macroeconomic dynamics of the US during the sixth Kondratieff cycle (2018-2050). For this forecast we use fixed production capital functional dependence on a long-term Kondratieff cycle and medium-term Juglar and Kuznets cycles. More accurate estimations of the time of crisis and recession are based on the model of accelerating log-periodic oscillations. The explosive growth of the prices of highly liquid commodities such as gold and oil is taken as real predictors of the global financial crisis. The second wave of crisis is expected to come in June 2011.展开更多
Let f be an entire solution of the Tumura-Clunie type non-linear delay differential equation.We mainly investigate the dynamical properties of Julia sets of f,and the lower bound estimates of the measure of related li...Let f be an entire solution of the Tumura-Clunie type non-linear delay differential equation.We mainly investigate the dynamical properties of Julia sets of f,and the lower bound estimates of the measure of related limiting directions is verified.展开更多
In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Trans...In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.展开更多
The dot product of the bases vectors on the super-surface of the non-linear nonholonomic constraints with one order, expressed by quasi-coorfinates, and Mishirskiiequalions are regarded as the fundamental equations of...The dot product of the bases vectors on the super-surface of the non-linear nonholonomic constraints with one order, expressed by quasi-coorfinates, and Mishirskiiequalions are regarded as the fundamental equations of dynamics with non-linear andnon-holononlic constraints in one order for the system of the variable mass. From thesethe variant ddferential-equations of dynamics expressed by quasi-coordinates arederived. The fundamental equations of dynamics are compatible with the principle ofJourdain. A case is cited.展开更多
In this paper we are concerned with the oscillation criteria of second order non-linear homogeneous differential equation. Example have been given to illustrate the results.
By the use of the Liapunov functional approach, a new result is obtained to ascertain the asymptotic stability of zero solution of a certain fourth-order non-linear differential equation with delay. The established re...By the use of the Liapunov functional approach, a new result is obtained to ascertain the asymptotic stability of zero solution of a certain fourth-order non-linear differential equation with delay. The established result is less restrictive than those reported in the literature.展开更多
In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarant...In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.展开更多
The basic objects of investigation in this article are nonlinear impulsive dif- ferential equations. The impulsive moments coincide with the moments of meeting of the integral curve and some of the so-called barrier c...The basic objects of investigation in this article are nonlinear impulsive dif- ferential equations. The impulsive moments coincide with the moments of meeting of the integral curve and some of the so-called barrier curves. For such type of equations, suf- ficient conditions are found under which the solutions are continuously dependent on the perturbations with respect to the initial conditions and barrier curves. The results are applied to a mathematical model of population dynamics.展开更多
We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by sol...We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.展开更多
The numerical analysis of the approximate inertial manifold in,weakly damped forced KdV equation is given. The results of numerical analysis under five models is the same as that of nonlinear spectral analysis.
There are given sufficient conditions for the ultimate boundedness of solutions and for the existence of periodic solutions of a certain vector differential equation of third-order.
Several nonlinear three-dimensional systems of ordinary differential equations are studied analytically and numerically in this paper in accordance with universal bifurcation theory of Feigenbaum-Sharkovskii-Magnitsky...Several nonlinear three-dimensional systems of ordinary differential equations are studied analytically and numerically in this paper in accordance with universal bifurcation theory of Feigenbaum-Sharkovskii-Magnitsky [1] [2]. All systems are autonomous and dissipative and display chaotic behaviour. The analysis confirms that transition to chaos in such systems is performed through cascades of bifurcations of regular attractors.展开更多
In this paper the authors study a class of non-linear singular partial differential equation in complex domain C-t x C-x(n). Under certain assumptions, they prove the existence and uniqueness of holomorphic solution n...In this paper the authors study a class of non-linear singular partial differential equation in complex domain C-t x C-x(n). Under certain assumptions, they prove the existence and uniqueness of holomorphic solution near origin of C-t x C-x(n).展开更多
We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatia...We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.展开更多
The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic ...The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic flow analysis,simulation,autonomous vehicle development,etc.Two-dimensional(2D)vehicular movement is basically stochastic and is the result of interactions between a driver's behavior and a vehicle's characteristics.Current microscopic models either neglect 2D noise,or overlook vehicle dynamics.The modeling capabilities,thus,are limited,so that stochastic lateral movement cannot be reproduced.The present research extends an intelligent driver model(IDM)by explicitly considering both vehicle dynamics and 2D noises to formulate a stochastic 2D IDM model,with vehicle dynamics based on the stochastic differential equation(SDE)theory.Control inputs from the vehicle include the steer rate and longitudinal acceleration,both of which are developed based on an idea from a traditional intelligent driver model.The stochastic stability condition is analyzed on the basis of Lyapunov theory.Numerical analysis is used to assess the two cases:(i)when a vehicle accelerates from a standstill and(ii)when a platoon of vehicles follow a leader with a stop-and-go speed profile,the formation of congestion and subsequent dispersion are simulated.The results show that the model can reproduce the stochastic 2D trajectories of the vehicle and the marginal distribution of lateral movement.The proposed model can be used in both a simulation platform and a behavioral analysis of a human driver in traffic flow.展开更多
There are vast constraint equations in conventional dynamics analysis of deployable structures,which lead to differential-algebraic equations(DAEs)solved hard.To reduce the difficulty of solving and the amount of equa...There are vast constraint equations in conventional dynamics analysis of deployable structures,which lead to differential-algebraic equations(DAEs)solved hard.To reduce the difficulty of solving and the amount of equations,a new flexible multibody dynamics analysis methodology of deployable structures with scissor-like elements(SLEs)is presented.Firstly,a precise model of a flexible bar of SLE is established by the higher order shear deformable beam element based on the absolute nodal coordinate formulation(ANCF),and the master/slave freedom method is used to obtain the dynamics equations of SLEs without constraint equations.Secondly,according to features of deployable structures,the specification matrix method(SMM)is proposed to eliminate the constraint equations among SLEs in the frame of ANCF.With this method,the inner and the boundary nodal coordinates of element characteristic matrices can be separated simply and efficiently,especially on condition that there are vast nodal coordinates.So the element characteristic matrices can be added end to end circularly.Thus,the dynamic model of deployable structure reduces dimension and can be assembled without any constraint equation.Next,a new iteration procedure for the generalized-a algorithm is presented to solve the ordinary differential equations(ODEs)of deployable structure.Finally,the proposed methodology is used to analyze the flexible multi-body dynamics of a planar linear array deployable structure based on three scissor-like elements.The simulation results show that flexibility has a significant influence on the deployment motion of the deployable structure.The proposed methodology indeed reduce the difficulty of solving and the amount of equations by eliminating redundant degrees of freedom and the constraint equations in scissor-like elements and among scissor-like elements.展开更多
The problem of preserving fidelity in numerical computation of nonlinear or-dinary differential equations is studied in terms of preserving local differential structure and approximating global integration structure o...The problem of preserving fidelity in numerical computation of nonlinear or-dinary differential equations is studied in terms of preserving local differential structure and approximating global integration structure of the dynamical system. The ordinary differen-tial equations are lifted to the corresponding partial differential equations in the framework of algebraic dynamics,and a new algorithm——algebraic dynamics algorithm is proposed based on the exact analytical solutions of the ordinary differential equations by the alge-braic dynamics method. In the new algorithm,the time evolution of the ordinary differential system is described locally by the time translation operator and globally by the time evo-lution operator. The exact analytical piece-like solution of the ordinary differential equa-tions is expressd in terms of Taylor series with a local convergent radius,and its finite order truncation leads to the new numerical algorithm with a controllable precision better than Runge Kutta Algorithm and Symplectic Geometric Algorithm.展开更多
文摘The new independent solutions of the nonlinear differential equation with time-dependent coefficients (NDE-TC) are discussed, for the first time, by employing experimental device called a drinking bird whose simple back-and-forth motion develops into water drinking motion. The solution to a drinking bird equation of motion manifests itself the transition from thermodynamic equilibrium to nonequilibrium irreversible states. The independent solution signifying a nonequilibrium thermal state seems to be constructed as if two independent bifurcation solutions are synthesized, and so, the solution is tentatively termed as the bifurcation-integration solution. The bifurcation-integration solution expresses the transition from mechanical and thermodynamic equilibrium to a nonequilibrium irreversible state, which is explicitly shown by the nonlinear differential equation with time-dependent coefficients (NDE-TC). The analysis established a new theoretical approach to nonequilibrium irreversible states, thermomechanical dynamics (TMD). The TMD method enables one to obtain thermodynamically consistent and time-dependent progresses of thermodynamic quantities, by employing the bifurcation-integration solutions of NDE-TC. We hope that the basic properties of bifurcation-integration solutions will be studied and investigated further in mathematics, physics, chemistry and nonlinear sciences in general.
文摘In this paper, we study the existence of the transcendental meromorphic solution of the delay differential equations , where a(z) is a rational function, and are polynomials in w(z) with rational coefficients, k is a positive integer. Under the assumption when above equations own transcendental meromorphic solutions with minimal hyper-type, we derive the concrete conditions on the degree of the right side of them. Specially, when w(z)=0 is a root of , its multiplicity is at most k. Some examples are given here to illustrate that our results are accurate.
文摘In this article we derive a general differential equation that describes long-term economic growth in terms of cyclical and trend components. Equation is based on the model of non-linear accelerator of induced investment. A scheme is proposed for obtaining approximate solutions of nonlinear differential equation by splitting solution into the rapidly oscillating business cycles and slowly varying trend using Krylov-Bogoliubov-Mitropolsky averaging. Simplest modes of the economic system are described. Characteristics of the bifurcation point are found and bifurcation phenomenon is interpreted as loss of stability making the economic system available to structural change and accepting innovations. System being in a nonequilibrium state has a dynamics with self-sustained undamped oscillations. The model is verified with economic development of the US during the fifth Kondratieff cycle (1982-2010). Model adequately describes real process of economic growth in both quantitative and qualitative aspects. It is one of major results that the model gives a rough estimation of critical points of system stability loss and falling into a crisis recession. The model is used to forecast the macroeconomic dynamics of the US during the sixth Kondratieff cycle (2018-2050). For this forecast we use fixed production capital functional dependence on a long-term Kondratieff cycle and medium-term Juglar and Kuznets cycles. More accurate estimations of the time of crisis and recession are based on the model of accelerating log-periodic oscillations. The explosive growth of the prices of highly liquid commodities such as gold and oil is taken as real predictors of the global financial crisis. The second wave of crisis is expected to come in June 2011.
基金supported by the National Natural Science Foundation of China(12171050,12071047)the Fundamental Research Funds for the Central Universities(500421126)。
文摘Let f be an entire solution of the Tumura-Clunie type non-linear delay differential equation.We mainly investigate the dynamical properties of Julia sets of f,and the lower bound estimates of the measure of related limiting directions is verified.
文摘In this paper, we discuss a new method employed to tackle non-linear partial differential equations, namely Double Elzaki Transform Decomposition Method (DETDM). This method is a combination of the Double ELzaki Transform and Adomian Decomposition Method. This technique is hereafter provided and supported with necessary illustrations, together with some attached examples. The results reveal that the new method is very efficient, simple and can be applied to other non-linear problems.
文摘The dot product of the bases vectors on the super-surface of the non-linear nonholonomic constraints with one order, expressed by quasi-coorfinates, and Mishirskiiequalions are regarded as the fundamental equations of dynamics with non-linear andnon-holononlic constraints in one order for the system of the variable mass. From thesethe variant ddferential-equations of dynamics expressed by quasi-coordinates arederived. The fundamental equations of dynamics are compatible with the principle ofJourdain. A case is cited.
文摘In this paper we are concerned with the oscillation criteria of second order non-linear homogeneous differential equation. Example have been given to illustrate the results.
文摘By the use of the Liapunov functional approach, a new result is obtained to ascertain the asymptotic stability of zero solution of a certain fourth-order non-linear differential equation with delay. The established result is less restrictive than those reported in the literature.
基金Supported by Russian Fund of Fund amental Investigations(Pr.990101064)and Russian Minister of Educatin
文摘In this paper the method of design of kinematical and dynamical equations of mechanical systems, applied to numerical ealization, is proposed. The corresponding difference equations, which are obtained, give a guarantee of computations with a given precision. The equations of programmed constraints and those of constraint perturbations are defined. The stability of the programmed manifold for numerical solutions of the kinematical and dynamical equations is obtained by corresponding construction of the constraint perturbation equations. The dynamical equations of system with programmed constraints are set up in the form of Lagrange’s equations in generalized coordinates. Certain inverse problems of rigid body dynamics are examined.
文摘The basic objects of investigation in this article are nonlinear impulsive dif- ferential equations. The impulsive moments coincide with the moments of meeting of the integral curve and some of the so-called barrier curves. For such type of equations, suf- ficient conditions are found under which the solutions are continuously dependent on the perturbations with respect to the initial conditions and barrier curves. The results are applied to a mathematical model of population dynamics.
基金The work of A.Kurganov was supported in part by the National Natural Science Foundation of China grant 11771201by the fund of the Guangdong Provincial Key Laboratory of Computational Science and Material Design(No.2019B030301001).
文摘We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.
文摘The numerical analysis of the approximate inertial manifold in,weakly damped forced KdV equation is given. The results of numerical analysis under five models is the same as that of nonlinear spectral analysis.
文摘There are given sufficient conditions for the ultimate boundedness of solutions and for the existence of periodic solutions of a certain vector differential equation of third-order.
文摘Several nonlinear three-dimensional systems of ordinary differential equations are studied analytically and numerically in this paper in accordance with universal bifurcation theory of Feigenbaum-Sharkovskii-Magnitsky [1] [2]. All systems are autonomous and dissipative and display chaotic behaviour. The analysis confirms that transition to chaos in such systems is performed through cascades of bifurcations of regular attractors.
文摘In this paper the authors study a class of non-linear singular partial differential equation in complex domain C-t x C-x(n). Under certain assumptions, they prove the existence and uniqueness of holomorphic solution near origin of C-t x C-x(n).
文摘We propose a novel computational framework that is capable of employing different time integration algorithms and different space discretized methods such as the Finite Element Method,particle methods,and other spatial methods on a single body sub-dividedintomultiple subdomains.This is in conjunctionwithimplementing thewell known Generalized Single Step Single Solve(GS4)family of algorithms which encompass the entire scope of Linear Multistep algorithms that have been developed over the past 50 years or so and are second order accurate into the Differential Algebraic Equation framework.In the current state of technology,the coupling of altogether different time integration algorithms has been limited to the same family of algorithms such as theNewmarkmethods and the coupling of different algorithms usually has resulted in reduced accuracy in one or more variables including the Lagrange multiplier.However,the robustness and versatility of the GS4 with its ability to accurately account for the numerical shifts in various time schemes it encompasses,overcomes such barriers and allows a wide variety of arbitrary implicit-implicit,implicit-explicit,and explicit-explicit pairing of the various time schemes while maintaining the second order accuracy in time for not only all primary variables such as displacement,velocity and acceleration but also the Lagrange multipliers used for coupling the subdomains.By selecting an appropriate spatialmethod and time scheme on the area with localized phenomena contrary to utilizing a single process on the entire body,the proposed work has the potential to better capture the physics of a given simulation.The method is validated by solving 2D problems for the linear second order systems with various combination of spatial methods and time schemes with great flexibility.The accuracy and efficacy of the present work have not yet been seen in the current field,and it has shown significant promise in its capabilities and effectiveness for general linear dynamics through numerical examples.
基金Project supported by the National Key Research and Development Program of China(Grant No.2021YFE0194400)the National Natural Science Foundation of China(Grant Nos.52272314 and 52131202)+1 种基金the Fund for Humanities and Social Science from the Ministry of Education of China(Grant No.21YJCZH116)the Public Welfare Scientific Research Project(Grant No.LGF22E080007)。
文摘The law of vehicle movement has long been studied under the umbrella of microscopic traffic flow models,especially the car-following(CF)models.These models of the movement of vehicles serve as the backbone of traffic flow analysis,simulation,autonomous vehicle development,etc.Two-dimensional(2D)vehicular movement is basically stochastic and is the result of interactions between a driver's behavior and a vehicle's characteristics.Current microscopic models either neglect 2D noise,or overlook vehicle dynamics.The modeling capabilities,thus,are limited,so that stochastic lateral movement cannot be reproduced.The present research extends an intelligent driver model(IDM)by explicitly considering both vehicle dynamics and 2D noises to formulate a stochastic 2D IDM model,with vehicle dynamics based on the stochastic differential equation(SDE)theory.Control inputs from the vehicle include the steer rate and longitudinal acceleration,both of which are developed based on an idea from a traditional intelligent driver model.The stochastic stability condition is analyzed on the basis of Lyapunov theory.Numerical analysis is used to assess the two cases:(i)when a vehicle accelerates from a standstill and(ii)when a platoon of vehicles follow a leader with a stop-and-go speed profile,the formation of congestion and subsequent dispersion are simulated.The results show that the model can reproduce the stochastic 2D trajectories of the vehicle and the marginal distribution of lateral movement.The proposed model can be used in both a simulation platform and a behavioral analysis of a human driver in traffic flow.
基金Supported by National Natural Science Foundation of China(Grant No.51175422)
文摘There are vast constraint equations in conventional dynamics analysis of deployable structures,which lead to differential-algebraic equations(DAEs)solved hard.To reduce the difficulty of solving and the amount of equations,a new flexible multibody dynamics analysis methodology of deployable structures with scissor-like elements(SLEs)is presented.Firstly,a precise model of a flexible bar of SLE is established by the higher order shear deformable beam element based on the absolute nodal coordinate formulation(ANCF),and the master/slave freedom method is used to obtain the dynamics equations of SLEs without constraint equations.Secondly,according to features of deployable structures,the specification matrix method(SMM)is proposed to eliminate the constraint equations among SLEs in the frame of ANCF.With this method,the inner and the boundary nodal coordinates of element characteristic matrices can be separated simply and efficiently,especially on condition that there are vast nodal coordinates.So the element characteristic matrices can be added end to end circularly.Thus,the dynamic model of deployable structure reduces dimension and can be assembled without any constraint equation.Next,a new iteration procedure for the generalized-a algorithm is presented to solve the ordinary differential equations(ODEs)of deployable structure.Finally,the proposed methodology is used to analyze the flexible multi-body dynamics of a planar linear array deployable structure based on three scissor-like elements.The simulation results show that flexibility has a significant influence on the deployment motion of the deployable structure.The proposed methodology indeed reduce the difficulty of solving and the amount of equations by eliminating redundant degrees of freedom and the constraint equations in scissor-like elements and among scissor-like elements.
基金This work was supported in part by the National Natural Science Foundation of China (Grant Nos.10375039 and 90503008)the Doctoral Program Foundation from the Ministry of Education of China,and the Center of Nuclear Physics of HIRFL of China.
文摘The problem of preserving fidelity in numerical computation of nonlinear or-dinary differential equations is studied in terms of preserving local differential structure and approximating global integration structure of the dynamical system. The ordinary differen-tial equations are lifted to the corresponding partial differential equations in the framework of algebraic dynamics,and a new algorithm——algebraic dynamics algorithm is proposed based on the exact analytical solutions of the ordinary differential equations by the alge-braic dynamics method. In the new algorithm,the time evolution of the ordinary differential system is described locally by the time translation operator and globally by the time evo-lution operator. The exact analytical piece-like solution of the ordinary differential equa-tions is expressd in terms of Taylor series with a local convergent radius,and its finite order truncation leads to the new numerical algorithm with a controllable precision better than Runge Kutta Algorithm and Symplectic Geometric Algorithm.