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.展开更多
In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the ...In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.展开更多
This study aims to optimize the inbound traffic flow on on-ramps by considering low time costs,good speed stability,and high driving safety for mixed traffic flow.The optimal inlet gap is identified in advance,and tra...This study aims to optimize the inbound traffic flow on on-ramps by considering low time costs,good speed stability,and high driving safety for mixed traffic flow.The optimal inlet gap is identified in advance,and trajectory guidance for vehicles entering the gap is determined under safety constraints.Based on the initial state and sequence of vehicles entering the merging area,individual vehicle trajectories are optimized sequentially.An optimization model and method for ramp entry trajectories in mixed traffic flow are developed,incorporating on-ramp vehicle entry sequencing and ordinary vehicle trajectory prediction.Key performance indicators,including driving safety,total travel time,parking wait probability,and trajectory smoothness,are compared and analyzed to evaluate the proposed approach.展开更多
Response analysis of structures involving non-probabilistic uncertain parameters can be closely related to optimization.This paper provides a review on optimization-based methods for uncertainty analysis,with focusing...Response analysis of structures involving non-probabilistic uncertain parameters can be closely related to optimization.This paper provides a review on optimization-based methods for uncertainty analysis,with focusing attention on specific properties of adopted numerical optimization approaches.We collect and discuss the methods based on nonlinear programming,semidefinite programming,mixed-integer programming,mathematical programming with complementarity constraints,difference-of-convex programming,optimization methods using surrogate models and machine learning techniques,and metaheuristics.As a closely related topic,we also overview the methods for assessing structural robustness using non-probabilistic uncertainty modeling.We conclude the paper by drawing several remarks through this review.展开更多
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.展开更多
Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patte...Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.展开更多
Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. ...Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.展开更多
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.展开更多
Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary a...Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms （EAs） and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.展开更多
In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model f...In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.展开更多
The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimiz...The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimization model for multibody system dynamics with multi-point contact and collision is established.The paper presents the study of the numerical solution scheme,in which particle swarm optimization method is used to deal with the corresponding optimization model.The article also presents the comparison of the Gauss optimization method(GOM)and the hybrid linear complementarity method(i.e.combining differential algebraic equations(DAEs)and linear complementarity problems(LCP)),commonly used to solve the dynamic contact problem of multibody systems with bilateral constraints.The results illustrate that,the GOM has the same advantage of dynamical modelling with LCP and when the redundant constraint exists,the GOM always has a unique solution and so no additional processing is needed,whereas the corresponding DAE-LCP method may have singular cases with multiple solutions or no solutions.Using numerical examples,the GOM is verified to effectively solve the dynamics of multibody systems with redundant unilateral and bilateral constraints without additional redundancy processing.The GOM can also be applied to the optimal control of systems in the future and combined with the parameter optimization of systems to handle dynamic problems.The work given provides the dynamics and control of the complex system with a new train of thought and method.展开更多
In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is...In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.展开更多
The present work dealt with the preconcentration of rare earth elements in Saghand ore(Yazd province,Iran)which was achieved by Humphrey spiral using orthogonal optimization method after scrubbing the sample at 45%sol...The present work dealt with the preconcentration of rare earth elements in Saghand ore(Yazd province,Iran)which was achieved by Humphrey spiral using orthogonal optimization method after scrubbing the sample at 45%solid pulp density for 30 min.The pulp was diluted and was fed to a Humphrey spiral for upgrading.The process parameters considered were feed size,feed solids and feed rate,and Taguchi’s L9(34)orthogonal array(OA)was selected for optimization of the process.The results show that the feed rate and feed size were more significant than the other operation parameters of the process.It was also found that under optimal conditions,the concentrate grade of rare earth elements increased from2860 10 6to 6050 10 6and recovery reached to 58%.展开更多
The optimization of metamorphic mechanisms is different from that of the conventional mechanisms for its characteristics of multi-configuration. There exist complex coupled design variables and constraints in its mult...The optimization of metamorphic mechanisms is different from that of the conventional mechanisms for its characteristics of multi-configuration. There exist complex coupled design variables and constraints in its multiple different configuration optimization models. To achieve the compatible optimized results of these coupled design variables, an optimization method for metamorphic mechanisms is developed in the paper based on the principle of multidisciplinary design optimization（MDO）. Firstly, the optimization characteristics of the metamorphic mechanism are summarized distinctly by proposing the classification of design variables and constraints as well as coupling interactions among its different configuration optimization models. Further, collaborative optimization technique which is used in MDO is adopted for achieving the overall optimization performance. The whole optimization process is then proposed by constructing a two-level hierarchical scheme with global optimizer and configuration optimizer loops. The method is demonstrated by optimizing a planar five-bar metamorphic mechanism which has two configurations,and results show that it can achieve coordinated optimization results for the same parameters in different configuration optimization models.展开更多
Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and ther...Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.展开更多
In this paper,a passive muzzle arc control device(PMACD)of the augmented railguns is studied.By discussing its performance at different numbers of extra rails,a parameter optimization model is proposed.Through the cal...In this paper,a passive muzzle arc control device(PMACD)of the augmented railguns is studied.By discussing its performance at different numbers of extra rails,a parameter optimization model is proposed.Through the calculation model,it is found that the PMACD works well in the simple railgun,which refers to the gun that there is only one pair of rails in the inner bore.The PMACD may decrease the simple railgun’s armature peak current and muzzle arc,but affect its muzzle velocity not much.However,in the augmented railguns it has different characteristics.If the parameters of the PMACD are not selected suitable.It may increase the armature peak current and muzzle arc,but greatly decrease the velocity.The reason for this problem is that the extra rails generate a strong magnetic field in front of the armature,which induces a large current to change the armature current.It is also found that when the resistance and inductance parameters of the PMACD satisfy with the optimization formula,the PMACD can also play a good role in arc suppression in the augmented railguns.Experiments of an augmented railgun with a stainless steel PMACD are carried out to verify this optimization method.Results show that the muzzle arc is obviously controlled.This work may provide a reference for the design of the muzzle arc control device.展开更多
The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization pro...The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization problem and then the optimization method is adopted to solve the nonlinear least squares problem with the analytic form of new Jacobian matrix.In this way,joint variables of the general serial robot can be searched out quickly under the desired precision when positions of the three non-collinear end effector points are given.Compared with the general Newton iterative method,the proposed algorithm can search out the solution when the robot is at the singular configuration and the initial configuration used in the optimization method may also be the singular configuration.So the convergence domain is bigger than that of the general Newton iterative method.Another advantage of the proposed algorithm is that positions of the three non-collinear end effector points are usually much easier to be measured than the orientation of the end effector.The inverse displacement analysis of the general 6R(six-revolute-joint) serial robot is illustrated as an example and the simulation results verify the efficiency of the proposed algorithm.Because the three non-collinear points can be selected at random,the method can be applied to any other types of serial robots.展开更多
To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectr...To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectral method（GPM）. Different from the traditional trajectory optimization problem which generally considers the RLV as a point mass, the coupling between translational dynamics and rotational dynamics is taken into account. An optimization problem is formulated to minimize a performance index subject to 6-DOF equations of motion, including translational and rotational dynamics. A two-step optimal strategy is then introduced to reduce the large calculations caused by multiple variables and convergence confinement in 6-DOF trajectory optimization. The simulation results demonstrate that the 6-DOF trajectory optimal strategy for RLV is feasible.展开更多
In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional co...In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional conforming mesh in the finite element method with structured cells that are fixed during the whole design process with a view to avoid the sophisticated re-meshing and eventual mesh distortion.Then, B-spline shape functions are further implemented to yield a high-order continuity field along the cell boundary in stress analysis. By means of the implicit description of the shape boundary, stress sensitivity is analytically derived with respect to shape design variables. Finally, we illustrate the efficiency and accuracy of the proposed protocol by several numerical test cases as well as a whole design procedure carried out on an aeronautic turbine disk.展开更多
Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this ...Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this paper, we present numerical simulation results of the thermal radiative properties of a silicon carbide（Si C） thermal emitter/absorber composed of periodic microstructures. We illustrate different electromagnetic resonant modes which can be excited with the structure,such as surface phonon polaritons, magnetic polaritons and photonic crystal modes, and the process of radiation spectrum optimization based on a non-linear optimization algorithm. We show that the spectral and directional control of thermal emission/absorption can be efficiently achieved by adjusting the geometrical parameters of the structure. Moreover, the optimized spectrum is insensitive to 3% dimension modification.展开更多
基金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 the Beijing Municipal Science&Technology Commission(Z211100004421012),the Key Reaserch and Development Pro⁃gram of China(2022YFF0605902)。
文摘In this paper,a linear optimization method(LOM)for the design of terahertz circuits is presented,aimed at enhancing the simulation efficacy and reducing the time of the circuit design workflow.This method enables the rapid determination of optimal embedding impedance for diodes across a specific bandwidth to achieve maximum efficiency through harmonic balance simulations.By optimizing the linear matching circuit with the optimal embedding impedance,the method effectively segregates the simulation of the linear segments from the nonlinear segments in the frequency multiplier circuit,substantially improving the speed of simulations.The design of on-chip linear matching circuits adopts a modular circuit design strategy,incorporating fixed load resistors to simplify the matching challenge.Utilizing this approach,a 340 GHz frequency doubler was developed and measured.The results demonstrate that,across a bandwidth of 330 GHz to 342 GHz,the efficiency of the doubler remains above 10%,with an input power ranging from 98 mW to 141mW and an output power exceeding 13 mW.Notably,at an input power of 141 mW,a peak output power of 21.8 mW was achieved at 334 GHz,corresponding to an efficiency of 15.8%.
文摘This study aims to optimize the inbound traffic flow on on-ramps by considering low time costs,good speed stability,and high driving safety for mixed traffic flow.The optimal inlet gap is identified in advance,and trajectory guidance for vehicles entering the gap is determined under safety constraints.Based on the initial state and sequence of vehicles entering the merging area,individual vehicle trajectories are optimized sequentially.An optimization model and method for ramp entry trajectories in mixed traffic flow are developed,incorporating on-ramp vehicle entry sequencing and ordinary vehicle trajectory prediction.Key performance indicators,including driving safety,total travel time,parking wait probability,and trajectory smoothness,are compared and analyzed to evaluate the proposed approach.
文摘Response analysis of structures involving non-probabilistic uncertain parameters can be closely related to optimization.This paper provides a review on optimization-based methods for uncertainty analysis,with focusing attention on specific properties of adopted numerical optimization approaches.We collect and discuss the methods based on nonlinear programming,semidefinite programming,mixed-integer programming,mathematical programming with complementarity constraints,difference-of-convex programming,optimization methods using surrogate models and machine learning techniques,and metaheuristics.As a closely related topic,we also overview the methods for assessing structural robustness using non-probabilistic uncertainty modeling.We conclude the paper by drawing several remarks through this review.
基金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.
基金the project PID2022-139202OB-I00Neural Networks and Optimization Techniques for the Design and Safe Maintenance of Transportation Infrastructures:Volcanic Rock Geotechnics and Slope Stability(IA-Pyroslope),funded by the Spanish State Research Agency of the Ministry of Science,Innovation and Universities of Spain and the European Regional Development Fund,MCIN/AEI/10.13039/501100011033/FEDER,EU。
文摘Volcanic terrains exhibit a complex structure of pyroclastic deposits interspersed with sedimentary processes,resulting in irregular lithological sequences that lack lateral continuity and distinct stratigraphic patterns.This complexity poses significant challenges for slope stability analysis,requiring the development of specialized techniques to address these issues.This research presents a numerical methodology that incorporates spatial variability,nonlinear material characterization,and probabilistic analysis using a Monte Carlo framework to address this issue.The heterogeneous structure is represented by randomly assigning different lithotypes across the slope,while maintaining predefined global proportions.This contrasts with the more common approach of applying probabilistic variability to mechanical parameters within a homogeneous slope model.The material behavior is defined using complex nonlinear failure criteria,such as the Hoek-Brown model and a parabolic model with collapse,both implemented through linearization techniques.The Discontinuity Layout Optimization(DLO)method,a novel numerical approach based on limit analysis,is employed to efficiently incorporate these advances and compute the factor of safety of the slope.Within this framework,the Monte Carlo procedure is used to assess slope stability by conducting a large number of simulations,each with a different lithotype distribution.Based on the results,a hybrid method is proposed that combines probabilistic modeling with deterministic design principles for the slope stability assessment.As a case study,the methodology is applied to a 20-m-high vertical slope composed of three lithotypes(altered scoria,welded scoria,and basalt)randomly distributed in proportions of 15%,60%,and 25%,respectively.The results show convergence of mean values after approximately 400 simulations and highlight the significant influence of spatial heterogeneity,with variations of the factor of safety between 5 and 12 in 85%of cases.They also reveal non-circular and mid-slope failure wedges not captured by traditional stability methods.Finally,an equivalent normal probability distribution is proposed as a reliable approximation of the factor of safety for use in risk analysis and engineering decision-making.
基金supported by the National Natural Science Foundation of China (Grant Nos.40334040 and 40974033)the Promoting Foundation for Advanced Persons of Talent of NCWU
文摘Local and global optimization methods are widely used in geophysical inversion but each has its own advantages and disadvantages. The combination of the two methods will make it possible to overcome their weaknesses. Based on the simulated annealing genetic algorithm （SAGA） and the simplex algorithm, an efficient and robust 2-D nonlinear method for seismic travel-time inversion is presented in this paper. First we do a global search over a large range by SAGA and then do a rapid local search using the simplex method. A multi-scale tomography method is adopted in order to reduce non-uniqueness. The velocity field is divided into different spatial scales and velocities at the grid nodes are taken as unknown parameters. The model is parameterized by a bi-cubic spline function. The finite-difference method is used to solve the forward problem while the hybrid method combining multi-scale SAGA and simplex algorithms is applied to the inverse problem. The algorithm has been applied to a numerical test and a travel-time perturbation test using an anomalous low-velocity body. For a practical example, it is used in the study of upper crustal velocity structure of the A＇nyemaqen suture zone at the north-east edge of the Qinghai-Tibet Plateau. The model test and practical application both prove that the method is effective and robust.
文摘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.
文摘Deterministic optimization methods are combined with the Pareto front concept to solve multi-criterion design problems. The algorithm and the numerical implementation are applied to aerodynamic designs. Evolutionary algorithms （EAs） and the Pareto front concept are used to solve practical design problems in industry for its robustness in capturing convex, concave, discrete or discontinuous Pareto fronts of multi-objective optimization problems. However, the process is time-consuming. Therefore, deterministic optimization methods are introduced to capture the Pareto front, and the types of the captured Pareto front are explained. Numerical experiments show that the deterministic optimization method is a good alternative to EAs for capturing any convex and some concave Pareto fronts in multi-criterion aerodynamic optimization problems due to its efficiency.
基金This work was supported by the National Natural Science Foundation of China(10071037)
文摘In this paper, a new derivative free trust region method is developed based on the conic interpolation model for the unconstrained optimization. The conic interpolation model is built by means of the quadratic model function, the collinear scaling formula, quadratic approximation and interpolation. All the parameters in this model are determined by objective function interpolation condition. A new derivative free method is developed based upon this model and the global convergence of this new method is proved without any information on gradient.
基金This study was funded by the National Natural Science Foundation of China(Grant 11272167).
文摘The discontinuous dynamical problem of multi-point contact and collision in multi-body system has always been a hot and difficult issue in this field.Based on the Gauss’principle of least constraint,a unified optimization model for multibody system dynamics with multi-point contact and collision is established.The paper presents the study of the numerical solution scheme,in which particle swarm optimization method is used to deal with the corresponding optimization model.The article also presents the comparison of the Gauss optimization method(GOM)and the hybrid linear complementarity method(i.e.combining differential algebraic equations(DAEs)and linear complementarity problems(LCP)),commonly used to solve the dynamic contact problem of multibody systems with bilateral constraints.The results illustrate that,the GOM has the same advantage of dynamical modelling with LCP and when the redundant constraint exists,the GOM always has a unique solution and so no additional processing is needed,whereas the corresponding DAE-LCP method may have singular cases with multiple solutions or no solutions.Using numerical examples,the GOM is verified to effectively solve the dynamics of multibody systems with redundant unilateral and bilateral constraints without additional redundancy processing.The GOM can also be applied to the optimal control of systems in the future and combined with the parameter optimization of systems to handle dynamic problems.The work given provides the dynamics and control of the complex system with a new train of thought and method.
基金supported in part by the CNRST Morocco,the Volkswagen Foundation:Grant number I/79315Hydromed project
文摘In the present work, we investigate the inverse problem of reconstructing the parameter of an integro-differential parabolic equation, which comes from pollution problems in porous media, when the final observation is given. We use the optimal control framework to establish both the existence and necessary condition of the minimizer for the cost func- tional. Furthermore, we prove the stability and the local uniqueness of the minimizer. Some numerical results will be presented and discussed.
基金the deputy director of Research and Development in Atomic Energy of Iran for financial support as well as Nuclear Science and Technology Research Institute for technical support
文摘The present work dealt with the preconcentration of rare earth elements in Saghand ore(Yazd province,Iran)which was achieved by Humphrey spiral using orthogonal optimization method after scrubbing the sample at 45%solid pulp density for 30 min.The pulp was diluted and was fed to a Humphrey spiral for upgrading.The process parameters considered were feed size,feed solids and feed rate,and Taguchi’s L9(34)orthogonal array(OA)was selected for optimization of the process.The results show that the feed rate and feed size were more significant than the other operation parameters of the process.It was also found that under optimal conditions,the concentrate grade of rare earth elements increased from2860 10 6to 6050 10 6and recovery reached to 58%.
基金the fundamental support of the National Natural Science Foundation of China (Nos. 51105013, 51125020)the Beijing Natural Science Foundation of China (No. 3133042)the fundamental support provided by the China Scholarship Council and the State Key Laboratory of Robotics and System (HIT)
文摘The optimization of metamorphic mechanisms is different from that of the conventional mechanisms for its characteristics of multi-configuration. There exist complex coupled design variables and constraints in its multiple different configuration optimization models. To achieve the compatible optimized results of these coupled design variables, an optimization method for metamorphic mechanisms is developed in the paper based on the principle of multidisciplinary design optimization（MDO）. Firstly, the optimization characteristics of the metamorphic mechanism are summarized distinctly by proposing the classification of design variables and constraints as well as coupling interactions among its different configuration optimization models. Further, collaborative optimization technique which is used in MDO is adopted for achieving the overall optimization performance. The whole optimization process is then proposed by constructing a two-level hierarchical scheme with global optimizer and configuration optimizer loops. The method is demonstrated by optimizing a planar five-bar metamorphic mechanism which has two configurations,and results show that it can achieve coordinated optimization results for the same parameters in different configuration optimization models.
基金University Putra Malaysia under Putra Grant No.9531200。
文摘Previous studies about optimizing earthquake structural energy dissipation systems indicated that most existing techniques employ merely one or a few parameters as design variables in the optimization process,and thereby are only applicable only to simple,single,or multiple degree-of-freedom structures.The current approaches to optimization procedures take a specific damper with its properties and observe the effect of applying time history data to the building;however,there are many different dampers and isolators that can be used.Furthermore,there is a lack of studies regarding the optimum location for various viscous and wall dampers.The main aim of this study is hybridization of the particle swarm optimization(PSO) and gravitational search algorithm(GSA) to optimize the performance of earthquake energy dissipation systems(i.e.,damper devices) simultaneously with optimizing the characteristics of the structure.Four types of structural dampers device are considered in this study:(ⅰ) variable stiffness bracing(VSB) system,(ⅱ) rubber wall damper(RWD),(ⅲ) nonlinear conical spring bracing(NCSB) device,(iv) and multi-action stiffener(MAS) device.Since many parameters may affect the design of seismic resistant structures,this study proposes a hybrid of PSO and GSA to develop a hybrid,multi-objective optimization method to resolve the aforementioned problems.The characteristics of the above-mentioned damper devices as well as the section size for structural beams and columns are considered as variables for development of the PSO-GSA optimization algorithm to minimize structural seismic response in terms of nodal displacement(in three directions) as well as plastic hinge formation in structural members simultaneously with the weight of the structure.After that,the optimization algorithm is implemented to identify the best position of the damper device in the structural frame to have the maximum effect and minimize the seismic structure response.To examine the performance of the proposed PSO-GSA optimization method,it has been applied to a three-story reinforced structure equipped with a seismic damper device.The results revealed that the method successfully optimized the earthquake energy dissipation systems and reduced the effects of earthquakes on structures,which significantly increase the building’s stability and safety during seismic excitation.The analysis results showed a reduction in the seismic response of the structure regarding the formation of plastic hinges in structural members as well as the displacement of each story to approximately 99.63%,60.5%,79.13% and 57.42% for the VSB device,RWD,NCSB device,and MAS device,respectively.This shows that using the PSO-GSA optimization algorithm and optimized damper devices in the structure resulted in no structural damage due to earthquake vibration.
基金acknowledge the Fundamental Research Funds for the Central Universities(Grants No 309190112102)the Natural Science Foundation of Jiangsu Province(Grants No BK20200493).
文摘In this paper,a passive muzzle arc control device(PMACD)of the augmented railguns is studied.By discussing its performance at different numbers of extra rails,a parameter optimization model is proposed.Through the calculation model,it is found that the PMACD works well in the simple railgun,which refers to the gun that there is only one pair of rails in the inner bore.The PMACD may decrease the simple railgun’s armature peak current and muzzle arc,but affect its muzzle velocity not much.However,in the augmented railguns it has different characteristics.If the parameters of the PMACD are not selected suitable.It may increase the armature peak current and muzzle arc,but greatly decrease the velocity.The reason for this problem is that the extra rails generate a strong magnetic field in front of the armature,which induces a large current to change the armature current.It is also found that when the resistance and inductance parameters of the PMACD satisfy with the optimization formula,the PMACD can also play a good role in arc suppression in the augmented railguns.Experiments of an augmented railgun with a stainless steel PMACD are carried out to verify this optimization method.Results show that the muzzle arc is obviously controlled.This work may provide a reference for the design of the muzzle arc control device.
基金Funded by National Natural Science Foundation of China (No. 50905102)the Natural Science Foundation of Guangdong Province (Nos. 10151503101000033 and 8351503101000001)the Building Fund for the Academic Innovation Team of Shantou University (No. ITC10003)
文摘The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization problem and then the optimization method is adopted to solve the nonlinear least squares problem with the analytic form of new Jacobian matrix.In this way,joint variables of the general serial robot can be searched out quickly under the desired precision when positions of the three non-collinear end effector points are given.Compared with the general Newton iterative method,the proposed algorithm can search out the solution when the robot is at the singular configuration and the initial configuration used in the optimization method may also be the singular configuration.So the convergence domain is bigger than that of the general Newton iterative method.Another advantage of the proposed algorithm is that positions of the three non-collinear end effector points are usually much easier to be measured than the orientation of the end effector.The inverse displacement analysis of the general 6R(six-revolute-joint) serial robot is illustrated as an example and the simulation results verify the efficiency of the proposed algorithm.Because the three non-collinear points can be selected at random,the method can be applied to any other types of serial robots.
基金supported by the National Basic Research Program of China(973 Program)(2012CB720003)the National Natural Science Foundation of China(10772011)
文摘To be close to the practical flight process and increase the precision of optimal trajectory, a six-degree-offreedom（6-DOF） trajectory is optimized for the reusable launch vehicle（RLV） using the Gauss pseudospectral method（GPM）. Different from the traditional trajectory optimization problem which generally considers the RLV as a point mass, the coupling between translational dynamics and rotational dynamics is taken into account. An optimization problem is formulated to minimize a performance index subject to 6-DOF equations of motion, including translational and rotational dynamics. A two-step optimal strategy is then introduced to reduce the large calculations caused by multiple variables and convergence confinement in 6-DOF trajectory optimization. The simulation results demonstrate that the 6-DOF trajectory optimal strategy for RLV is feasible.
基金supported by the National Natura Science Foundation of China (Grant 51275424)973 Program (Gran2011CB610304)+1 种基金Research Fund for the Doctoral Program of Higher Education of China (Grant 20126102130003)the opening project (Grant KFJJ13-6M) of the State Key Laboratory of Explosion Science and Technology (Beijing Institute of Technology)
文摘In this work, a design procedure extending the B-spline based finite cell method into shape optimization is developed for axisymmetric solids involving the centrifugal force effect. We first replace the traditional conforming mesh in the finite element method with structured cells that are fixed during the whole design process with a view to avoid the sophisticated re-meshing and eventual mesh distortion.Then, B-spline shape functions are further implemented to yield a high-order continuity field along the cell boundary in stress analysis. By means of the implicit description of the shape boundary, stress sensitivity is analytically derived with respect to shape design variables. Finally, we illustrate the efficiency and accuracy of the proposed protocol by several numerical test cases as well as a whole design procedure carried out on an aeronautic turbine disk.
基金Project supported by the National Natural Science Foundation of China(Grant No.51076002)the National Basis Research Program of China(Grant No.2013CA328900)the Key Project of Complicated Electromagnetic Environment Laboratory of CAEP,China(Grant No.2015E0-01-1)
文摘Spectral and directional control of thermal emission based on excitation of confined electromagnetic resonant modes paves a viable way for the design and construction of microscale thermal emitters/absorbers. In this paper, we present numerical simulation results of the thermal radiative properties of a silicon carbide（Si C） thermal emitter/absorber composed of periodic microstructures. We illustrate different electromagnetic resonant modes which can be excited with the structure,such as surface phonon polaritons, magnetic polaritons and photonic crystal modes, and the process of radiation spectrum optimization based on a non-linear optimization algorithm. We show that the spectral and directional control of thermal emission/absorption can be efficiently achieved by adjusting the geometrical parameters of the structure. Moreover, the optimized spectrum is insensitive to 3% dimension modification.
