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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
As a basic mathematical structure,the system of inequalities over symmetric cones and its solution can provide an effective method for solving the startup problem of interior point method which is used to solve many o...As a basic mathematical structure,the system of inequalities over symmetric cones and its solution can provide an effective method for solving the startup problem of interior point method which is used to solve many optimization problems.In this paper,a non-interior continuation algorithm is proposed for solving the system of inequalities under the order induced by a symmetric cone.It is shown that the proposed algorithm is globally convergent and well-defined.Moreover,it can start from any point and only needs to solve one system of linear equations at most at each iteration.Under suitable assumptions,global linear and local quadratic convergence is established with Euclidean Jordan algebras.Numerical results indicate that the algorithm is efficient.The systems of random linear inequalities were tested over the second-order cones with sizes of 10,100,,1 000 respectively and the problems of each size were generated randomly for 10 times.The average iterative numbers show that the proposed algorithm can generate a solution at one step for solving the given linear class of problems with random initializations.It seems possible that the continuation algorithm can solve larger scale systems of linear inequalities over the secondorder cones quickly.Moreover,a system of nonlinear inequalities was also tested over Cartesian product of two simple second-order cones,and numerical results indicate that the proposed algorithm can deal with the nonlinear cases.展开更多
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 exam...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.展开更多
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)-solve...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.展开更多
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 provi...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.展开更多
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 co...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.展开更多
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.展开更多
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 influenc...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.展开更多
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 c...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.展开更多
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 adopt...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.展开更多
On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these meth...On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming （bi）linear macroelements or （bi）quadratic elements, and the rotation by conforming （bi）linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.展开更多
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 th...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.展开更多
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 componen...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.展开更多
基金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 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.
文摘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 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.
基金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.
基金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.
文摘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.
基金Supported by National Natural Science Foundation of China (No.10871144)the Seed Foundation of Tianjin University (No.60302023)
文摘As a basic mathematical structure,the system of inequalities over symmetric cones and its solution can provide an effective method for solving the startup problem of interior point method which is used to solve many optimization problems.In this paper,a non-interior continuation algorithm is proposed for solving the system of inequalities under the order induced by a symmetric cone.It is shown that the proposed algorithm is globally convergent and well-defined.Moreover,it can start from any point and only needs to solve one system of linear equations at most at each iteration.Under suitable assumptions,global linear and local quadratic convergence is established with Euclidean Jordan algebras.Numerical results indicate that the algorithm is efficient.The systems of random linear inequalities were tested over the second-order cones with sizes of 10,100,,1 000 respectively and the problems of each size were generated randomly for 10 times.The average iterative numbers show that the proposed algorithm can generate a solution at one step for solving the given linear class of problems with random initializations.It seems possible that the continuation algorithm can solve larger scale systems of linear inequalities over the secondorder cones quickly.Moreover,a system of nonlinear inequalities was also tested over Cartesian product of two simple second-order cones,and numerical results indicate that the proposed algorithm can deal with the nonlinear cases.
文摘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.
文摘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.
文摘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.
基金the National Natural Science Foundation of China (Grants 11072009 and 11872080)Beijing Education Committee Development Project (Grant SQKM201610005001).
文摘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.
基金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.
文摘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.
基金Project(50634030) supported by the National Natural Science Foundation of ChinaProject(20090042120005) supported by the Doctorate Foundation of the Ministry of Education of ChinaProject(2006CB605208-1) supported by the State Basic Research Program of China
文摘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.
基金supported by the Shanghai Municipal Natural Science Foundation(No.13ZR1415700)
文摘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.
基金supported by NSFC(11571266,91430106,11171168,11071132)NSFC-RGC(China-Hong Kong)(11661161017)
文摘On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming （bi）linear macroelements or （bi）quadratic elements, and the rotation by conforming （bi）linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.
文摘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.
基金project supported by the National High-Technology Research and Development Program of China(Grant No.8632005AA642010)
文摘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.
