INVESTIGATION OF ASSESSMENT METHODS FOR RAILWAY MASONRY ARCH BRIDGES

Similar documents
Chapter 5 Bending Moments and Shear Force Diagrams for Beams

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

The Wave Equation I. MA 436 Kurt Bryan

Shear and torsion interaction of hollow core slabs

63. Representation of functions as power series Consider a power series. ( 1) n x 2n for all 1 < x < 1

Module 1. Energy Methods in Structural Analysis

Design Data 1M. Highway Live Loads on Concrete Pipe

Module 2. Analysis of Statically Indeterminate Structures by the Matrix Force Method. Version 2 CE IIT, Kharagpur

Solution Manual. for. Fracture Mechanics. C.T. Sun and Z.-H. Jin

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

Review of Calculus, cont d

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!

99/105 Comparison of OrcaFlex with standard theoretical results

Chapter 6 Notes, Larson/Hostetler 3e

1.2. Linear Variable Coefficient Equations. y + b "! = a y + b " Remark: The case b = 0 and a non-constant can be solved with the same idea as above.

V. DEMENKO MECHANICS OF MATERIALS LECTURE 6 Plane Bending Deformation. Diagrams of Internal Forces (Continued)

CE 160 Lab 2 Notes: Shear and Moment Diagrams for Beams

The Algebra (al-jabr) of Matrices

Comparison of the Design of Flexural Reinforced Concrete Elements According to Albanian Normative

approaches as n becomes larger and larger. Since e > 1, the graph of the natural exponential function is as below

Part I: Basic Concepts of Thermodynamics

Improved Results on Stability of Time-delay Systems using Wirtinger-based Inequality

A. Limits - L Hopital s Rule. x c. x c. f x. g x. x c 0 6 = 1 6. D. -1 E. nonexistent. ln ( x 1 ) 1 x 2 1. ( x 2 1) 2. 2x x 1.

Advanced Calculus: MATH 410 Notes on Integrals and Integrability Professor David Levermore 17 October 2004

MECHANICS OF MATERIALS

APPROXIMATE INTEGRATION

1 Bending of a beam with a rectangular section

DETERMINATION OF MECHANICAL PROPERTIES OF NANOSTRUCTURES WITH COMPLEX CRYSTAL LATTICE USING MOMENT INTERACTION AT MICROSCALE

Unit #9 : Definite Integral Properties; Fundamental Theorem of Calculus

New Expansion and Infinite Series

Conservation Law. Chapter Goal. 5.2 Theory

Plates on elastic foundation

CBE 291b - Computation And Optimization For Engineers

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

Pre-Session Review. Part 1: Basic Algebra; Linear Functions and Graphs

Solutions to Supplementary Problems

Math 1B, lecture 4: Error bounds for numerical methods

Physics 116C Solution of inhomogeneous ordinary differential equations using Green s functions

ELE B7 Power Systems Engineering. Power System Components Modeling

Explain shortly the meaning of the following eight words in relation to shells structures.

DIRECT CURRENT CIRCUITS

13: Diffusion in 2 Energy Groups

Math 8 Winter 2015 Applications of Integration

Effects of peripheral drilling moment on delamination using special drill bits

MAC-solutions of the nonexistent solutions of mathematical physics

Operations with Polynomials

Mathcad Lecture #1 In-class Worksheet Mathcad Basics

Advanced Computational Analysis

Strong acids and bases. Strong acids and bases. Systematic Treatment of Equilibrium & Monoprotic Acid-base Equilibrium.

The steps of the hypothesis test

5.7 Improper Integrals

THE INTERVAL LATTICE BOLTZMANN METHOD FOR TRANSIENT HEAT TRANSFER IN A SILICON THIN FILM

Theme 8 Stability and buckling of members

Properties of Integrals, Indefinite Integrals. Goals: Definition of the Definite Integral Integral Calculations using Antiderivatives

Math& 152 Section Integration by Parts

ODE: Existence and Uniqueness of a Solution

Intro to Nuclear and Particle Physics (5110)

