This paper presents two complementary mathematical programming based approaches for the accurate safety assessment of semirigid elastoplastic frames under quasistatic loads. The inelastic behavior of the flexible connections and material plasticity are accommodated through piecewise linearized nonlinear yield surfaces. As is necessary for this class of structures, geometric nonlinearity is taken into account. Moreover, only a 2nd-order geometric approximation is included as this is sufficiently accurate for practical structures. The work described has a twofold contribution. First, we develop an algorithm that can robustly and efficiently process the complete (path-dependent) nonholonomic response of the structure in a stepwise (path-independent) holonomic fashion. The governing formulation is cast in mixed statickinematic variables and leads naturally to what is known in the mathematical programming literature as a mixed complementarity problem (MCP). The novelty of the proposed algorithm is that it processes the MCP directly without using some iterative (and often cumbersome) predictor-corrector procedure. Second, in the spirit of simplified analyses, the classical limit analysis approach is extended to compute the limit load multiplier under the simultaneous influence of joint flexibility, material and geometric nonlinearities, and limited ductility. Our formulation is an instance of the challenging class of optimization problems known as a mathematical program with equilibrium constraints (MPEC). Various nonlinear programming based algorithms are proposed to solve the MPEC. Finally, four numerical examples, concerning practical structures and benchmark cases, are provided to illustrate application of the analyses as well as to validate the accuracy and robustness of the proposed schemes.