Trigonometric Regularization and Continuation Method Based Time-Optimal Control of Hypersonic Vehicles

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 APPLIED IN KINEMATICS OF PARALLEL ROBOT

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 GLOBAL LINEAR AND LOCAL QUADRATIC SINGLE-STEP NONINTERIOR CONTINUATION METHOD FOR MONOTONE SEMIDEFINITE COMPLEMENTARITY PROBLEMS

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.

INVESTIGATION OF FLIGHT DYNAMICS OF THRUST VECTORING AIRCRAFT USING EXTENDED CONTINUATION METHODS

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.

Predicting wing-pylon-nacelle configuration flutter characteristics using adaptive continuation method

In this study,the aeroelastic response of a wing-pylon-nacelle system in subsonic and low supersonic flow regimes is analyzed using the continuation method in conjunction with an adaptive step size control algorithm.Idealizing the pylon and nacelle as a point mass,the computed effects of a standard structural analysis of the wing together with the pylon and nacelle are compared with those of a clean wing to build a reduced-order model for analysis.The aerodynamic forces relating to different reduced frequencies are assessed using the Doublet Lattice Method(DLM)in the subsonic flow regime and supersonic lifting surface theory relying on the unsteady linearized small-disturbance potential flow model in the low supersonic flow regime.The Rational Function Approximation(RFA)method is then utilized for the state-space formulation of the system equations,appended with the continuation method for flutter prediction.Thereafter,the linearized aeroelastic equations are resolved using the continuation method with adaptive step size,the results of which are matched with those obtained from the traditional p-k method to emphasize that the continuation method exhibits a distinct advantage in achieving better accuracy in estimating the flutter speed and identifying the“mode switching”phenomenon.

The non-interior continuation methods for solving the function nonlinear complementarity problem

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.

A Combined Feasible-Nonfeasible Point Continuation Method for Strongly Monotone Variational Inequalities

In this paper, we discuss the variational inequality problems VIP(X, F), where Fis a strongly monotone function and the convex feasible set X is described by some inequality eonstraints. We present a continuation method for VIP(X. F). which solves a sequence ofperturbed variational inequality problems PVIP(X. F, ε. μ) depending on two parameters ε≥ 0and μ＞0. It is worthy to point out that the method will be a feasible point type whenε = 0 and a nonfeasible point type when ε＞0, i.e., it is a combined feasible-nonfeasible point(CFNFP for short) method. We analyse the existence, uniqueness and continuity of the solutionto PVIP(X, F, ε,μ), and prove that any sequence generated by this method converges to theunique solution of VIP(X, F).

Solving a Class of Nonlinear Programming Problems via a Homotopy Continuation Method

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.

Downward continuation of airborne geomagnetic data based on two iterative regularization methods in the frequency domain

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.

An efficient method for tracing planar implicit curves

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.

Existence and Uniqueness of 2π-Periodic Solution About Duffing Equation and Its Numerical Method

A new provement of the existence and uniqueness about periodic boundary value Duffing equation is established by using global inverse function theorem. An algorithm for solving differential equation that has a large convergence domain is given. Finally, a numerical example is given.

Establishment of antenna membrane states FEM on analogue method

On the basis of choosing the basic element as the bar and choosing the basic mesh as the triangle as well as supposing the conditions of the element, the membrane states of an antenna reflector were researched by the analogue method, because the membrane effect was not omitted during the ending deployment process of the radial rib antenna. The expressions of the bar element’s section area and density were obtained, while the expression of the stress state during the ending deployment process of antenna was attained. During the establishment process of the analogue method, the analysis method of the net shell structure was employed. Moreover, during the backward deduction of membrane stress, the continuation method was adopted. Because the expression of the membrane stress state can realize the analysis on the antenna membrane state, this research has great significance of theoretical direction to the normal operation of the space deployable antenna.

Fatigue Topology Optimization Design Based on Distortion Energy Theory and Independent Continuous Mapping Method

Fatigue failure is a common failure mode under the action of cyclic loads in engineering applications,which often occurs with no obvious signal.The maximum structural stress is far below the allowable stress when the structures are damaged.Aiming at the lightweight structure,fatigue topology optimization design is investigated to avoid the occurrence of fatigue failure in the structural conceptual design beforehand.Firstly,the fatigue life is expressed by topology variables and the fatigue life filter function.The continuum fatigue optimization model is established with the independent continuous mapping(ICM)method.Secondly,fatigue life constraints are transformed to distortion energy constraints explicitly by taking advantage of the distortion energy theory.Thirdly,the optimization formulation is solved by the dual sequence quadratic programming(DSQP).And the design scheme of lightweight structure considering the fatigue characteristics is obtained.Finally,numerical examples illustrate the practicality and effectiveness of the fatigue optimization method.This method further expands the theoretical application of the ICM method and provides a novel approach for the fatigue optimization problem.

