This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended met...This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended methods are capable of calculating the continuation curves of the equilibrium points for the particular type of trimming flight. Therefore, these methods can not only give the performance measures of aircraft, but also determine the stability of trimming points. In this paper, the methods are used to verify the effectiveness of the thrust vectoring control law, to define the flight envelope boundary, to analyze the stability and controllability of trimming flight, and to predict the departures of the instable flight. The result shows that the extended methods provide more flight dynamic information and are useful in preliminary design of the thrust vectoring aircraft.展开更多
In this paper, we propose a new smooth function that possesses a property not satisfied by the existing smooth functions. Based on this smooth function, we discuss the existence and continuity of the smoothing path fo...In this paper, we propose a new smooth function that possesses a property not satisfied by the existing smooth functions. Based on this smooth function, we discuss the existence and continuity of the smoothing path for solving the P0 function nonlinear complementarity problem (NCP). Using the characteristics of the new smooth function, we investigate the boundedness of the iteration sequence generated by the non-interior continuation methods for solving the P0 function NCP under the assumption that the solution set of the NCP is nonempty and bounded. We show that the assumption that the solution set of the NCP is nonempty and bounded is weaker than those required by a few existing continuation methods for solving the NCP.展开更多
Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed ...Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform（FFT）algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.展开更多
Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analy...Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.展开更多
Continuation method solving forward kinematics problem of parallel robot was discussed. And through a coefficient-parameter continuation method the efficiency and feasibility of continuation method were improved. Usin...Continuation method solving forward kinematics problem of parallel robot was discussed. And through a coefficient-parameter continuation method the efficiency and feasibility of continuation method were improved. Using this method all forward solutions of a new parallel robot model which was put forward lately by Robot Open Laboratory of Science Institute of China were obtained. Therefore it provided the basis of mechanism analysis and real-time control for new model.展开更多
A noninterior continuation method is proposed for semidefinite complementarity problem （SDCP）. This method improves the noninterior continuation methods recently developed for SDCP by Chen and Tseng. The main proper...A noninterior continuation method is proposed for semidefinite complementarity problem （SDCP）. This method improves the noninterior continuation methods recently developed for SDCP by Chen and Tseng. The main properties of our method are： （i） it is well d.efined for the monotones SDCP; （ii） it has to solve just one linear system of equations at each step; （iii） it is shown to be both globally linearly convergent and locally quadratically convergent under suitable assumptions.展开更多
An idea of relaxing the effect of delay when computing the Runge-Kutta stages in the current step and a class of two-step continuity Runge-Kutta methods （TSCRK） is presented. Their construction, their order conditio...An idea of relaxing the effect of delay when computing the Runge-Kutta stages in the current step and a class of two-step continuity Runge-Kutta methods （TSCRK） is presented. Their construction, their order conditions and their convergence are studied. The two-step continuity Runge-Kutta methods possess good numerical stability properties and higher stage-order, and keep the explicit process of computing the Runge-Kutta stages. The numerical experiments show that the TSCRK methods are efficient.展开更多
In this paper we present a homotopy continuation method for finding the Karush-Kuhn-Tucker point of a class of nonlinear non-convex programming problems. Two numerical examples are given to show that this method is ef...In this paper we present a homotopy continuation method for finding the Karush-Kuhn-Tucker point of a class of nonlinear non-convex programming problems. Two numerical examples are given to show that this method is effective. It should be pointed out that we extend the results of Lin et al. （see Appl. Math. Comput., 80（1996）, 209-224） to a broader class of non-convex programming problems.展开更多
By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved hav...By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.展开更多
In this paper,according to the Fujian Seismic Network earthquake catalog records,the T now method and the Four Stations Continuous Location method( hereinafter called FSCL)put forward by Jin Xing are inspected by usin...In this paper,according to the Fujian Seismic Network earthquake catalog records,the T now method and the Four Stations Continuous Location method( hereinafter called FSCL)put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations of each seismic event. Results show that for earthquakes within the network,both methods can obtain similar location results and location deviations are small for the majority of the events. For earthquakes outside the network,the location deviation may be amplified as the epicentral distance increases,owing to the seismic station distribution which spread toward the side of the epicenter and the small opening angle between seismic stations used for locating and epicenter. For the FSCL method,the impacts of the wave velocity on the location results may be significant for earthquakes outside the network.Thus,selecting a velocity model which is similar to the actual structure of the wave velocity will contribute to improving location results of earthquakes. The FSCL method can locate more seismic events than the T now method. It concludes that the T now method makes use of mistake information from some non-triggering stations in earthquake catalog,and some P-wave arrivals are not included in the earthquake catalog due to discontinuous records or unclear records of the seismic phase,which induces incorrect location.展开更多
A sufficient condition is given to assert that a continuous mapping between Rm and Rn has a zero. The constructive proof of the result is based upon continuation methods and supplies the existence of a path leading to...A sufficient condition is given to assert that a continuous mapping between Rm and Rn has a zero. The constructive proof of the result is based upon continuation methods and supplies the existence of a path leading to the zero point.展开更多
It is proven that an autonomous system verifying some conditions has at least one stable stationary trajectory and it is also given a lower bound to the number of unstable stationary trajectorlies.
Sufficient conditions are given to assert that two differentiable mappings between Banach spaces have common values. The proof is essentially based upon continuation methods.
Sufficient conditions are given to assert that two C1-mappings share only one value in a connected compact Banach manifold modelled over Rn. The proof of the result, which is based upon continuation methods, is constr...Sufficient conditions are given to assert that two C1-mappings share only one value in a connected compact Banach manifold modelled over Rn. The proof of the result, which is based upon continuation methods, is constructive.展开更多
Sufficient conditions were given to assert that between any two Banach spaces over K, Fredholm mappings share at least one .value in a specific open ball. The proof of the result is constructive and based upon continu...Sufficient conditions were given to assert that between any two Banach spaces over K, Fredholm mappings share at least one .value in a specific open ball. The proof of the result is constructive and based upon continuation methods.展开更多
The present work introduces a novel concurrent optimization formulation to meet the requirements of lightweight design and various constraints simultaneously.Nodal displacement of macrostructure and effective thermal ...The present work introduces a novel concurrent optimization formulation to meet the requirements of lightweight design and various constraints simultaneously.Nodal displacement of macrostructure and effective thermal conductivity of microstructure are regarded as the constraint functions, which means taking into account both the loadcarrying capabilities and the thermal insulation properties.The effective properties of porous material derived from numerical homogenization are used for macrostructural analysis. Meanwhile, displacement vectors of macrostructures from original and adjoint load cases are used for sensitivity analysis of the microstructure. Design variables in the form of reciprocal functions of relative densities are introduced and used for linearization of the constraint function. The objective function of total mass is approximately expressed by the second order Taylor series expansion. Then, the proposed concurrent optimization problem is solved using a sequential quadratic programming algorithm, by splitting into a series of sub-problems in the form of the quadratic program. Finally, several numerical examples are presented to validate the effectiveness of the proposed optimization method. The various effects including initial designs, prescribed limits of nodal displacement, and effective thermal conductivity on optimized designs are also investigated. An amount of optimized macrostructures and their corresponding microstructures are achieved.展开更多
It is an important topic to improve the redundancy of optimized configuration to resist the local failure in topology optimization of continuum structures.Such a fail-safe topology optimization problem has been solved...It is an important topic to improve the redundancy of optimized configuration to resist the local failure in topology optimization of continuum structures.Such a fail-safe topology optimization problem has been solved effectively in the ficld of statics.In this paper,the fail-safe topology optimization problem is extended to the field of frequency topology optimization.Based on the independent continuous mapping(ICM)method,the model of fail-safe topology optimization is established with the objective of minimal weight integrating with the discrete condition of topological variables and the constraint of the fundamental frequency.The fail-safe optimization model established above is substituted by a sequence of subproblems in the form of the quadratic program with exact second-order information and solved efficiently by the dual sequence quadratic programming(DSQP)algorithm.The numerical result reveals that the optimized fail-safe structure has more complex configuration and preserved materials than the structure obtained from the traditional frequency topology optimization,which means that the optimized fail-safe structure has higher redundancy.Moreover,the optimized fail-safe structure guarantees that the natural frequency meets the constraint of fundamental frequency when the local failure ocurs,which can avoid the structural frequency to be sensitive to local failure.The fail-safe optimirzation topology model is proved effective and feasible by four numerical examples.展开更多
To develop an effective numerical method for the cable sliding problem in cable structures, two-node catenary cable element was built to model the cables based on analytical solution of elastic catenary. Cooperated wi...To develop an effective numerical method for the cable sliding problem in cable structures, two-node catenary cable element was built to model the cables based on analytical solution of elastic catenary. Cooperated with Newton method, continuation method was used to solve the nonlinear equations. This approach is more efficient than using Newton method only and has a wider range to select initial values for the process to converge. The relationship between the tension on a cable segment and its unstrained length was derived and used to calculate the unbalanced cable tensions at the supports. An example is presented to show the correctness and efficiency of the proposed method.展开更多
This paper presents a method for tracing a planar implicit curve f（x, y）=0 on a rectangular region based on continuation scheme. First, according to the starting track-point and the starting track-direction of the c...This paper presents a method for tracing a planar implicit curve f（x, y）=0 on a rectangular region based on continuation scheme. First, according to the starting track-point and the starting track-direction of the curve, make a new fimction F（x, y）=0 where the same curve withf（x, y）=0 is defined. Then we trace the curve between the two domains where F（x, y）〉0 and F（x, y）〈0 alternately, according to the two rules presented in this paper. Equal step size or adaptive step size can be used, when we trace the curve. An irregular planar implicit curve （such as the curve with large curvatures at some points on the curve）, can be plotted if an adaptive step size is used. Moreover, this paper presents a scheme to search for the multiple points on the curve. Our method has the following advantages： （1） it can plot Co planar implicit curves; （2） it can plot the planar implicit curves with multiple points; （3） by the help of using the two rules, our method does not need to compute the tangent vector at the points on the curve, and directly searches for the direction of the tracing curve; （4） the tracing procedure costs only one of two evaluations of function f（x, y）=0 per moving step, while most existing similar methods cost more evaluations of the function.展开更多
文摘This paper presents the flight dynamical behavior of the thrust vectoring aircraft with extended bifurcation and continuation methods. In contrast to the standard bifurcation and continuation methods, the extended methods are capable of calculating the continuation curves of the equilibrium points for the particular type of trimming flight. Therefore, these methods can not only give the performance measures of aircraft, but also determine the stability of trimming points. In this paper, the methods are used to verify the effectiveness of the thrust vectoring control law, to define the flight envelope boundary, to analyze the stability and controllability of trimming flight, and to predict the departures of the instable flight. The result shows that the extended methods provide more flight dynamic information and are useful in preliminary design of the thrust vectoring aircraft.
基金the National Natural Science Foundation of China (Grant Nos. 19871016 and 19731001).
文摘In this paper, we propose a new smooth function that possesses a property not satisfied by the existing smooth functions. Based on this smooth function, we discuss the existence and continuity of the smoothing path for solving the P0 function nonlinear complementarity problem (NCP). Using the characteristics of the new smooth function, we investigate the boundedness of the iteration sequence generated by the non-interior continuation methods for solving the P0 function NCP under the assumption that the solution set of the NCP is nonempty and bounded. We show that the assumption that the solution set of the NCP is nonempty and bounded is weaker than those required by a few existing continuation methods for solving the NCP.
基金supported by the National Natural Science Foundation of China(41304022,41174026,41104047)the National 973 Foundation(61322201,2013CB733303)+1 种基金the Key laboratory Foundation of Geo-space Environment and Geodesy of the Ministry of Education(13-01-08)the Youth Innovation Foundation of High Resolution Earth Observation(GFZX04060103-5-12)
文摘Downward continuation is a key step in processing airborne geomagnetic data. However,downward continuation is a typically ill-posed problem because its computation is unstable; thus, regularization methods are needed to realize effective continuation. According to the Poisson integral plane approximate relationship between observation and continuation data, the computation formulae combined with the fast Fourier transform（FFT）algorithm are transformed to a frequency domain for accelerating the computational speed. The iterative Tikhonov regularization method and the iterative Landweber regularization method are used in this paper to overcome instability and improve the precision of the results. The availability of these two iterative regularization methods in the frequency domain is validated by simulated geomagnetic data, and the continuation results show good precision.
基金supported by the Na-tional Natural Science Foundation of China(No.52272369).
文摘Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.
文摘Continuation method solving forward kinematics problem of parallel robot was discussed. And through a coefficient-parameter continuation method the efficiency and feasibility of continuation method were improved. Using this method all forward solutions of a new parallel robot model which was put forward lately by Robot Open Laboratory of Science Institute of China were obtained. Therefore it provided the basis of mechanism analysis and real-time control for new model.
基金This work was supported by the National Natural Science Foundation of China (10201001, 70471008)
文摘A noninterior continuation method is proposed for semidefinite complementarity problem （SDCP）. This method improves the noninterior continuation methods recently developed for SDCP by Chen and Tseng. The main properties of our method are： （i） it is well d.efined for the monotones SDCP; （ii） it has to solve just one linear system of equations at each step; （iii） it is shown to be both globally linearly convergent and locally quadratically convergent under suitable assumptions.
文摘An idea of relaxing the effect of delay when computing the Runge-Kutta stages in the current step and a class of two-step continuity Runge-Kutta methods （TSCRK） is presented. Their construction, their order conditions and their convergence are studied. The two-step continuity Runge-Kutta methods possess good numerical stability properties and higher stage-order, and keep the explicit process of computing the Runge-Kutta stages. The numerical experiments show that the TSCRK methods are efficient.
文摘In this paper we present a homotopy continuation method for finding the Karush-Kuhn-Tucker point of a class of nonlinear non-convex programming problems. Two numerical examples are given to show that this method is effective. It should be pointed out that we extend the results of Lin et al. （see Appl. Math. Comput., 80（1996）, 209-224） to a broader class of non-convex programming problems.
基金Project supported by the National Natural Science Foundation of China (No.10471038)
文摘By applying the continuous finite element methods of ordinary differential equations, the linear element methods are proved having second-order pseudo-symplectic scheme and the quadratic element methods are proved having third-order pseudo- symplectic scheme respectively for general Hamiltonian systems, and they both keep energy conservative. The finite element methods are proved to be symplectic as well as energy conservative for linear Hamiltonian systems. The numerical results are in agree-ment with theory.
基金funded by the National Key Technology R&D Program of China(2009BAK55B02)
文摘In this paper,according to the Fujian Seismic Network earthquake catalog records,the T now method and the Four Stations Continuous Location method( hereinafter called FSCL)put forward by Jin Xing are inspected by using P-wave arrival information of the first four stations of each seismic event. Results show that for earthquakes within the network,both methods can obtain similar location results and location deviations are small for the majority of the events. For earthquakes outside the network,the location deviation may be amplified as the epicentral distance increases,owing to the seismic station distribution which spread toward the side of the epicenter and the small opening angle between seismic stations used for locating and epicenter. For the FSCL method,the impacts of the wave velocity on the location results may be significant for earthquakes outside the network.Thus,selecting a velocity model which is similar to the actual structure of the wave velocity will contribute to improving location results of earthquakes. The FSCL method can locate more seismic events than the T now method. It concludes that the T now method makes use of mistake information from some non-triggering stations in earthquake catalog,and some P-wave arrivals are not included in the earthquake catalog due to discontinuous records or unclear records of the seismic phase,which induces incorrect location.
基金This work is partially supported by D.G.E.S. PB 96-1338-CO2-01 and the Junta de Andalucla.
文摘A sufficient condition is given to assert that a continuous mapping between Rm and Rn has a zero. The constructive proof of the result is based upon continuation methods and supplies the existence of a path leading to the zero point.
基金This work is partially supported by D.G.Y.C.T.PB 96-1338-CO 2-01 and the Junta de Andalucía.
文摘It is proven that an autonomous system verifying some conditions has at least one stable stationary trajectory and it is also given a lower bound to the number of unstable stationary trajectorlies.
文摘Sufficient conditions are given to assert that two differentiable mappings between Banach spaces have common values. The proof is essentially based upon continuation methods.
基金partially supported by D.G.E.S.Pb96-1338-CO 2-01 and the Junta de Andalucia
文摘Sufficient conditions are given to assert that two C1-mappings share only one value in a connected compact Banach manifold modelled over Rn. The proof of the result, which is based upon continuation methods, is constructive.
基金Project supported by D.G.E.S. Pb 96-1338-CO 2-01 and the Junta de Andalucia
文摘Sufficient conditions were given to assert that between any two Banach spaces over K, Fredholm mappings share at least one .value in a specific open ball. The proof of the result is constructive and based upon continuation methods.
基金supported by the National Natural Science Foundation of China (Grants 11202078, 51405123)the Fundamental Research Funds for the Central Universities (Grant 2017MS077)
文摘The present work introduces a novel concurrent optimization formulation to meet the requirements of lightweight design and various constraints simultaneously.Nodal displacement of macrostructure and effective thermal conductivity of microstructure are regarded as the constraint functions, which means taking into account both the loadcarrying capabilities and the thermal insulation properties.The effective properties of porous material derived from numerical homogenization are used for macrostructural analysis. Meanwhile, displacement vectors of macrostructures from original and adjoint load cases are used for sensitivity analysis of the microstructure. Design variables in the form of reciprocal functions of relative densities are introduced and used for linearization of the constraint function. The objective function of total mass is approximately expressed by the second order Taylor series expansion. Then, the proposed concurrent optimization problem is solved using a sequential quadratic programming algorithm, by splitting into a series of sub-problems in the form of the quadratic program. Finally, several numerical examples are presented to validate the effectiveness of the proposed optimization method. The various effects including initial designs, prescribed limits of nodal displacement, and effective thermal conductivity on optimized designs are also investigated. An amount of optimized macrostructures and their corresponding microstructures are achieved.
基金the National Natural Science Foundation of China(Grant 11872080).
文摘It is an important topic to improve the redundancy of optimized configuration to resist the local failure in topology optimization of continuum structures.Such a fail-safe topology optimization problem has been solved effectively in the ficld of statics.In this paper,the fail-safe topology optimization problem is extended to the field of frequency topology optimization.Based on the independent continuous mapping(ICM)method,the model of fail-safe topology optimization is established with the objective of minimal weight integrating with the discrete condition of topological variables and the constraint of the fundamental frequency.The fail-safe optimization model established above is substituted by a sequence of subproblems in the form of the quadratic program with exact second-order information and solved efficiently by the dual sequence quadratic programming(DSQP)algorithm.The numerical result reveals that the optimized fail-safe structure has more complex configuration and preserved materials than the structure obtained from the traditional frequency topology optimization,which means that the optimized fail-safe structure has higher redundancy.Moreover,the optimized fail-safe structure guarantees that the natural frequency meets the constraint of fundamental frequency when the local failure ocurs,which can avoid the structural frequency to be sensitive to local failure.The fail-safe optimirzation topology model is proved effective and feasible by four numerical examples.
文摘To develop an effective numerical method for the cable sliding problem in cable structures, two-node catenary cable element was built to model the cables based on analytical solution of elastic catenary. Cooperated with Newton method, continuation method was used to solve the nonlinear equations. This approach is more efficient than using Newton method only and has a wider range to select initial values for the process to converge. The relationship between the tension on a cable segment and its unstrained length was derived and used to calculate the unbalanced cable tensions at the supports. An example is presented to show the correctness and efficiency of the proposed method.
文摘This paper presents a method for tracing a planar implicit curve f（x, y）=0 on a rectangular region based on continuation scheme. First, according to the starting track-point and the starting track-direction of the curve, make a new fimction F（x, y）=0 where the same curve withf（x, y）=0 is defined. Then we trace the curve between the two domains where F（x, y）〉0 and F（x, y）〈0 alternately, according to the two rules presented in this paper. Equal step size or adaptive step size can be used, when we trace the curve. An irregular planar implicit curve （such as the curve with large curvatures at some points on the curve）, can be plotted if an adaptive step size is used. Moreover, this paper presents a scheme to search for the multiple points on the curve. Our method has the following advantages： （1） it can plot Co planar implicit curves; （2） it can plot the planar implicit curves with multiple points; （3） by the help of using the two rules, our method does not need to compute the tangent vector at the points on the curve, and directly searches for the direction of the tracing curve; （4） the tracing procedure costs only one of two evaluations of function f（x, y）=0 per moving step, while most existing similar methods cost more evaluations of the function.
