In this paper, we present a nonmonotone smoothing Newton algorithm for solving the circular cone programming(CCP) problem in which a linear function is minimized or maximized over the intersection of an affine space w...In this paper, we present a nonmonotone smoothing Newton algorithm for solving the circular cone programming(CCP) problem in which a linear function is minimized or maximized over the intersection of an affine space with the circular cone. Based on the relationship between the circular cone and the second-order cone(SOC), we reformulate the CCP problem as the second-order cone problem(SOCP). By extending the nonmonotone line search for unconstrained optimization to the CCP, a nonmonotone smoothing Newton method is proposed for solving the CCP. Under suitable assumptions, the proposed algorithm is shown to be globally and locally quadratically convergent. Some preliminary numerical results indicate the effectiveness of the proposed algorithm for solving the CCP.展开更多
Based on the differential properties of the smoothing metric projector onto the second-order cone,we prove that,for a locally optimal solution to a nonlinear second-order cone programming problem,the nonsingularity of...Based on the differential properties of the smoothing metric projector onto the second-order cone,we prove that,for a locally optimal solution to a nonlinear second-order cone programming problem,the nonsingularity of the Clarke's generalized Jacobian of the smoothing Karush-Kuhn-Tucker system,constructed by the smoothing metric projector,is equivalent to the strong second-order sufficient condition and constraint nondegeneracy,which is in turn equivalent to the strong regularity of the Karush-Kuhn-Tucker point.Moreover,this nonsingularity property guarantees the quadratic convergence of the corresponding smoothing Newton method for solving a Karush-Kuhn-Tucker point.Interestingly,the analysis does not need the strict complementarity condition.展开更多
In this paper,we accomplish the unified convergence analysis of a second-order method of multipliers(i.e.,a second-order augmented Lagrangian method)for solving the conventional nonlinear conic optimization problems.S...In this paper,we accomplish the unified convergence analysis of a second-order method of multipliers(i.e.,a second-order augmented Lagrangian method)for solving the conventional nonlinear conic optimization problems.Specifically,the algorithm that we investigate incorporates a specially designed nonsmooth(generalized)Newton step to furnish a second-order update rule for the multipliers.We first show in a unified fashion that under a few abstract assumptions,the proposed method is locally convergent and possesses a(nonasymptotic)superlinear convergence rate,even though the penalty parameter is fixed and/or the strict complementarity fails.Subsequently,we demonstrate that for the three typical scenarios,i.e.,the classic nonlinear programming,the nonlinear second-order cone programming and the nonlinear semidefinite programming,these abstract assumptions are nothing but exactly the implications of the iconic sufficient conditions that are assumed for establishing the Q-linear convergence rates of the method of multipliers without assuming the strict complementarity.展开更多
The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess v...The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.展开更多
In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a gene...In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.展开更多
In high-speed multiuser Time Reversal(TR)downlink systems,the transmission rate is degraded due to the presence of severe inter-user and inter-symbol interference.Moreover,maximizing the weighted sum rate in such syst...In high-speed multiuser Time Reversal(TR)downlink systems,the transmission rate is degraded due to the presence of severe inter-user and inter-symbol interference.Moreover,maximizing the weighted sum rate in such systems is a critical objective,since the weighting factors represent the priority of different users in different applications.However,it faces significant challenges as it is an NP-hard and non-convex problem.In order to suppress these interferences and maximize the weighted sum rate,in this paper we present a novel approach for the joint design of the pre-filters.The proposed method applies successive convex approximation to transform the original problem into a Second-Order Cone Programming(SOCP)problem.Then,a low-complexity iterative algorithm is provided to effectively solve the resulting SOCP problem.According to the simulation results,the proposed method reaches a local optimum within a few iterations and demonstrates superior performance in terms of weighted sum rate compared to the current algorithm.展开更多
An optimal operation scheme is of great significance in islanded distribution networks to restore critical loads and has recently attracted considerable attention.In this paper,an optimal power flow(OPF)model for isla...An optimal operation scheme is of great significance in islanded distribution networks to restore critical loads and has recently attracted considerable attention.In this paper,an optimal power flow(OPF)model for islanded distribution networks equipped with soft open points(SOPs)is proposed.Unlike in the grid-connected mode,the adequacy of local power generation in distribution networks is critical for islanded systems.The proposed approach utilizes the power output of local distributed generations(DGs)and the benefits of reactive power compensation provided by SOPs to allow maximum loadability.To exploit the available resources,an optimal secondary droop control strategy is introduced for the islanded distribution networks,thereby minimizing load shedding.The formulated OPF problem is essentially a mixed-integer nonlinear programming(MINLP)model.To guarantee the computation efficiency and accuracy.A successive mixed-integer second-order cone programming(SMISOCP)algorithm is proposed for handling the nonlinear islanded power flow formulations.Two case studies,incorporating a modified IEEE 33-bus system and IEEE 123-bus system,are performed to test the effectiveness of the proposed approach.展开更多
A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The ...A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.展开更多
More demand-side flexible resources(DFRs)are participating in the frequency regulation of renewable power systems,whose heterogeneous characteristics have a significant impact on the system frequency response.Conseque...More demand-side flexible resources(DFRs)are participating in the frequency regulation of renewable power systems,whose heterogeneous characteristics have a significant impact on the system frequency response.Consequently,selecting suitable DFRs poses a formidable challenge for independent system operators(ISO).In this paper,a reserve allocation methodology for heterogeneous DFRs is proposed to manage the risk of power system frequency.Firstly,a performance curve is developed to describe the cost,capacity,and response speed of DFRs.Moreover,a clustering method for multiple distributed DFRs is conducted to calculate the aggregated performance curves and uncertainty coefficients.Then,the frequency security criterion considering DFRs’performance is constructed,whose linearity makes it can be easily coupled into the system scheduling model and solved.Furthermore,a risk management model for DFRs considering frequency-chance-constraint is proposed to make a trade-off between cost and frequency security.Finally,the model is transformed into mixed integer second-order cone programming(MISOCP)and solved by the commercial solver.The proposed model is validated by the IEEE 30 and IEEE 118 bus systems.展开更多
基金supported by the National Natural Science Foundation of China(11401126,71471140 and 11361018)Guangxi Natural Science Foundation(2016GXNSFBA380102 and 2014GXNSFFA118001)+2 种基金Guangxi Key Laboratory of Cryptography and Information Security(GCIS201618)Guangxi Key Laboratory of Automatic Detecting Technology and Instruments(YQ15112 and YQ16112)China
文摘In this paper, we present a nonmonotone smoothing Newton algorithm for solving the circular cone programming(CCP) problem in which a linear function is minimized or maximized over the intersection of an affine space with the circular cone. Based on the relationship between the circular cone and the second-order cone(SOC), we reformulate the CCP problem as the second-order cone problem(SOCP). By extending the nonmonotone line search for unconstrained optimization to the CCP, a nonmonotone smoothing Newton method is proposed for solving the CCP. Under suitable assumptions, the proposed algorithm is shown to be globally and locally quadratically convergent. Some preliminary numerical results indicate the effectiveness of the proposed algorithm for solving the CCP.
基金supported by National Natural Science Foundation of China (Grant Nos.10771026,10901094)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry of China
文摘Based on the differential properties of the smoothing metric projector onto the second-order cone,we prove that,for a locally optimal solution to a nonlinear second-order cone programming problem,the nonsingularity of the Clarke's generalized Jacobian of the smoothing Karush-Kuhn-Tucker system,constructed by the smoothing metric projector,is equivalent to the strong second-order sufficient condition and constraint nondegeneracy,which is in turn equivalent to the strong regularity of the Karush-Kuhn-Tucker point.Moreover,this nonsingularity property guarantees the quadratic convergence of the corresponding smoothing Newton method for solving a Karush-Kuhn-Tucker point.Interestingly,the analysis does not need the strict complementarity condition.
基金supported by National Natural Science Foundation of China (Grant No. 11801158)the Hunan Provincial Natural Science Foundation of China (Grant No. 2019JJ50040)+2 种基金the Fundamental Research Funds for the Central Universities in Chinasupported by National Natural Science Foundation of China (Grant No. 11871002)the General Program of Science and Technology of Beijing Municipal Education Commission (Grant No. KM201810005004)
文摘In this paper,we accomplish the unified convergence analysis of a second-order method of multipliers(i.e.,a second-order augmented Lagrangian method)for solving the conventional nonlinear conic optimization problems.Specifically,the algorithm that we investigate incorporates a specially designed nonsmooth(generalized)Newton step to furnish a second-order update rule for the multipliers.We first show in a unified fashion that under a few abstract assumptions,the proposed method is locally convergent and possesses a(nonasymptotic)superlinear convergence rate,even though the penalty parameter is fixed and/or the strict complementarity fails.Subsequently,we demonstrate that for the three typical scenarios,i.e.,the classic nonlinear programming,the nonlinear second-order cone programming and the nonlinear semidefinite programming,these abstract assumptions are nothing but exactly the implications of the iconic sufficient conditions that are assumed for establishing the Q-linear convergence rates of the method of multipliers without assuming the strict complementarity.
基金part of the TPS projecta Vied-Newton PhD scholarship+1 种基金a Dixon scholarship from Imperial College London,UKthe Dean’s Fund from Imperial College London for financial support(2017-2020)。
文摘The node-based smoothed finite element method(NS-FEM)is shortly presented for calculations of the static and seismic bearing capacities of shallow strip footings.A series of computations has been performed to assess variations in seismic bearing capacity factors with both horizontal and vertical seismic accelerations.Numerical results obtained agree very well with those using the slip-line method,revealing that the magnitude of the seismic bearing capacity is highly dependent upon the combinations of various directions of both components of the seismic acceleration.An upward vertical seismic acceleration reduces the seismic bearing capacity compared to the downward vertical seismic acceleration in calculations.In addition,particular emphasis is placed on a separate estimation of the effects of soil and superstructure inertia on each seismic bearing capacity component.While the effect of inertia forces arising in the soil on the seismic bearing capacity is non-trivial,and the superstructure inertia is the major contributor to reductions in the seismic bearing capacity.Both tables and charts are given for practical application to the seismic design of the foundations.
基金supported by the Swiss National Science Foundation(Grant No.189882)the National Natural Science Foundation of China(Grant No.41961134032)support provided by the New Investigator Award grant from the UK Engineering and Physical Sciences Research Council(Grant No.EP/V012169/1).
文摘In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.
基金partially supported by the following funding sources:The National Natural Science Foundation of China(No.61771084)the Chongqing Graduate Scientific Research Innovation Project(No.CYB21200)。
文摘In high-speed multiuser Time Reversal(TR)downlink systems,the transmission rate is degraded due to the presence of severe inter-user and inter-symbol interference.Moreover,maximizing the weighted sum rate in such systems is a critical objective,since the weighting factors represent the priority of different users in different applications.However,it faces significant challenges as it is an NP-hard and non-convex problem.In order to suppress these interferences and maximize the weighted sum rate,in this paper we present a novel approach for the joint design of the pre-filters.The proposed method applies successive convex approximation to transform the original problem into a Second-Order Cone Programming(SOCP)problem.Then,a low-complexity iterative algorithm is provided to effectively solve the resulting SOCP problem.According to the simulation results,the proposed method reaches a local optimum within a few iterations and demonstrates superior performance in terms of weighted sum rate compared to the current algorithm.
基金This work was supported in part by the science and technology project of State Grid Corporation of China under Grant 5400-201955369A-0-0-00。
文摘An optimal operation scheme is of great significance in islanded distribution networks to restore critical loads and has recently attracted considerable attention.In this paper,an optimal power flow(OPF)model for islanded distribution networks equipped with soft open points(SOPs)is proposed.Unlike in the grid-connected mode,the adequacy of local power generation in distribution networks is critical for islanded systems.The proposed approach utilizes the power output of local distributed generations(DGs)and the benefits of reactive power compensation provided by SOPs to allow maximum loadability.To exploit the available resources,an optimal secondary droop control strategy is introduced for the islanded distribution networks,thereby minimizing load shedding.The formulated OPF problem is essentially a mixed-integer nonlinear programming(MINLP)model.To guarantee the computation efficiency and accuracy.A successive mixed-integer second-order cone programming(SMISOCP)algorithm is proposed for handling the nonlinear islanded power flow formulations.Two case studies,incorporating a modified IEEE 33-bus system and IEEE 123-bus system,are performed to test the effectiveness of the proposed approach.
基金This is part of the TPS projecta Vied-Newton PhD scholarship and a Dixon scholarship from Imperial College London, UK, for supporting his studies at Imperial College Londonthe Dean’s Fund from Imperial College London for financial support (2017-2020).
文摘A numerical procedure using a stable cell-based smoothed finite element method(CS-FEM)is presented for estimation of stability of a square tunnel in the soil where the shear strength increases linearly with depth.The kinematically admissible displacement fields are approximated by uniform quadrilateral elements in conjunction with the strain smoothing technique,eliminating volumetric locking issues and the singularity associated with the MohreCoulomb model.First,a rich set of simulations was performed to compute the static stability of a square tunnel with different geometries and soil conditions.The presented results are in excellent agreement with the upper and lower bound solutions using the standard finite element method(FEM).The stability charts and tables are given for practical use in the tunnel design,along with a newly proposed formulation for predicting the undrained stability of a single square tunnel.Second,the seismic stability number was computed using the present numerical approach.Numerical results reveal that the seismic stability number reduces with an increasing value of the horizontal seismic acceleration(a_(h)),for both cases of the weightless soil and the soil with unit weight.Third,the link between the static and seismic stability numbers is described using corrective factors that represent reductions in the tunnel stability due to seismic loadings.It is shown from the numerical results that the corrective factor becomes larger as the unit weight of soil mass increases;however,the degree of the reduction in seismic stability number tends to reduce for the case of the homogeneous soil.Furthermore,this advanced numerical procedure is straightforward to extend to three-dimensional(3D)limit analysis and is readily applicable for the calculation of the stability of tunnels in highly anisotropic and heterogeneous soils which are often encountered in practice.
基金supported by the Key Science and Technology Project of China Southern Power Grid Corporation(Grant No.090000KK52220020)。
文摘More demand-side flexible resources(DFRs)are participating in the frequency regulation of renewable power systems,whose heterogeneous characteristics have a significant impact on the system frequency response.Consequently,selecting suitable DFRs poses a formidable challenge for independent system operators(ISO).In this paper,a reserve allocation methodology for heterogeneous DFRs is proposed to manage the risk of power system frequency.Firstly,a performance curve is developed to describe the cost,capacity,and response speed of DFRs.Moreover,a clustering method for multiple distributed DFRs is conducted to calculate the aggregated performance curves and uncertainty coefficients.Then,the frequency security criterion considering DFRs’performance is constructed,whose linearity makes it can be easily coupled into the system scheduling model and solved.Furthermore,a risk management model for DFRs considering frequency-chance-constraint is proposed to make a trade-off between cost and frequency security.Finally,the model is transformed into mixed integer second-order cone programming(MISOCP)and solved by the commercial solver.The proposed model is validated by the IEEE 30 and IEEE 118 bus systems.
