Computational & Technology Resources
an online resource for computational,
engineering & technology publications 

CivilComp Proceedings
ISSN 17593433 CCP: 77
PROCEEDINGS OF THE NINTH INTERNATIONAL CONFERENCE ON CIVIL AND STRUCTURAL ENGINEERING COMPUTING Edited by: B.H.V. Topping
Paper 104
The Computational Efficiency of Two Rigid Block Analysis Formulations for Application to Masonry Structures H.M. Ahmed and M. Gilbert
Department of Civil and Structural Engineering, University of Sheffield, United Kingdom H.M. Ahmed, M. Gilbert, "The Computational Efficiency of Two Rigid Block Analysis Formulations for Application to Masonry Structures", in B.H.V. Topping, (Editor), "Proceedings of the Ninth International Conference on Civil and Structural Engineering Computing", CivilComp Press, Stirlingshire, UK, Paper 104, 2003. doi:10.4203/ccp.77.104
Keywords: limit analysis, masonry arch bridges, rigid block, linear programming, joint equilibrium, redundant forces.
Summary
Rigid block analysis is a computational limit analysis method which is now
widely applied to the analysis of masonry gravity structures such as arch bridges.
Various formulations of the rigid block analysis method have been proposed, with
linear programming (LP) generally being used in the solution process. However,
there appears to be little information in the literature on the relative computational
efficiencies of the various formulations (when used either with traditional simplex or
newer interior point LP solvers). This paper aims to at least partially address this, by
considering two alternative formulations.
According to the classification of LP formulations presented by Tam and Jennings for plastic frame design [1], Livesley's original rigid block analysis formulation [2] may be classified as a `joint equilibrium' approach. In this approach the problem variables comprise both the interface forces and the resultant forces acting on each and every block. Conversely, Gilbert and Melbourne's formulation for single and multispan arch bridges [3] may be classified as a `redundant forces' approach (though the dual form of this was actually presented). In this approach the problem variables are both the interface forces and (only) the redundant forces. In general the redundant forces method will produce LP tableaux containing fewer variables and constraints than the joint equilibrium method. However, in the case of multiring masonry arch analysis problems, a problem with the redundant forces method is that the presence of interring contacts leads to numerous redundant forces (i.e. additional forces at every contact point on the interface between rings). Hence the resulting tableaux, although compact, will be very heavily populated with nonzero elements; this can be quite computationally expensive to solve. A recent reappraisal of the situation has also indicated that modern LP methods (i.e. interiorpoint based LP methods) may in fact be well suited to the solving the type of large sparse tableaux produced by the joint equilibrium method. This is investigated in the paper. Thus `redundant forces' and `joint equilibrium' formulations are described and then tested by applying them to a number of case study problems (these comprise: multiring arch, multispan arch and block wall problems). It is found that whilst the redundant forces formulation is the most computationally efficient formulation for simple single and multispan arch problems, the joint equilibrium formulation is the most computationally efficient formulation when applied to more complex geometrical arrangements of blocks, such as those found in multiring arches. This is because of the low proportion of nonzero elements present in the albeit large LP tableau generated when using the joint equilibrium formulation. This sparseness greatly speeds the solution process when using most modern linear programming solvers, particularly those based on interior point methods. Furthermore, using the joint equilibrium formulation it is also relatively easy to automatically set up and then solve problems involving arbitrary two and three dimensional assemblages of rigid blocks. However, it is well known that adoption of the standard LP limit analysis assumption of normality in the presence of sliding at frictional interfaces [i.e. using an associated flow rule] can potentially lead to unsafe estimates of the load factor. Whilst in many cases the overestimate of the computed load factor when friction is modelled in this way will be small, unfortunately this cannot be guaranteed, with the magnitude of the overestimate generally being highly problem dependent. Thus when certain geometrical configurations of blocks are considered the overestimate may be quite high. For example, in the case of a wall problem documented in the paper, a load factor was obtained which was some 20 percent higher than that obtained when using a nonassociative constitutive model [4]. Unfortunately modelling nonassociative friction in rigid block analysis is difficult and the methods proposed to date have tended to rely on the use of specialised solution algorithms (not LP), with only relatively small problems being tractable at the present time. The authors are therefore currently working on the development of more efficient solution methods for such problems. References
purchase the fulltext of this paper (price £20)
go to the previous paper 
