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.

Fourier Continuation Discontinuous Galerkin Methods for Linear Hyperbolic Problems

Fourier continuation(FC)is an approach used to create periodic extensions of non-periodic functions to obtain highly-accurate Fourier expansions.These methods have been used in partial differential equation(PDE)-solvers and have demonstrated high-order convergence and spectrally accurate dispersion relations in numerical experiments.Discontinuous Galerkin(DG)methods are increasingly used for solving PDEs and,as all Galerkin formulations,come with a strong framework for proving the stability and the convergence.Here we propose the use of FC in forming a new basis for the DG framework.

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.

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.

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.

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.

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.

Fatigue Safety Assessment of Concrete Continuous Rigid Frame Bridge Based on Rain Flow Counting Method and Health Monitoring Data

The fatigue of concrete structures will gradually appear after being subjected to alternating loads for a long time,and the accidents caused by fatigue failure of bridge structures also appear from time to time.Aiming at the problem of degradation of long-span continuous rigid frame bridges due to fatigue and environmental effects,this paper suggests a method to analyze the fatigue degradation mechanism of this type of bridge,which combines long-term in-site monitoring data collected by the health monitoring system(HMS)and fatigue theory.In the paper,the authors mainly carry out the research work in the following aspects:First of all,a long-span continuous rigid frame bridge installed with HMS is used as an example,and a large amount of health monitoring data have been acquired,which can provide efficient information for fatigue in terms of equivalent stress range and cumulative number of stress cycles;next,for calculating the cumulative fatigue damage of the bridge structure,fatigue stress spectrum got by rain flow counting method,S-N curves and damage criteria are used for fatigue damage analysis.Moreover,it was considered a linear accumulation damage through the Palmgren-Miner rule for the counting of stress cycles.The health monitoring data are adopted to obtain fatigue stress data and the rain flow counting method is used to count the amplitude varying fatigue stress.The proposed fatigue reliability approach in the paper can estimate the fatigue damage degree and its evolution law of bridge structures well,and also can help bridge engineers do the assessment of future service duration.

A Differential Continuation－Regularization Method for Solving Inverse Problems of Acoustic Wave Equations

ＡＤｉｆｆｅｒｅｎｔｉａｌＣｏｎｔｉｎｕａｔｉｏｎ－ＲｅｇｕｌａｒｉｚａｔｉｏｎＭｅｔｈｏｄｆｏｒＳｏｌｖｉｎｇＩｎｖｅｒｓｅＰｒｏｂｌｅｍｓｏｆＡｃｏｕｓｔｉｃＷａｖｅＥｑｕａｔｉｏｎｓ￥（张大力）（韩波）（刘家琦）（姜立功）ＺＨＡＮＧＤａｌｉ；Ｈ...

Solving Load Flow Problems of Power System by Explicit Pseudo-Transient Continuation (E-ψtc) Method

This paper is a further study of two papers [1] and [2], which were related to Ill-Conditioned Load Flow Problems and were published by IEEE Trans. PAS. The authors of this paper have some different opinions, for example, the 11-bus system is not an ill-conditioned system. In addition, a new approach to solve Load Flow Problems, E-ψtc, is introduced. It is an explicit method;solving linear equations is not needed. It can handle very tough and very large systems. The advantage of this method has been fully proved by two examples. The authors give this new method a detailed description of how to use it to solve Load Flow Problems and successfully apply it to the 43-bus and the 11-bus systems. The authors also propose a strategy to test the reliability, and by solving gradient equations, this new method can answer if the solution exists or not.

A Two-step Inverse Method for Quantitative Damage Evaluation of Plate-like Structures Using Vibration Approach

This study presents a novel two-step approach to assess plate-like structural laminar damages,particularly for delamination damage detection of composite structures.Firstly,a 2-D continuous wavelet transform is employed to identify the damage location and sizes from vibration curvature data.An inverse method is subsequently then used to determine the bending stiffness reduction ratio along a specified direction,enabling the quantification of the delamination severity.The method employed in this study is an extension of the one-dimensional inverse method developed in a previous work of the authors.The applicability of the two-step inverse approach is demonstrated in a simulation analysis and by an experimental study on a cantilever composite plate containing a single delamination.The inverse method is shown to have the capacity to reveal the detailed damage information of delamination within a constrained searching space and can be used to determine the effective flexural stiffness of composite plate structures,even in cases of complex delamination damage.

Analysis of the Causes of Continuous Cropping Obstacles for Atractylodes macrocephala Koidz in Pingjiang County and Its Control Methods

