A method as claimed in claim 15, additionally comprising the step of maximising an applied load factor, ??, subject to BTt+??fL???q=???fD and NTq???g, where B is a suitable compatibility matrix, t represents a vector of equivalent internal nodal actions for all n nodes, where q is a vector of actions acting on the potential discontinuities, fL and fD are vectors representing live and dead loads, N