Computation of Critical Boundaries on Equilibrium Manifolds

The study of various equilibrium phenomena leads to nonlinear equations (1) $F(y,u) = 0$, where $y \in R^n $ is a vector of behavior or state variables, $u \in R^p $ a vector of $p \geqq 2$ parameters or controls and $F:D \subset R^n \times R^p \to R^n $ a sufficiently differentiable map. The solution set of (1) in $R^n \times R^p $ is often called the equilibrium surface of the problem, and in many applications it is of interest to determine the critical boundary, that is, the set of solutions of (1) where the derivative $D_y F(x)$ is singular. For example, in structural problems these points may represent buckling points. After characterizing the various properties of the problem, we present three different numerical methods which allow for a computational trace of paths in the critical boundary. These methods represent extensions of the earlier-developed locally parametrized continuation method for tracing regular paths on the equilibrium surface. They permit, for instance, a direct computational deter...