First at all, it introduced the concept and the damages of continuous cropping obstacle. Then, it analyzed the causes of continuous cropping obstacles for Atractylodes macrocephala Koidz. In the end, in order to provide guidance for pro- moting sustainable development of Atractylodes macrocephala Koidz industry in Pingjiang County, it put forward some control methods for eliminating continuous cropping obstacles of Atractylodes macrocephala Koidz, including breeding varieties with high resistance; applying rotation cropping and intercropping reasonable; rational fertilization and soil disinfection; introducing antagonistic bacterial and eliminating au- tointoxication.

ICM method for topology optimization of multimaterial continuum structure with displacement constraint

A new topology optimization method is formulated for lightweight design of multimaterial structures, using the independent continuous mapping (ICM) method to minimize the weight with a prescribed nodal displacement constraint. Two types of independent topological variable are used to identify the presence of elements and select the material for each phase, to realize the interpolations of the element stiffness matrix and total weight. Furthermore, an explicit expression for the optimized formulation is derived, using approximations of the displacement and weight given by first- and second-order Taylor expansions. The optimization problem is thereby transformed into a standard quadratic programming problem that can be solved using a sequential quadratic programming approach. The feasibility and effectiveness of the proposed multimaterial topology optimization method are demonstrated by determining the best load transfer path for four numerical examples. The results reveal that the topologically optimized configuration of the multimaterial structure varies with the material properties, load conditions, and constraint. Firstly, the weight of the optimized multimaterial structure is found to be lower than that composed of a single material. Secondly, under the precondition of a displacement constraint, the weight of the topologically optimized multimaterial structure decreases as the displacement constraint value is increased. Finally, the topologically optimized multimaterial structures differ depending on the elastic modulus of the materials. Besides, the established optimization formulation is more reliable and suitable for use in practical engineering applications with structural performance parameters as constraint.

Energy method and numerical simulation of critical backfillheight in non-pillar continuous mining

Non-pillar continuous mining(NPCM) is regarded as a high-efficient, high-level and one-step mining technology, which can be divided into two substopes. Back fill stability status in substope I, which directly influence the loss rate and dilution rate, etc, will determine whether the experimental research is successful or not. By employing energy method of limit analysis and finite element numerical simulation method, the critical backfill height was determined under the prerequisite condition of its stability, which put forward theoretical basis for reasonable and correct selection of backfill’s parameters. The result showed that the first backfill could not keep stable for NPCM, while the other was able to.

Longitudinal crack on slab surface at straightening stage during continuous casting using finite element method

Deformation behavior of slab at the straightening stage during continuous casting was simulated by the explicit dynamic finite element method,and the stress distribution along the width direction of the slab and its change regularity at slab center during continuous casting were obtained.The influence of distribution and change of stress on the propagation of longitudinal cracks on slab surface was discussed.The results show that the tensional stress appears on slab surface at the inner arc side and the compressive stress appears on slab surface at the outer arc side at stages 6-8 in straightening zone during continuous casting.Longitudinal cracks generally appear on slab top surface and do not appear on slab bottom surface,which are also observed in industry.

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.

Particle path tracking method in two- and three-dimensional continuously rotating detonation engines

The particle path tracking method is proposed and used in two-dimensional（2D） and three-dimensional（3D） numerical simulations of continuously rotating detonation engines（CRDEs）. This method is used to analyze the combustion and expansion processes of the fresh particles, and the thermodynamic cycle process of CRDE. In a 3D CRDE flow field, as the radius of the annulus increases, the no-injection area proportion increases, the non-detonation proportion decreases, and the detonation height decreases. The flow field parameters on the 3D mid annulus are different from in the 2D flow field under the same chamber size. The non-detonation proportion in the 3D flow field is less than in the 2D flow field. In the 2D and 3D CRDE, the paths of the flow particles have only a small fluctuation in the circumferential direction. The numerical thermodynamic cycle processes are qualitatively consistent with the three ideal cycle models, and they are right in between the ideal F–J cycle and ideal ZND cycle. The net mechanical work and thermal efficiency are slightly smaller in the 2D simulation than in the 3D simulation. In the 3D CRDE, as the radius of the annulus increases, the net mechanical work is almost constant, and the thermal efficiency increases. The numerical thermal efficiencies are larger than F–J cycle, and much smaller than ZND cycle.

Ant colony algorithm based on genetic method for continuous optimization problem

A new algorithm is presented by using the ant colony algorithm based on genetic method （ACG） to solve the continuous optimization problem. Each component has a seed set. The seed in the set has the value of component, trail information and fitness. The ant chooses a seed from the seed set with the possibility determined by trail information and fitness of the seed. The genetic method is used to form new solutions from the solutions got by the ants. Best solutions are selected to update the seeds in the sets and trail information of the seeds. In updating the trail information, a diffusion function is used to achieve the diffuseness of trail information. The new algorithm is tested with 8 different benchmark functions.

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.

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.