Fraddosio, Aguinaldo

Piccioni, Mario Daniele

Ricci, Eleonora

Sacco, Elio

Archivio della ricerca - Università degli studi di Napoli Federico II

AIMETA 2017 XXIII Conference The Italian Association of Theoretical and Applied Mechanics Luigi Ascione, Valentino Berardi, Luciano Feo, Fernando Fraternali and Antonio Michele Tralli (eds.) Salerno, Italy, 4-7September  2017 NUMERICAL METHODS FOR THE LOWER BOUND LIMIT ANALYSIS OF MASONRY ARCHES Eleonora Ricci1, Aguinaldo Fraddosio2, Mario Daniele Piccioni2, and Elio Sacco1 1 Department of Mechanical and Civil Engineering, University of Cassino and SL Via di Biasio, 43, Cassino, Italy  {e.ricci, sacco}@unicas.it 2Dipartimento di Scienze dell'Ingegneria Civile e dell'Architettura, Politecnico di Bari Via Orabona, 4, Bari, Italy {aguinaldo.fraddosio, mariodaniele.piccioni}@poliba.it Keywords: masonry, arches, limit analysis, numerical methods, point collocation method. Abstract. The present paper deals with the problem of stability of masonry arches. In particular, the problem is approached invoking the lower bound theorem of Limit Analysis; thus, the exis-tence of a thrust-line entirely contained in the thickness of the arch ensures that the arch does not collapse under the assigned load. With this aim, the Milankovitch theory [8]of the prob-lem of the equilibrium of the arches is provided in a general framework, regardless of the shape of the arch and of the nature of the applied loads. Here, in order to formulate the lower bound limit analysis problem in general context, the Milankovitch’s theory is reviewed, for-mulating the problem of the determination of the thrust-line in a form suitable for the imple-mentation in numerical procedures. In particular, the thrust curve is approximated by polynomial functions that are solved employing the Point Collocation Method [10]. Moreover, an optimization procedure is formulated for determining admissible equilibrium minimum and maximum thrust solutions. For the special case of a circular arch subjected to vertical load, the numerical procedure is assessed comparing the results obtained by the Collocation technique with the corresponding closed form solutions of the equilibrium problem. 15261 INTRODUCTION The growing interest in the preservation of heritage and historic constructions justifies the current research efforts on the development of efficient methods for structural analysis of arches and vaults, which are one of the fundamental and most fascinating parts of these con-structions [1, 2]. In particular, although the studies on the stability of masonry arches have a long story [6, 7], the problem of determining the load bearing capacity of these structures, es-pecially in the case of seismic loads, still is an open research topic. In this vein, Limit Analy-sis approaches are considered among the most effective tools for the structural assessment of masonry arches and vaults. Indeed, the failure of these structures is mainly related to their shape and not to the crushing of the masonry material in compression. The current formula-tion of the Limit Analysis for masonry constructions is mainly due to the fundamental contri-butions by Heyman [1]; the theory is based on the following three hypotheses: masonry has no tensile strength, infinite compressive strength, and sliding between masonry blocks is not possible. Under these hypotheses, the lower bound theorem of the Limit Analysis for masonry arches can be stated as follows: a masonry arch is safe if it is possible to find a line of thrust in equilibrium with external loads and completely lying within the thickness of the arch [6, 7]. In this paper, a numerical lower bound Limit Analysis approach for the study of masonry arches is proposed, starting from the Milankovitch’s general equilibrium theory [8]. We recall that Milankovitch provided a rigorous treatment of the equilibrium problem for masonry arches from both a mechanical and a mathematical point of view [9]. One of the fundamental results of this theory is the differential equation of the line of thrust. Here, we review the Mi-lankovitch’s theory of the line of thrust; since, depending on the geometry of the arch and on the loading conditions, closed form solutions can be prohibitive to find, we formulate the problem in a form suitable for the implementation in numerical procedures. In particular, we tackle numerically the problem by employing the Point Collocation Method (PCM) [10], a very handy and stable technique characterized by low computational costs. Finally, since it is possible to find infinite equilibrium lower bound solution, depending on the considered boun-dary conditions, we couple the PCM procedure with a constrained optimization routine, aimed at forcing the thrust line to be contained within the thickness of the arch, and at finding suita-ble lower bound solutions like, e.g., minimum and maximum thrust solutions. The effective-ness of the implemented method is discussed by numerical examples concerning the case of a circular arch subjected to the self-weight, for which analytical closed form solutions can also be found. 2 EQUILIBRIUM EQUATIONS Let us consider a masonry arch of generic shape and of constant depth t, whose middle plane is represented in Figure 1. In polar coordinates  ,O r   the geometry of the arch is as-signed by the mid-curve  mr  and by the thickness d , which in general can be variable. For each transversal section of the arch we denote by C the center of curvature of the mid-curve and by mR its radii of curvature; then, we name eR  and iR  the distances of the extrados line and the intrados line from C, with  e iR R d .Moreover, for the considered transversal sec-tion we denote by R  the distance of the thrust line from C and by   the angle between the transversal section and the vertical-axis. We restrict our attention on an arbitrary small part of the arch defined by an angular neighborhood d  of the transversal section; for this part, we consider a fully general loading condition, including: - the self-weight dg , applied in the center of mass G; 1527Eleonora Ricci, Aguinaldo Fraddosio, Mario D. Piccioni and Elio Sacco - horizontal inertial forces dp , also applied in G; - vertical and horizontal distributed load on the extrados, ef  and ep , respectively; - vertical and horizontal distributed load on the intrados, if  and ip , respectively. Let T  be the resultant of the forces transmitted by the rest of the arch to the right face of the considered part, and let be V  and H the vertical and horizontal components of T , respective-ly. Finally, let dT T  be the resultant of the forces transmitted by the rest of the arch to the left face of the considered part, with vertical and horizontal componentsV dV  and H dH , respectively.  Figure 1: Geometry and loading conditions of a small part of the arch.  In a Cartesian coordinate system O(x,y), by following an approach similar to that of Milan-kovitch [8], it is possible to show that the translational equilibrium in the horizontal and ver-tical directions and the rotational equilibrium around the point E of the thrust line (see Figure 1) are described by the following system of three ordinary differential equations (ODE):             2 22 222 222 21021021sin sin2 121sin cos2 12cos cosh e i e e i ie i e e i ie i e e e mmi i m i h e ime e e m i i m idH t R R p R p RddV t R R f R f Rddx dy dV H t R R f R R Rd d Rdf R R R t R RRp R R R p R R R                                               0, (1) where is the weight per unit volume of the material, h is an apparent density of the horizon-1528tal inertial forces, which may represent a peak seismic acceleration, and  is the distance be-tween the middle line and the thrust line.  In equations (1) the thrust line ( )R   with the vertical and horizontal thrust forces ( )V   and ( )H  , respectively, are the unknown functions to be determined; it is apparent that closed form solutions can be very hard to find, even in case of simplified problems. For this reason, numerical methods are needed. 3 THE CIRCULAR ARCH SUBJECTED TO SELF-WEIGHT For the purpose of the present paper, we specialize the analysis to the simple case of a cir-cular arch, which not only is relevant for applications but allows the evaluation of the analyti-cal solutions that represent a reference for the assessment of the numerical procedure. We consider the arch subjected only to the self-weight and to vertical loads acting on the extrados, thus, we assume:  0, 0, 0, 0.h i e i     f     p     p      (2) In this case, it is convenient to rewrite equations (1) in polar coordinates  ,R   whose origin is coincident with the (fixed) center of curvature C. After some simple manipulations, we get:      2 222 20102sin cos cos sin1sin sin 02 12e i e ee i m e e emdHddV t R R f RddR dRV V R H H Rd ddt R R R R f R R RR                          (3) representing a system of 3 ODE where the thrust line  R   and the vertical and horizontal thrusts  V   and  H   are the unknown functions. From equation (3)1 we see that the hori-zontal thrust  H   is constant, and we denote by 0H  this constant value. Integrating equa-tions (3)2 and (3)3, the vertical thrust  V   and the thrust line  R   can be determined; 0H  and the two integration constants in the solutions of (3)2 and (3)3 have to be determined by setting suitable boundary conditions. Let now further specialize the problem by considering a circular arch of constant thickness, subjected only to its self-weight, that is:         cost , 0, cost , cost .e i i e e        f         R R          R R             (4) By integrating equation (3)2, which now becomes uncoupled from equation (3)3, it results:   V q k    (5) where    2 212e iq t R R   (6) 1529Eleonora Ricci, Aguinaldo Fraddosio, Mario D. Piccioni and Elio Sacco and k is an integration constant. Given the symmetry of the problem with respect to the vertic-al axis, we can assume in the key section of the arch that  0 0V  ,which implies 0k  ; con-sequently, formula (5) becomes:  .V q  (7) Thus, equation (3)3 reduces to the following linear ODE in  R  :     20 0sin cos cos sin sin sin 0.12mmdR dq H q H q R q Rd R                 (8) Setting   00R R , the differential equation (8) admits the following closed form solution:    20 001 cos12.sin cosmmdH R q RRRq H        (9) which describes the thrust line and depends upon two integration constants 0R and 0H . Then, as it is expected, the equilibrium problem admits infinite solutions, depending on the assigned boundary conditions. 4 THE NUMERICAL APPROACH For more general loading conditions, and/or for arches with more complex shapes and eventually with variable thickness, analytical solutions of the ODE system (1) is prohibitive, and numerical methods have to be employed. In particular, we adopt the Point Collocation Method (PCM) [10], which is very effective for the numerical solution of boundary value problems like that under examination, and is characterized by a reduced computational cost.  PCM consists in approximating the unknown function  R   with a polynomial of degree v , and in exactly satisfying the differential equation in a discrete number m of points, called collocation points.  In particular, let us consider the differential equation (8) in the unknown function  R  , defined on the domain ;2 2      , function also of the unknown constant 0H . The domain is divided in en intervals by introducing 1d en n   points. For the i th interval, the function    iR  is approximated as follows:       41i ij jjR r   (10) where  ijr  are the coefficients of the polynomial and j are interpolating functions, satisfying the following conditions: 1530                                         1 1 1 21 1 1 22 1 2 22 1 2 23 1 3 23 1 3 24 1 4 24 1 4 21, 0, 0, 00, 1, 0, 00, 0, 1, 00, 0, 0, 1i ii ii ii ii ii ii ii id dj=1      d dd dj=2      d dd dj=3      d dd dj=4      d d                                         (11) where 1i  and  2i are the initial and final angles of the i th interval. For each element, expression (10)is substituted into the ODE (8)that have to be exactly sa-tisfied at the following three collocation points:                          1 1 2 12 1 2 13 1 2 11 15,2 101,21 15.2 10i i i ii i i ii i i icpcpcp                       (12) Therefore, starting from the differential equation (8) a set of 3en  linear algebraic equations in the polynomial coefficients is obtained; moreover, further 1en   algebraic equations coming from the requirement of the continuity of the approximating functions when passing from an element to the next are considered:         1 12 1 .i ii iR R        (13) Finally, if we also assign boundary conditions on  R  ,we obtain a system of 4en  linear algebraic equations in 4en   unknowns (the coefficients in the functions   iR  ); notice that at this stage the thrust 0H  is undetermined. Now, in the spirit of the lower bound theorem of Limit Analysis, the thrust line must be contained in the thickness of the arch. Moreover, if the boundary conditions are not a priori prescribed (indeed, in practical problems, they are unknown), we can search between the infi-nite possible equilibrium solutions, those that are relevant: usually, the minimum horizontal thrust and the maximum horizontal thrust solutions. To do this, we add to the numerical pro-cedure an optimization routine, which solve the following constrained optimization problem:                 0 0min max.i iiei iefind  H   or  HR Rsuch that         with i=1,...,nR R    (14) The above approach has been applied for finding the thrust line of a circular arch subjected 1531Eleonora Ricci, Aguinaldo Fraddosio, Mario D. Piccioni and Elio Sacco to the self-weight; the arch has external radius eR =1250 mm, internal radius iR =950 mm, width t =500 mm and specific weight =2000 daN/m3. The domain  / 2, / 2   has been subdivided in 8en   intervals (angular amplitude / 8   ). For what concerns the maxi-mum thrust solution, Figure 2 shows a comparison between the analytical results (dotted line), obtained from formula (9) by setting 0 iR R  and  / 2 eR R  , and the numerical results (solid line). Notice that the two curves are practically identical; moreover, also the analytical-ly determined horizontal thrust on the spring, equal to 297.58 daN, is the same of that numeri-cally determined, resulting equal to 297.06 daN. The same applies for the minimum thrust solution: by formula (9) we obtained a horizontal thrust on the springing equal to 131.10 daN, practically the same of that numerically determined, equal to 130.57 daN. The analytical (dot-ted line) and the numerical (solid line) thrust lines are compared in Figure 3.  Figure 2: Maximum thrust solution: analytical (dotted line) and numerical (solid line) results.   Figure 3: Minimum thrust solution: analytical (dotted line) and numerical (solid line) results.  5 CONCLUSIONS   We revised the Milankovitch's theory of the equilibrium of arches.  We proposed a numerical procedure for solving the differential equilibrium equations of a masonry arch, based on the Point Collocation Method. By coupling this procedure with an optimization routine, it is possible to find equilibrium lower bound solutions relevant for practical applications.  The method has been successfully applied for determining the maximum and minimum thrust of a circular arch subjected to the self-weight.  1532 The method is simple, stable and requires low computational costs. REFERENCES  [1] M. Como, Statics of Historic Masonry Constructions, 2nd ed..Springer Series in Solid and Structural Mechanics, 2016. [2] P. Roca, M. Cervera, G. Gariup, L. Pelà, Structural analysis of masonry historical con-structions, Classical and Advanced Approaches, Archives of Computational Methods in Engineering, 17(3), 299-325, 2010. [3] E. Benvenuto, An Introduction to the History of Structural Mechanics, Part II: Vaulted Structures and Elastic Systems. Springer-Verlag, 1991. [4] P. Block, M. De Jong, J. Ochsendorf, As Hangs the Flexible Line: Equilibrium of Ma-sonry Arches, Nexus Network Journal, 8(2), 13-24, 2006. [5] J. Heyman, Arches, vaults and buttress: masonry structures and their engineering. Variorum, 1996. [6] E. Sacco, Some Aspects on the Statics of Masonry Arches. Masonry Structures: Be-tween Mechanics and Architecture (D. Aita, O. Pedemonte, K. Williams (Editors)). Springer International Publishing Switzerland, 2015, pp. 265-290. [7] E. Sacco, Stress approaches for the analysis of masonry arches. Proceedings of 8th In-ternational Conference on Structural Analysis of Historical Constructions (SAHC 2012) (Editor: Jerzy Jasienko), Wroclaw (Polonia), 15-17 October 2012. [8] M. Milankovitch, Theorie der Drukccurven, Zeitschriftfür Mathematik und Physik, 55, 1-27, 1907. [9] F. Foce, Milankovitch’s Theorie der Druckkurven: Good mechanics for masonry archi-tecture, Nexus Network Journal, 9(2), 185-210, 2007. [10] A. Iserles, A first course in the numerical analysis of differential. Cambridge University Press, 1996.  1533

Numerical methods for the lower bound limit analysis of masonry arches

https://www.iris.unina.it/retrieve/handle/11588/718154/196066/Ricci.Fraddosio.Piccioni.Sacco.pdf

Numerical methods for the lower bound limit analysis of masonry arches

Abstract

Similar works

Full text

Available Versions

Archivio della ricerca - Università degli studi di Napoli Federico II