The Fundamental Theorem of Calculus

We know that if f is a continuous nonnegative function on the interval [a, b], then b

Minimum Energy State of Plasmas with an Internal Transport Barrier

12. DYNAMIC ANALYSIS. Force Equilibrium is Fundamental in the Dynamic Analysis of Structures 12.1 INTRODUCTION

Goals: Determine how to calculate the area described by a function. Define the definite integral. Explore the relationship between the definite

8.2 ESTIMATING DETENTION VOLUMES

Job No. Sheet 1 of 8 Rev B. Made by IR Date Aug Checked by FH/NB Date Oct Revised by MEB Date April 2006

STEP FUNCTIONS, DELTA FUNCTIONS, AND THE VARIATION OF PARAMETERS FORMULA. 0 if t < 0, 1 if t > 0.

Kirchhoff and Mindlin Plates

CHM Physical Chemistry I Chapter 1 - Supplementary Material

MATH , Calculus 2, Fall 2018

Modification Adomian Decomposition Method for solving Seventh OrderIntegro-Differential Equations

PhysicsAndMathsTutor.com

13.3. The Area Bounded by a Curve. Introduction. Prerequisites. Learning Outcomes

Numerical Analysis: Trapezoidal and Simpson s Rule

Quantum Physics II (8.05) Fall 2013 Assignment 2

ESCI 343 Atmospheric Dynamics II Lesson 14 Inertial/slantwise Instability

Evaluation of Allowable Hold Loading of, Hold No. 1 with Cargo Hold No. 1 Flooded, for Existing Bulk Carriers

Calculus - Activity 1 Rate of change of a function at a point.

Estimation of the particle concentration in hydraulic liquid by the in-line automatic particle counter based on the CMOS image sensor

Entropy ISSN

INVESTIGATION OF BURSA, ESKIKARAAGAC USING VERTICAL ELECTRICAL SOUNDING METHOD

Advanced Calculus: MATH 410 Uniform Convergence of Functions Professor David Levermore 11 December 2015

x = b a N. (13-1) The set of points used to subdivide the range [a, b] (see Fig. 13.1) is

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

The Properties of Stars

Dynamics of an Inertially Driven Robot

n f(x i ) x. i=1 In section 4.2, we defined the definite integral of f from x = a to x = b as n f(x i ) x; f(x) dx = lim i=1

INTRODUCTION. The three general approaches to the solution of kinetics problems are:

Lecture 23: Interpolatory Quadrature

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Chapter 4 Contravariance, Covariance, and Spacetime Diagrams

Fundamentals of Analytical Chemistry

A. Limits - L Hopital s Rule ( ) How to find it: Try and find limits by traditional methods (plugging in). If you get 0 0 or!!, apply C.! 1 6 C.

FUZZY HOMOTOPY CONTINUATION METHOD FOR SOLVING FUZZY NONLINEAR EQUATIONS

Math 124A October 04, 2011

New data structures to reduce data size and search time

Chapter 3 Exponential and Logarithmic Functions Section 3.1

Experimental Study, Stiffness of Semi-Rigid Beam-to-Column Connections Using Bolts and Angles

E S dition event Vector Mechanics for Engineers: Dynamics h Due, next Wednesday, 07/19/2006! 1-30

MA 124 January 18, Derivatives are. Integrals are.

MATH 144: Business Calculus Final Review

Transcription:

