In the independent electro-hydrogen system(IEHS)with hybrid energy storage(HESS),achieving optimal scheduling is crucial.Still,it presents a challenge due to the significant deviations in values ofmultiple optimizatio...In the independent electro-hydrogen system(IEHS)with hybrid energy storage(HESS),achieving optimal scheduling is crucial.Still,it presents a challenge due to the significant deviations in values ofmultiple optimization objective functions caused by their physical dimensions.These deviations seriously affect the scheduling process.A novel standardization fusion method has been established to address this issue by analyzing the variation process of each objective function’s values.The optimal scheduling results of IEHS with HESS indicate that the economy and overall energy loss can be improved 2–3 times under different optimization methods.The proposed method better balances all optimization objective functions and reduces the impact of their dimensionality.When the cost of BESS decreases by approximately 30%,its participation deepens by about 1 time.Moreover,if the price of the electrolyzer is less than 15¥/kWh or if the cost of the fuel cell drops below 4¥/kWh,their participation will increase substantially.This study aims to provide a more reasonable approach to solving multi-objective optimization problems.展开更多
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%.展开更多
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 proposes an optimal midcourse guidance method for dual pulse air-to-air missiles,which is based on the framework of the linear Gauss pseudospectral model predictive control method.Firstly,a multistage optim...This paper proposes an optimal midcourse guidance method for dual pulse air-to-air missiles,which is based on the framework of the linear Gauss pseudospectral model predictive control method.Firstly,a multistage optimal control problem with unspecified terminal time is formulated.Secondly,the control and terminal time update formulas are derived analytically.In contrast to previous work,the derivation process fully considers the Hamiltonian function corresponding to the unspecified terminal time,which is coupled with control,state,and costate.On the assumption of small perturbation,a special algebraic equation is provided to represent the equivalent optimal condition for the terminal time.Also,using Gauss pseudospectral collocation,error propagation dynamical equations involving the first-order correction term of the terminal time are transformed into a set of algebraic equations.Furthermore,analytical modification formulas can be derived by associating those equations and optimal conditions to eliminate terminal error and approach nonlinear optimal control.Even with their mathematical complexity,these formulas produce more accurate control and terminal time corrections and remove reliance on task-related parameters.Finally,several numerical simulations,comparisons with typical methods,and Monte Carlo simulations have been done to verify its optimality,high convergence rate,great stability and robustness.展开更多
In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are reg...In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are regulated through sensitivity filtering tomitigate numerical instabilities associatedwith stress concentrations.Ap-norm aggregation function is employed to globalize local stress constraints,and a normalization technique linearly weights strain energy and stress,transforming the multi-objective problem into a single-objective formulation.The sensitivity of the objective function with respect to design variables is rigorously derived.Three numerical examples are presented,comparing the optimized structures in terms of strain energy,mass,and stress across five different mathematical models with varying combinations of optimization objectives.The results validate the effectiveness and feasibility of the proposed method for achieving a balanced design between structural stiffness and strength.This approach offers a new perspective for future research on stiffness-strength coordinated structural optimization.展开更多
Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method ge...Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method generated a specific trajectory for the UAV to effectively induce the proportional navigation missile to successfully intercept the obstacle,thereby accomplishing the evasive maneuver.The evasive maneuver was divided into two distinct stages,namely the collision-inducing phase and the fast departure phase.The obstacle potential field-based target selection algorithm was employed to identify the most appropriate target obstacle,while the induced trajectory was determined through a combination of receding horizon optimization and the hp-adaptive pseudo-spectral method.Simulation experiments were carried out under three different types of obstacle environments and one multiobstacle environment,and the simulation results show that the method proposed in this paper greatly improves the success rate of UAV evasive maneuvers,proving the effectiveness of this method.展开更多
Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation o...Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.展开更多
The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through...The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through post-processing,potentially altering the mechanical properties of the optimized structure.A topology optimization method of Movable Morphable Smooth Boundary(MMSB)is proposed based on the idea of mesh adaptation to solve the problem of jagged boundaries and the influence of post-processing.Based on the ICM method,the rational fraction function is introduced as the filtering function,and a topology optimization model with the minimum weight as the objective and the displacement as the constraint is established.A triangular mesh is utilized as the base mesh in this method.The mesh is re-divided in the optimization process based on the contour line,and a smooth boundary parallel to the contour line is obtained.Numerical examples demonstrate that the MMSB method effectively resolves the jagged boundary issues,leading to enhanced structural performance.展开更多
The quasi-rectangular tunnel represents a novel cross-section design,intended to supersede the traditional circular and rectangular tunnel formats.Due to the limited capacity of the tunnel vault to withstand vertical ...The quasi-rectangular tunnel represents a novel cross-section design,intended to supersede the traditional circular and rectangular tunnel formats.Due to the limited capacity of the tunnel vault to withstand vertical loads,an interior column is often installed at the center to enhance its load-bearing capacity.This study aims to develop a hyperstatic reaction method(HRM)for the analysis of deformation and structural integrity in this specific tunnel type.The computational model is validated through comparison with the corresponding finite element method(FEM)analysis.Following comprehensive validation,an ensemble machine learning(ML)model is proposed,using numerical benchmark data,to facilitate real-time design and optimization.Subsequently,three widely used ensemble models,i.e.random forest(RF),gradient boosting decision tree(GBDT),and extreme gradient boosting(XGBoost)are compared to identify the most efficient ML model for replacing the HRM model in the design optimization process.The performance metrics,such as the coefficient of determination R2 of about 0.999 and the mean absolute percentage error(MAPE)of about 1%,indicate that XGBoost outperforms the others,exhibiting excellent agreement with the HRM analysis.Additionally,the model demonstrates high computational efficiency,with prediction times measured in seconds.Finally,the HRM-XGBoost model is integrated with the well-known particle swarm optimization(PSO)for the real-time design optimization of quasi-rectangular tunnels,both with and without the interior column.A feature importance assessment is conducted to evaluate the sensitivity of design input features,enabling the selection of the most critical features for the optimization task.展开更多
To quantify the seismic effectiveness of the most commonly used fishing line tie up method for securing museum collections and optimize fixed strategies for exhibitions,shaking table tests of the seismic systems used ...To quantify the seismic effectiveness of the most commonly used fishing line tie up method for securing museum collections and optimize fixed strategies for exhibitions,shaking table tests of the seismic systems used for typical museum collection replicas have been carried out.The influence of body shape and fixed measure parameters on the seismic responses of replicas and the interaction behavior between replicas and fixed measures have been explored.Based on the results,seismic effectiveness evaluation indexes of the tie up method are proposed.Reasonable suggestions for fixed strategies are given,which provide a basis for the exhibition of delicate museum collections considering the principle of minimizing seismic responses and intervention.The analysis results show that a larger ratio of height of mass center to bottom diameter led to more intense rocking responses.Increasing the initial pretension of fishing lines was conducive to reducing the seismic responses and stress variation of the lines.Through comprehensive consideration of the interaction forces and effective securement,it is recommended to apply 20%of breaking stress as the initial pretension.For specific museum collections that cannot be effectively protected by the independent tie up method,an optimized strategy of a combination of fishing lines and fasteners is recommended.展开更多
Refined 3D modeling of mine slopes is pivotal for precise prediction of geological hazards.Aiming at the inadequacy of existing single modeling methods in comprehensively representing the overall and localized charact...Refined 3D modeling of mine slopes is pivotal for precise prediction of geological hazards.Aiming at the inadequacy of existing single modeling methods in comprehensively representing the overall and localized characteristics of mining slopes,this study introduces a new method that fuses model data from Unmanned aerial vehicles(UAV)tilt photogrammetry and 3D laser scanning through a data alignment algorithm based on control points.First,the mini batch K-Medoids algorithm is utilized to cluster the point cloud data from ground 3D laser scanning.Then,the elbow rule is applied to determine the optimal cluster number(K0),and the feature points are extracted.Next,the nearest neighbor point algorithm is employed to match the feature points obtained from UAV tilt photogrammetry,and the internal point coordinates are adjusted through the distanceweighted average to construct a 3D model.Finally,by integrating an engineering case study,the K0 value is determined to be 8,with a matching accuracy between the two model datasets ranging from 0.0669 to 1.0373 mm.Therefore,compared with the modeling method utilizing K-medoids clustering algorithm,the new modeling method significantly enhances the computational efficiency,the accuracy of selecting the optimal number of feature points in 3D laser scanning,and the precision of the 3D model derived from UAV tilt photogrammetry.This method provides a research foundation for constructing mine slope model.展开更多
The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the effi...The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the efficiency of RBDO algorithm,which hinders their application to high-dimensional engineering problems.To address these issues,this paper proposes an efficient decoupled RBDO method combining high dimensional model representation(HDMR)and the weight-point estimation method(WPEM).First,we decouple the RBDO model using HDMR and WPEM.Second,Lagrange interpolation is used to approximate a univariate function.Finally,based on the results of the first two steps,the original nested loop reliability optimization model is completely transformed into a deterministic design optimization model that can be solved by a series of mature constrained optimization methods without any additional calculations.Two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system with 28 design variables are analyzed to illustrate the performance and practicability of the proposed method.展开更多
This paper focuses on the online distributed optimization problem based on multi-agent systems. In this problem, each agent can only access its own cost function and a convex set, and can only exchange local state inf...This paper focuses on the online distributed optimization problem based on multi-agent systems. In this problem, each agent can only access its own cost function and a convex set, and can only exchange local state information with its current neighbors through a time-varying digraph. In addition, the agents do not have access to the information about the current cost functions until decisions are made. Different from most existing works on online distributed optimization, here we consider the case where the cost functions are strongly pseudoconvex and real gradients of the cost functions are not available. To handle this problem, a random gradient-free online distributed algorithm involving the multi-point gradient estimator is proposed. Of particular interest is that under the proposed algorithm, each agent only uses the estimation information of gradients instead of the real gradient information to make decisions. The dynamic regret is employed to measure the proposed algorithm. We prove that if the cumulative deviation of the minimizer sequence grows within a certain rate, then the expectation of dynamic regret increases sublinearly. Finally, a simulation example is given to corroborate the validity of our results.展开更多
Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coeff...Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.展开更多
Component sequence preservation is an intrinsic requirement in typical engineering applications, such as deployable chain-likestructures, 3D printing structures with contour-parallel toolpaths, additive manufacturing ...Component sequence preservation is an intrinsic requirement in typical engineering applications, such as deployable chain-likestructures, 3D printing structures with contour-parallel toolpaths, additive manufacturing of continuous fibre-reinforcedpolymer structures, customized stents, and soft robotics parts. This study presents a feature-driven method that preservescomponent sequences accounting for engineering requirements. The chain-of-bars design variables setting scheme is developedto realize the sequential component’s layout, which sets the current bar’s end point as the next bar’s start point. The total lengthof the printing path is constrained to reduce the consumption of material accurately. Also, the angle between adjacent bars isconstrained to avoid sharp angles at the turning point of the 3D printing path. Next, the sensitivity analysis considering theinter-dependence of substructures is performed. Several numerical examples are given to demonstrate the validity and merits ofthe proposed method in designing structures preserving component sequences.展开更多
A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization...A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.展开更多
Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analy...Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.展开更多
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.展开更多
To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivi...To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivity and thickness of hardened concrete layers on the detection of target objects are explored.Both numerical simulations and indoor tank tests indicate that when the ratio method is employed to correct the original collected data,the maximum allowable error for the isolation coefficient should not exceed 1%.Notably,when the ratio of hardened layer thickness to electrode spacing does not exceed 1,correction through this method significantly enhances the recognition capability of target objects.However,when the hardened layer thickness is greater than the electrode spacing by a factor of 2 or more,the ratio method cannot achieve satisfactory results.The case study of flood control engineering detection in the Zhangxi section of the Huangpen River in Dongzhi County demonstrates that the detection effect after correction by the ratio method is comparable to that for the adjacent unhardened pavement,and the influence of the hardened layer is obviously weakened,resulting in more reliable results.展开更多
A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines...A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.展开更多
基金sponsored by R&D Program of Beijing Municipal Education Commission(KM202410009013).
文摘In the independent electro-hydrogen system(IEHS)with hybrid energy storage(HESS),achieving optimal scheduling is crucial.Still,it presents a challenge due to the significant deviations in values ofmultiple optimization objective functions caused by their physical dimensions.These deviations seriously affect the scheduling process.A novel standardization fusion method has been established to address this issue by analyzing the variation process of each objective function’s values.The optimal scheduling results of IEHS with HESS indicate that the economy and overall energy loss can be improved 2–3 times under different optimization methods.The proposed method better balances all optimization objective functions and reduces the impact of their dimensionality.When the cost of BESS decreases by approximately 30%,its participation deepens by about 1 time.Moreover,if the price of the electrolyzer is less than 15¥/kWh or if the cost of the fuel cell drops below 4¥/kWh,their participation will increase substantially.This study aims to provide a more reasonable approach to solving multi-objective optimization problems.
基金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%.
基金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 the National Natural Science Foundation of China(No.62003019)the Young Talents Support Program of Beihang University,China(No.YWF-21-BJ-J-1180).
文摘This paper proposes an optimal midcourse guidance method for dual pulse air-to-air missiles,which is based on the framework of the linear Gauss pseudospectral model predictive control method.Firstly,a multistage optimal control problem with unspecified terminal time is formulated.Secondly,the control and terminal time update formulas are derived analytically.In contrast to previous work,the derivation process fully considers the Hamiltonian function corresponding to the unspecified terminal time,which is coupled with control,state,and costate.On the assumption of small perturbation,a special algebraic equation is provided to represent the equivalent optimal condition for the terminal time.Also,using Gauss pseudospectral collocation,error propagation dynamical equations involving the first-order correction term of the terminal time are transformed into a set of algebraic equations.Furthermore,analytical modification formulas can be derived by associating those equations and optimal conditions to eliminate terminal error and approach nonlinear optimal control.Even with their mathematical complexity,these formulas produce more accurate control and terminal time corrections and remove reliance on task-related parameters.Finally,several numerical simulations,comparisons with typical methods,and Monte Carlo simulations have been done to verify its optimality,high convergence rate,great stability and robustness.
基金funded by National Nature Science Foundation of China(92266203)National Nature Science Foundation of China(52205278)+1 种基金Key Projects of Shijiazhuang Basic Research Program(241791077A)Central Guide Local Science and Technology Development Fund Project of Hebei Province(246Z1022G).
文摘In this paper,a topology optimization method for coordinated stiffness and strength design is proposed under mass constraints,utilizing the Solid Isotropic Material with Penalization approach.Element densities are regulated through sensitivity filtering tomitigate numerical instabilities associatedwith stress concentrations.Ap-norm aggregation function is employed to globalize local stress constraints,and a normalization technique linearly weights strain energy and stress,transforming the multi-objective problem into a single-objective formulation.The sensitivity of the objective function with respect to design variables is rigorously derived.Three numerical examples are presented,comparing the optimized structures in terms of strain energy,mass,and stress across five different mathematical models with varying combinations of optimization objectives.The results validate the effectiveness and feasibility of the proposed method for achieving a balanced design between structural stiffness and strength.This approach offers a new perspective for future research on stiffness-strength coordinated structural optimization.
基金Natural Science Foundation of Heilongjiang Province of China(Grant No.YQ2022F012)the Fundamental Research Funds for the Central Universities(Grant No.HIT.OCEF.2023010)to provide fund for conducting experiments.
文摘Aiming at the missile avoidance problem of the unmanned aerial vehicle(UAV)in complex obstacle environments,this work proposes a collision-avoidance method based on receding horizon optimization.The proposed method generated a specific trajectory for the UAV to effectively induce the proportional navigation missile to successfully intercept the obstacle,thereby accomplishing the evasive maneuver.The evasive maneuver was divided into two distinct stages,namely the collision-inducing phase and the fast departure phase.The obstacle potential field-based target selection algorithm was employed to identify the most appropriate target obstacle,while the induced trajectory was determined through a combination of receding horizon optimization and the hp-adaptive pseudo-spectral method.Simulation experiments were carried out under three different types of obstacle environments and one multiobstacle environment,and the simulation results show that the method proposed in this paper greatly improves the success rate of UAV evasive maneuvers,proving the effectiveness of this method.
基金supported by the National Natural Science Foundation of China(Grant Nos.12372200 and 12072242).
文摘Parameterized level-set method(PLSM)has been proposed and developed for many years,and is renowned for its efficacy in ad-dressing topology optimization challenges associated with intricate boundaries and nucleation of new holes.However,most pertinent investigations in the field rely predominantly on fixed background mesh,which is never remeshed.Consequently,the mesh element partitioned by material interface during the optimization process necessitates approximation by using artificial interpolation models to obtain its element stiffness or other properties.This paper introduces a novel approach to topology op-timization by integrating the PLSM with body-fitted adaptive mesh and Helmholtz-type filter.Primarily,combining the PLSM with body-fitted adaptive mesh enables the regeneration of mesh based on the zero level-set interface.This not only precludes the direct traversal of the material interface through the mesh element during the topology optimization process,but also improves the accuracy of calculation.Additionally,the incorporation of a Helmholtz-type partial differential equation filter,relying solely on mesh information essential for finite element discretization,serves to regulate the topological complexity and the minimum feature size of the optimized structure.Leveraging these advantages,the topology optimization program demonstrates its versa-tility by successfully addressing various design problems,encompassing the minimum mean compliance problem and minimum energy dissipation problem.Ultimately,the result of numerical example indicates that the optimized structure exhibits a dis-tinct and smooth boundary,affirming the effective control over both topological complexity and the minimum feature size of the optimized structure.
基金supported by the National Natural Science Foundation of China(Grant 12472113).
文摘The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through post-processing,potentially altering the mechanical properties of the optimized structure.A topology optimization method of Movable Morphable Smooth Boundary(MMSB)is proposed based on the idea of mesh adaptation to solve the problem of jagged boundaries and the influence of post-processing.Based on the ICM method,the rational fraction function is introduced as the filtering function,and a topology optimization model with the minimum weight as the objective and the displacement as the constraint is established.A triangular mesh is utilized as the base mesh in this method.The mesh is re-divided in the optimization process based on the contour line,and a smooth boundary parallel to the contour line is obtained.Numerical examples demonstrate that the MMSB method effectively resolves the jagged boundary issues,leading to enhanced structural performance.
基金funded by the Hanoi University of Mining and Geology(Grant No.T23-44)The research is also funded by the German Research Foundation(DFG e Project number 518862444)in collaboration with the Vietnam National Foundation for Science and Technology Development(NAFOSTED)under grant number DFG.105e2022.03The third author was funded by the Postdoctoral Scholarship Program of the Vingroup Innovation Foundation(VINIF)(VINIF.2023.STS.15).
文摘The quasi-rectangular tunnel represents a novel cross-section design,intended to supersede the traditional circular and rectangular tunnel formats.Due to the limited capacity of the tunnel vault to withstand vertical loads,an interior column is often installed at the center to enhance its load-bearing capacity.This study aims to develop a hyperstatic reaction method(HRM)for the analysis of deformation and structural integrity in this specific tunnel type.The computational model is validated through comparison with the corresponding finite element method(FEM)analysis.Following comprehensive validation,an ensemble machine learning(ML)model is proposed,using numerical benchmark data,to facilitate real-time design and optimization.Subsequently,three widely used ensemble models,i.e.random forest(RF),gradient boosting decision tree(GBDT),and extreme gradient boosting(XGBoost)are compared to identify the most efficient ML model for replacing the HRM model in the design optimization process.The performance metrics,such as the coefficient of determination R2 of about 0.999 and the mean absolute percentage error(MAPE)of about 1%,indicate that XGBoost outperforms the others,exhibiting excellent agreement with the HRM analysis.Additionally,the model demonstrates high computational efficiency,with prediction times measured in seconds.Finally,the HRM-XGBoost model is integrated with the well-known particle swarm optimization(PSO)for the real-time design optimization of quasi-rectangular tunnels,both with and without the interior column.A feature importance assessment is conducted to evaluate the sensitivity of design input features,enabling the selection of the most critical features for the optimization task.
基金Beijing Nova Program under Grant No.2022036National Key Research and Development Program under Grant No.2019YFC1521000。
文摘To quantify the seismic effectiveness of the most commonly used fishing line tie up method for securing museum collections and optimize fixed strategies for exhibitions,shaking table tests of the seismic systems used for typical museum collection replicas have been carried out.The influence of body shape and fixed measure parameters on the seismic responses of replicas and the interaction behavior between replicas and fixed measures have been explored.Based on the results,seismic effectiveness evaluation indexes of the tie up method are proposed.Reasonable suggestions for fixed strategies are given,which provide a basis for the exhibition of delicate museum collections considering the principle of minimizing seismic responses and intervention.The analysis results show that a larger ratio of height of mass center to bottom diameter led to more intense rocking responses.Increasing the initial pretension of fishing lines was conducive to reducing the seismic responses and stress variation of the lines.Through comprehensive consideration of the interaction forces and effective securement,it is recommended to apply 20%of breaking stress as the initial pretension.For specific museum collections that cannot be effectively protected by the independent tie up method,an optimized strategy of a combination of fishing lines and fasteners is recommended.
基金funded by National Natural Science Foundation of China(Grant Nos.42272333,42277147).
文摘Refined 3D modeling of mine slopes is pivotal for precise prediction of geological hazards.Aiming at the inadequacy of existing single modeling methods in comprehensively representing the overall and localized characteristics of mining slopes,this study introduces a new method that fuses model data from Unmanned aerial vehicles(UAV)tilt photogrammetry and 3D laser scanning through a data alignment algorithm based on control points.First,the mini batch K-Medoids algorithm is utilized to cluster the point cloud data from ground 3D laser scanning.Then,the elbow rule is applied to determine the optimal cluster number(K0),and the feature points are extracted.Next,the nearest neighbor point algorithm is employed to match the feature points obtained from UAV tilt photogrammetry,and the internal point coordinates are adjusted through the distanceweighted average to construct a 3D model.Finally,by integrating an engineering case study,the K0 value is determined to be 8,with a matching accuracy between the two model datasets ranging from 0.0669 to 1.0373 mm.Therefore,compared with the modeling method utilizing K-medoids clustering algorithm,the new modeling method significantly enhances the computational efficiency,the accuracy of selecting the optimal number of feature points in 3D laser scanning,and the precision of the 3D model derived from UAV tilt photogrammetry.This method provides a research foundation for constructing mine slope model.
基金supported by the Innovation Fund Project of the Gansu Education Department(Grant No.2021B-099).
文摘The objective of reliability-based design optimization(RBDO)is to minimize the optimization objective while satisfying the corresponding reliability requirements.However,the nested loop characteristic reduces the efficiency of RBDO algorithm,which hinders their application to high-dimensional engineering problems.To address these issues,this paper proposes an efficient decoupled RBDO method combining high dimensional model representation(HDMR)and the weight-point estimation method(WPEM).First,we decouple the RBDO model using HDMR and WPEM.Second,Lagrange interpolation is used to approximate a univariate function.Finally,based on the results of the first two steps,the original nested loop reliability optimization model is completely transformed into a deterministic design optimization model that can be solved by a series of mature constrained optimization methods without any additional calculations.Two numerical examples of a planar 10-bar structure and an aviation hydraulic piping system with 28 design variables are analyzed to illustrate the performance and practicability of the proposed method.
基金supported by the National Natural Science Foundation of China(Nos.62103169,51875380)the China Postdoctoral Science Foundation(No.2021M691313).
文摘This paper focuses on the online distributed optimization problem based on multi-agent systems. In this problem, each agent can only access its own cost function and a convex set, and can only exchange local state information with its current neighbors through a time-varying digraph. In addition, the agents do not have access to the information about the current cost functions until decisions are made. Different from most existing works on online distributed optimization, here we consider the case where the cost functions are strongly pseudoconvex and real gradients of the cost functions are not available. To handle this problem, a random gradient-free online distributed algorithm involving the multi-point gradient estimator is proposed. Of particular interest is that under the proposed algorithm, each agent only uses the estimation information of gradients instead of the real gradient information to make decisions. The dynamic regret is employed to measure the proposed algorithm. We prove that if the cumulative deviation of the minimizer sequence grows within a certain rate, then the expectation of dynamic regret increases sublinearly. Finally, a simulation example is given to corroborate the validity of our results.
文摘Equivalent staggered-grid(ESG) as a new family of schemes has been utilized in seismic modeling,imaging,and inversion.Traditionally,the Taylor series expansion is often applied to calculate finite-difference(FD) coefficients on spatial derivatives,but the simulation results suffer serious numerical dispersion on a large frequency zone.We develop an optimized equivalent staggered-grid(OESG) FD method that can simultaneously suppress temporal and spatial dispersion for solving the second-order system of the 3 D elastic wave equation.On the one hand,we consider the coupling relations between wave speeds and spatial derivatives in the elastic wave equation and give three sets of FD coefficients with respect to the P-wave,S-wave,and converted-wave(C-wave) terms.On the other hand,a novel plane wave solution for the 3 D elastic wave equation is derived from the matrix decomposition method to construct the time-space dispersion relations.FD coefficients of the OESG method can be acquired by solving the new dispersion equations based on the Newton iteration method.Finally,we construct a new objective function to analyze P-wave,S-wave,and C-wave dispersion concerning frequencies.The dispersion analyses show that the presented method produces less modeling errors than the traditional ESG method.The synthetic examples demonstrate the effectiveness and superiority of the presented method.
基金supported by the Chinese Studentship Council(Grant No.201908060224)the Young Talent Fund of Association for Science and Technology in Shaanxi,China(Grant No.20230240)+1 种基金the National Natural Science Foundation of China(Grant No.11972308)Queen Mary University of London with the PhD fee waiver.
文摘Component sequence preservation is an intrinsic requirement in typical engineering applications, such as deployable chain-likestructures, 3D printing structures with contour-parallel toolpaths, additive manufacturing of continuous fibre-reinforcedpolymer structures, customized stents, and soft robotics parts. This study presents a feature-driven method that preservescomponent sequences accounting for engineering requirements. The chain-of-bars design variables setting scheme is developedto realize the sequential component’s layout, which sets the current bar’s end point as the next bar’s start point. The total lengthof the printing path is constrained to reduce the consumption of material accurately. Also, the angle between adjacent bars isconstrained to avoid sharp angles at the turning point of the 3D printing path. Next, the sensitivity analysis considering theinter-dependence of substructures is performed. Several numerical examples are given to demonstrate the validity and merits ofthe proposed method in designing structures preserving component sequences.
基金supported by a Major Research Project in Higher Education Institutions in Henan Province,with Project Number 23A560015.
文摘A new approach for flexoelectricmaterial shape optimization is proposed in this study.In this work,a proxymodel based on artificial neural network(ANN)is used to solve the parameter optimization and shape optimization problems.To improve the fitting ability of the neural network,we use the idea of pre-training to determine the structure of the neural network and combine different optimizers for training.The isogeometric analysis-finite element method(IGA-FEM)is used to discretize the flexural theoretical formulas and obtain samples,which helps ANN to build a proxy model from the model shape to the target value.The effectiveness of the proposed method is verified through two numerical examples of parameter optimization and one numerical example of shape optimization.
基金supported by the Na-tional Natural Science Foundation of China(No.52272369).
文摘Aiming at the time-optimal control problem of hypersonic vehicles(HSV)in ascending stage,a trigonometric regularization method(TRM)is introduced based on the indirect method of optimal control.This method avoids analyzing the switching function and distinguishing between singular control and bang-bang control,where the singular control problem is more complicated.While in bang-bang control,the costate variables are unsmooth due to the control jumping,resulting in difficulty in solving the two-point boundary value problem(TPBVP)induced by the indirect method.Aiming at the easy divergence when solving the TPBVP,the continuation method is introduced.This method uses the solution of the simplified problem as the initial value of the iteration.Then through solving a series of TPBVP,it approximates to the solution of the original complex problem.The calculation results show that through the above two methods,the time-optimal control problem of HSV in ascending stage under the complex model can be solved conveniently.
基金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.
基金National Key Research and Development Program of China(No.2021YFC3000103).
文摘To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivity and thickness of hardened concrete layers on the detection of target objects are explored.Both numerical simulations and indoor tank tests indicate that when the ratio method is employed to correct the original collected data,the maximum allowable error for the isolation coefficient should not exceed 1%.Notably,when the ratio of hardened layer thickness to electrode spacing does not exceed 1,correction through this method significantly enhances the recognition capability of target objects.However,when the hardened layer thickness is greater than the electrode spacing by a factor of 2 or more,the ratio method cannot achieve satisfactory results.The case study of flood control engineering detection in the Zhangxi section of the Huangpen River in Dongzhi County demonstrates that the detection effect after correction by the ratio method is comparable to that for the adjacent unhardened pavement,and the influence of the hardened layer is obviously weakened,resulting in more reliable results.
基金supported by the National Natural Science Foundation of China(11871312,12131014)the Natural Science Foundation of Shandong Province,China(ZR2023MA086)。
文摘A bicubic B-spline finite element method is proposed to solve optimal control problems governed by fourth-order semilinear parabolic partial differential equations.Its key feature is the selection of bicubic B-splines as trial functions to approximate the state and costate variables in two space dimensions.A Crank-Nicolson difference scheme is constructed for time discretization.The resulting numerical solutions belong to C2in space,and the order of the coefficient matrix is low.Moreover,the Bogner-Fox-Schmit element is considered for comparison.Two numerical experiments demonstrate the feasibility and effectiveness of the proposed method.
