In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
Computing electrostatic interaction on non-cooperative targets with unknown meshes is crucial for electrostatic-based space on-orbit services.Although meshes for electrostatic interaction computations can be reconstru...Computing electrostatic interaction on non-cooperative targets with unknown meshes is crucial for electrostatic-based space on-orbit services.Although meshes for electrostatic interaction computations can be reconstructed from point clouds,they are usually too dense,leading to high computational costs.This paper presents an optimization method for converting dense meshes into optimal meshes,enabling fast and accurate computation of the electrostatic interaction by point clouds.First,the dense mesh reconstructed from point clouds is simplified into a coarse mesh using local operators.Second,the simplified mesh is refined by an iterative strategy that integrates a lightweight method of moments and an impedance matrix inheritance technique,ultimately yielding an optimal mesh for computing the electrostatic interaction.Simulation results show that our method effectively optimizes dense meshes,making electrostatic interaction computations using point clouds approximately 63.4 times more efficient than the previous method.展开更多
In two-scale topology optimization,enhancing the connectivity between adjacent microstructures is crucial for achieving the collaborative optimization of micro-scale performance and macro-scale manufacturability.This ...In two-scale topology optimization,enhancing the connectivity between adjacent microstructures is crucial for achieving the collaborative optimization of micro-scale performance and macro-scale manufacturability.This paper proposes a two-scale concurrent topology optimization strategy aimed at improving the interface connection strength.This method employs a parametric approach to explicitly divide the micro-design domain into a“boundary connection region”and a“free design domain”at the initial stage of optimization.The boundary connection region is used to generate a connection layer that enhances the interface strength,while the free design domain is not constrained by this layer,thus fully exploiting the design potential of the material layout.During the optimization process,the solid isotropic material with penalization(SIMP)method is first used to optimize the material distribution in the free design domain,and filtering and projection techniques are employed to alleviate numerical instability and obtain a clear topological structure.Subsequently,the effective performance of the microstructure is calculated through homogenization and transferred to the macro-scale for global response analysis.Throughout the iterative process,the geometry of the connection layer remains unchanged,and only the free design domain is optimized,thereby achieving a balance between high performance and good manufacturability.The effectiveness of the proposed method is verified through numerical examples.展开更多
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.展开更多
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.展开更多
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.展开更多
This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are glob...This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions.展开更多
Sensitivity analysis and topology optimization of microstructures using strain energy-based method is presented. Compared with homogenization method, the strain energy-based method has advantages of higher computing e...Sensitivity analysis and topology optimization of microstructures using strain energy-based method is presented. Compared with homogenization method, the strain energy-based method has advantages of higher computing efficiency and simplified programming. Both the dual convex programming method and perimeter constraint scheme are used to optimize the 2D and 3D microstructures. Numerical results indicate that the strain energy-based method has the same effectiveness as that of homogenization method for orthotropic materials.展开更多
This paper presents a novel experimental design to greatly improve the calibration accuracy of the acceleration-insensitive bias and the acceleration-sensitive bias of the dynamically tuned gyroscopes（DTGs）.In order...This paper presents a novel experimental design to greatly improve the calibration accuracy of the acceleration-insensitive bias and the acceleration-sensitive bias of the dynamically tuned gyroscopes（DTGs）.In order to reduce experimental cost,the D-optimal criteria with constraints are constructed.The turntable positions and the number of test points are chosen to build D-optimal experimental designs.The D-optimal experimental designs are tested by multi-position calibration experiment for tactical-grade DTGs.Test results show that,with the same cost,the fit uncertainty is reduced by about 50% by using the D-optimal 8-position experimental procedure,compared to using a defacto standard experimental procedure in ANSI/IEEE Std 813-1988.Furthermore,the new experimental procedure almost achieves optimal accuracy with only 12-position which is half the cost of the widely adopted 24-position experimental procedure for achieving optimal accuracy.展开更多
The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practi...The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practice.An improved optimal elemental method is presented that defines a new objective function,and as a byproduct,circumvents the need for mass normalized modal shapes,which are also not readily available in practice.To solve the group of nonlinear equations created by the improved optimal method,the Lagrange multiplier method and Matlab function fmincon are employed.To deal with actual complex structures, the float-encoding genetic algorithm(FGA)is introduced to enhance the capability of the improved method.Two examples,a 7- degree of freedom(DOF)mass-spring system and a 53-DOF planar frame,respectively,are updated using the improved method. The example results demonstrate the advantages of the improved method over existing optimal methods,and show that the genetic algorithm is an effective way to update the models used for actual complex structures.展开更多
The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine...The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine. However, unreasonable design in the control schedule causes not only performance deterioration but also serious aerodynamic stability problems. Thus, in this work,a hybrid optimization method that automatically chooses the working modes and identifies the optimal and smooth control schedules is proposed, by combining the differential evolution algorithm and the Latin hypercube sampling method. The control schedule architecture does not only optimize the engine steady-state performance under different working modes but also solves the control-schedule discontinuity problem, especially during mode transition. The optimal control schedules are continuous and almost monotonic, and hence are strongly suitable for a control system, and are designed for two different working conditions, i.e., supersonic and subsonic throttling,which proves that the proposed hybrid method applies to various working conditions. The evaluation demonstrates that the proposed control method optimizes the engine performance, the surge margin of the compression components, and the range of the thrust during throttling.展开更多
A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to ...A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.展开更多
Target assembly is a key consumable material for producing thin fi lm used in the electronic packaging and devices. The residual stresses induced during the process of soldering are detrimental to the performance of t...Target assembly is a key consumable material for producing thin fi lm used in the electronic packaging and devices. The residual stresses induced during the process of soldering are detrimental to the performance of target assembly. In this work, the intensity and distribution of the soldering residual stress of Co/In/Cu target assembly subjected to a 20 W/(m^2 K) cooling condition corresponding to the actual air cooling process were studied, based on fi nite element simulation and Taguchi method, to optimize the sputtering target assembly. Effects of different control factors, including solder material, thickness of solder layer, target and backing plate, on the soldering residual stress of target assembly are investigated. The maximum residual stress is calculated as 9.28 MPa in the target located at 0.16 mm from target–solder layer interface and at a distance of 0.78 mm from symmetry axis. The optimal design in target assembly has the combination of indium solder material, cobalt target at 12 mm thick, solder layer at 0.8 mm thick, copper backing plate at 15 mm thick. Moreover, solder material is the most important factor among control factors in the target assembly.展开更多
The capacity and size of hydro-generator units are increasing with the rapid development of hydroelectric enterprises, and the vibration of the powerhouse structure has increasingly become a major problem. Field testi...The capacity and size of hydro-generator units are increasing with the rapid development of hydroelectric enterprises, and the vibration of the powerhouse structure has increasingly become a major problem. Field testing is an important method for research on dynamic identification and vibration mechanisms. Research on optimal sensor placement has become a very important topic due to the need to obtain effective testing information from limited test resources. To overcome inadequacies of the present methods, this paper puts forward the triaxial effective independence driving-point residue （EfI3-DPR3） method for optimal sensor placement. The Efl3-DPR3 method can incorporate both the maximum triaxial modal kinetic energy and linear independence of the triaxial target modes at the selected nodes. It was applied to the optimal placement oftriaxial sensors for vibration testing in a hydropower house, and satisfactory results were obtained. This method can provide some guidance for optimal placement of triaxial sensors of underground powerhouses.展开更多
The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the ...The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.展开更多
Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite ele...Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.展开更多
The classical iterative methods for finding roots of nonlinear equations,like the Newton method,Halley method,and Chebyshev method,have been modified previously to achieve optimal convergence order.However,the Househo...The classical iterative methods for finding roots of nonlinear equations,like the Newton method,Halley method,and Chebyshev method,have been modified previously to achieve optimal convergence order.However,the Householder method has so far not been modified to become optimal.In this study,we shall develop two new optimal Newton-Householder methods without memory.The key idea in the development of the new methods is the avoidance of the need to evaluate the second derivative.The methods fulfill the Kung-Traub conjecture by achieving optimal convergence order four with three functional evaluations and order eight with four functional evaluations.The efficiency indices of the methods show that methods perform better than the classical Householder’s method.With the aid of convergence analysis and numerical analysis,the efficiency of the schemes formulated in this paper has been demonstrated.The dynamical analysis exhibits the stability of the schemes in solving nonlinear equations.Some comparisons with other optimal methods have been conducted to verify the effectiveness,convergence speed,and capability of the suggested methods.展开更多
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.展开更多
A new multi-level analysis method of introducing the super-element modeling method, derived from the multi-level analysis method first proposed by O. F. Hughes, has been proposed in this paper to solve the problem of ...A new multi-level analysis method of introducing the super-element modeling method, derived from the multi-level analysis method first proposed by O. F. Hughes, has been proposed in this paper to solve the problem of high time cost in adopting a rational-based optimal design method for ship structural design. Furthermore,the method was verified by its effective application in optimization of the mid-ship section of a container ship. A full 3-D FEM model of a ship,suffering static and quasi-static loads, was used as the analyzing object for evaluating the structural performance of the mid-ship module, including static strength and buckling performance. Research results reveal that this new method could substantially reduce the computational cost of the rational-based optimization problem without decreasing its accuracy, which increases the feasibility and economic efficiency of using a rational-based optimal design method in ship structural design.展开更多
Since the connection of small-scale wind farms to distribution networks,power grid voltage stability has been reduced with increasing wind penetration in recent years,owing to the variable reactive power consumption o...Since the connection of small-scale wind farms to distribution networks,power grid voltage stability has been reduced with increasing wind penetration in recent years,owing to the variable reactive power consumption of wind generators.In this study,a two-stage reactive power optimization method based on the alternating direction method of multipliers(ADMM)algorithm is proposed for achieving optimal reactive power dispatch in wind farm-integrated distribution systems.Unlike existing optimal reactive power control methods,the proposed method enables distributed reactive power flow optimization with a two-stage optimization structure.Furthermore,under the partition concept,the consensus protocol is not needed to solve the optimization problems.In this method,the influence of the wake effect of each wind turbine is also considered in the control design.Simulation results for a mid-voltage distribution system based on MATLAB verified the effectiveness of the proposed method.展开更多
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金supported by the National Natural Science Foundation of China(No.62003269).
文摘Computing electrostatic interaction on non-cooperative targets with unknown meshes is crucial for electrostatic-based space on-orbit services.Although meshes for electrostatic interaction computations can be reconstructed from point clouds,they are usually too dense,leading to high computational costs.This paper presents an optimization method for converting dense meshes into optimal meshes,enabling fast and accurate computation of the electrostatic interaction by point clouds.First,the dense mesh reconstructed from point clouds is simplified into a coarse mesh using local operators.Second,the simplified mesh is refined by an iterative strategy that integrates a lightweight method of moments and an impedance matrix inheritance technique,ultimately yielding an optimal mesh for computing the electrostatic interaction.Simulation results show that our method effectively optimizes dense meshes,making electrostatic interaction computations using point clouds approximately 63.4 times more efficient than the previous method.
基金supported by the Science and Technology Research Project of Henan Province(242102241055)the Industry-University-Research Collaborative Innovation Base on Automobile Lightweight of“Science and Technology Innovation in Central Plains”(2024KCZY315)the Opening Fund of State Key Laboratory of Structural Analysis,Optimization and CAE Software for Industrial Equipment(GZ2024A03-ZZU).
文摘In two-scale topology optimization,enhancing the connectivity between adjacent microstructures is crucial for achieving the collaborative optimization of micro-scale performance and macro-scale manufacturability.This paper proposes a two-scale concurrent topology optimization strategy aimed at improving the interface connection strength.This method employs a parametric approach to explicitly divide the micro-design domain into a“boundary connection region”and a“free design domain”at the initial stage of optimization.The boundary connection region is used to generate a connection layer that enhances the interface strength,while the free design domain is not constrained by this layer,thus fully exploiting the design potential of the material layout.During the optimization process,the solid isotropic material with penalization(SIMP)method is first used to optimize the material distribution in the free design domain,and filtering and projection techniques are employed to alleviate numerical instability and obtain a clear topological structure.Subsequently,the effective performance of the microstructure is calculated through homogenization and transferred to the macro-scale for global response analysis.Throughout the iterative process,the geometry of the connection layer remains unchanged,and only the free design domain is optimized,thereby achieving a balance between high performance and good manufacturability.The effectiveness of the proposed method is verified through numerical examples.
基金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.
基金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 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.
文摘This paper explores the convergence of a class of optimally conditioned self scaling variable metric (OCSSVM) methods for unconstrained optimization. We show that this class of methods with Wolfe line search are globally convergent for general convex functions.
基金National Natural Science Foundation of China （90405016, 10676028） 973 Program （2006CB601205）+1 种基金 863 Project （2006AA04Z 122） Aeronautical Science Foundation （04B53080, 2006ZA 53006） and 111 Project （B07050）
文摘Sensitivity analysis and topology optimization of microstructures using strain energy-based method is presented. Compared with homogenization method, the strain energy-based method has advantages of higher computing efficiency and simplified programming. Both the dual convex programming method and perimeter constraint scheme are used to optimize the 2D and 3D microstructures. Numerical results indicate that the strain energy-based method has the same effectiveness as that of homogenization method for orthotropic materials.
基金National Natural Science Foundation of China (61071014)National Basic Research Program of China(2009CB72400201)
文摘This paper presents a novel experimental design to greatly improve the calibration accuracy of the acceleration-insensitive bias and the acceleration-sensitive bias of the dynamically tuned gyroscopes（DTGs）.In order to reduce experimental cost,the D-optimal criteria with constraints are constructed.The turntable positions and the number of test points are chosen to build D-optimal experimental designs.The D-optimal experimental designs are tested by multi-position calibration experiment for tactical-grade DTGs.Test results show that,with the same cost,the fit uncertainty is reduced by about 50% by using the D-optimal 8-position experimental procedure,compared to using a defacto standard experimental procedure in ANSI/IEEE Std 813-1988.Furthermore,the new experimental procedure almost achieves optimal accuracy with only 12-position which is half the cost of the widely adopted 24-position experimental procedure for achieving optimal accuracy.
基金The China Hi-Tech R&D Program(863 Program) Project Number 2001AA602023
文摘The optimal matrix method and optimal elemental method used to update finite element models may not provide accurate results.This situation occurs when the test modal model is incomplete,as is often the case in practice.An improved optimal elemental method is presented that defines a new objective function,and as a byproduct,circumvents the need for mass normalized modal shapes,which are also not readily available in practice.To solve the group of nonlinear equations created by the improved optimal method,the Lagrange multiplier method and Matlab function fmincon are employed.To deal with actual complex structures, the float-encoding genetic algorithm(FGA)is introduced to enhance the capability of the improved method.Two examples,a 7- degree of freedom(DOF)mass-spring system and a 53-DOF planar frame,respectively,are updated using the improved method. The example results demonstrate the advantages of the improved method over existing optimal methods,and show that the genetic algorithm is an effective way to update the models used for actual complex structures.
基金funded by National Nature Science Foundation of China(Nos.51776010 and 91860205)supported by the Academic Excellence Foundation of BUAA for PhD Students,China。
文摘The alternative working modes and flexible working states are the outstanding features of an adaptive cycle engine, with a proper control schedule design being the only way to exploit the performance of such an engine. However, unreasonable design in the control schedule causes not only performance deterioration but also serious aerodynamic stability problems. Thus, in this work,a hybrid optimization method that automatically chooses the working modes and identifies the optimal and smooth control schedules is proposed, by combining the differential evolution algorithm and the Latin hypercube sampling method. The control schedule architecture does not only optimize the engine steady-state performance under different working modes but also solves the control-schedule discontinuity problem, especially during mode transition. The optimal control schedules are continuous and almost monotonic, and hence are strongly suitable for a control system, and are designed for two different working conditions, i.e., supersonic and subsonic throttling,which proves that the proposed hybrid method applies to various working conditions. The evaluation demonstrates that the proposed control method optimizes the engine performance, the surge margin of the compression components, and the range of the thrust during throttling.
基金Projects(51275235, 50975135) supported by the National Natural Science Foundation of ChinaProject(U0934004) supported by the Natural Science Foundation of Guangdong Province, ChinaProject(2011CB707602) supported by the National Basic Research Program of China
文摘A new method for optimizing a butterfly-shaped linear ultrasonic motor was proposed to maximize its mechanical output. The finite element analysis technology and response surface methodology were combined together to realize the optimal design of the butterfly-shaped linear ultrasonic motor. First, the operation principle of the motor was introduced. Second, the finite element parameterized model of the stator of the motor was built using ANSYS parametric design language and some structure parameters of the stator were selected as design variables. Third, the sample points were selected in design variable space using latin hypercube Design. Through modal analysis and harmonic response analysis of the stator based on these sample points, the target responses were obtained. These sample points and response values were combined together to build a response surface model. Finally, the simplex method was used to find the optimal solution. The experimental results showed that many aspects of the design requirements of the butterfly-shaped linear ultrasonic motor have been fulfilled. The prototype motor fabricated based on the optimal design result exhibited considerably high dynamic performance, such as no-load speed of 873 ram/s, maximal thrust of 27.5 N, maximal efficiency of 43%, and thrust-weight ratio of 45.8.
基金supported financially by the National Key R&D Program of China (No. 2017YFB0305501)the National Natural Science Foundation of China (No. 51475220)the China Postdoctoral Science Foundation Funded Project (No. 2016M591464)
文摘Target assembly is a key consumable material for producing thin fi lm used in the electronic packaging and devices. The residual stresses induced during the process of soldering are detrimental to the performance of target assembly. In this work, the intensity and distribution of the soldering residual stress of Co/In/Cu target assembly subjected to a 20 W/(m^2 K) cooling condition corresponding to the actual air cooling process were studied, based on fi nite element simulation and Taguchi method, to optimize the sputtering target assembly. Effects of different control factors, including solder material, thickness of solder layer, target and backing plate, on the soldering residual stress of target assembly are investigated. The maximum residual stress is calculated as 9.28 MPa in the target located at 0.16 mm from target–solder layer interface and at a distance of 0.78 mm from symmetry axis. The optimal design in target assembly has the combination of indium solder material, cobalt target at 12 mm thick, solder layer at 0.8 mm thick, copper backing plate at 15 mm thick. Moreover, solder material is the most important factor among control factors in the target assembly.
基金supported by the National Natural Science Foundation of China (Grant No. 50909072)the New Teachers' Fund for Doctor Station, the Ministry of Education of China (Grant No. 20090032120082)the Communication Research Item for the West Area, the Ministry of Communications of China (Grant No. 2009328000084)
文摘The capacity and size of hydro-generator units are increasing with the rapid development of hydroelectric enterprises, and the vibration of the powerhouse structure has increasingly become a major problem. Field testing is an important method for research on dynamic identification and vibration mechanisms. Research on optimal sensor placement has become a very important topic due to the need to obtain effective testing information from limited test resources. To overcome inadequacies of the present methods, this paper puts forward the triaxial effective independence driving-point residue （EfI3-DPR3） method for optimal sensor placement. The Efl3-DPR3 method can incorporate both the maximum triaxial modal kinetic energy and linear independence of the triaxial target modes at the selected nodes. It was applied to the optimal placement oftriaxial sensors for vibration testing in a hydropower house, and satisfactory results were obtained. This method can provide some guidance for optimal placement of triaxial sensors of underground powerhouses.
基金supported by the National Natural Science Foundation of China(No.11472058)
文摘The attitude optimal control problem （OCP） of a two-rigid-body space- craft with two rigid bodies coupled by a ball-in-socket joint is considered. Based on conservation of angular momentum of the system without the external torque, a dynamic equation of three-dimensional attitude motion of the system is formulated. The attitude motion planning problem of the coupled-rigid-body spacecraft can be converted to a dis- crete nonlinear programming （NLP） problem using the Chebyshev-Gauss pseudospectral method （CGPM）. Solutions of the NLP problem can be obtained using the sequential quadratic programming （SQP） algorithm. Since the collocation points of the CGPM are Chebyshev-Gauss （CG） points, the integration of cost function can be approximated by the Clenshaw-Curtis quadrature, and the corresponding quadrature weights can be calculated efficiently using the fast Fourier transform （FFT）. To improve computational efficiency and numerical stability, the barycentric Lagrange interpolation is presented to substitute for the classic Lagrange interpolation in the approximation of state and con- trol variables. Furthermore, numerical float errors of the state differential matrix and barycentric weights can be alleviated using trigonometric identity especially when the number of CG points is large. A simple yet efficient method is used to avoid sensitivity to the initial values for the SQP algorithm using a layered optimization strategy from a feasible solution to an optimal solution. Effectiveness of the proposed algorithm is perfect for attitude motion planning of a two-rigid-body spacecraft coupled by a ball-in-socket joint through numerical simulation.
基金the National Natural Science Foundation of China(No.50678093)Program for Changjiang Scholars and Innovative Research Team in University(No.IRT00736)
文摘Based on the newly-developed element energy projection （EEP） method with optimal super-convergence order for computation of super-convergent results, an improved self-adaptive strategy for one-dimensional finite element method （FEM） is proposed. In the strategy, a posteriori errors are estimated by comparing FEM solutions to EEP super-convergent solutions with optimal order of super-convergence, meshes are refined by using the error-averaging method. Quasi-FEM solutions are used to replace the true FEM solutions in the adaptive process. This strategy has been found to be simple, clear, efficient and reliable. For most problems, only one adaptive step is needed to produce the required FEM solutions which pointwise satisfy the user specified error tolerances in the max-norm. Taking the elliptical ordinary differential equation of the second order as the model problem, this paper describes the fundamental idea, implementation strategy and computational algorithm and representative numerical examples are given to show the effectiveness and reliability of the proposed approach.
基金This research was supported by Universiti Kebangsaan Malaysia under research grant GUP-2019-033.
文摘The classical iterative methods for finding roots of nonlinear equations,like the Newton method,Halley method,and Chebyshev method,have been modified previously to achieve optimal convergence order.However,the Householder method has so far not been modified to become optimal.In this study,we shall develop two new optimal Newton-Householder methods without memory.The key idea in the development of the new methods is the avoidance of the need to evaluate the second derivative.The methods fulfill the Kung-Traub conjecture by achieving optimal convergence order four with three functional evaluations and order eight with four functional evaluations.The efficiency indices of the methods show that methods perform better than the classical Householder’s method.With the aid of convergence analysis and numerical analysis,the efficiency of the schemes formulated in this paper has been demonstrated.The dynamical analysis exhibits the stability of the schemes in solving nonlinear equations.Some comparisons with other optimal methods have been conducted to verify the effectiveness,convergence speed,and capability of the suggested methods.
基金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 Project of Ministry of Education and Finance(No.200512)the Project of the State Key Laboratory of ocean engineering(GKZD010053-10)
文摘A new multi-level analysis method of introducing the super-element modeling method, derived from the multi-level analysis method first proposed by O. F. Hughes, has been proposed in this paper to solve the problem of high time cost in adopting a rational-based optimal design method for ship structural design. Furthermore,the method was verified by its effective application in optimization of the mid-ship section of a container ship. A full 3-D FEM model of a ship,suffering static and quasi-static loads, was used as the analyzing object for evaluating the structural performance of the mid-ship module, including static strength and buckling performance. Research results reveal that this new method could substantially reduce the computational cost of the rational-based optimization problem without decreasing its accuracy, which increases the feasibility and economic efficiency of using a rational-based optimal design method in ship structural design.
基金support of The National Key Research and Development Program of China(Basic Research Class)(No.2017YFB0903000)the National Natural Science Foundation of China(No.U1909201)。
文摘Since the connection of small-scale wind farms to distribution networks,power grid voltage stability has been reduced with increasing wind penetration in recent years,owing to the variable reactive power consumption of wind generators.In this study,a two-stage reactive power optimization method based on the alternating direction method of multipliers(ADMM)algorithm is proposed for achieving optimal reactive power dispatch in wind farm-integrated distribution systems.Unlike existing optimal reactive power control methods,the proposed method enables distributed reactive power flow optimization with a two-stage optimization structure.Furthermore,under the partition concept,the consensus protocol is not needed to solve the optimization problems.In this method,the influence of the wake effect of each wind turbine is also considered in the control design.Simulation results for a mid-voltage distribution system based on MATLAB verified the effectiveness of the proposed method.