th Interntionl onference on Arch Bridges October -7, 06, Wrocłw, Polnd INVESTIGATION OF ASSESSMENT METODS FOR RAIWAY MASONRY AR BRIDGES J. Wng, J. ynes,. Melbourne University of Slford, Directorte of ivil Engineering, Slford, UNITED KINGDOM. e-mils: j.wng@slford.c.uk, b.j.hynes@slford.c.uk, c.melbourne@slford.c.uk SUMMARY Although msonry rch bridge reserch over the lst few decdes hs spwned severl dvnced cpcity ssessment tools, the mintennce nd ssessment of these bridges is still constnt concern for the bridge owners. There hve been mny ppels for the MEXE method to be replced but no greed lterntive hs yet been found. Though the development of new, simple nd robust cpcity ssessment is importnt, there is n immedite need to identify criticl issues in the current ssessment methods which my led to unsfe cpcity ssessments. This pper provides criticl pprisl of the MEXE method proposed in the current Rilwy odes. Detils of the specific concerns re estblished, justified nd ssessed. onclusions re drwn which re potentilly of significnt prcticl interest. Keywords: Msonry, rch bridges, MEXE, rilwy, ssessment method.. INTRODUTION In the UK, the current method of determining the lod crrying cpcity of rilwy msonry rch bridge is embodied in the guidnce note for the structurl ssessment of underbridges NR/GN/IV/0 [] herefter referred to s the Rilwy code. The first level in the ssessment is to pply the so-clled MEXE method. The MEXE method evolved from the work undertken by Pipprd in the 90 s which included both field nd lbortory tests to clibrte theoreticl work []. The full review of Pipprd s equtions is presented by Wng nd Melbourne []. According to the Rilwy code, the MEXE method for rilwy structures hs been developed differently to tht for highwy structures s defined in BD []. The Provisionl Ale pcity for rilwy msonry rch bridge cn be determined using the envelope curves or clculted using the formule presented in the Rilwy code. This pper presents derivtion of the bsic equtions used for the envelope curves in the Rilwy code. Anomlies re identified nd discussed. 0

Msonry rch bridges. ONSIDERING ONGITUDINA IVE OAD DISTRIBUTION Pipprd s originl equtions re bsed on the nlysis of two-pinned prbolic rch with spn/rise rtio of four, subjected to point live lod pplied t the crown, i.e. no live lod distribution is considered. When live lod distribution through the fill in the longitudinl direction is considered, s shown in Fig., the following equtions re derived. 0 ls A q B V A VB Fig.. Unit live lod distribution through fill in longitudinl direction... orizontl thrust nd bending moment ls P= q Y A X V A B A B = - = V B V A = q ls V B = q ls b Fig.. Two-hinged prbolic rch I =I 0 sec. For concentrted unit lod, s shown in Fig., the totl horizontl thrust nd bending moment t the crown re: M 0 For uniform distributed lod q, s shown in Fig. b, thrust nd bending moment due to the increse of length d re given by: 0

d q d q d d q d q dm 0 0 Using symmetry, where l s Then, 0. q d q 6 0 0. q d q M 7 For uniformly distributed lod where s ql, the totl horizontl thrust nd bending moment t the crown re given by: l s l M s 9.. orizontl thrust stress nd bending stress If PU represents the horizontl thrust stress nd PUB represents the bending stress t the crown due to the effect of the unit live lod s shown in Fig., then: bd l bd PU s 0 6 6 bd l bd M PUB s.. The provisionl le cpcity If the vilble live lod stress P represents subtrcting the ded lod stress from the limiting compressive stress f c t the crown, then the provisionl le cpcity is given by: 0 th Interntionl onference on Arch Bridges October -7, 06, Wrocłw, Polnd

Msonry rch bridges P PA PU PUB The full derivtion of the provisionl le cpcity for uniformly distributed lod cross the crown is presented elsewhere [].. MEXE METOD IN TE RAIWAY ODE.. Avilble live lod stress The mthemticl epression of vilble live lod stress presented in the Rilwy code is s follows: d h d d d P 00 kn/m The derivtion of the bove eqution is presented elsewhere [],[]... ive lod effects... ongitudinl loded length od distribution is considered in the longitudinl direction nd the Rilwy code offers equtions nd to obtin the longitudinl loded length. m m l s h 0. loded length due to single lelod l b h. loded length due to bogie lod These equtions suggest tht lods re distributed from the bse of the sleeper through the fill onto the rch brrel t slope of horizontlly to verticlly in the longitudinl direction, s shown in Fig.. Fig.. ongitudinl live lod distribution through fill. 06

