Numerical analysis of block stone displacements in ancient masonry structures: A new method to estimate historic ground motions
|
|
- Christina Bryan
- 5 years ago
- Views:
Transcription
1 INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Int. J. Numer. Anal. Meth. Geomech. 2008; 32: Published online 15 November 2007 in Wiley InterScience ( Numerical analysis of block stone displacements in ancient masonry structures: A new method to estimate historic ground motions Ronnie Kamai 1,2 and Yossef H. Hatzor 1,, 1 Department of Geological and Environmental Sciences, Ben-Gurion University of the Negev, Beer-Sheva, Israel 2 Department of Civil and Environmental Engineering, University of California, Davis, CA, U.S.A. SUMMARY An innovative approach is presented, in which the discontinuous deformation analysis (DDA) method is used to estimate historic ground motions by back analysis of unique structural failures in archaeological sites. Two archaeological sites in Israel are investigated using this new approach and results are presented in terms of displacement evolution of selected structural elements in the studied masonry structure. The response of the structure is studied up to the point of incipient failure, in a mechanism similar to the one observed in the field. Structural response is found to be very sensitive to dynamic parameters of the loading function such as amplitude and frequency. Prior to back analysis of case studies, two validations are presented. Both compare the performance of DDA with analytical solutions and present strong agreement between the two. Using comprehensive sensitivity analyses, the most likely peak ground acceleration (PGA) and frequency that must have driven the observed block displacements are found for the two case studies the Nabatean city of Mamshit and the medieval fortress of Nimrod in southern and northern Israel, respectively. It is found that horizontal peak ground accelerations (HPGA) of 0.5g and 1g were required to generate the observed deformations in Mamshit and Nimrod, respectively. Although these might seem too high, considering structural and topographic amplifications it is concluded that the analyses suggest ground motions of 0.2g at a frequency of 1.5 Hz for Mamshit and up to 0.4g at a frequency of 1 Hz for Nimrod. These values provide constraints on the seismic risk associated with these regions as appears in the local building code using a completely independent approach. Copyright q 2007 John Wiley & Sons, Ltd. Received 31 August 2006; Revised 15 August 2007; Accepted 15 October 2007 KEY WORDS: ground motions; seismic risk; PGA; discontinuous deformation analysis; masonry structures; validation Correspondence to: Yossef H. Hatzor, Department of Geological and Environmental Sciences, Ben-Gurion University of the Negev, Beer-Sheva 84105, Israel. hatzor@bgu.ac.il Copyright q 2007 John Wiley & Sons, Ltd.
2 1322 R. KAMAI AND Y. H. HATZOR 1. INTRODUCTION The expected ground motions at a site, a critical component in seismic risk assessments, are generally evaluated using empirical relationships between magnitude, distance, and site effects [1, 2]. Here, we present a quantitative approach to estimate historic ground motions using back analysis of block displacements in archaeological masonry structures. Reliable estimates of paleoseismic ground motions at a given site can provide constraints on ground motion predictions which are typically obtained using empirical methods. The analysis focuses on man-made masonry structures such as arches, where the hewn stones that form the building provide a well-defined initial geometrical reference. Where failure is confined to displaced blocks within an otherwise intact structure, block displacement is measurable and a mechanical analysis is possible; this cannot be achieved in completely collapsed structures. Additionally, sound bedrock foundation and high-quality masonry workmanship enhance the reliability of the results by eliminating other failure modes associated with poor construction methods or weak foundation material [3]. The suggested method is tested using two archeological sites in Israel: The Nabatean city of Mamshit in the Negev desert and the Ayyûbid-Crusader Fortress of Nimrod, at the foot of Mt. Hermon (Figure 1). For each site, a sensitivity analysis is performed to determine the most likely amplitude and frequency of motion that must have generated the observed failure in the field. The numerical analysis is performed using the discontinuous deformation analysis (DDA) method, a member of the family of implicit DEM methods [5] developed by Shi during the late 1980s [6]. Here, we use a research version of DDA software modified by DDA author Dr Gen-hua Shi for the specific targets of this research program. The DDA method solves a finite element type of mesh, where all elements are real isolated blocks and the unknowns of the equations are the displacements and deformations of the blocks. The blocks are simply deformable, namely stresses and strains are constant throughout the block. When the blocks are in contact, Coulomb s friction law applies to the contact interface, and the simultaneous equilibrium equations are formulated and solved for each loading or time increment. The formulation is based on the minimization of the system potential energy, following the second law of thermodynamics. One of the strengths of the DDA method, which is extremely important for this research, is that the mode of failure of the block system is a result of the analysis and not an assumption. This feature enables us to study the response of the structure to different loading functions and to isolate the parameters that produce the unique failure that is observed in the field. Over the last decade, researchers in the DDA community have dedicated a great deal of effort to prove the accuracy of the method by performing theoretical validation studies. MacLaughlin and Doolin [7] review more than 100 validation studies with respect to analytical solutions, laboratory and field data, and other numerical techniques. DDA performance was found to be more than adequate for engineering applications. Up to date, DDA has been used for numerous rock engineering projects. In Israel alone, DDA has been applied to solve dynamic stability issues in the highly discontinuous Masada rock slopes [8] to stabilize a discontinuous overhanging cliff [9] and to determine the stability of historic underground excavations in discontinuous rock [10, 11]. Unlike these applications, which involve potentially hazardous rock masses, this work uses the DDA method to examine the stability of man-made masonry structures.
3 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1323 Figure 1. Location map of studied sites (after [4]) shown on two sides of the seismically active Dead Sea transform (DST).
4 1324 R. KAMAI AND Y. H. HATZOR 2. NUMERICAL CODE VALIDATION Validation studies in this section were performed for calibration purposes only, using both existing and newly developed analytical solutions. The validations are fully dynamic, meaning that the velocity of each block at the end of a time step is fully transferred to the next time step with no energy dissipation (in DDA notation k01= 1) Dynamic sliding of a block on an incline The case of a sliding block on an incline subjected to gravitational loads only has been successfully analyzed and discussed in detail in previous papers [12], and hence will not be repeated here. The case of a sliding block subjected to dynamic loading, however, is addressed here in detail. The analysis presents higher accuracy than previously described [13] due to a correction that has been made here to the analytical solution. A displacement-based sliding block model was first proposed by Newmark [14] as well as Goodman and Seed [15], now largely referred to as Newmark -type analysis. Determination of the amount of displacement during an earthquake necessarily involves two steps [15]: (1) determination of horizontal acceleration required to initiate down slope motion, also known as yield acceleration (a y ), which can be found by pseudo-static analysis, and (2) Evaluation of the displacement developed during time intervals when yield acceleration is exceeded, by double integration of the acceleration time history, with the yield acceleration used as reference datum. Figure 2 displays a single block resting on a plane inclined at an angle α with friction along the interface, and subjected to gravitational acceleration g and a horizontal, time-based, sinusoidal acceleration input a as driving forces. Goodman and Seed [15] showed that for this case the yield acceleration is given by a y =tan( α)g (1) g(cosα-ksin( ωt)sin α)*tg φ a=kgsin( ωt) gsin α+kgsin( ωt)cosα g Figure 2. Schematic presentation of the forces acting on a single block on an incline. The driving forces are gravity and a horizontal, time-based, sinusoidal acceleration. The resisting frictional force is constant and depends only on the friction angle (μ=tan ). α
5 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1325 For an acceleration record of the form a =kg sin(ϖt), where ω corresponds to the frequency of the function and k calibrates the proportion between a and g, the corresponding time interval θ until yield acceleration is attained is given by θ= sin 1 (a y /kg) ϖ (2) The down slope acceleration of the sliding block can be determined by subtracting the resisting forces from the driving forces, as described in Figure 2 a t =[kg sin(ϖt) cos α+g sin α] [g cos α kg sin(ϖt) sin α]tan (3) Similarly, the displacement of the block at any time is determined by double integration on the acceleration, with θ as reference datum, in a conditional manner: if a>a y or v>0 t d t = v = a = g[(sinα cosαtan )(t 2 /2 θt)] (4) θ θ + ag [(cosα+sinαtan )(ϖ cos(ϖθ)(t θ) sin(ϖt)+sin(ϖθ))] ϖ2 otherwise d t =d t 1 Equation (4) provides the analytical solution for the dynamic displacement of a block on an inclined plane with inclination α and friction angle, starting from rest and subjected to a sinusoidal loading function with frequency ω. Figure 3 displays a comparison between the analytical and DDA solutions for three friction angles ( =22, 30, and 35 ) with θ=0.035,0.18, and 0.28 s, respectively. The accumulated displacements are computed for a plane inclination of 20. The numerical error E N showninthelowerpanelof Figure 3 is defined as E N = d d N d 100% (5) where d and d N are the analytical and the numeric displacement vectors, respectively. is the norm operator, which for a 2D displacement vector is d = u 2 +v 2 (6) Other than for the first second, in which the relative error reaches 50%, due to a very small absolute error of only 2.6E10 5 m, the relative error is less than 1% for all three friction angles for all accumulated displacements. Time-step size is kept constant in all DDA runs and is s.
6 1326 R. KAMAI AND Y. H. HATZOR Figure 3. Validation of the dynamic case of a block on an incline for three different interface friction angles: top input acceleration function; center comparison between analytical (solid line) and DDA (symbols) solutions; and bottom relative error for each simulation Block response to induced displacements in the foundation In the previous section, a validation was performed for a block on an incline subjected to dynamic input applied to its center of mass. We now present a validation for the case of a block resting on a horizontal interface and responding to dynamic input applied to the underlying block. Such an application would be suitable for dynamic analysis of structures which respond to earthquake motions at their foundations. This approach was not adopted in this current paper because application of this loading mechanism for cases that involve multiple blocks seems to require further research. Nevertheless, we feel it is instructive to present this original validation here because this is the first time such a validation is performed and the obtained results are promising. The studied block system consists of three blocks: a fixed foundation block (no. 0), the induced block (no. 1), and the responding block (no. 2) (Figure 4). Block 1 is subjected to a horizontal displacement input function in the form of a cosine, starting from 0 (shown in the upper panel of Figure 5) d(t)= D(1 cos(2πωt)) (6 )
7 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1327 y Block 2 x Block 1 d(t) Block 0 Figure 4. The DDA block system and sign convention used for the case of a single block resting on an underlying block which is subjected to time-dependant displacement. Input Displacement (m) D=0.3 m D=0.5 m D=1 m 2D DDA DDA DDA Analytic Analytic Analytic Input Motion Displacement of upper block (m) relative error (%) Time (sec) 2 Figure 5. The responding block motion (Figure 4) for three different input displacement amplitudes, all at f =1Hz: top input displacement function of the form d(t)= D(1 cos(2πωt)) applied to Block 1; center comparison between analytical (solid line) and DDA (symbols) solutions; and bottom relative error for each simulation. The dynamic response of Block 2 is investigated and DDA results are compared with the analytical solution, the essentials of which are provided below.
8 1328 R. KAMAI AND Y. H. HATZOR The only force acting on Block 2 other than gravity is the frictional force, which immediately determines the acceleration of Block 2: m 2 a 2 = F friction (7) m 2 a 2 =μ m 2 g (8) a 2 =μ g (9) where μ is the friction coefficient. The direction of the driving force is determined by the direction of the relative velocity between Block 1 and 2 (v1 ). When Block 1 moves to the right relative to Block 2, the frictional force pulls Block 2 in the same direction, and determines the sign of a 2. When Block 2 is at rest in relation to the Block 1, the frictional force is determined by the acceleration of the bottom block (a 1 ). The threshold acceleration, under which the two blocks move in harmony, is equal to the friction coefficient multiplied by the gravitation acceleration (μg). When the acceleration of Block 1 passes the threshold value, the frictional forces act in the same direction as a 1. The positive direction is determined by the sign convention in Figure 4, and the relative velocity of Block 1 is given by v 1 =v 1 v 2 (10) The direction of the acceleration of Block 2 is set by the following boundary conditions: if v1 =0 and a 1 <μg, a 2 =a 1 and a 1 >μg and a 1 >0, a 2 =μg and a 1 <0, a 2 = μg if v1 =0 and v 1 >0, a 2 =μg and v1 <0, a 2 = μg (11) In the solution of Equation (11) we employed Matlab 7.0 [16] software package because the analytical solution must be computed iteratively as the relative velocity and the direction of the force are dependent on each other. In order to compare DDA and the analytical solution, the mode of failure of the analyzed block in DDA has to be constrained to horizontal sliding only, so that it would simulate a single degree of freedom system as does the analytical solution. Therefore, Block 2 of the DDA block system was made very flat, to eliminate possible rotational motions (Figure 4). Figures 5 and 6 present the accumulated displacement of Block 2 and the relative error between analytical (Matlab) and DDA solutions, using a set of sensitivity analyses for amplitude and friction coefficient. Figure 5 presents the response of Block 2 to changing displacement amplitudes (D),
9 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1329 Figure 6. The responding block motion (Figure 4) for three different friction coefficients, all at D = 0.5m, f =1Hz: top the input displacement function applied to Block 1; center comparison between analytical (solid line) and DDA (symbols) solutions; and bottom relative error for each simulation. with a constant input frequency of 1 Hz and a friction coefficient of 0.6. The obtained magnitude of accumulated displacement is directly proportional to the input amplitude, as expected. Note that the three displacement curves follow the periodic behavior of the input displacement function (T =1s), and that divergence between curves starts after 0.25s, where the input displacement function has an inflection point. The relative error, as shown in the upper panel of Figure 5, is mostly between 1 and 2%, reaching a maximum of 3.5% for D =0.3m, which actually corresponds to a difference of 1 cm only in the displacements. Figure 6 presents the response of Block 2 to changing friction coefficients across the interface, with a constant displacement function of D =0.5m, f =1Hz. Note that the accumulated displacement is directly proportional to the friction coefficient up to 0.5 s, after which the input displacement function changes direction. After that point, the accumulated displacement for μ=0.6 is larger than for μ=1, since the higher friction works in both forward and backward directions. It is interesting to note that the relative errors (lower panel of Figure 6) show large sensitivity to the friction coefficient. For a reasonable friction coefficient of say 0.6, the relative error remains within 2%. However, when the friction coefficient is very low (μ=0.1) or relatively high (μ=1), the errors increase up to 12%. Note that the errors for the high friction follow the displacement
10 1330 R. KAMAI AND Y. H. HATZOR pattern, with the maximum errors closely corresponding to change in the direction of motion, where the system experiences high accelerations and inertia forces. In general, it is shown that DDA follows the analytical results in both cases with a remarkable agreement, both for changing friction coefficients and for changing motion amplitudes. It is important to note that while application of this loading mechanism yielded good agreement for a single responding block as shown above, when we experimented with problems involving multiple blocks we obtained a complex structural response not supported by field observations. We feel that while this loading mechanism is more appropriate for the problem at hand, further research of its numerical implementation is required. 3. RESULTS We will present our results in the form of sensitivity analyses for amplitude and frequency of the input loading function. As mentioned previously, the input acceleration function is applied simultaneously at the centroid of each block in the mesh, simulating therefore structural accelerations. The relevant ground motion corresponding to this excitation mode is discussed in the next section. Since the final objective of the analysis is to successfully mimic an existing failure in the field, the sensitivity analyses must be performed in trial and error: an artificial, sinusoidal acceleration record is applied and the responses of the masonry structures are studied up to the point of incipient failure, in a mode similar to the one observed in the field. Note that the magnitude of the measured displacement in the field could be a result of structural response to several, consecutive, large earthquakes over the life span of the studied structure, which in our case is in the order of years. The analysis herein assumes that the observed failure is a result of a single shaking event, in which all shaking characteristics have come together to create the unique pattern of deformation. Therefore, in addition to mimicking the mode of failure, the results are also discussed in terms of the magnitude of displacement, where applicable. In order for the analysis to be as accurate as possible, genuine mechanical properties of the analyzed structures are used as input for DDA. Original building stones were taken from the archaeological sites, under supervision of the Israel Nature and Parks Authority, to the Rock Mechanics Laboratory of the Negev at Ben-Gurion University. Lab tests were performed in order to obtain physical and mechanical properties of intact rock samples and results are summarized in Table I. The input stiffness of the embedding wall is reduced by four orders of magnitude in order to simulate the soft filling material between stones in the wall. Higher stiffness values for the confining wall were investigated, and it was found that the arch rocking behavior is extremely limited at high confining wall stiffness, restricting the initiation of the unique failure mechanisms observed in the field. Choosing a low value, in the order of soil stiffness, allows for large deformations under low stresses. Consequently, the ground motion that is claimed to have caused the deformations is a lower bound on the expected ground motions at that site. A friction angle of 35 for the stone interfaces was chosen as a representative value for all analyses. It was found that the response is not very sensitive to the friction angle, and that friction angles between 25 and 40 yielded keystone displacements of cm, respectively, with the same unique failure mechanism. All numerical analyses are performed with kinetic damping of 1% (k01= 0.99), namely, the initial velocity in each time step is 1% lower than the terminal velocity of the previous step. The introduction of kinetic damping in the analysis stems from previous work that suggests that with an increasing number of blocks in the mesh, a certain degree of energy dissipation is required for
11 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1331 Table I. Mechanical properties obtained in the lab from the original building blocks taken from the studied sites. Mechanical property Mamshit Nimrod Fortress Lithology and formation Hazera limestone Hermon limestone Density (kg/m 3 ) Porosity (%) Dynamic Young s modulus (GPa) Dynamic Poisson s ratio Dynamic shear modulus (GPa) Point load index (MPa) Uniaxial compressive strength (MPa) Peak interface friction angle ( ) Owing to sampling difficulties at NF, these parameters are derived from tests performed on blocks sampled at the Avdat site [17] with similar rock properties. Figure 7. The deformed arch at Mamshit: (a) general view of the arch which is embedded in a very heterogenic wall and (b) the keystone has slid 4 cm downwards while the rest of the arch remained intact. realistic modeling [8, 18], because DDA does not consider inelastic energy dissipation mechanisms such as crushing at block ends, out-of-plane rotations, heat generation, etc The case studies Two archaeological sites in Israel, where a confined structural failure is identified, measured and mapped are selected for analysis (Figure 1). Mamshit is the last Nabatean city built in the Negev on the trade route between Petra, Hebron, and Jerusalem [19, 20]. A unique structural failure is noticed in a tower at the corner of the Eastern Church, where a key stone has slid approximately 4 cm downwards out of a still standing semi-circular arch (Figure 7). The numerical analysis of the arch at Mamshit is performed on a block system (Figure 8) that contains the arch, with its accurate measurements from the field, confined by a uniform masonry wall. The arch is intersected by material lines which enable us to assign separate sets of mechanical parameters for the wall and the arch. These simulate the great difference between the
12 1332 R. KAMAI AND Y. H. HATZOR Figure 8. The DDA block system for the embedded arch at Mamshit. The modeled masonry wall rests on two fixed blocks. The lines intersecting the arch blocks represent material lines and a measurement point (circle) is assigned for the keystone. hewn stones of the arch itself and the heterogenic confining wall material (Figure 7) by assigning intact rock stiffness values to the arch stones (E arch =17GPa) and soil-like stiffness values to the wall stones (E wall =1MPa). A measurement point is assigned at the keystone, and its vertical displacement versus time is plotted in the relevant figures. In Nimrod fortress, the largest medieval fortress in Israel [21, 22], evidence of destruction caused by a severe earthquake can be seen throughout the fortress. The damage is dated to one of two large earthquakes of the year 1759 (Ms=6.6 and 7.4) [23]. The most impressive evidence from the earthquake that has damaged the fortress can be observed in the gate tower. There, three parallel arches as well as an arched passageway display the same mode of failure, in which a single block, not necessarily the keystone, has slid downwards while the rest of the arch remained intact (Figure 9). In the modeled block system (Figure 10), material lines intersect the arch blocks and allow us to assign separate sets of mechanical parameters to the arch and the wall, as before. The wall in which the arch is embedded is confined by a boundary block on its right-hand side and its base, added in order to simulate the geographical constraints at the site: the fortress is built on an elongated range, cut off by deep canyons from three sides. As a result, the gate tower is supported by a 30 m high retaining wall on its western side, while the eastern side rests directly on bedrock. This setting may certainly cause asymmetric structural response and possibly seismic wave amplifications in preferred orientations. Five measurement points are assigned to five blocks in the arch, the label and position of which are shown in Figure 11. During shaking, the entire modeled structure is gradually displaced to the left, most probably due to the confining wall on the right-hand side which represents a fixed boundary. In order to represent correctly the computed displacement of block A R (Figure 9(a)), the block that exhibits the greatest displacement in the field in relation to the computed displacements of the entire block system, the selected parameter for sensitivity analyses was the relative displacement of A R with respect to the average displacement of the other four keystones. The following equations describe the derivation of that parameter: ui u avg = 4, v vi avg = (12) 4 u = u(a R ) u avg, v =v(a R ) v avg (13)
13 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1333 Figure 9. Evidence for a destructive earthquake that hit the Nimrod Fortress most probably at 1759 (see text for explanation): (a) the same stone has slid downwards out of three parallel arches at the main gate arch; (b) a similar deformation is noticed at a passage-way parallel to the gate arch; (c) an adjacent arch facing east was not deformed; and (d) location of the passage-way, 10 m west of main gate. where u versus v is the relative spatial displacement of block A R, plotted in the relevant figures Sensitivity analyses A very interesting and unique behavior of the two structures is revealed through the sensitivity analyses of structural response to changing amplitude of motion. Instead of an intuitive proportion between induced motion amplitude and degree of damage, measured here by the magnitude of the displacement vector, there is apparently a structural preference to a certain range of amplitudes:
14 1334 R. KAMAI AND Y. H. HATZOR Figure 10. The DDA block system for the studied arch at Nimrod Fortress. Four fixed points (squares) are assigned to the confining block and five measurement points (circles) are assigned to the top arch blocks. The lines intersecting the arch blocks represent material lines. Figure 11. Sign convention for the Nimrod arch blocks: K =keystone, A and B for the first and second block from the keystone respectively, and R or L to indicate right or left. in both Mamshit and Nimrod, the studied block exhibits the greatest displacement under specific amplitudes, not necessarily the largest. Figure 12 shows that the keystone of the arch at Mamshit exhibits the greatest downward displacement under an amplitude of 0.5g, when everything else is kept equal. Downward displacement increases with acceleration amplitude up to an amplitude of 0.5g. When the acceleration amplitude is greater than 0.5g, the keystone response exhibits strong fluctuations and even a shift in displacement direction. A similar phenomenon is observed in Nimrod. Figure 13 shows the spatial displacement of block A R in response to changing amplitudes. Under the smallest amplitude presented (A =0.8g), a very minor response is detected where the block displaces inwards and later outwards, resulting in almost no accumulated displacement. The largest amplitude (A = 1.5g), however, has a strong affect on the entire structure leading to partial destruction, but consequently block A R is locked
15 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1335 Time (sec) Key Block Displacement (cm) g 0.32 g 0.5 g 1g 0.8 g Figure 12. Influence of input motion amplitude on vertical keystone displacement at Mamshit ( f = 1Hz). Figure 13. Influence of input motion amplitude on relative spatial displacement of Block A R at Nimrod Fortress ( f =2Hz). after 16 s of shaking. A continuous, accumulating inward displacement at a constant rate and direction is achieved only when A =1g. The inward displacement of block A R decelerates only after 70 s when further displacement is locked, as shown in Figure 14. Note that also in the real event the structure must have locked at a certain point, and the measured displacement that was preserved matches the modeled displacement with A =1g. Structural response to frequency of motion is usually discussed in terms of the natural period of the structure, at which the structure will develop a resonance mode and collapse. Since the studied failures are local and not complete, each mode of failure will have its own natural period which can be different than that of the whole structure. The structural sensitivity to frequency, revealed in our sensitivity analyses, is significant considering that the common terminology for seismic risk evaluation uses mainly PGA (peak ground acceleration) and largely ignores frequency.
16 1336 R. KAMAI AND Y. H. HATZOR Figure 14. Normalized displacement of stone A R in the direction of the displacement vector over time (d = u 2 +v 2 ). After approximately 70 s, the displacement stops and the structural deformation is locked. Time (sec) Key Block Displacement (cm) Hz 5Hz 1Hz 15*Nuweiba 1.5 Hz 10 Hz Figure 15. Influence of input motion frequency on vertical displacement of keystone at Mamshit ( A = 0.5g). In the case of the keystone at Mamshit (Figure 15), a clear preference is detected for frequencies in the range of Hz; only under those input frequencies the downward displacement of the keystone is continuous, and accumulates more than 3 cm of displacement, similar to the amount of displacement measured in the field. Lower or higher frequencies result in other modes of failure such as oscillations in the case of low frequencies (e.g. 0.5 Hz) and locking of the structure in the case of higher frequencies (e.g Hz). A very important observation is made when a real earthquake record, that of Nuweiba 1995 (Figure 16), is used as input motion. The original record, de-convoluted to rock response [8], is amplified by 15 in order for its PGA to reach the same amplitude of the synthetic motions used for the results that are plotted in Figure 15, namely 0.5g. Although the Nuweiba quake loads the structure with a wide range of frequencies and with two simultaneous components of motion (horizontal and vertical), the structural response to the natural
17 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS g(z) g(ew) g(ns) a (g) time (sec) Figure 16. The Nuweiba 1995 record after de-convolution to rock response. The rectangle marks the 10 s that were used for the analysis of the Mamshit block system. Figure 17. Influence of frequency on relative spatial displacement of Block A R at Nimrod Fortress (A =1g). quake is very similar to that of the simple sinusoidal ones, though more moderate. This finding strongly suggests that the results of the sensitivity analysis, using synthetic records of horizontal motion only, are valid enough to be further discussed. A similar observation is made in Figure 17, which shows that the mode of failure of the arch at Nimrod Fortress has a clear preference to a frequency of 2 Hz over both 1 and 3 Hz. 4. DISCUSSION 4.1. Limitations of our approach The suggested method uses a 2D numerical model with simple deformable blocks, meaning that stresses and strains are identical everywhere in the block. Although the mode of failure in DDA is a result of the analysis and not an assumption, the 2D constraint imposes limitations on the mode of
18 1338 R. KAMAI AND Y. H. HATZOR failure, for example, it does not allow for out-of-plane rotations, etc. To minimize this constraint, the analysis was performed on planar structures such as arches, where most of the deformation is concentrated in the plane of analysis. An analysis of a rectangular structure, such as a tower, could have been significantly biased by the 2D constraint. Moreover, the simple deformable block assumption may cause small deviations in the energy involved in the deformation, since in reality stress concentrations at block corners may result in block cracking or chipping, a mechanical process that cannot be modeled in this version of DDA. We assume, though, that an energy dissipation of 1%, which is implemented through the numerical kinetic damping, is enough to compensate for such losses [18]. Another limitation of this study is the use of a synthetic sinusoidal loading function as direct input, consisting of a single amplitude and frequency value. This is obviously not the case in a real earthquake which consists of a train of frequencies at different amplitudes. Nevertheless, the results presented in Figure 15 indicate that the structural response in this case does not vary all that much when a real earthquake record is used as input Applicability of our results The obtained PGA and frequency values in this research can be compared with Israel Building Code # 413 [2], generated largely on empirical grounds by the Geophysical Institute of Israel. IBC 413 predicts PGA values of 0.11g for Mamshit and 0.25g for Nimrod Fortress. Although our results are significantly higher, it can be argued that at least two amplification factors can explain the discrepancies: structural and topographic. The loading function in the numerical analysis is inserted into all the blocks simultaneously, namely, the analyzed block is subjected to a uniform structural input motion. The building code however relates to bedrock PGA, which may be amplified significantly at the top of the structure due to inherent structural amplifications as well as rock structure interactions. Common response spectra predict amplification of up to 2.5 [24], due to wave propagation phenomena that take place when the motion is transferred from bedrock to the structure. Therefore, the PGA values found by inversion in our study must be reduced by a factor of 2.5, resulting in bedrock PGA predictions of g and 0.4 1g for the Mamshit and Nimrod sites, respectively. No further modifications should be made for the PGA value obtained for Mamshit, because topographic effects are not expected in this site. In Kalaat Nimrod, a further modification should be made due to the possible topographic effect at that site, the precise magnitude of which is not known and must be measured in the field. In the absence of such data we can only suggest that the actual PGA value for the rock foundations at Kalaat Nimrod should be lower than 0.4g by up to a factor of 4, if we adopt the results of the topographic site effect measurements at four different sites in Israel [25]. 5. SUMMARY AND CONCLUSIONS DDA method was validated with respect to analytical solutions for two separate cases of a Newmark sliding block type of analysis: in one, a block is subjected to time-dependent input acceleration applied at its center of mass, and in the other, the block responds to a time-dependent input displacement applied to an underlying block. In both simulations, DDA results matched the analytical solutions very well.
19 NUMERICAL ANALYSIS OF BLOCK STONE DISPLACEMENTS 1339 The structural deformation pattern observed in the two case studies was duplicated very nicely with DDA, with localized failure of particular blocks in the mesh-matching field observations. Both deformation pattern and displacement magnitude were duplicated accurately in two different structural geometries and in two different sites. For the semi-circular arch in Mamshit, DDA results imply that most of the damage took place during the first two seconds of motion. Our best estimate for the structural acceleration parameters for Mamshit is f =1Hz, A =0.5g. Considering structural amplifications these results can be downscaled to f =1Hz, A =0.2g. The obtained amplitude is similar to PGA=0.11g predicted by Israel seismic building code for that region. Interestingly, introduction of vertical motions using a real earthquake record does not change the obtained deformation picture very significantly, justifying our simplified methodology of using synthetic sinusoidal input motions which are restricted to the horizontal axis only. For the asymmetric arch at Nimrod our best estimate for the structural acceleration parameters is f =2Hz, A =1g. Considering both structural and topographic amplifications, the actual range of the paleoseismic PGA values can be between 0.1g and 0.4g. Assuming structural amplification of 2.5 and topographic amplification for this elevated site of 2.0, the expected paleoseismic PGA is 0.2g. This estimation can be compared with Israel seismic building code prediction of 0.25g for this region. Further research is required concerning failure modes associated with application of the loading function at the foundation of a multi-block structure. The application of this loading mechanism will generate true structural amplifications at the top of the structures and will enable more rigorous conclusions regarding paleoseismic PGA values. ACKNOWLEDGEMENTS We thank Gen-Hua Shi for modifying his latest version of the DDA code for the purpose of this research, Yuli Zaslavski for the de-convoluted Nuweiba earthquake record, Shmuel Marco and Itai Leviathan for discussions, and two anonymous reviewers for critical reading and instructive comments that greatly improved this paper. REFERENCES 1. Boore DM, Joyner WB, Fumal TE. Equations for estimating horizontal response spectra and peak acceleration from western North American earthquakes: a summary of recent work. Seismological Research Letters 1997; 68: Shapira A. Explanations to the PGA Map for the Israeli Building Code 413. Available from: Mazor E, Korjenkov A. In Applied Archeoseismology, Mazor E, Krasnov B (eds), vol. 4. The Makhteshim Country Laboratory of Nature, Pensoft: Moscow, 2001; Hall JK. Digital Shaded-relief Map of Israel and Environs 1:500,000. I. G. Survey, Jing L. A review of techniques, advances and outstanding issues in numerical modelling for rock mechanics and rock engineering. International Journal of Rock Mechanics and Mining Sciences 2003; 40: Shi G-H. Block System Modeling by Discontinuous Deformation Analysis. Computational Mechanics Publication: Southampton, U.K., MacLaughlin MM, Doolin DM. Review of validation of the discontinuous deformation analysis (DDA) method. International Journal for Numerical and Analytical Methods in Geomechanics 2006; 30(4): Hatzor YH, Arzi AA, Zaslavsky Y, Shapira A. Dynamic stability analysis of jointed rock slopes using the DDA method: King Herod s Palace, Masada, Israel. International Journal of Rock Mechanics and Mining Sciences 2004; 41(5):
20 1340 R. KAMAI AND Y. H. HATZOR 9. Tsesarsky M, Hatzor YH, Levithan I, Saltzman U, Sokolowsky M. Structural control on the stability of overhanging, discontinuous rock slopes. Fortieth U.S. Symposium on Rock Mechanics (USRMS): Rock Mechanics for Energy, Mineral and Infrastructure Development in the Northern Regions, Anchorage, AK, Hatzor YH, Benary R. The stability of a laminated Voussoir beam: back analysis of a historical roof collapse using DDA. International Journal of Rock Mechanics and Mining Sciences 1998; 35(2): Hatzor YH, Tsesarsky M, Eimermacher RC. In Structural Stability of Historic Underground Openings in Rocks: Two Case Studies from Israel, Kouroulis SK (ed.). Springer: Berlin, 2006; Sitar N, MacLaughlin MM, Doolin DM. Influence of kinematics on landslide mobility and failure mode. Journal of Geotechnical and Geoenvironmental Engineering 2005; 131(6): Hatzor YH, Feintuch A. The validity of dynamic block displacement prediction using DDA. International Journal of Rock Mechanics and Mining Sciences 2001; 38(4): Newmark N. Effects of earthquakes on dams and embankments. Geotechnique 1965; 15(2): Goodman RE, Seed HB. Earthquake-induced displacements in sand embankments. Journal of Soil Mechanics and Foundation Division (ASCE) 1966; 90(SM2): MATLAB. MATLAB, version 7. The MathWorks, Inc., Natick, MA. 17. Kamai R. Estimation of historical seismic ground-motions using back analysis of structural failures in archaeological sites. M.Sc. Thesis, The Department of Geological and Environmental Sciences, Ben-Gurion University of the Negev, Beer-Sheva, 2006; 127. Available from: Ronnie MSc thesis.pdf. 18. Tsesarsky M, Hatzor YH, Sitar N. Dynamic displacement of a block on an inclined plane: analytical, experimental and DDA results. Rock Mechanics and Rock Engineering 2005; 38(2): Negev A. The Architecture of Mampsis. Quedem, Monographs of the Institute of Archaeology, The Hebrew University of Jerusalem, Jerusalem, 1988; Negev A. The Nabatean Cities in the Negev. Ariel: Jerusalem, Hartal M. The Al-Subayba (Nimrod) Fortress, Towers 11 and 9. IAA, Israel Antiquities Authority: Jerusalem, Ellenblum R. Who Built Qalat Al-Subayba? Dumbarton Oaks Papers 1989; 41: Ambraseys NN, Barazangi M. The 1759 Earthquake in the Bekka valley: implications for earthquake hazard assessment in the eastern Mediterranean region. Journal of Geophysical Research 1989; 94(B4): Stewart JP, Kim S, Bielak J, Dobry R, Power MS. Revisions to soil structure interaction procedures in NEHRP design provisions. Earthquake Spectra 2003; 19(3): Zaslavsky Y, Shapira A. Experimental study of topographic amplification using the Israel seismic network. Journal of Earthquake Engineering 2000; 4(1):43 65.
Site Response Analysis with 2D-DDA
Site Response Analysis with 2D-DDA Yossef H. Hatzor Sam and Edna Lemkin Professor of Rock Mechanics Dept. of Geological and Environmental Sciences Ben-Gurion University of the Negev, Beer-Sheva, Israel
More informationTHE VOUSSOIR BEAM REACTION CURVE
THE VOUSSOIR BEAM REACTION CURVE Yossef H. Hatzor Ben-Gurion University, Department of Geological and Environmental Sciences Beer-Sheva, Israel, 84105 ABSTRACT: The influence of joint spacing (s) on the
More informationEarthquake response analysis of rock-fall models by discontinuous deformation analysis
c Earthquake response analysis of rock-fall models by discontinuous deformation analysis T. Sasaki, I. Hagiwara & K. Sasaki Rock Engineering Lab., Suncoh Consultants Co. Ltd., Tokyo, Japan R. Yoshinaka
More informationKINETIC EEFCT ON FLEXIBLE BODIES BEHAVIOR
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2 Paper No. 99 KINETIC EEFCT ON FLEXIBLE BODIES BEHAVIOR Eduardo BOTERO 1 and Miguel P. ROMO 2 SUMMARY When a flexible
More informationRisk assessment of collapse in shallow caverns using numerical modeling of block interactions with DDA: Suggested approach and case studies
Risk assessment of collapse in shallow caverns using numerical modeling of block interactions with DDA: Suggested approach and case studies Yossef H. Hatzor Professor and Chair in Rock Mechanics, BGU,
More informationDynamic modelling in slopes using finite difference program
Bulletin of the Department of Geology Bulletin of the Department of Geology, Tribhuvan University, Kathmandu, Nepal, Vol. 12, 2009, pp. 89 94 Dynamic modelling in slopes using finite difference program
More informationHatzor, Y.H. and Glaser, S., Environmentally controlled, multi scale, dynamic behavior of rock masse SUMMARY... 3
Hatzor, Y.H. and Glaser, S., Environmentally controlled, multi scale, dynamic behavior of rock masse TABLE OF CONTENT SUMMARY... 3 CHAPTER 1 - A HYBRID geodfn-dda PREPROCESSOR FOR SIMULATING MECHANICALLY
More informationDynamic Soil Pressures on Embedded Retaining Walls: Predictive Capacity Under Varying Loading Frequencies
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 2015 Christchurch, New Zealand Dynamic Soil Pressures on Embedded Retaining Walls: Predictive Capacity Under Varying Loading
More informationKinematics of Overhanging Slopes in Discontinuous Rock
Kinematics of Overhanging Slopes in Discontinuous Rock Michael Tsesarsky, Ph.D. 1 ; and Yossef H. Hatzor, Ph.D., M.ASCE 2 Abstract: The kinematics of overhanging rock slopes and the mechanical constraints
More informationTechnical Note Dynamic Displacement of a Block on an Inclined Plane: Analytical, Experimental and DDA Results
Rock Mech. Rock Engng. (2005) 38 (2), 153 167 DOI 10.1007/s00603-004-0043-2 Technical Note Dynamic Displacement of a Block on an Inclined Plane: Analytical, Experimental and DDA Results By M. Tsesarsky
More informationLimit analysis of brick masonry shear walls with openings under later loads by rigid block modeling
Limit analysis of brick masonry shear walls with openings under later loads by rigid block modeling F. Portioli, L. Cascini, R. Landolfo University of Naples Federico II, Italy P. Foraboschi IUAV University,
More informationAPPENDIX J. Dynamic Response Analysis
APPENDIX J Dynamic Response Analysis August 25, 216 Appendix J Dynamic Response Analysis TABLE OF CONTENTS J1 INTRODUCTION... 1 J2 EARTHQUAKE TIME HISTORIES... 1 J3 MODEL AND INPUT DATA FOR SITE RESPONSE
More informationSHAKE TABLE STUDY OF SOIL STRUCTURE INTERACTION EFFECTS ON SEISMIC RESPONSE OF SINGLE AND ADJACENT BUILDINGS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1918 SHAKE TABLE STUDY OF SOIL STRUCTURE INTERACTION EFFECTS ON SEISMIC RESPONSE OF SINGLE AND ADJACENT
More informationCOMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM
COMPARISON BETWEEN 2D AND 3D ANALYSES OF SEISMIC STABILITY OF DETACHED BLOCKS IN AN ARCH DAM Sujan MALLA 1 ABSTRACT The seismic safety of the 147 m high Gigerwald arch dam in Switzerland was assessed for
More informationAN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS
AN IMPORTANT PITFALL OF PSEUDO-STATIC FINITE ELEMENT ANALYSIS S. Kontoe, L. Pelecanos & D.M. Potts ABSTRACT: Finite Element (FE) pseudo-static analysis can provide a good compromise between simplified
More informationRealistic dynamic analysis of jointed rock slopes using DDA
Plenary lectures Stability of Rock Structures, Hatzor (ed) 22 Swets & Zeitlinger, Lisse, ISBN 9 589 519 3 Realistic dynamic analysis of jointed rock slopes using DDA Y.H. Hatzor, A.A. Arzi, and M. Tsesarsky
More informationRESPONSE SPECTRA RECOMMENDED FOR AUSTRALIA
RESPONSE SPECTRA RECOMMENDED FOR AUSTRALIA Malcolm Somerville, Kevin McCue and Cvetan Sinadinovski Australian Geological Survey Organisation, Canberra SUMMARY Response spectra suitable for intraplate regions
More informationDynamic Soil Structure Interaction
Dynamic Soil Structure Interaction Kenji MIURA, Dr. Eng. Professor Graduate School of Engineering Hiroshima University Dynamic Soil Structure Interaction Chapter 1 : Introduction Kenji MIURA, Dr. Eng.
More informationThe Transmission Gully Project: Method for Assessing Seismic Stability of Large Rock Cuts In Steep Terrain, Considering Topographic Amplification
6 th International Conference on Earthquake Geotechnical Engineering 1-4 November 2015 Christchurch, New Zealand The Transmission Gully Project: Method for Assessing Seismic Stability of Large Rock Cuts
More informationProceedings oh the 18th IAHR International Symposium on Ice (2006) DISCRETE ELEMENT SIMULATION OF ICE PILE-UP AGAINST AN INCLINED STRUCTURE
DISCRETE ELEMENT SIMULATION OF ICE PILE-UP AGAINST AN INCLINED STRUCTURE Jani Paavilainen, Jukka Tuhkuri and Arttu Polojärvi Helsinki University of Technology, Laboratory for Mechanics of Materials, Finland
More informationGravity dam and earthquake
Gravity dam and earthquake Tardieu s Dynamic simplified method Patrick LIGNIER, Tractebel Engineering Coyne et Bellier Château des Comtes de Challes 9 octobre 2014 CONTENTS 2 Vulnerability of gravity dam
More informationSEISMIC HAZARD ANALYSIS. Instructional Material Complementing FEMA 451, Design Examples Seismic Hazard Analysis 5a - 1
SEISMIC HAZARD ANALYSIS Instructional Material Complementing FEMA 451, Design Examples Seismic Hazard Analysis 5a - 1 Seismic Hazard Analysis Deterministic procedures Probabilistic procedures USGS hazard
More informationEvaluation of dynamic behavior of culverts and embankments through centrifuge model tests and a numerical analysis
Computer Methods and Recent Advances in Geomechanics Oka, Murakami, Uzuoka & Kimoto (Eds.) 2015 Taylor & Francis Group, London, ISBN 978-1-138-00148-0 Evaluation of dynamic behavior of culverts and embankments
More informationDynamic Analysis Contents - 1
Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2
More informationDynamic Response of EPS Blocks /soil Sandwiched Wall/embankment
Proc. of Second China-Japan Joint Symposium on Recent Development of Theory and Practice in Geotechnology, Hong Kong, China Dynamic Response of EPS Blocks /soil Sandwiched Wall/embankment J. C. Chai 1
More informationEvaluating the Seismic Coefficient for Slope Stability Analyses
Evaluating the Seismic Coefficient for Slope Stability Analyses by Edward Kavazanjian, Jr., Ph.D., P.E.,D.GE., NAE Ira A. Fulton Professor of Geotechnical Engineering School of Sustainable Engineering
More informationInfluences of material dilatancy and pore water pressure on stability factor of shallow tunnels
Influences of material dilatancy and pore water pressure on stability factor of shallow tunnels YANG Xiao-li( ), HUANG Fu( ) School of Civil and Architectural Engineering, Central South University, Changsha
More informationEmbedment Depth Effect on the Shallow Foundation - Fault Rupture Interaction
Embedment Depth Effect on the Shallow Foundation - Fault Rupture Interaction M. Ashtiani & A. Ghalandarzadeh Faculty of Civil Engineering, University of Tehran, Iran SUMMARY: The 1999 earthquakes in Turkey
More informationEffect of buttress on reduction of rock slope sliding along geological boundary
Paper No. 20 ISMS 2016 Effect of buttress on reduction of rock slope sliding along geological boundary Ryota MORIYA *, Daisuke FUKUDA, Jun-ichi KODAMA, Yoshiaki FUJII Faculty of Engineering, Hokkaido University,
More informationAnalytical and Numerical Investigations on the Vertical Seismic Site Response
Analytical and Numerical Investigations on the Vertical Seismic Site Response Bo Han, Lidija Zdravković, Stavroula Kontoe Department of Civil and Environmental Engineering, Imperial College, London SW7
More informationBack Analysis of the Lower San Fernando Dam Slide Using a Multi-block Model
Proceedings Geohazards Engineering Conferences International Year 2006 Back Analysis of the Lower San Fernando Dam Slide Using a Multi-block Model C. A. Stamatopoulos P. Petridis Stamatopoulos and Associates
More informationEmbedded Foundation with Different Parameters under Dynamic Excitations
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 2287 Embedded Foundation with Different Parameters under Dynamic Excitations Jaya K P 1 and Meher Prasad
More informationINELASTIC RESPONSES OF LONG BRIDGES TO ASYNCHRONOUS SEISMIC INPUTS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 24 Paper No. 638 INELASTIC RESPONSES OF LONG BRIDGES TO ASYNCHRONOUS SEISMIC INPUTS Jiachen WANG 1, Athol CARR 1, Nigel
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More informationOn seismic landslide hazard assessment: Reply. Citation Geotechnique, 2008, v. 58 n. 10, p
Title On seismic landslide hazard assessment: Reply Author(s) Yang, J; Sparks, A.D.W. Citation Geotechnique, 28, v. 58 n. 1, p. 831-834 Issued Date 28 URL http://hdl.handle.net/1722/58519 Rights Geotechnique.
More informationTunnelling and Underground Space Technology 21 (2006) 29 45
Tunnelling and Underground Space Technology 21 (26) 29 45 Tunnelling and Underground Space Technology incorporating Trenchless Technology Research www.elsevier.com/locate/tust Tunnel roof deflection in
More informationSeismic Slope Stability
ISSN (e): 2250 3005 Volume, 06 Issue, 04 April 2016 International Journal of Computational Engineering Research (IJCER) Seismic Slope Stability Mohammad Anis 1, S. M. Ali Jawaid 2 1 Civil Engineering,
More informationCOMPARATIVE STUDY OF LINEAR-ELASTIC AND NONLINEAR- INELASTIC SEISMIC RESPONSES OF FLUID-TANK SYSTEMS
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1127 COMPARATIVE STUDY OF LINEAR-ELASTIC AND NONLINEAR- INELASTIC SEISMIC RESPONSES OF FLUID-TANK SYSTEMS
More informationNumerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation
Numerical Modeling of Interface Between Soil and Pile to Account for Loss of Contact during Seismic Excitation P. Sushma Ph D Scholar, Earthquake Engineering Research Center, IIIT Hyderabad, Gachbowli,
More informationShake Table Study of Soil Structure Interaction Effects in Surface and Embedded Foundations
Shake Table Study of Soil Structure Interaction Effects in Surface and Embedded Foundations Naghdali Hosseinzadeh Structural Engineering Research Center, International Institute of Earthquake Engineering
More informationNUMERICAL ANALYSIS OF PASSIVE EARTH PRESSURES WITH INTERFACES
III European Conference on Computational Mechanics Solids, Structures and Coupled Problems in Engineering C.A. Mota Soares et.al. (eds.) Lisbon, Portugal, 5-8 June 2006 NUMERICAL ANALYSIS OF PASSIVE EARTH
More informationLandslide FE Stability Analysis
Landslide FE Stability Analysis L. Kellezi Dept. of Geotechnical Engineering, GEO-Danish Geotechnical Institute, Denmark S. Allkja Altea & Geostudio 2000, Albania P. B. Hansen Dept. of Geotechnical Engineering,
More informationStructural control on the stability of overhanging, discontinuous rock slopes
ARMA/USRMS 5-77 Structural control on the stability of overhanging, discontinuous rock slopes Tsesarsky, M. and Hatzor Y. H Department of Geological and Environmental Sciences. Ben Gurion University of
More informationSeismic stability analysis of quay walls: Effect of vertical motion
Proc. 18 th NZGS Geotechnical Symposium on Soil-Structure Interaction. Ed. CY Chin, Auckland J. Yang Department of Civil Engineering, The University of Hong Kong, Hong Kong. Keywords: earthquakes; earth
More informationSLOPE STABILITY EVALUATION AND ACCEPTANCE STANDARDS
INFORMATION BULLETIN / PUBLIC - BUILDING CODE REFERENCE NO.: LABC 7006.3, 7014.1 Effective: 01-01-2017 DOCUMENT NO.: P/BC 2017-049 Revised: 12-21-2016 Previously Issued As: P/BC 2014-049 SLOPE STABILITY
More informationNumerical Modelling of Dynamic Earth Force Transmission to Underground Structures
Numerical Modelling of Dynamic Earth Force Transmission to Underground Structures N. Kodama Waseda Institute for Advanced Study, Waseda University, Japan K. Komiya Chiba Institute of Technology, Japan
More informationA modified model of a single rock joint s shear behavior in
This paper is accepted for publication in the International Journal of Mining Science and Technology A modified model of a single rock joint s shear behavior in limestone specimens Dindarloo Saeid R a*,
More informationWeak Rock - Controlling Ground Deformations
EOSC 547: Tunnelling & Underground Design Topic 7: Ground Characteristic & Support Reaction Curves 1 of 35 Tunnelling Grad Class (2014) Dr. Erik Eberhardt Weak Rock - Controlling Ground Deformations To
More informationThe Dynamic Response Analysis of Concrete Gravity Dam under the Earthquake
Copyright 2013 Tech Science Press SL, vol.9, no.1, pp.23-36, 2013 The Dynamic Response Analysis of Concrete Gravity Dam under the Earthquake Yang Lu 1, Li Shi-Min 1, Cao Peng 2 and Shen Xin-Pu 1 Abstract:
More informationSeismic Design of a Hydraulic Fill Dam by Nonlinear Time History Method
Seismic Design of a Hydraulic Fill Dam by Nonlinear Time History Method E. Yıldız & A.F. Gürdil Temelsu International Engineering Services Inc., Ankara, Turkey SUMMARY: Time history analyses conducted
More informationNON-LINEAR ANALYSIS OF SOIL-PILE-STRUCTURE INTERACTION UNDER SEISMIC LOADS
NON-LINEAR ANALYSIS OF SOIL-PILE-STRUCTURE INTERACTION UNDER SEISMIC LOADS Yingcai Han 1 and Shin-Tower Wang 2 1 Fluor Canada Ltd., Calgary AB, Canada Email: yingcai.han@fluor.com 2 Ensoft, Inc. Austin,
More informationNONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS
October 12-17, 28, Beijing, China NONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS J. Huang 1 and A. Zerva 2 1 Dept. of Civil, Architectural and Environmental
More informationON THE PREDICTION OF EXPERIMENTAL RESULTS FROM TWO PILE TESTS UNDER FORCED VIBRATIONS
Transactions, SMiRT-24 ON THE PREDICTION OF EXPERIMENTAL RESULTS FROM TWO PILE TESTS UNDER FORCED VIBRATIONS 1 Principal Engineer, MTR & Associates, USA INTRODUCTION Mansour Tabatabaie 1 Dynamic response
More informationMEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 1 : THEORY SUBJECT CODE: COMRMC MODERATOR: H YILMAZ EXAMINATION DATE: OCTOBER 2017 TIME:
MEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 1 : THEORY EXAMINER: WM BESTER SUBJECT CODE: COMRMC EXAMINATION DATE: OCTOBER 2017 TIME: MODERATOR: H YILMAZ TOTAL MARKS: [100] PASS MARK: (60%)
More informationNONLINEAR SEISMIC SOIL-STRUCTURE (SSI) ANALYSIS USING AN EFFICIENT COMPLEX FREQUENCY APPROACH
NONLINEAR SEISMIC SOIL-STRUCTURE (SSI) ANALYSIS USING AN EFFICIENT COMPLEX FREQUENCY APPROACH Dan M. GHIOCEL 1 ABSTRACT The paper introduces a novel approach for modeling nonlinear hysteretic behavior
More informationDynamic sliding of tetrahedral wedge: The role of interface friction
INTERNATIONAL JOURNAL FOR NUMERICAL AND ANALYTICAL METHODS IN GEOMECHANICS Published online 19 January 211 in Wiley Online Library (wileyonlinelibrary.com). DOI: 1.12/nag.19 Dynamic sliding of tetrahedral
More informationEarthquakes and Earth s Interior
- What are Earthquakes? Earthquakes and Earth s Interior - The shaking or trembling caused by the sudden release of energy - Usually associated with faulting or breaking of rocks - Continuing adjustment
More informationNumerical analysis of effect of mitigation measures on seismic performance of a liquefiable tailings dam foundation
Numerical analysis of effect of mitigation measures on seismic performance of a liquefiable tailings dam foundation Yong-Beom Lee, Jorge Castillo Ausenco, USA Aurelian C. Trandafir Fugro GeoConsulting
More informationThe Seismic Performance of Tousheh Dam During the Chi-Chi Earthquake
( C023) Proceedings of 9 th Conference on Current Researches in Geotechnical Engineering, Shihman Reservoir, Tai-Yuan, Taiwan, R.O.C. August 30-3 and September, 200 92 () (the semi-analysis-testing method)(2)
More informationStructural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake
Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake Seung Keun Park and Hae Sung Lee ABSTRACT This paper presents a system identification (SI) scheme
More informationModel tests and FE-modelling of dynamic soil-structure interaction
Shock and Vibration 19 (2012) 1061 1069 1061 DOI 10.3233/SAV-2012-0712 IOS Press Model tests and FE-modelling of dynamic soil-structure interaction N. Kodama a, * and K. Komiya b a Waseda Institute for
More informationMethods of Interpreting Ground Stress Based on Underground Stress Measurements and Numerical Modelling
University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2006 Methods of Interpreting Ground Stress Based on Underground Stress Measurements and
More informationApplication of pseudo-symmetric technique in dynamic analysis of concrete gravity dams
Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams V. Lotfi Department of Civil and Environmental Engineering, Amirkabir University, Iran Abstract A new approach is
More informationDeterministic Seismic Hazard Assessment of Quetta, Pakistan
Deterministic Seismic Hazard Assessment of Quetta, Pakistan M.A. Shah Micro Seismic Studies Programme, Islamabad, Pakistan Pakistan Institute of Engineering and Applied Sciences, Islamabad, Pakistan M.
More informationON THE FACE STABILITY OF TUNNELS IN WEAK ROCKS
33 rd 33 Annual rd Annual General General Conference conference of the Canadian of the Canadian Society for Society Civil Engineering for Civil Engineering 33 e Congrès général annuel de la Société canadienne
More informationWave Propagation through Discontinuous Media in Rock Engineering
Wave Propagation through Discontinuous Media in Rock Engineering THESIS SUMMARY Andrea Perino Department of Structural and Geotechnical Engineering Politecnico di Torino, Italy 2 Thesis summary Abstract
More informationAn Overview of Geotechnical Earthquake Engineering
An Overview of Geotechnical Earthquake Engineering Sudhir K Jain Slide 1 Outline Introduction to Seismic Design Principle Dynamic Soil Properties Site Effects Soil Structure Interaction Issues for Foundation
More informationSeismic Analysis of Retaining Structures. Nanjundaswamy P. Department of Civil Engineering S J College of Engineering, Mysore
Seismic Analysis of Retaining Structures Nanjundaswamy P. Department of Civil Engineering S J College of Engineering, Mysore pnswamy@yahoo.com Retaining Walls Retaining Walls. Where? Retaining Walls. Road
More informationA METHODOLOGY FOR ASSESSING EARTHQUAKE-INDUCED LANDSLIDE RISK. Agency for the Environmental Protection, ITALY (
A METHODOLOGY FOR ASSESSING EARTHQUAKE-INDUCED LANDSLIDE RISK Roberto W. Romeo 1, Randall W. Jibson 2 & Antonio Pugliese 3 1 University of Urbino, ITALY (e-mail: rwromeo@uniurb.it) 2 U.S. Geological Survey
More informationPROTECTING MONUMENTS AND HISTORICAL SETTINGS FROM THE NEXT EARTHQUAKE
PROTECTING MONUMENTS AND HISTORICAL SETTINGS FROM THE NEXT EARTHQUAKE R.PAPADHMHTRIOU, L.PELLI EUROPEAN CENTER OF PREVENTING & FORECASTING OF EARTHQUAKES Confronting the problem SEISMIC RISK R SEISMIC
More informationWhere αh and αv are the horizontal and vertical pseudostatic k h, k v = coefficient of horizontal and vertical pseudostatic
International Journal of Scientific & Engineering Research, Volume 5, Issue 12, December-2014 75 Seismic Active Earth Pressure behind Retaining Wall Roshni John 1, K. Preethakumari 2, Pankaj Sethi 3 1
More informationEMPIRICAL EVIDENCE FROM THE NORTHRIDGE EARTHQUAKE FOR SITE- SPECIFIC AMPLIFICATION FACTORS USED IN US BUILDING CODES
EMPIRICAL EVIDENCE FROM THE NORTHRIDGE EARTHQUAKE FOR SITE- SPECIFIC AMPLIFICATION FACTORS USED IN US BUILDING CODES Roger D BORCHERDT And Thomas E FUMAL SUMMARY Site-specific amplification factors, F
More informationBehavior of Concrete Dam under Seismic Load
Behavior of Concrete Dam under Seismic Load A. Hebbouche Ecole Nationale Supérieure d Hydraulique, Blida-Algérie Departement of civil engineering, Saad dahlab university, Blida, Algeria M. Bensaibi Departement
More informationModelling Progressive Failure with MPM
Modelling Progressive Failure with MPM A. Yerro, E. Alonso & N. Pinyol Department of Geotechnical Engineering and Geosciences, UPC, Barcelona, Spain ABSTRACT: In this work, the progressive failure phenomenon
More informationDeformation And Stability Analysis Of A Cut Slope
Deformation And Stability Analysis Of A Cut Slope Masyitah Binti Md Nujid 1 1 Faculty of Civil Engineering, University of Technology MARA (Perlis), 02600 Arau PERLIS e-mail:masyitahmn@perlis.uitm.edu.my
More information2C09 Design for seismic and climate changes
2C09 Design for seismic and climate changes Lecture 10: Characterisation of seismic motion Aurel Stratan, Politehnica University of Timisoara 07/04/2017 European Erasmus Mundus Master Course Sustainable
More informationSeismic Analysis of Concrete Dams Workshop Field investigations and foundation material properties USSD Annual Conference April 6-7, 2017
Seismic Analysis of Concrete Dams Workshop Field investigations and foundation material properties Foundation deformation modulus For analysis purposes, approach used at Reclamation is to reduce the laboratory
More informationEARTHQUAKE SAFETY OF AN ARCH-GRAVITY DAM WITH A HORIZONTAL CRACK IN THE UPPER PORTION OF THE DAM
EARTHQUAKE SAFETY OF AN ARCH-GRAVITY DAM WITH A HORIZONTAL CRACK IN THE UPPER PORTION OF THE DAM Martin WIELAND 1 And Sujan MALLA 2 SUMMARY A horizontal crack first appeared along the downstream wall of
More informationLecture # 6. Geological Structures
1 Lecture # 6 Geological Structures ( Folds, Faults and Joints) Instructor: Dr. Attaullah Shah Department of Civil Engineering Swedish College of Engineering and Technology-Wah Cantt. 2 The wavy undulations
More informationAnalysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis
Analysis of Blocky Rock Slopes with Finite Element Shear Strength Reduction Analysis R.E. Hammah, T. Yacoub, B. Corkum & F. Wibowo Rocscience Inc., Toronto, Canada J.H. Curran Department of Civil Engineering
More informationDUCTILITY BEHAVIOR OF A STEEL PLATE SHEAR WALL BY EXPLICIT DYNAMIC ANALYZING
The 4 th World Conference on arthquake ngineering October -7, 008, Beijing, China ABSTRACT : DCTILITY BHAVIOR OF A STL PLAT SHAR WALL BY XPLICIT DYNAMIC ANALYZING P. Memarzadeh Faculty of Civil ngineering,
More informationStatic & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering
Static & Dynamic Analysis of Structures A Physical Approach With Emphasis on Earthquake Engineering Edward LWilson Professor Emeritus of Civil Engineering University of California, Berkeley Fourth Edition
More informationEvaluation of Fault Foundation Interaction, Using Numerical Studies
Evaluation of Fault Foundation Interaction, Using Numerical Studies Jabbary, M. Msc Student, School of Civil Engineering, Iran University of Science and Technology, Tehran, Iran, Nabizadeh, A. PhD Candidate,
More informationRole of hysteretic damping in the earthquake response of ground
Earthquake Resistant Engineering Structures VIII 123 Role of hysteretic damping in the earthquake response of ground N. Yoshida Tohoku Gakuin University, Japan Abstract Parametric studies are carried out
More informationEvaluating the effects of near-field earthquakes on the behavior of moment resisting frames
Comp. Meth. Civil Eng., Vol. 3, 2 (2012) 79-91 Copyright by the University of Guilan, Printed in I.R. Iran CMCE Computational Methods in Civil Engineering Evaluating the effects of near-field earthquakes
More informationA study on seismic vulnerability of Inca constructions
A study on seismic vulnerability of Inca constructions J. Ogawa 1, C. Cuadra 1, M. B. Karkee 1 & J. Rojas 2 1 Graduate School of System Science & Technology, Akita Prefectural University, Tsuchiya, Honjo,
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationSoil-Structure Interaction in Nonlinear Pushover Analysis of Frame RC Structures: Nonhomogeneous Soil Condition
ABSTRACT: Soil-Structure Interaction in Nonlinear Pushover Analysis of Frame RC Structures: Nonhomogeneous Soil Condition G. Dok ve O. Kırtel Res. Assist., Department of Civil Engineering, Sakarya University,
More informationThe effect of discontinuities on stability of rock blocks in tunnel
International Journal of the Physical Sciences Vol. 6(31), pp. 7132-7138, 30 November, 2011 Available online at http://www.academicjournals.org/ijps DOI: 10.5897/IJPS11.777 ISSN 1992-1950 2011 Academic
More informationThe Stability of a Laminated Voussoir Beam: Back Analysis of a Historic Roof Collapse Using DDA
PII: S0148-9062(97)00309-4 Int. J. Rock Mech. Min. Sci. Vol. 35, No. 2, pp. 165±181, 1998 # 1998 Elsevier Science Ltd. All rights reserved Printed in Great Britain 0148-9062/98 $19.00 + 0.00 The Stability
More informationCHAPTER 8 SUMMARY, CONCLUSIONS AND RECOMMENDATIONS
CHAPTER 8 SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 8.1 SUMMARY This thesis aimed to investigate the mechanisms behind valley closure and upsidence over unmined coal and old longwall panels using UDEC.
More informationEffect of structural design on fundamental frequency of reinforced-soil retaining walls
Soil Dynamics and Earthquake Engineering 19 (2000) 137 157 www.elsevier.com/locate/soildyn Effect of structural design on fundamental frequency of reinforced-soil retaining walls K. Hatami*, R.J. Bathurst
More informationRock slope failure along non persistent joints insights from fracture mechanics approach
Rock slope failure along non persistent joints insights from fracture mechanics approach Louis N.Y. Wong PhD(MIT), BSc(HKU) Assistant Professor and Assistant Chair (Academic) Nanyang Technological University,
More informationEFFECT OF SEISMIC WAVE INCLINATION ON STRUCTURAL RESPONSE
EFFECT OF SEISMIC WAVE INCLINATION ON STRUCTURAL RESPONSE t I M. Tabatabaie,l A.M. ASCEj N. Abrahamson,2 and J. P. Singh, M. ASCE ABSTRACT A new procedure is developed using the computer program SASSI
More informationComparative study between the push-over analysis and the method proposed by the RPA for the evaluation of seismic reduction coefficient
33, Issue (27) 5-23 Journal of Advanced Research in Materials Science Journal homepage: www.akademiabaru.com/arms.html ISSN: 2289-7992 Comparative study between the push-over analysis and the method proposed
More informationDynamics of structures
Dynamics of structures 2.Vibrations: single degree of freedom system Arnaud Deraemaeker (aderaema@ulb.ac.be) 1 Outline of the chapter *One degree of freedom systems in real life Hypothesis Examples *Response
More informationRecent Research on EPS Geofoam Seismic Buffers. Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada
Recent Research on EPS Geofoam Seismic Buffers Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada What is a wall (SEISMIC) buffer? A compressible inclusion placed between
More informationSeismic Analysis of Soil-pile Interaction under Various Soil Conditions
Seismic Analysis of Soil-pile Interaction under Various Soil Conditions Preeti Codoori Assistant Professor, Department of Civil Engineering, Gokaraju Rangaraju Institute of Engineering and Technology,
More information(Refer Slide Time: 1: 19)
Mechanical Measurements and Metrology Prof. S. P. Venkateshan Department of Mechanical Engineering Indian Institute of Technology, Madras Module - 4 Lecture - 46 Force Measurement So this will be lecture
More informationALASKA ENERGY AUTHORITY AEA ENGINEERING FEASIBILITY REPORT. Appendix B8. Finite Element Analysis
ALASKA ENERGY AUTHORITY AEA11-022 ENGINEERING FEASIBILITY REPORT Appendix B8 Finite Element Analysis Susitna-Watana Hydroelectric Project Alaska Energy Authority FERC Project No. 14241 December 2014 Seismic
More information