The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And...The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.展开更多
A novel control method for a general class of nonlinear systems using fuzzy logic systems （FLSs） is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a ...A novel control method for a general class of nonlinear systems using fuzzy logic systems （FLSs） is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a high-gain observer is used to estimate the derivatives of the system output. The closed-loop system is proven to be semiglobally uniformly ultimately bounded. In addition, it is shown that if the approximation accuracy of the fuzzy logic system is high enough and the observer gain is chosen sufficiently large, an arbitrarily small tracking error can be achieved. Simulation results verify the effectiveness of the newly designed scheme and the theoretical discussion.展开更多
A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, th...A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, the search space of free parameters has been decreased. Then, in order to obtain the global optimal solution effectively and robustly, the simulated annealing and penalty function strategies were used to handle the constraints, and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem, in which, a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement. Comparing to the classical direct method, this novel method has fewer free parameters, needs not initial guesses, and has higher computation precision. An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach. The results of simulation indicate that our approach is available to solve the problem of optimal muhi-revolution transfer between Earth-orbits.展开更多
In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices ou...In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices outside of the active set, while the projected gradient method is used to update the active variables. At each iterative level, the search direction consists of three parts, one of which is a subspace truncated Newton direction, the other two are subspace gradient and modified gradient directions. The subspace truncated Newton direction is obtained by solving a sparse system of linear equations. The global convergence and quadratic convergence rate of the algorithm are proved and some numerical tests are given.展开更多
This paper presents a trust region two phase model algorithm for solving the equality and bound constrained nonlinear optimization problem. A concept of substationary point is given. Under suitable assumptions,the gl...This paper presents a trust region two phase model algorithm for solving the equality and bound constrained nonlinear optimization problem. A concept of substationary point is given. Under suitable assumptions,the global convergence of this algorithm is proved without assuming the linear independence of the gradient of active constraints. A numerical example is also presented.展开更多
The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolat...The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolation is not only reliable but also can save the amount of calculation by nearly 36%. Large amount of calculation and lacking strict theoretical bas is has been th e two disadvantage of direct method by new. If this defect is not overcome, they will not only s eriously affect the application of this method, but also hinder its further rese arch. Based on sufficient calculation practice, this article has made a primary discussion about the theory and method of reducing the amount of calculation, an d has achieved some satisfactory results.展开更多
Current topology optimization methods for nonlinear continuum structures often suffer from low computational efficiency and limited applicability to complex nonlinear problems.To address these issues,this paper propos...Current topology optimization methods for nonlinear continuum structures often suffer from low computational efficiency and limited applicability to complex nonlinear problems.To address these issues,this paper proposes an improved bi-directional evolutionary structural optimization(BESO)method tailored for maximizing stiffness in nonlinear structures.The optimization program is developed in Python and can be combined with Abaqus software to facilitate finite element analysis(FEA).To accelerate the speed of optimization,a novel adaptive evolutionary ratio(ER)strategy based on the BESO method is introduced,with four distinct adaptive ER functions proposed.The Newton-Raphson method is utilized for iteratively solving nonlinear equilibrium equations,and the sensitivity information for updating design variables is derived using the adjoint method.Additionally,this study extends topology optimization to account for both material nonlinearity and geometric nonlinearity,analyzing the effects of various nonlinearities.A series of comparative studies are conducted using benchmark cases to validate the effectiveness of the proposed method.The results show that the BESO method with adaptive ER significantly improves the optimization efficiency.Compared to the BESO method with a fixed ER,the convergence speed of the four adaptive ER BESO methods is increased by 37.3%,26.7%,12%and 18.7%,respectively.Given that Abaqus is a powerful FEA platform,this method has the potential to be extended to large-scale engineering structures and to address more complex optimization problems.This research proposes an improved BESO method with novel adaptive ER,which significantly accelerates the optimization process and enables its application to topology optimization of nonlinear structures.展开更多
A new two-level subspace method is proposed for solving the general unconstrained minimization formulations discretized from infinite-dimensional optimization problems. At each iteration, the algorithm executes either...A new two-level subspace method is proposed for solving the general unconstrained minimization formulations discretized from infinite-dimensional optimization problems. At each iteration, the algorithm executes either a direct step on the current level or a coarse subspace correction step. In the coarse subspace correction step, we augment the traditional coarse grid space by a two-dimensional subspace spanned by the coordinate direction and the gradient direction at the current point. Global convergence is proved and convergence rate is studied under some mild conditions on the discretized functions. Preliminary numerical experiments on a few variational problems show that our two-level subspace method is promising.展开更多
A numerical embedding method was proposed for solving the nonlinear optimization problem. By using the nonsmooth theory, the existence and the continuation of the following path for the corresponding homotopy equation...A numerical embedding method was proposed for solving the nonlinear optimization problem. By using the nonsmooth theory, the existence and the continuation of the following path for the corresponding homotopy equations were proved. Therefore the basic theory for the algorithm of the numerical embedding method for solving the non-linear optimization problem was established. Based on the theoretical results, a numerical embedding algorithm was designed for solving the nonlinear optimization problem, and prove its convergence carefully. Numerical experiments show that the algorithm is effective.展开更多
A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structur...A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structural strength in engineering applications.First,a topology optimization model is established for a lightweight structure with element stress as constraints.Second,the stress globalization method is adopted to convert local stress constraints into strain energy constraints,which overcomes the difficulties caused by local stress constraints,such as model establishment,sensitivity analysis,and massive solution calculations.Third,the sensitivity of the objective function and constraint function is analyzed,and the method of moving asymptotes is employed to solve the optimization model.In addition,the additive hyperelasticity technique is utilized to solve the numerical instability induced by structures undergoing large deformation.Numerical examples are given to validate the feasibility of the proposed method.The method provides a significant reference for geometrically nonlinear optimization design.展开更多
Trust region methods are powerful and effective optimization methods. The conic model method is a new type of method with more information available at each iteration than standard quadratic-based methods. The adva...Trust region methods are powerful and effective optimization methods. The conic model method is a new type of method with more information available at each iteration than standard quadratic-based methods. The advantages of the above two methods can be combined to form a more powerful method for constrained optimization. The trust region subproblem of our method is to minimize a conic function subject to the linearized constraints and trust region bound. At the same time, the new algorithm still possesses robust global properties. The global convergence of the new algorithm under standard conditions is established.展开更多
In this paper, a distributed algorithm is proposed to solve a kind of multi-objective optimization problem based on the alternating direction method of multipliers. Compared with the centralized algorithms, this algor...In this paper, a distributed algorithm is proposed to solve a kind of multi-objective optimization problem based on the alternating direction method of multipliers. Compared with the centralized algorithms, this algorithm does not need a central node. Therefore, it has the characteristics of low communication burden and high privacy. In addition, numerical experiments are provided to validate the effectiveness of the proposed algorithm.展开更多
Weight reduction has attracted much attention among ship designers and ship owners.In the present work,based on an improved bi-directional evolutionary structural optimization(BESO) method and surrogate model method,w...Weight reduction has attracted much attention among ship designers and ship owners.In the present work,based on an improved bi-directional evolutionary structural optimization(BESO) method and surrogate model method,we propose a hybrid optimization method for the structural design optimization of beam-plate structures,which covers three optimization levels:dimension optimization,topology optimization and section optimization.The objective of the proposed optimization method is to minimize the weight of design object under a group of constraints.The kernel optimization procedure(KOP) uses BESO to obtain the optimal topology from a ground structure.To deal with beam-plate structures,the traditional BESO method is improved by using cubic box as the unit cell instead of solid unit to construct periodic lattice structure.In the first optimization level,a series of ground structures are generated based on different dimensional parameter combinations,the KOP is performed to all the ground structures,the response surface model of optimal objective values and dimension parameters is created,and then the optimal dimension parameters can be obtained.In the second optimization level,the optimal topology is obtained by using the KOP according to the optimal dimension parameters.In the third optimization level,response surface method(RSM) is used to determine the section parameters.The proposed method is applied to a hatch cover structure design.The locations and shapes of all the structural members are determined from an oversized ground structure.The results show that the proposed method leads to a greater weight saving,compared with the original design and genetic algorithm(GA) based optimization results.展开更多
Orthogonal conditional nonlinear optimal perturbations(O-CNOPs)have been used to generate ensemble forecasting members for achieving high forecasting skill of high-impact weather and climate events.However,highly effi...Orthogonal conditional nonlinear optimal perturbations(O-CNOPs)have been used to generate ensemble forecasting members for achieving high forecasting skill of high-impact weather and climate events.However,highly efficient calculations for O-CNOPs are still challenging in the field of ensemble forecasting.In this study,we combine a gradient-based iterative idea with the Gram‒Schmidt orthogonalization,and propose an iterative optimization method to compute O-CNOPs.This method is different from the original sequential optimization method,and allows parallel computations of O-CNOPs,thus saving a large amount of computational time.We evaluate this method by using the Lorenz-96 model on the basis of the ensemble forecasting ability achieved and on the time consumed for computing O-CNOPs.The results demonstrate that the parallel iterative method causes O-CNOPs to yield reliable ensemble members and to achieve ensemble forecasting skills similar to or even slightly higher than those produced by the sequential method.Moreover,the parallel method significantly reduces the computational time for O-CNOPs.Therefore,the parallel iterative method provides a highly effective and efficient approach for calculating O-CNOPs for ensemble forecasts.Expectedly,it can play an important role in the application of the O-CNOPs to realistic ensemble forecasts for high-impact weather and climate events.展开更多
This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element couplin...This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.展开更多
In this work, we present the direct discontinuous Galerkin （DDG） method for the one-dimensional coupled nonlinear Schr5dinger （CNLS） equation. We prove that the new discontinuous Galerkin method preserves the disc...In this work, we present the direct discontinuous Galerkin （DDG） method for the one-dimensional coupled nonlinear Schr5dinger （CNLS） equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.展开更多
Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock ma...Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock mass,the implementation of this technology often encounters design challenges,leading to suboptimal results and increased costs.This paper establishes a structural analysis model of the goaf working face roof,revealing the failure mechanism of DRC,and clarifies the positive role of DRC in improving the stress of the roadway surrounding rock and reducing the subsidence of the roof through numerical simulation experiments.On this basis,the paper further analyses the roadway pressure and roof settlement under different DRC design heights,and ultimately proposes an optimized design method for the DRC height.The results indicate that the implementation of DRC can significantly optimize the stress environment of the working face roadway surrounding rock.At the same time,during the application of DRC,three scenarios may arise:insufficient,reasonable,and excessive DRC height.Insufficient height will significantly reduce the effectiveness of the technology,while excessive height has little impact on the implementation effect but will greatly increase construction costs and difficulty.Engineering verification shows that the optimized DRC design method proposed in this paper reduces the peak stress of the protective coal pillar in the roadway by 27.2%and the central subsidence of the roof by 41.8%,demonstrating excellent application results.This method provides technical support for the further promotion of NCMSE mining method.展开更多
We extend Lou＇s direct perturbation method for solving the nonlinear Schrodinger equation to the case of the derivative nonlinear Schrodinger equation （DNLSE）. By applying this method, different types of perturbati...We extend Lou＇s direct perturbation method for solving the nonlinear Schrodinger equation to the case of the derivative nonlinear Schrodinger equation （DNLSE）. By applying this method, different types of perturbation solutions are obtained. Based on these approximate solutions, the analytical forms of soliton parameters, such as the velocity, the width and the initial position, are carried out and the effects of perturbation on solitons are analyzed at the same time. A numerical simulation of perturbed DNLSE finally verifies the results of the perturbation method.展开更多
In this work,we present the direct discontinuous Galerkin(DDG) method for the one-dimensional coupled nonlinear Schrdinger(CNLS) equation.We prove that the new discontinuous Galerkin method preserves the discrete mass...In this work,we present the direct discontinuous Galerkin(DDG) method for the one-dimensional coupled nonlinear Schrdinger(CNLS) equation.We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system.The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method.Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.展开更多
Low-thrust Earth-orbit transfers with 10^- 5-order thrust-to-weight ratios involve a large number of orbital revolutions which poses a real challenge to trajectory optimization. This article develops a direct method t...Low-thrust Earth-orbit transfers with 10^- 5-order thrust-to-weight ratios involve a large number of orbital revolutions which poses a real challenge to trajectory optimization. This article develops a direct method to optimize minimum-time low-thrust many-revolution Earth-orbit transfers. A parameterized control law in each orbit, in the form of the true optimal control, is proposed, and the time history of the parameters governing the control law is interpolated through a finite number of nodal values. The orbital averaging method is used to significantly reduce the computational workload and the trajectory optimization is conducted based on the orbital averaging dynamics expressed by nonsingular equinoctial elements. Furthermore, Earth＇s shadowing and perturbation effects are taken into account. The optimal transfer problem is thus converted to the parameter optimization problem that can be solved by nonlinear programming. Taking advantage of the mapping between the parameterized control law and the Lyapunov control law, a technique is proposed to acquire good initial guesses for optimization variables, which results in enlarged convergence domain of the direct optimization method. Numerical examples of optimal Earth-orbit transfers are presented.展开更多
文摘The principle of direct method used in optimal control problem is introduced. Details of applying this method to flight trajectory generation are presented including calculation of velocity and controls histories. And capabilities of flight and propulsion systems are considered also. Combined with digital terrain map technique, the direct method is applied to the three dimensional trajectory optimization for low altitude penetration, and simplex algorithm is used to solve the parameters in optimization. For the small number of parameters, the trajectory can be optimized in real time on board.
基金This project was supported by the National Natural Science Foundation of China (90405011).
文摘A novel control method for a general class of nonlinear systems using fuzzy logic systems （FLSs） is presertted. Indirect and direct methods are combined to design the adaptive fuzzy output feedback controller and a high-gain observer is used to estimate the derivatives of the system output. The closed-loop system is proven to be semiglobally uniformly ultimately bounded. In addition, it is shown that if the approximation accuracy of the fuzzy logic system is high enough and the observer gain is chosen sufficiently large, an arbitrarily small tracking error can be achieved. Simulation results verify the effectiveness of the newly designed scheme and the theoretical discussion.
基金Sponsored by the National Natural Science Foundation of China(Grant No.10672044)
文摘A modified direct optimization method is proposed to solve the optimal multi-revolution transfer with low-thrust between Earth-orbits. First, through parameterizing the control steering angles by costate variables, the search space of free parameters has been decreased. Then, in order to obtain the global optimal solution effectively and robustly, the simulated annealing and penalty function strategies were used to handle the constraints, and a GA/SQP hybrid optimization algorithm was utilized to solve the parameter optimization problem, in which, a feasible suboptimal solution obtained by GA was submitted as an initial parameter set to SQP for refinement. Comparing to the classical direct method, this novel method has fewer free parameters, needs not initial guesses, and has higher computation precision. An optimal-fuel transfer problem from LEO to GEO was taken as an example to validate the proposed approach. The results of simulation indicate that our approach is available to solve the problem of optimal muhi-revolution transfer between Earth-orbits.
基金The research was supported by the State Education Grant for Retumed Scholars
文摘In this paper we report a sparse truncated Newton algorithm for handling large-scale simple bound nonlinear constrained minimixation problem. The truncated Newton method is used to update the variables with indices outside of the active set, while the projected gradient method is used to update the active variables. At each iterative level, the search direction consists of three parts, one of which is a subspace truncated Newton direction, the other two are subspace gradient and modified gradient directions. The subspace truncated Newton direction is obtained by solving a sparse system of linear equations. The global convergence and quadratic convergence rate of the algorithm are proved and some numerical tests are given.
文摘This paper presents a trust region two phase model algorithm for solving the equality and bound constrained nonlinear optimization problem. A concept of substationary point is given. Under suitable assumptions,the global convergence of this algorithm is proved without assuming the linear independence of the gradient of active constraints. A numerical example is also presented.
文摘The general interpolation mentioned in this a rticle provides an effective way for reducing the amount of calculation of direc t optimal exploration. It has been testified by real case calculations that the interpolation is not only reliable but also can save the amount of calculation by nearly 36%. Large amount of calculation and lacking strict theoretical bas is has been th e two disadvantage of direct method by new. If this defect is not overcome, they will not only s eriously affect the application of this method, but also hinder its further rese arch. Based on sufficient calculation practice, this article has made a primary discussion about the theory and method of reducing the amount of calculation, an d has achieved some satisfactory results.
基金Supported by National Natural Science Foundation of China(Grant No.52105271).
文摘Current topology optimization methods for nonlinear continuum structures often suffer from low computational efficiency and limited applicability to complex nonlinear problems.To address these issues,this paper proposes an improved bi-directional evolutionary structural optimization(BESO)method tailored for maximizing stiffness in nonlinear structures.The optimization program is developed in Python and can be combined with Abaqus software to facilitate finite element analysis(FEA).To accelerate the speed of optimization,a novel adaptive evolutionary ratio(ER)strategy based on the BESO method is introduced,with four distinct adaptive ER functions proposed.The Newton-Raphson method is utilized for iteratively solving nonlinear equilibrium equations,and the sensitivity information for updating design variables is derived using the adjoint method.Additionally,this study extends topology optimization to account for both material nonlinearity and geometric nonlinearity,analyzing the effects of various nonlinearities.A series of comparative studies are conducted using benchmark cases to validate the effectiveness of the proposed method.The results show that the BESO method with adaptive ER significantly improves the optimization efficiency.Compared to the BESO method with a fixed ER,the convergence speed of the four adaptive ER BESO methods is increased by 37.3%,26.7%,12%and 18.7%,respectively.Given that Abaqus is a powerful FEA platform,this method has the potential to be extended to large-scale engineering structures and to address more complex optimization problems.This research proposes an improved BESO method with novel adaptive ER,which significantly accelerates the optimization process and enables its application to topology optimization of nonlinear structures.
文摘A new two-level subspace method is proposed for solving the general unconstrained minimization formulations discretized from infinite-dimensional optimization problems. At each iteration, the algorithm executes either a direct step on the current level or a coarse subspace correction step. In the coarse subspace correction step, we augment the traditional coarse grid space by a two-dimensional subspace spanned by the coordinate direction and the gradient direction at the current point. Global convergence is proved and convergence rate is studied under some mild conditions on the discretized functions. Preliminary numerical experiments on a few variational problems show that our two-level subspace method is promising.
文摘A numerical embedding method was proposed for solving the nonlinear optimization problem. By using the nonsmooth theory, the existence and the continuation of the following path for the corresponding homotopy equations were proved. Therefore the basic theory for the algorithm of the numerical embedding method for solving the non-linear optimization problem was established. Based on the theoretical results, a numerical embedding algorithm was designed for solving the nonlinear optimization problem, and prove its convergence carefully. Numerical experiments show that the algorithm is effective.
基金This work was supported by the National Natural Science Foundation of China(11872080)Beijing Natural Science Foundation(3192005)。
文摘A geometrically nonlinear topology optimization method for continuum structures is proposed based on the independent continuous mapping method.The stress constraint problem is studied due to the importance of structural strength in engineering applications.First,a topology optimization model is established for a lightweight structure with element stress as constraints.Second,the stress globalization method is adopted to convert local stress constraints into strain energy constraints,which overcomes the difficulties caused by local stress constraints,such as model establishment,sensitivity analysis,and massive solution calculations.Third,the sensitivity of the objective function and constraint function is analyzed,and the method of moving asymptotes is employed to solve the optimization model.In addition,the additive hyperelasticity technique is utilized to solve the numerical instability induced by structures undergoing large deformation.Numerical examples are given to validate the feasibility of the proposed method.The method provides a significant reference for geometrically nonlinear optimization design.
文摘Trust region methods are powerful and effective optimization methods. The conic model method is a new type of method with more information available at each iteration than standard quadratic-based methods. The advantages of the above two methods can be combined to form a more powerful method for constrained optimization. The trust region subproblem of our method is to minimize a conic function subject to the linearized constraints and trust region bound. At the same time, the new algorithm still possesses robust global properties. The global convergence of the new algorithm under standard conditions is established.
文摘In this paper, a distributed algorithm is proposed to solve a kind of multi-objective optimization problem based on the alternating direction method of multipliers. Compared with the centralized algorithms, this algorithm does not need a central node. Therefore, it has the characteristics of low communication burden and high privacy. In addition, numerical experiments are provided to validate the effectiveness of the proposed algorithm.
基金the National Natural Science Foundation of China(No.51509033)
文摘Weight reduction has attracted much attention among ship designers and ship owners.In the present work,based on an improved bi-directional evolutionary structural optimization(BESO) method and surrogate model method,we propose a hybrid optimization method for the structural design optimization of beam-plate structures,which covers three optimization levels:dimension optimization,topology optimization and section optimization.The objective of the proposed optimization method is to minimize the weight of design object under a group of constraints.The kernel optimization procedure(KOP) uses BESO to obtain the optimal topology from a ground structure.To deal with beam-plate structures,the traditional BESO method is improved by using cubic box as the unit cell instead of solid unit to construct periodic lattice structure.In the first optimization level,a series of ground structures are generated based on different dimensional parameter combinations,the KOP is performed to all the ground structures,the response surface model of optimal objective values and dimension parameters is created,and then the optimal dimension parameters can be obtained.In the second optimization level,the optimal topology is obtained by using the KOP according to the optimal dimension parameters.In the third optimization level,response surface method(RSM) is used to determine the section parameters.The proposed method is applied to a hatch cover structure design.The locations and shapes of all the structural members are determined from an oversized ground structure.The results show that the proposed method leads to a greater weight saving,compared with the original design and genetic algorithm(GA) based optimization results.
基金sponsored by the National Natural Science Foundation of China(Grant Nos.41930971,42330111,and 42405061)the National Key Scientific and Technological Infrastructure project“Earth System Numerical Simulation Facility”(Earth Lab).
文摘Orthogonal conditional nonlinear optimal perturbations(O-CNOPs)have been used to generate ensemble forecasting members for achieving high forecasting skill of high-impact weather and climate events.However,highly efficient calculations for O-CNOPs are still challenging in the field of ensemble forecasting.In this study,we combine a gradient-based iterative idea with the Gram‒Schmidt orthogonalization,and propose an iterative optimization method to compute O-CNOPs.This method is different from the original sequential optimization method,and allows parallel computations of O-CNOPs,thus saving a large amount of computational time.We evaluate this method by using the Lorenz-96 model on the basis of the ensemble forecasting ability achieved and on the time consumed for computing O-CNOPs.The results demonstrate that the parallel iterative method causes O-CNOPs to yield reliable ensemble members and to achieve ensemble forecasting skills similar to or even slightly higher than those produced by the sequential method.Moreover,the parallel method significantly reduces the computational time for O-CNOPs.Therefore,the parallel iterative method provides a highly effective and efficient approach for calculating O-CNOPs for ensemble forecasts.Expectedly,it can play an important role in the application of the O-CNOPs to realistic ensemble forecasts for high-impact weather and climate events.
基金supported by grants from the National Natural Science Foundation of China (51478130)the Guangzhou Municipal Education Bureau’s Scientific Research Project, China (2024312217)+1 种基金the China Scholarship Council (201808440070)the 111 Project of China (D21021).
文摘This paper presents an improved level set method for topology optimization of geometrically nonlinear structures accounting for the effect of thermo-mechanical couplings.It derives a new expression for element coupling stress resulting from the combination of mechanical and thermal loading,using geometric nonlinear finite element analysis.A topological model is then developed to minimize compliance while meeting displacement and frequency constraints to fulfill design requirements of structural members.Since the conventional Lagrange multiplier search method is unable to handle convergence instability arising from large deformation,a novel Lagrange multiplier search method is proposed.Additionally,the proposed method can be extended to multi-constrained geometrically nonlinear topology optimization,accommodating multiple physical field couplings.
基金Project supported by the National Natural Science Foundation of China （Grant No 11171038）.
文摘In this work, we present the direct discontinuous Galerkin （DDG） method for the one-dimensional coupled nonlinear Schr5dinger （CNLS） equation. We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system. The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge Kutta method. Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.
基金funded by the National Natural Science Foundation of China(52074298)Beijing Municipal Natural Science Foundation(8232056)+1 种基金Guizhou Province science and technology plan project([2020]3008)Liulin Energy and Environment Academician Workstation(2022XDHZ12).
文摘Directional roof cutting(DRC)is one of the key techniques in non-pillar coal mining with self-formed entries(NCMSE)mining method.Due to the inability to accurately measure the expansion coefficient of the goaf rock mass,the implementation of this technology often encounters design challenges,leading to suboptimal results and increased costs.This paper establishes a structural analysis model of the goaf working face roof,revealing the failure mechanism of DRC,and clarifies the positive role of DRC in improving the stress of the roadway surrounding rock and reducing the subsidence of the roof through numerical simulation experiments.On this basis,the paper further analyses the roadway pressure and roof settlement under different DRC design heights,and ultimately proposes an optimized design method for the DRC height.The results indicate that the implementation of DRC can significantly optimize the stress environment of the working face roadway surrounding rock.At the same time,during the application of DRC,three scenarios may arise:insufficient,reasonable,and excessive DRC height.Insufficient height will significantly reduce the effectiveness of the technology,while excessive height has little impact on the implementation effect but will greatly increase construction costs and difficulty.Engineering verification shows that the optimized DRC design method proposed in this paper reduces the peak stress of the protective coal pillar in the roadway by 27.2%and the central subsidence of the roof by 41.8%,demonstrating excellent application results.This method provides technical support for the further promotion of NCMSE mining method.
基金supported by National Natural Science Foundation of China under Grant No.10575087
文摘We extend Lou＇s direct perturbation method for solving the nonlinear Schrodinger equation to the case of the derivative nonlinear Schrodinger equation （DNLSE）. By applying this method, different types of perturbation solutions are obtained. Based on these approximate solutions, the analytical forms of soliton parameters, such as the velocity, the width and the initial position, are carried out and the effects of perturbation on solitons are analyzed at the same time. A numerical simulation of perturbed DNLSE finally verifies the results of the perturbation method.
基金Project supported by the National Natural Science Foundation of China (Grant No 11171038)
文摘In this work,we present the direct discontinuous Galerkin(DDG) method for the one-dimensional coupled nonlinear Schrdinger(CNLS) equation.We prove that the new discontinuous Galerkin method preserves the discrete mass conservations corresponding to the properties of the CNLS system.The ordinary differential equations obtained by the DDG space discretization is solved via a third-order stabilized Runge-Kutta method.Numerical experiments show that the new DDG scheme gives stable and less diffusive results and has excellent long-time numerical behaviors for the CNLS equations.
基金National Natural Science Foundation of China (10603005)
文摘Low-thrust Earth-orbit transfers with 10^- 5-order thrust-to-weight ratios involve a large number of orbital revolutions which poses a real challenge to trajectory optimization. This article develops a direct method to optimize minimum-time low-thrust many-revolution Earth-orbit transfers. A parameterized control law in each orbit, in the form of the true optimal control, is proposed, and the time history of the parameters governing the control law is interpolated through a finite number of nodal values. The orbital averaging method is used to significantly reduce the computational workload and the trajectory optimization is conducted based on the orbital averaging dynamics expressed by nonsingular equinoctial elements. Furthermore, Earth＇s shadowing and perturbation effects are taken into account. The optimal transfer problem is thus converted to the parameter optimization problem that can be solved by nonlinear programming. Taking advantage of the mapping between the parameterized control law and the Lyapunov control law, a technique is proposed to acquire good initial guesses for optimization variables, which results in enlarged convergence domain of the direct optimization method. Numerical examples of optimal Earth-orbit transfers are presented.