th Interntionl onference on Arch Bridges October -7, 06, Wrocłw, Polnd... Influence line for bending moment The mthemticl epression of n influence line for bending moment t the crown is presented in the Rilwy code s: F 0.7 0. 0. 6 According to the Rilwy code, the eqution is obtined from close fit to the influence line curve for bending moment in the previous report nd ssessment code. Fig. shows comprison of this influence line nd tht constructed using eqution. Fig.. omprison of the bending moment influence lines. By integrting eqution 6 which clcultes the re under the influence line over the loded length, the following eqution is given in the Rilwy code A 0. 0.09 0. 7 The eqution is lter used to clculte the bending stress s presented in Section...... Influence line for horizontl thrust The mthemticl epression of n influence line for horizontl thrust t the crown is presented in the rilwy code s: T 0.7 07

Msonry rch bridges According to the Rilwy code, the eqution is n pproimtion to the influence line for horizontl thrust given in the previous report nd ssessment code. Fig. shows comprison of this influence line nd tht constructed using eqution for n rch with spn/rise rtio of. For n rch with spn/rise rtio of, the vlue of T is multiplied by to obtin the ctul thrust. Fig.. omprison of the horizontl thrust influence lines.... Section properties Section properties re given in the Rilwy code by the following epressions, d h. m h 0. m, B h. d m, M 9 6 d h. m 0. h 0.9 m, B h. d m, M 0 6 Where, B is defined s the cross sectionl re, M is defined s the section modulus, d is the thickness of rch ring t the crown, nd h is depth of fill between sleeper soffit nd rch ring t the crown. This suggests tht the effective width of the rch brrel dopted by the Rilwy code is h 0. m, b h. m 0. h 0.9 m, b h. m 0

th Interntionl onference on Arch Bridges October -7, 06, Wrocłw, Polnd... Unit live lod bending stress nd horizontl thrust stress The unit live lod bending stress is presented for single le lod in the Rilwy code by the following epression: PUB S A l M S For bogie lod, the subscript chnges to b, i.e. PUB b nd the loded length chnges to l b correspondingly. Replcing A with eqution 7, the unit live lod bending stress PUB S becomes PUB S 6 ls bd 0. 0.09 0. This suggests tht the bending moment t the crown in the Rilwy code is: M 0. 0.09 0. l S omprison of the results between this eqution nd those from Eq. 9, for unit spn rch, is s shown in Fig. 6. They produce lmost identicl moments s the loded length increses. Fig. 6. omprison of the bending moment t the crown. 09

Msonry rch bridges..6. Unit live lod horizontl thrust stress The unit live lod horizontl thrust stress is presented in the Rilwy code by the following epression: 0.9 T PU 6 B Where B is the cross sectionl re. This suggests tht the horizontl thrust t the crown in the Rilwy code is: 0.9 T 7 It is not cler how this horizontl thrust is obtined. owever, if the influence line in the Rilwy code, i.e. eqution, is used nd integrted over the loded length, different horizontl thrust eqution cn be obtined, which does produce the sme results s those from the bove eqution, i.e. eqution 7. Substituting T from eqution into eqution 7, the following epressions re obtined: 0.9 T 0.9 0.7 omprison of horizontl thrust obtined from this eqution nd tht from eqution spn/rise= is s shown in Fig. 7. Fig. 7. omprison of the horizontl thrust. 0

th Interntionl onference on Arch Bridges October -7, 06, Wrocłw, Polnd.. Provisionl le cpcity The provisionl le cpcity in the Rilwy code is obtined by dividing the vilble live lod stress by the sum of the unit live lod horizontl thrust nd bending stresses nd epressed by equtions 9 nd 0. For single le lod, P S 9 PU S PUBS For bogie lod, b PU 0.P 0 b PUB b The provisionl le cpcity is tken s the minimum of s or b nd both re represented grphiclly in Fig.. Fig.. Provisionl le cpcity for n rch with d=00mm nd h= 00mm.. ANOMAIES IN TE RAIWAY ODE.. od distribution in the longitudinl direction luse 6.. in the Rilwy code sttes tht lods should be distributed from the bse of the sleeper through the fill onto the rch brrel t slope of horizontlly to verticlly in the longitudinl direction. owever, the loded length i.e. equtions nd, dopted in the clcultions for producing the envelope curves re bsed on the slope of horizontlly to verticlly. For this inconsistency, typicl m spn rch with 00 mm