Design optimization of unsteady airfoils with continuous adjoint method

In this paper, a new unsteady aerodynamic design method is presented based on the Navier-Stokes equations and a continuous adjoint approach. A basic framework of time-accurate unsteady airfoil optimization which adopts time-averaged aerodynamic coefficients as objective functions is presented. The time-accurate continuous adjoint equation and its boundary conditions are derived. The flow field and the adjoint equation are simulated numerically by the finite volume method （FVM）. Feasibility and accuracy of the approach are perfectly validated by the design optimization results of the plunging NACA0012 airfoil.

A class of twostep continuity Runge-Kutta methods for solving singular delay differential equations and its convergence

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.

Continuity control method of cutter posture vector for efficient five-axis machining

During five-axis machining of impeller, the excessive local interference avoidance leads to inconsistency of cutter posture, low quality of machined surface and increase of processing time. Therefore, in order to improve the efficiency of five-axis machining of impellers, it is necessary to minimize the cutter posture changes and create a continuous tool path while avoiding interference. By using an MC-space algorithm for interference avoidance, an MB-spline algorithm for continuous control was intended to create a five-axis machining tool path with excellent surface quality and economic feasibility. A five-axis cutting experiment was performed to verify the effectiveness of the continuity control. The result shows that the surface shape with continuous method is greatly improved, and the surface roughness is generally favorable. Consequently, the effectiveness of the suggested method is verified by identifying the improvement of efficiency of five-axis machining of an impeller in aspects of surface quality and machining time.

ICM Method Combined with Meshfree Approximation for Continuum Structure

The independent continuous mapping（ICM） method is integrated into element free Galerkin method and a new implementation of topology optimization for continuum structure is presented.To facilitate the enforcement of the essential boundary condition and derivative of various sensitivities,a singular weight function in element free Galerkin method is introduced.Material point variable is defined to illustrate the condition of material point and its vicinity instead of element or node.The topological variables field is constructed by moving least square approximation which inherits the continuity and smoothness of the weight function.Due to reciprocal relationships between the topological variables and design variables,various structural responses sensitivities are derived according to the method for calculating the partial derivatives of compound functions.Numerical examples indicate that checkerboard pattern and mesh-dependence phenomena are overcome without additional restriction methods.

EFFICIENT ANALYSIS OF ELECTROMAGNETIC RADIATION FROM MICROSTRIP ANTENNA USING VOLUME-SURFACE CURRENT CONTUNUITY METHOD

The Volume-Surface Current Continuity Method (VSCCM) is presented to analyze electromagnetic radiation from microstrip antenna. The microstrip antenna is discretized into small triangular patches on conducting surface and tetrahedral volume cells in dielectric region. The Method of Moments (MoM) is applied to solve the integral equation. An equation contains the restriction relation between the volume and surface current coefficient is derived from the current continuity equation at those parts where the conducting surface is in contact with the dielectric material. A simple equivalent strip model is introduced in the treatment of the feeding probe in VSCCM. The VSCCM can reduce the unknowns required to be solved in MoM, as well as the condition number of the matrix equation. Numerical results are given to validate the accuracy and efficiency of this method.

Continuous finite element methods for Hamiltonian systems

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.

Stepwise Method Based on Confidence Bound and Information Incorporation for Identifying the Maximum Tolerable Dose

The primary goal of a phase I clinical trial is to find the maximum tolerable dose of a treatment. In this paper, we propose a new stepwise method based on confidence bound and information incorporation to determine the maximum tolerable dose among given dose levels. On the one hand, in order to avoid severe even fatal toxicity to occur and reduce the experimental subjects, the new method is executed from the lowest dose level, and then goes on in a stepwise fashion. On the other hand, in order to improve the accuracy of the recommendation, the final recommendation of the maximum tolerable dose is accomplished through the information incorporation of an additional experimental cohort at the same dose level. Furthermore, empirical simulation results show that the new method has some real advantages in comparison with the modified continual reassessment method.