No Arabic abstract
We consider the problem of optimal placement of concentrated masses along a massless elastic column that is clamped at one end and loaded by a nonconservative follower force at the free end. The goal is to find the largest possible interval such that the variation in the loading parameter within this interval preserves stability of the structure. The stability constraint is nonconvex and nonsmooth, making the optimization problem quite challenging. We give a detailed analytical treatment for the case of two masses, arguing that the optimal parameter configuration approaches the flutter and divergence boundaries of the stability region simultaneously. Furthermore, we conjecture that this property holds for any number of masses, which in turn suggests a simple formula for the maximal load interval for $n$ masses. This conjecture is strongly supported by extensive computational results, obtained using the recently developed open-source software package GRANSO (GRadient-based Algorithm for Non-Smooth Optimization) to maximize the load interval subject to an appropriate formulation of the nonsmooth stability constraint. We hope that our work will provide a foundation for new approaches to classical long-standing problems of stability optimization for nonconservative elastic systems arising in civil and mechanical engineering.
The concept of leader--follower (or Stackelberg) equilibrium plays a central role in a number of real--world applications of game theory. While the case with a single follower has been thoroughly investigated, results with multiple followers are only sporadic and the problem of designing and evaluating computationally tractable equilibrium-finding algorithms is still largely open. In this work, we focus on the fundamental case where multiple followers play a Nash equilibrium once the leader has committed to a strategy---as we illustrate, the corresponding equilibrium finding problem can be easily shown to be $mathcal{FNP}$--hard and not in Poly--$mathcal{APX}$ unless $mathcal{P} = mathcal{NP}$ and therefore it is one among the hardest problems to solve and approximate. We propose nonconvex mathematical programming formulations and global optimization methods to find both exact and approximate equilibria, as well as a heuristic black box algorithm. All the methods and formulations that we introduce are thoroughly evaluated computationally.
We investigate the buckling under compression of a slender beam with a distributed lateral elastic support, for which there is an associated cost. For a given cost, we study the optimal choice of support to protect against Euler buckling. We show that with only weak lateral support, the optimum distribution is a delta-function at the centre of the beam. When more support is allowed, we find numerically that the optimal distribution undergoes a series of bifurcations. We obtain analytical expressions for the buckling load around the first bifurcation point and corresponding expansions for the optimal position of support. Our theoretical predictions, including the critical exponent of the bifurcation, are confirmed by computer simulations.
We consider the critical points (equilibria) of a planar potential generated by $n$ Newtonian point masses augmented with a quadratic term (such as arises from a centrifugal effect). Particular cases of this problem have been considered previously in studies of the circular restricted $n$-body problem. We show that the number of equilibria is finite for a generic set of parameters, and we establish estimates for the number of equilibria. We prove that the number of equilibria is bounded below by $n+1$, and we provide examples to show that this lower bound is sharp. We prove an upper bound on the number of equilibria that grows exponentially in $n$. In order to establish a lower bound on the maximum number of equilibria, we analyze a class of examples, referred to as ``ring configurations, consisting of $n-1$ equal masses positioned at vertices of a regular polygon with an additional mass located at the center. Previous numerical observations indicate that these configurations can produce as many as $5n-5$ equilibria. We verify analytically that the ring configuration has at least $5n-5$ equilibria when the central mass is sufficiently small. We conjecture that the maximum number of equilibria grows linearly with the number of point masses. We also discuss some mathematical similarities to other equilibrium problems in mathematical physics, namely, Maxwells problem from electrostatics and the image counting problem from gravitational lensing.
The load-flow equations are the main tool to operate and plan electrical networks. For transmission or distribution networks these equations can be simplified into a linear system involving the graph Laplacian and the power input vector. We show, using spectral graph theory, how to solve this system efficiently. This spectral approach gives a new geometric view of the network and power vector. This formulation yields a Parseval-like relation for the $L_2$ norm of the power in the lines. Using this relation as a guide, we show that a small number of eigenvector components of the power vector are enough to obtain an estimate of the solution. This would allow fast reconfiguration of networks and better planning.
We set a shortcut-to-adiabaticity strategy to design the trolley motion in a double-pendulum bridge crane. The trajectories found guarantee payload transport without residual excitation regardless of the initial conditions within the small oscillations regime. The results are compared with exact dynamics to set the working domain of the approach. The method is free from instabilities due boundary effects or to resonances with the two natural frequencies.