Msonry rch bridges crown ring thickness nd 00mm fill cover, hs provisionl le cpcity of 6. tonnes for : lod distribution, which is % higher thn tht for : lod distribution... Effective width The influence of effective width hs been presented elsewhere [7]. The rch cpcity envelope curves set out in the Rilwy code re bsed upon effective width ssumptions which re shown in Fig. 9, however, there ppers to be no logicl connection between these effective width equtions nd the descriptions given in Figure 6. of the Rilwy code. Fig. 9. Effective width of ril rch bridges... Provisionl le cpcity luse 6... in the Rilwy code sttes tht the provisionl le cpcities obtined from the envelope curves re likely to be conservtive becuse ssumptions bout lod dispersl lterlly nd longitudinlly re pessimistic. It goes on to suggest tht where higher cpcity is sought, there is scope to justify greter lod spred thn originlly ssumed. owever, when cpcity predictions using severl sources re compred, s shown in Fig., it is evident tht those using Pipprd s originl equtions suggest much lower cpcities thn the Rilwy code Single le nd Bogie. The cpcities from the Rilwy code drop to zero becuse the Rilwy MEXE method pplies Pipprd s compressive stress criteri i.e. 00 kn/m under the self-weight only.

th Interntionl onference on Arch Bridges October -7, 06, Wrocłw, Polnd Fig. 0. omprison of provisionl le cpcity d=00 mm nd h= 00 mm.. ONUSIONS.. onclusions The elstic cpcity ssessment produced by Pipprd is bsed upon well estblished nd ccepted ssumptions. To derive the envelope curves, the Rilwy code inconsistently pplies its own ssumptions bout: od distribution ngle through fill Effective width of the rch brrel Additionlly there is no check on the tensile stress level included in the Rilwy code. pcity predictions using the Rilwy code re significntly higher thn those obtined using Pipprd s equtions. This is most significnt t shorter spns... Recommendtions The equtions used in the Rilwy code re bsed upon best fit pproimtions, however these hve lost their physicl representtion nd it would be more intuitive to use the fundmentl reltionships upon which they re bsed. In order to estblish the rel dispersl of loding through fill, further physicl testing is required. This should ddress both longitudinl nd lterl lod spred through shllow nd deep bllst. In prticulr this must estblish useble dispersl pttern for deep fill where lod dispersl ptterns interfere, nd must be comptible with equilibrium. REFERENES [] NR/GN/IV/0. 06. The structurl ssessment of underbridges. Network Ril Guidnce Note. Network Ril.

Msonry rch bridges [] MKIBBINS,., MEBOURNE,., SAWAR, N. nd GAIARD,.S., Msonry rch bridges: condition pprisl nd remedil tretment. IRIA 66. 006 [] WANG, J. nd MEBOURNE,., Mechnics behind the MEXE method for msonry rch ssessment. Proceedings of the Institution of ivil Engineers, Engineering nd omputtionl Mechnics. Volume 6, EM, 00, pp.7 0. [] BD /0. The Assessment of highwy bridges nd structures. Design mnul for rods nd bridges, Volume, Section, Prt, ighwys Agency, 00 [] MEBOURNE,. nd WANG, J., Independent check on the proposed modifictions nd of limittions for the MEXE method, Improving Assessment, Optimistion of Mintennce nd Development of Dtbse for Msonry Arch Bridges UI Project I/0/U/, 00. [6] BA 6/97, The Assessment of highwy bridges nd structures. Design mnul for rods nd bridges, Volume, Section, Prt, ighwys Agency, 997 [7] WANG, J., AYNES, J. nd MEBOURNE,., A comprison between the MEXE nd Pipprd's methods of ssessing the lod crrying cpcity of msonry rch bridges, AR 7th Interntionl onference on Arch Bridges, roti, 0, pp.9 96