Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams

Size: px
Start display at page:

Download "Application of pseudo-symmetric technique in dynamic analysis of concrete gravity dams"

Transcription

1 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 presented for seismic analysis of concrete gravity dams in the time domain. The technique is based on symmetric formulation of dam-reservoir systems utilizing pressure degrees of freedom in the fluid domain. The method is explained initially in relation to other available techniques. Following from this, the analysis of a typical concrete dam is considered as a verification example by utilizing a modified program based on this technique. The methodology introduced is very convenient and can be easily implemented in general-purpose finite element programs with regard to their fluid-structure interaction modules. 1 Introduction There has been extensive research in relation to Fluid-Structure Interaction problem encountered in dynamic analysis of concrete dams. Some of these studies are formulated very efficiently in frequency domain [l, 21. However, these methods are limited to linear dynamic analysis. Others, which have studied techniques suitable for time domain, are actually encountered with the solution of an unsymmetric coupled equations. In general, algorithms relying on unsymmetric solvers can be applied as the first alternative. However, the efficiency would be very low. There is also a technique available which transforms the coupled Fluid-Structure equations to symmetric relations by formulating the fluid elements in terms of velocity potential degrees of freedom [31. In this paper, a symmetric approach based on nodal pressures degrees of freedom [4], is used as the basis of the study. The method is modified and a technique is proposed that is referred to as Pseudo-Symmetric Approach. The procedure is described initially in relation to other available techniques. Meanwhile, the methodology is implemented in a special computer program "MAP-76" [5]. Subsequently, the analysis of Pine Flat Dam is considered as a verification

2 208 Fluid Structure Interaction 11 example. The results are discussed and they are also compared with a previous study from accuracy point of view. 2 Formulation Let us consider a general two-dimensional dam-reservoir system interacting internally due to ground excitations. The dam and water domain, are assumed to be discretized by plane solid and fluid finite elements, respectively. It is easily shown that the coupled equations of dam-reservoir system can be written in combined form as [4]: M, C, K in this relation represent mass, damping and stiffness matrices of the dam body, respectively. r is the vector of dam nodal relative displacements. a, denotes the vector of ground accelerations, and J is a matrix with each two rows equal to a 2 X 2 identity matrix (its columns correspond to unit rigid body motion in horizontal and vertical directions) [l, 41. Furthermore, G, L, H in this relation, are assembled matrices of fluid domain, and the matrix B is usually referred to as the interaction matrix 14, 61. Meanwhile, p is the vector of nodal pressures, and p, c denote mass density and pressure wave velocity of fluid. It is clear that the overall mass and stiffness matrices of the system are unsymmetric in the above relation. Therefore, this coupled equation is usually solved by using the scattered approach [4, 71. That is separating dam and reservoir equations and solving the problem by an iterative scheme. Although there is no need for unsymmetric solvers in this approach, the major drawback is that it requires iterations even for linear problems. Meanwhile, the algorithm is not appealing to be coded in a general-purpose finite element program, which considers the fluid-structure interaction problem as one of its special cases. Another remedy for avoiding unsymmetric solvers is to use velocity potentials instead of pressures as degrees of freedom in the water domain. Because in that case, one would obtain the following coupled equations for the system [3, 61 that replaces equation (1): represents the vector of nodal velocity potentials. Meanwhile, v, is the vector of ground velocities. It is easy to see that this relation can become symmetrical by simply multiplying the lower partitions by a factor of -p. In this case, the mass, damping and stiffness matrices of the system are all symmetric. This is the formulation applied in some of the general-purpose finite

3 Fluid Structure Interaction I1 209 element programs [3]. The disadvantage of this method is that vector of ground velocities are required to be calculated which could also introduce some error into the solution. Furthermore, the calculation of nodal velocity potential quantities, are less attractive for most engineers in comparison to nodal pressures. 2.1 Symmetric approach based on nodal pressures It is well known that one can obtain symmetric coupled equations even in the case where nodal pressures are used as degrees of freedom in the fluid domain [4]. However, the method is presented in a slightly different technique in this section, which can be utilized as the base of the Pseudo-Symmetric approach. The method commences by modifying the relation (l), such that the lower partitions of that matrix equation, is multiplied by a factor of - (pao)-l. The parameter a, utilized in this factor, is one of the coefficients obtained in Newmark's algorithm [4], defined in terms of Newmark's parameter a (a =0.25 in constant average acceleration method) and time increment At : Therefore, the coupled equations of dam-reservoir system would be: Or alternatively in a more compact form: - - It is obvious that matrices M, C and K can be written as sum of the symmetric and unsymmetric parts as below: It is noted from equation (4) that the damping matrix (-6) is totally symmetric, and the following relation also holds:

4 21 0 Fluid Structure Interaction 11 Applying Newrnark7s algorithm on equation (5), the effective stiffness matrix (K) and the force vector at instant n+l (R,,, ) would be: + (Ms +Mu )(ao % +a, Gn +a, ;n ) + CS (a, F, +a, pn +a, ) Substituting (7) into (g), it yields, Although total mass and stiffness matrices of the system are unsymmetric, it is clear from relation (10) that the effective stiffness matrix is symmetric. This means that the need for unsymmetric solver is eliminated. 2.2 Pseudo-symmetric approach Considering the coupled equations (1) of dam-reservoir system, it is noticed that unsymmetric terms in the total mass and stiffness matrices are related to B matrix or its transpose. This matrix is usually obtained by assemblage of contributing submatrices of interface elements located at fluid-solid contact, or even surfaces where fluid elements are adjacent to rigid or absorbing boundaries. However, to make it more efficient from programming point of view, one can eliminate these interface elements and consider its effect as part of adjacent fluid element matrices. In that case, matrices of the i& fluid element which contribute to the corresponding total mass, damping and stiffness matrices of the system would be generally as follows, respectively: In the proposed method, interface elements are excluded and their effects are considered as part of fluid element matrices as explained above. Furthermore, everything is made symmetric from the very beginning. This means that fluid element matrices are made symmetric by transferring the unsymmetric part of the stiffness matrix to the mass matrix while multiplying it by a factor of ail. Meanwhile, rows corresponding to the pressure degrees of freedom are multiplied by a factor of - (pao)-l. Therefore, when the total mass, damping and stiffness matrices of the system are formed, one obtains:

5 Fluid Structure Interaction I The total matrices of the system are different from corresponding matrices obtained in the previous method (relation (6)), and in this case all matrices are symmetric. However, the effective stiffness matrix, which results from applying Newmark's method becomes, Which is exactly the same as effective stiffness matrix resulted for the previous method (relation (10)). Therefore, the need for unsymmetric solver is once more eliminated. The formulation of the effective force vector in this method is written as follows: - = + (Ms +MU +M; );(a, 5 +a, Pn +a, ;n ) + CS (a, i;, +a, kn +a, In ) In this relation, a special operator * is utilized, which emphasizes that in this matrix operation, the stored terms corresponding to non-zero elements of M; must not participate. Therefore, the resulting effective force vector would become exactly as in relation (9), obtained for previous method. The location of these storages would be very clear in the total mass matrix, if all of the pressure degrees of freedom are stored at the end. The unhown vector calculated at each step of the Pseudo-Symmetric algorithm would definitely be the same as in previous method, since the effective stiffness and force vector are exactly the same. Meanwhile, all the matrices of fluid elements are symmetric like any solid elements, even though this was forced artificially at the element level. This approach is very convenient for general-purpose programs, as far as symmetric upper skyline storage scheme that can be utilized for total mass, damping and stiffness matrices. Furthermore, it takes advantage of symmetric solvers to obtain the unknowns at each time step. 3 Application on concrete gravity dams In this section, the analysis of Pine Flat Dam is considered as a verification example. This particular dam was selected, because many researchers have analyzed it based on different approaches [l, 81. The dam is m high, with the crest length of m and it is located on the King's River near Fresno, California. A special computer program "MAP-76" is used as the basis of the present study [5]. The program was modified by introducing two-dimensional fluid elements with pressure degrees of freedom into the program, as well as implementing the other required changes in the code based on the Pseudo-Symmetric methodology. All fluid element matrices were defined symmetrical artificially, as explained in the previous section.

6 21 2 Fluid Structure Interaction Finite element idealization Selected models Most researchers have considered the analysis of Pine Flat Dam based on a two dimensional model, which is well accepted for a typical gravity dam. That means, a section of Dam-Reservoir system is considered with unit thickness. The two dimensional elements of dam body is usually assumed in a state of plane stress, while fluid elements formulation is presuming plane strain conditions for fluid domain. In the present study, a two-dimensional model of Pine Flat Dam is considered (Figure 1). The model is based on normal practice idealization as mentioned above. The 8-node plane solid and fluid finite elements are used for the dam and reservoir domains, respectively. The model consists of a total of 439 nodes and 568 degrees of freedom and the mesh includes 40, 90 plane solid and fluid elements, respectively. Figure 1. Finite Element mesh of Dam-Reservoir system Basic parameters The concrete is assumed to be homogeneous and isotropic with the following basic properties: Elastic modulus E, = GPa Poisson's ratio v, = 0.20 Unit weight y, = 24.8 w/m3 The water is taken as compressible, inviscid fluid, with weight density of 10.0 kn/m3 and pressure wave velocity of mls. The water level is considered at the height of m above the base, similar to the reference [l]. The length of the reservoir discretized is 200 m, and Sornmerfeld boundary condition is applied at the upstream boundary, while the reservoir bottom condition is assumed completely reflective. Moreover, in the analysis carried out, the Rayleigh damping matrix is applied and the corresponding coefficients are determined such that equivalent damping for frequencies close to the first and third modes of vibration would be 5% of the critical damping Loading It should be mentioned that static loads (weight, hydrostatic pressures) are each visualized as being applied in one separate increment of time. Therefore, the

7 Fluid Structure Interaction I same time step of 0.01 second, which is chosen in dynamic analysis, is also considered as time increment of static loads application. It is noted that time for static analysis is just a convenient tool for applying the load sequentially, but it is obvious that inertia and damping effects are disregarded in the process. In this respect, the dead load is applied in one increment and hydrostatic pressures thereafter in another increment at negative range of time. At time zero, the actual dynamic analysis begins with the static displacements and stresses being applied as initial conditions. The dynamic excitation considered, is the S69E component of Taft earthquake records, which is applied in the horizontal direction. The time duration utilized, is 13 seconds. 3.2 Analysis results As mentioned in the previous section, a two-dimensional model of Pine Flat Dam is considered. The analysis is carried out based on the Pseudo-Symmetric methodology by utilizing the MAP-76 program. The model is analyzed, and results corresponding to the envelope of maximum tensile and compressive stresses are illustrated in Figures 2 and 3. Furthermore, time histories of some important quantities are presented in Figures 4 and 5. As for verification, some of the available results taken from the work of Chopra et a1 [l] are illustrated in Figure 6 for comparison purposes. However, it should be noticed that even though similar properties are used in both studies, there are still slight differences. In this reference, the dam is discretized by incompatible quadratic finite elements. Meanwhile, closed form solution is utilized for the effects of impounded water, which is derived for a region extending to infinity in the upstream direction, and also assuming a vertical boundary at the damreservoir interface. Furthermore, the analysis is carried out in frequency domain and hysteretic type of damping is employed. While in the present study, the solution is obtained in time domain and Rayleigh type of damping is utilized. Having the above comments in perspective, one can compare the envelope of maximum tensile stresses between the present study (Figure 2) and the reference (Figure 6a). It is observed that very good agreement exist for the magnitude and distribution of maximum tensile stresses in different regions. However, there seems to be some difference in the maximum value of tensile stress, which occurs at the dam heel. This discrepancy can be explained by the fact that, stress concentration is present in that location, and the reference is not providing additional contour lines close to this corner. The interpolated value in the present study results, which correspond to the location of maximum contour shown in the reference, would be roughly o, =2.35 Mpa. This shows a reasonable agreement (less than 5% difference) as far as the envelope of maximum tensile stresses is concerned. Additional result of the reference that can be compared, is the time history of horizontal displacement at dam crest. This is shown in Figure 6b and can be compared with the result of present work depicted in Figure 4a. It is observed that the history of displacements follow very similar trends and the maximum values are different by less than 2%. Therefore, it can be concluded that results of the present study are in good agreement with the reference [l].

8 21 4 Fluid Structure Interaction 11 Figure 2. Envelope of maximum tensile principal stresses (MPa) Figure 3. Envelope of maximum compressive principal stresses (MPa) Figure 4. Displacement histories at dam crest

9 DISPL - INCHES r&lo-w

10 21 6 Fluid Structure Interaction 11 4 Conclusions A technique is proposed for the seismic analysis of concrete gravity dams. The method is described, and the procedure is implemented in a special computer program "MAP-76". Meanwhile, analysis of Pine Flat Dam was considered as a verifying example. A two-dimensional model of dam-reservoir system is considered and analyzed. The results are also compared with the reference. Overall, the main conclusions obtained can be listed as follows: The results of present study are in close agreement with the reference. The small differences, which exist with respect to the reference, are in acceptable range, having in mind that there were some differences in the idealizations utilized as mentioned in the discussion of the results. The Pseudo-Symmetric approach is proved to be an effective technique for dynamic analysis of concrete gravity dams. The methodology introduced is very convenient and can be easily implemented in general-purpose finite element programs in regard to their fluid-structure interaction modules. Although the procedure was explained for linear dynamic analysis, it is easily expandable to nonlinear dynamic analysis of concrete dams, since it is carried out in time domain. References 1. Chopra, A. K., Chakrabarti, P., and Gupta, S., "Earthquake response of concrete gravity dams including hydrodynamic and foundation interaction effects", Report No. EERC-80101, University of California, Berkeley, Jan Fok, K.-L. and Chopra, A. K., "Earthquake analysis and response of concrete arch dams", Report No. EERC-85/07, University of California, Berkeley, July Olson, L. G. and Bathe, K. J., "Analysis of fluid-structure interactions. A direct symmetric coupled formulation based on the fluid velocity potential", Intern. Journal of Computers & Structures, Vol. 21, Nos. 112, Special Issue: Nonlinear Finite Element Analysis and ADINA, Zienkiewicz, O.C. and Taylor, R.L., "The Finite Element Method", 5th Edition, Butterworth-Heinemann, Oxford, UK, Lotfi, V., "MAP-76: A program for analysis of concrete dams", Amirkabir University of Technology, Tehran, Iran, Lotfi, V., "Dynamic analysis of concrete dams", Mahab Ghodss Consulting Engineers, Tehran, Iran, Cervera, M., Oliver, J. and Faria, R. "Seismic evaluation of concrete dams via continuum damage models", Journal of Earthquake Engineering & Structural Dynamics, Vol. 24, No. 9, Sept Feltrin, G., Wepf, D. and Bachmann, H., "Seismic cracking of concrete gravity dams", Journal of Dam Engineering, Vol. I, Issue 4, Oct

EFFECTS OF RESERVOIR LENGTH ON DYNAMIC ANALYSIS OF CONCRETE GRAVITY DAMS

EFFECTS OF RESERVOIR LENGTH ON DYNAMIC ANALYSIS OF CONCRETE GRAVITY DAMS The th October -,, Beijing, China EFFECTS OF RESERVOIR LENGTH ON DYNAMIC ANALYSIS OF CONCRETE GRAVITY DAMS A. Fathi and V. Lotfi M.Sc. Student, Dept. of Civil and Environmental Engineering, Amirkabir University

More information

COMPARISON 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 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 information

Modelling of Concrete Gravity Dam Including Dam-Water-Foundation Rock Interaction

Modelling of Concrete Gravity Dam Including Dam-Water-Foundation Rock Interaction World Applied Sciences Journal 22 (4): 538-546, 2013 ISSN 1818-4952 IDOSI Publications, 2013 DOI: 10.5829/idosi.wasj.2013.22.04.551 Modelling of Concrete Gravity Dam Including Dam-Water-Foundation Rock

More information

NONLINEAR ANALYSIS OF A DAM-RESERVOIR-FOUNDATION SYSTEM UNDER SPATIALLY VARIABLE SEISMIC EXCITATIONS

NONLINEAR 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 information

Behavior of Concrete Dam under Seismic Load

Behavior 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 information

Effects of perimetrical and vertical joint opening on the seismic response of a concrete arch dam

Effects of perimetrical and vertical joint opening on the seismic response of a concrete arch dam Effects of perimetrical and vertical joint opening on the seismic response of a concrete arch dam Soheil RAZAVI DARBAR', Hiroyuki WATANABE' *Ph.D. Candidate, Earthquake Engineering Laboratory, Department

More information

STUDY OF DYNAMIC SOIL-STRUCTURE INTERACTION OF CONCRETE GRAVITY DAMS

STUDY OF DYNAMIC SOIL-STRUCTURE INTERACTION OF CONCRETE GRAVITY DAMS STUDY OF DYNAMIC SOIL-STRUCTURE INTERACTION OF CONCRETE GRAVITY DAMS Djamel OUZANDJA 1, Fatiha BENKECHIDA 2, Toufiq OUZANDJA 3, Hamza BELHADED 4 ABSTRACT The safety evaluation of the dams subjected to

More information

Analytical Solution for a Fluid-Structure Interaction Problem in Comparison with Finite Element Solution

Analytical Solution for a Fluid-Structure Interaction Problem in Comparison with Finite Element Solution Analytical Solution for a Fluid-Structure Interaction Problem in Comparison with Finite Element Solution Amirhossein Keivani & Ahmad Shooshtari Ferdowsi University of Mashhad, Mashhad, Iran. Ahmad Aftabi

More information

Evaluation of the behaviour of an arch-gravity dam. featuring a pre-existing crack

Evaluation of the behaviour of an arch-gravity dam. featuring a pre-existing crack Evaluation of the behaviour of an arch-gravity dam featuring a pre-existing crack Dr Aïssa Mellal, Civil Engineer STUCKY SA, Switzerland NUMERICS IN GEOTECHNICS AND STRUCTURES - ZSoil Days - 1-2 September

More information

ALASKA ENERGY AUTHORITY AEA ENGINEERING FEASIBILITY REPORT. Appendix B8. Finite Element Analysis

ALASKA 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

Static & Dynamic. Analysis of Structures. Edward L.Wilson. University of California, Berkeley. Fourth Edition. Professor Emeritus of Civil Engineering

Static & 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 information

International Journal Of Scientific & Engineering Research Volume 4, Issue 5, May-2013 ISSN

International Journal Of Scientific & Engineering Research Volume 4, Issue 5, May-2013 ISSN 180 SEISMIC FRACTURE ANALYSIS IN CONCRETE GRAVITY DAMS Deeja Alora, Indrani Gogoi Abstract Modelling of crack propagation in solids has been a major area of focus both in industry and research communities.

More information

Dynamic Analysis Contents - 1

Dynamic 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 information

Seismic Analyses of Concrete Gravity Dam with 3D Full Dam Model

Seismic Analyses of Concrete Gravity Dam with 3D Full Dam Model Seismic Analyses of Concrete Gravity Dam with 3D Full Dam Model Haibo Wang, Deyu Li & Huichen Yang China Institute of Water Resources and Hydropower Research, Beijing, China SUMMARY: Seismic analyses of

More information

A.0 IDUKKI ARCH DAM - ANALYSIS FOR VARIOUS LOAD CASES AND DISCRETISATIONS

A.0 IDUKKI ARCH DAM - ANALYSIS FOR VARIOUS LOAD CASES AND DISCRETISATIONS Annexure A.0 IDUKKI ARCH DAM - ANALYSIS FOR VARIOUS LOAD CASES AND DISCRETISATIONS In this Annexure, Idukki arch dam chosen for case study is analyzed with the developed program for various load conditions

More information

A simple lumped mass-damper idealization for dam-reservoir- foundation system for seismic analysis

A simple lumped mass-damper idealization for dam-reservoir- foundation system for seismic analysis ICCM2015, 14-17 th July, Auckland, NZ A simple lumped mass-damper idealization for dam-reservoir- foundation system for seismic analysis *D.K. Paul¹, A. Banerjee 2, 3, R.N. Dubey 1, M.K.Alam-Chowdhury

More information

The Dynamic Response Analysis of Concrete Gravity Dam under the Earthquake

The 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 information

THEME A. Analysis of the elastic behaviour of La Aceña arch-gravity dam

THEME A. Analysis of the elastic behaviour of La Aceña arch-gravity dam THEME A Analysis of the elastic behaviour of La Aceña arch-gravity dam Gjorgi KOKALANOV, Professor, Faculty of Civil Eng., Skopje, Republic of Macedonia Ljubomir TANČEV, Professor, Faculty of Civil Eng.,

More information

Seismic stability safety evaluation of gravity dam with shear strength reduction method

Seismic stability safety evaluation of gravity dam with shear strength reduction method Water Science and Engineering, 2009, 2(2): 52-60 doi:10.3882/j.issn.1674-2370.2009.02.006 http://kkb.hhu.edu.cn e-mail: wse@hhu.edu.cn Seismic stability safety evaluation of gravity dam with shear strength

More information

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma).

Structural Dynamics. Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). Structural Dynamics Spring mass system. The spring force is given by and F(t) is the driving force. Start by applying Newton s second law (F=ma). We will now look at free vibrations. Considering the free

More information

VULNERABILITY OF AGED CONCRETE GRAVITY DAMS

VULNERABILITY OF AGED CONCRETE GRAVITY DAMS 13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 004 Paper No. 1839 VULNERABILIY OF AGED CONCREE GRAVIY DAMS INDRANI GOGOI 1 and DAMODAR MAIY SUMMARY he ageing of a

More information

EARTHQUAKE 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 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 information

Chapter 2 Finite Element Formulations

Chapter 2 Finite Element Formulations Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are

More information

Earthquake Analysis of Arch Dams

Earthquake Analysis of Arch Dams Earthquake Analysis of Arch Dams Anil K. Chopra Electricity Generating Authority of Thailand Bangkok, Thailand December 7-8, 2010 Earthquake Analysis of Arch Dams: Factors To Be Considered Focus on linear

More information

CURRENT METHODOLOGY AT THE BUREAU OF RECLAMATION FOR THE NONLINEAR ANALYSES OF ARCH DAMS USING EXPLICIT FINITE ELEMENT TECHNIQUES

CURRENT METHODOLOGY AT THE BUREAU OF RECLAMATION FOR THE NONLINEAR ANALYSES OF ARCH DAMS USING EXPLICIT FINITE ELEMENT TECHNIQUES CURRENT METHODOLOGY AT THE BUREAU OF RECLAMATION FOR THE NONLINEAR ANALYSES OF ARCH DAMS USING EXPLICIT FINITE ELEMENT TECHNIQUES Barbara Mills-Bria, P.E., 1 Larry Nuss, P.E. 2 and Dr. Anil Chopra 2 1

More information

Investigation of crack development in concrete dams due to water level fluctuation during earthquakes

Investigation of crack development in concrete dams due to water level fluctuation during earthquakes River Basin Management VIII 271 Investigation of crack development in concrete dams due to water level fluctuation during earthquakes S. A. Neshaei 1, M. K. Abadi 2 & R. Vahedzadegan 3 1 Faculty of Engineering,

More information

Analysis of the Body of the Double Concrete Arch Dam Bakhtiari by Three Dimensional Finite Element Encoding

Analysis of the Body of the Double Concrete Arch Dam Bakhtiari by Three Dimensional Finite Element Encoding World Applied Sciences Journal 19 (8): 1194-101, 01 ISSN 1818-495 IDOSI Publications, 01 DOI: 10.589/idosi.wasj.01.19.08.313 Analysis of the Body of the Double Concrete Arch Dam Bakhtiari by Three Dimensional

More information

Investigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis

Investigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis Print ISSN: 2322-2093; Online ISSN: 2423-6691 DOI: 10.7508/ceij.2016.02.011 Technical Note Investigation of Utilizing a Secant Stiffness Matrix for 2D Nonlinear Shape Optimization and Sensitivity Analysis

More information

* ABAQUS DDA DEM Yu Houhe et al. BT Email: vahlotfi@aut.ac.ir * ABAQUS c. n tan 2C cos ft 2Rt 1 sin f t Max C.cot 2C cos f c 2Rc 1 sin f c f t Parameter E C Definition Modulus of Elasticity (GPa) Poisson's

More information

Plastic Damage Analysis of Koyna Dam in Different Damping Mechanisms with Dam Water Interaction

Plastic Damage Analysis of Koyna Dam in Different Damping Mechanisms with Dam Water Interaction Plastic Damage Analysis of Koyna Dam in Different Damping Mechanisms with Dam Water Interaction O. Omidi & V. Lotfi Amirkabir University of Technology, Iran S. Valliappan University of New South Wales,

More information

The Finite Element Method for Solid and Structural Mechanics

The Finite Element Method for Solid and Structural Mechanics The Finite Element Method for Solid and Structural Mechanics Sixth edition O.C. Zienkiewicz, CBE, FRS UNESCO Professor of Numerical Methods in Engineering International Centre for Numerical Methods in

More information

Numerical 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 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 information

Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams

Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Practical methodology for inclusion of uplift and pore pressures in analysis of concrete dams Michael McKay 1 and Francisco Lopez 2 1 Dams Engineer, GHD Pty 2 Principal Dams/Structural Engineer, GHD Pty

More information

Codal Provisions IS 1893 (Part 1) 2002

Codal Provisions IS 1893 (Part 1) 2002 Abstract Codal Provisions IS 1893 (Part 1) 00 Paresh V. Patel Assistant Professor, Civil Engineering Department, Nirma Institute of Technology, Ahmedabad 38481 In this article codal provisions of IS 1893

More information

Analysis of Planar Truss

Analysis of Planar Truss Analysis of Planar Truss Although the APES computer program is not a specific matrix structural code, it can none the less be used to analyze simple structures. In this example, the following statically

More information

Dynamic Response of EPS Blocks /soil Sandwiched Wall/embankment

Dynamic 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 information

DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION

DYNAMIC 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 information

DISPENSA FEM in MSC. Nastran

DISPENSA FEM in MSC. Nastran DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation

More information

SEISMIC RESPONSE EVALUATION OF AN RC BEARING WALL BY DISPLACEMENT-BASED APPROACH

SEISMIC RESPONSE EVALUATION OF AN RC BEARING WALL BY DISPLACEMENT-BASED APPROACH 3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -, 4 Paper No. 49 SEISMIC RESPONSE EVALUATION OF AN RC BEARING WALL BY DISPLACEMENT-BASED APPROACH Chang-Hun HYUN, Sanghyun

More information

Gravity dam and earthquake

Gravity 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 information

Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground

Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground IOP Conference Series: Earth and Environmental Science PAPER OPEN ACCESS Computational Simulation of Dynamic Response of Vehicle Tatra T815 and the Ground To cite this article: Jozef Vlek and Veronika

More information

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS

VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS The 4 th World Conference on Earthquake Engineering October -7, 008, Beijing, China VORONOI APPLIED ELEMENT METHOD FOR STRUCTURAL ANALYSIS: THEORY AND APPLICATION FOR LINEAR AND NON-LINEAR MATERIALS K.

More information

Selection of Rayleigh Damping Coefficients for Seismic Response Analysis of Soil Layers

Selection of Rayleigh Damping Coefficients for Seismic Response Analysis of Soil Layers Selection of Rayleigh Damping Coefficients for Seismic Response Analysis of Soil Layers Huai-Feng Wang, Meng-Lin Lou, Ru-Lin Zhang Abstract One good analysis method in seismic response analysis is direct

More information

RESPONSE-SPECTRUM-BASED ESTIMATES OF MOHR S CIRCLE

RESPONSE-SPECTRUM-BASED ESTIMATES OF MOHR S CIRCLE 3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -, 004 Paper No. 5 RESPONSE-SPECTRUM-BASED ESTIMATES OF MOHR S CIRCLE Charles MENUN SUMMARY The utility of Mohr s circle as

More information

2D Liquefaction Analysis for Bridge Abutment

2D Liquefaction Analysis for Bridge Abutment D Liquefaction Analysis for Bridge Abutment Tutorial by Angel Francisco Martinez Integrated Solver Optimized for the next generation 64-bit platform Finite Element Solutions for Geotechnical Engineering

More information

Behaviour of Earth Dam under Seismic Load Considering Nonlinearity of the Soil

Behaviour of Earth Dam under Seismic Load Considering Nonlinearity of the Soil Open Journal of Civil Engineering, 216, 6, 75-83 Published Online March 216 in SciRes. http://www.scirp.org/journal/ojce http://dx.doi.org/1.4236/ojce.216.627 Behaviour of Earth Dam under Seismic Load

More information

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS

NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS NUMERICAL SIMULATION OF THE INELASTIC SEISMIC RESPONSE OF RC STRUCTURES WITH ENERGY DISSIPATORS ABSTRACT : P Mata1, AH Barbat1, S Oller1, R Boroschek2 1 Technical University of Catalonia, Civil Engineering

More information

midas Civil Dynamic Analysis

midas Civil Dynamic Analysis Edgar De Los Santos Midas IT August 23 rd 2017 Contents: Introduction Eigen Value Analysis Response Spectrum Analysis Pushover Analysis Time History Analysis Seismic Analysis Seismic Analysis The seismic

More information

Principles of Finite Element for Design Engineers and Analysts. Ayman Shama, Ph.D., P.E., F.ASCE

Principles of Finite Element for Design Engineers and Analysts. Ayman Shama, Ph.D., P.E., F.ASCE Principles of Finite Element for Design Engineers and Analysts Ayman Shama, Ph.D., P.E., F.ASCE Outline Principles of Engineering Analysis The development of the finite element method Types of elements

More information

The Influence of Initial Cracks on the Crack Propagation Process of Concrete Gravity Dam-Reservoir-Foundation Systems

The Influence of Initial Cracks on the Crack Propagation Process of Concrete Gravity Dam-Reservoir-Foundation Systems See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/281382524 The Influence of Initial Cracks on the Crack Propagation Process of Concrete Gravity

More information

2D & 3D Nonlinear Dynamic Analysis of an Asphaltic Concrete Core Rockfill Dam (a Case Study)

2D & 3D Nonlinear Dynamic Analysis of an Asphaltic Concrete Core Rockfill Dam (a Case Study) 1 2D & 3D Nonlinear Dynamic Analysis of an Asphaltic Concrete Core Rockfill Dam (a Case Study) A. Akhtarpour Ph.D., Maharab Consulting Engineers Co., Mashhad, Iran. A_akhtarpur@aut.ac.ir A. Khodaii Ph.D.,

More information

APPLICATION OF RESPONSE SPECTRUM METHOD TO PASSIVELY DAMPED DOME STRUCTURE WITH HIGH DAMPING AND HIGH FREQUENCY MODES

APPLICATION OF RESPONSE SPECTRUM METHOD TO PASSIVELY DAMPED DOME STRUCTURE WITH HIGH DAMPING AND HIGH FREQUENCY MODES 3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 5 APPLICATION OF RESPONSE SPECTRUM METHOD TO PASSIVELY DAMPED DOME STRUCTURE WITH HIGH DAMPING AND HIGH FREQUENCY

More information

Analysis of Concrete Walls under Earthquake Action Influence from different types of Foundation

Analysis of Concrete Walls under Earthquake Action Influence from different types of Foundation Analysis of Concrete Walls under Earthquake Action Influence from different types of Foundation E.R. Thorhallsson & H.S. Olafsson Reykjavik University, Iceland ABSTRACT In these paper concrete walls are

More information

International Journal of Advanced Engineering Technology E-ISSN

International Journal of Advanced Engineering Technology E-ISSN Research Article INTEGRATED FORCE METHOD FOR FIBER REINFORCED COMPOSITE PLATE BENDING PROBLEMS Doiphode G. S., Patodi S. C.* Address for Correspondence Assistant Professor, Applied Mechanics Department,

More information

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d

EDEM DISCRETIZATION (Phase II) Normal Direction Structure Idealization Tangential Direction Pore spring Contact spring SPRING TYPES Inner edge Inner d Institute of Industrial Science, University of Tokyo Bulletin of ERS, No. 48 (5) A TWO-PHASE SIMPLIFIED COLLAPSE ANALYSIS OF RC BUILDINGS PHASE : SPRING NETWORK PHASE Shanthanu RAJASEKHARAN, Muneyoshi

More information

DUCTILITY BEHAVIOR OF A STEEL PLATE SHEAR WALL BY EXPLICIT DYNAMIC ANALYZING

DUCTILITY 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 information

2D Embankment and Slope Analysis (Numerical)

2D Embankment and Slope Analysis (Numerical) 2D Embankment and Slope Analysis (Numerical) Page 1 2D Embankment and Slope Analysis (Numerical) Sunday, August 14, 2011 Reading Assignment Lecture Notes Other Materials FLAC Manual 1. 2. Homework Assignment

More information

Multi-Point Constraints

Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Multi-Point Constraints Single point constraint examples Multi-Point constraint examples linear, homogeneous linear, non-homogeneous linear, homogeneous

More information

Modal Analysis: What it is and is not Gerrit Visser

Modal Analysis: What it is and is not Gerrit Visser Modal Analysis: What it is and is not Gerrit Visser What is a Modal Analysis? What answers do we get out of it? How is it useful? What does it not tell us? In this article, we ll discuss where a modal

More information

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method

Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method 9210-220 Post Graduate Diploma in Mechanical Engineering Computational mechanics using finite element method You should have the following for this examination one answer book scientific calculator No

More information

MODELLING DAMAGE AND PROGRESSIVE COLLAPSE OF FRAMES USING A GAUSSIAN SPRINGS BASED APPLIED ELEMENT METHOD

MODELLING DAMAGE AND PROGRESSIVE COLLAPSE OF FRAMES USING A GAUSSIAN SPRINGS BASED APPLIED ELEMENT METHOD 6th European Conference on Computational Mechanics (ECCM 6) 7th European Conference on Computational Fluid Dynamics (ECFD 7) 1115 June 018, Glasgow, UK MODELLING DAMAGE AND PROGRESSIVE COLLAPSE OF FRAMES

More information

The Seismic Performance of Tousheh Dam During the Chi-Chi Earthquake

The 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 information

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models

Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models Dynamic Analysis of a Reinforced Concrete Structure Using Plasticity and Interface Damage Models I. Rhee, K.J. Willam, B.P. Shing, University of Colorado at Boulder ABSTRACT: This paper examines the global

More information

EXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro

EXTENDED ABSTRACT. Dynamic analysis of elastic solids by the finite element method. Vítor Hugo Amaral Carreiro EXTENDED ABSTRACT Dynamic analysis of elastic solids by the finite element method Vítor Hugo Amaral Carreiro Supervisor: Professor Fernando Manuel Fernandes Simões June 2009 Summary The finite element

More information

Seismic Design of a Hydraulic Fill Dam by Nonlinear Time History Method

Seismic 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 information

Comparison of Structural Models for Seismic Analysis of Multi-Storey Frame Buildings

Comparison of Structural Models for Seismic Analysis of Multi-Storey Frame Buildings Comparison of Structural Models for Seismic Analysis of Multi-Storey Frame Buildings Dj. Ladjinovic, A. Raseta, A. Radujkovic & R. Folic University of Novi Sad, Faculty of Technical Sciences, Novi Sad,

More information

Seismic analysis of horseshoe tunnels under dynamic loads due to earthquakes

Seismic analysis of horseshoe tunnels under dynamic loads due to earthquakes University of Wollongong Research Online Coal Operators' Conference Faculty of Engineering and Information Sciences 2010 Seismic analysis of horseshoe tunnels under dynamic loads due to earthquakes Navid

More information

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS

ALGORITHM FOR NON-PROPORTIONAL LOADING IN SEQUENTIALLY LINEAR ANALYSIS 9th International Conference on Fracture Mechanics of Concrete and Concrete Structures FraMCoS-9 Chenjie Yu, P.C.J. Hoogenboom and J.G. Rots DOI 10.21012/FC9.288 ALGORITHM FOR NON-PROPORTIONAL LOADING

More information

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in

Sabah Shawkat Cabinet of Structural Engineering Walls carrying vertical loads should be designed as columns. Basically walls are designed in Sabah Shawkat Cabinet of Structural Engineering 17 3.6 Shear walls Walls carrying vertical loads should be designed as columns. Basically walls are designed in the same manner as columns, but there are

More information

BEM-FEM DYNAMIC ANALYSES OF ARCH DAM-FLUID SYSTEM

BEM-FEM DYNAMIC ANALYSES OF ARCH DAM-FLUID SYSTEM BEM-FEM DYNAMIC ANALYSES OF ARCH DAM-FLUID SYSTEM Violeta MIRCEVSKA 1 and Miroslav NASTEV 2 ABSTRACT The dam-fluid interaction is significantly affected by the irregularity of the terrain in the vicinity

More information

Numerical Solution of Equation of Motion

Numerical Solution of Equation of Motion Class Notes: Earthquake Engineering, Ahmed Elgamal, September 25, 2001 (DRAFT) Numerical Solution of Equation of Motion Average Acceleration Method (Trapezoidal method) m a + c v + k d = f (t) In the above

More information

Effect of Mass Matrix Formulation Schemes on Dynamics of Structures

Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Effect of Mass Matrix Formulation Schemes on Dynamics of Structures Swapan Kumar Nandi Tata Consultancy Services GEDC, 185 LR, Chennai 600086, India Sudeep Bosu Tata Consultancy Services GEDC, 185 LR,

More information

Theoretical Manual Theoretical background to the Strand7 finite element analysis system

Theoretical Manual Theoretical background to the Strand7 finite element analysis system Theoretical Manual Theoretical background to the Strand7 finite element analysis system Edition 1 January 2005 Strand7 Release 2.3 2004-2005 Strand7 Pty Limited All rights reserved Contents Preface Chapter

More information

Dynamic Analysis Using Response Spectrum Seismic Loading

Dynamic Analysis Using Response Spectrum Seismic Loading Dynamic Analysis Using Response Spectrum Seismic Loading Paleti Teja M.Tech (Structural Engineering) Department of Civil Engineering Jogaiah Institute of Technology & Sciences College of Engineering Kalagampudi,

More information

To initiate a dynamic analysis in FormWorks, the Analysis Mode in the JOB CONTROL tab of DEFINE JOB window has to be set to one of the two Dynamic

To initiate a dynamic analysis in FormWorks, the Analysis Mode in the JOB CONTROL tab of DEFINE JOB window has to be set to one of the two Dynamic Dynamic Analysis with VecTor2 This bulletin describes the dynamic analysis facilities within VecTor2. Details of VecTor2 and FormWorks can be found in the VecTor2 and FormWorks Manual (Vecchio and Wong,

More information

Finite Element Method in Geotechnical Engineering

Finite Element Method in Geotechnical Engineering Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps

More information

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06

Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras. Lecture - 06 Finite Element Analysis Prof. Dr. B. N. Rao Department of Civil Engineering Indian Institute of Technology, Madras Lecture - 06 In the last lecture, we have seen a boundary value problem, using the formal

More information

ON THE PREDICTION OF EXPERIMENTAL RESULTS FROM TWO PILE TESTS UNDER FORCED VIBRATIONS

ON 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 information

Earthquake response analysis of rock-fall models by discontinuous deformation analysis

Earthquake 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 information

Abstract. 1 Introduction

Abstract. 1 Introduction Contact analysis for the modelling of anchors in concrete structures H. Walter*, L. Baillet** & M. Brunet* *Laboratoire de Mecanique des Solides **Laboratoire de Mecanique des Contacts-CNRS UMR 5514 Institut

More information

Comparison of 2D vs. 3D Modeling Approaches for the Analyses of Concrete Faced Rockfill Dams

Comparison of 2D vs. 3D Modeling Approaches for the Analyses of Concrete Faced Rockfill Dams Comparison of 2D vs. 3D Modeling Approaches for the Analyses of Concrete Faced Rockfill Dams H. F. Özel & Y. Arici Dept. of Civil Engineering, Structural Mechanics Laboratory,K2-303, Middle East Technical

More information

Dynamic Soil Structure Interaction

Dynamic 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 information

Comparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties

Comparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties Comparison between the visco-elastic dampers And Magnetorheological dampers and study the Effect of temperature on the damping properties A.Q. Bhatti National University of Sciences and Technology (NUST),

More information

ANALYSIS OF HIGHRISE BUILDING STRUCTURE WITH SETBACK SUBJECT TO EARTHQUAKE GROUND MOTIONS

ANALYSIS OF HIGHRISE BUILDING STRUCTURE WITH SETBACK SUBJECT TO EARTHQUAKE GROUND MOTIONS ANALYSIS OF HIGHRISE BUILDING SRUCURE WIH SEBACK SUBJEC O EARHQUAKE GROUND MOIONS 157 Xiaojun ZHANG 1 And John L MEEK SUMMARY he earthquake response behaviour of unframed highrise buildings with setbacks

More information

NON-LINEAR RESPONSE OF A CONCRETE GRAVITY DAM SUBJECTED TO NEAR-FIELD GROUND MOTIONS

NON-LINEAR RESPONSE OF A CONCRETE GRAVITY DAM SUBJECTED TO NEAR-FIELD GROUND MOTIONS NON-LINEAR RESPONSE OF A CONCRETE GRAVITY DAM SUBJECTED TO NEAR-FIELD GROUND MOTIONS by K. Jagan Mohan, Pradeep Kumar Ramancharla in 15th Symposium on Earthquake Engineering (15SEE) Report No: IIIT/TR/2014/-1

More information

Iterative coupling in fluid-structure interaction: a BEM-FEM based approach

Iterative coupling in fluid-structure interaction: a BEM-FEM based approach oundary Elements and Other Mesh Reduction Methods XXV 205 terative coupling in fluid-structure interaction: a EM-EM based approach D. Soares Jr 1,2, W. J. Mansur 1 & O. von Estorff 3 1 Department of Civil

More information

INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA)

INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA) INTRODUCCION AL ANALISIS DE ELEMENTO FINITO (CAE / FEA) Title 3 Column (full page) 2 Column What is Finite Element Analysis? 1 Column Half page The Finite Element Method The Finite Element Method (FEM)

More information

Centrifuge Shaking Table Tests and FEM Analyses of RC Pile Foundation and Underground Structure

Centrifuge Shaking Table Tests and FEM Analyses of RC Pile Foundation and Underground Structure Centrifuge Shaking Table s and FEM Analyses of RC Pile Foundation and Underground Structure Kenji Yonezawa Obayashi Corporation, Tokyo, Japan. Takuya Anabuki Obayashi Corporation, Tokyo, Japan. Shunichi

More information

INELASTIC SEISMIC DISPLACEMENT RESPONSE PREDICTION OF MDOF SYSTEMS BY EQUIVALENT LINEARIZATION

INELASTIC SEISMIC DISPLACEMENT RESPONSE PREDICTION OF MDOF SYSTEMS BY EQUIVALENT LINEARIZATION INEASTIC SEISMIC DISPACEMENT RESPONSE PREDICTION OF MDOF SYSTEMS BY EQUIVAENT INEARIZATION M. S. Günay 1 and H. Sucuoğlu 1 Research Assistant, Dept. of Civil Engineering, Middle East Technical University,

More information

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method

Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method 9210-203 Level 7 Postgraduate Diploma in Engineering Computational mechanics using finite element method You should have the following for this examination one answer book No additional data is attached

More information

ANALYSIS OF EARTHQUAKE SAFETY OF A LARGE ARCH DAM

ANALYSIS OF EARTHQUAKE SAFETY OF A LARGE ARCH DAM ANALYSIS OF EARTHQUAKE SAFETY OF A LARGE ARCH DAM Miodrag SEKULOVIC 1, Radenko PEJOVIC 2 And Radivoje MRDAK 3 SUMMARY In this paper, safety assessment of high reinforce concrete arch dam, built in seismic

More information

Effects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device

Effects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device Effects of Damping Ratio of Restoring force Device on Response of a Structure Resting on Sliding Supports with Restoring Force Device A. Krishnamoorthy Professor, Department of Civil Engineering Manipal

More information

Title. Author(s)DONG, Q.; OKAZAKI, T.; MIDORIKAWA, M.; RYAN, K.; SAT. Issue Date Doc URL. Type. Note. File Information BEARINGS

Title. Author(s)DONG, Q.; OKAZAKI, T.; MIDORIKAWA, M.; RYAN, K.; SAT. Issue Date Doc URL. Type. Note. File Information BEARINGS Title ANALYSIS OF SHAKE-TABLE TESTS OF A FULL-SCALE BUILDI BEAINGS Author(s)DONG, Q.; OKAZAKI, T.; MIDOIKAWA, M.; YAN, K.; SAT Issue Date -9- Doc UL http://hdl.handle.net// Type proceedings Note The Thirteenth

More information

Fujinuma Dam Performance during 2011 Tohoku Earthquake, Japan and Failure Mechanism by FEM

Fujinuma Dam Performance during 2011 Tohoku Earthquake, Japan and Failure Mechanism by FEM Fujinuma Dam Performance during 2011 Tohoku Earthquake, Japan and Failure Mechanism by FEM Mahdavian Abbas Powue and Water University of Technology, Tehran, Iran Shiro Takada Tehran University, Tehran,

More information

Structural Damage Detection Using Time Windowing Technique from Measured Acceleration during Earthquake

Structural 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 information

CIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass

CIVL 8/7117 Chapter 12 - Structural Dynamics 1/75. To discuss the dynamics of a single-degree-of freedom springmass CIV 8/77 Chapter - /75 Introduction To discuss the dynamics of a single-degree-of freedom springmass system. To derive the finite element equations for the time-dependent stress analysis of the one-dimensional

More information

Stretching of a Prismatic Bar by its Own Weight

Stretching of a Prismatic Bar by its Own Weight 1 APES documentation (revision date: 12.03.10) Stretching of a Prismatic Bar by its Own Weight. This sample analysis is provided in order to illustrate the correct specification of the gravitational acceleration

More information

Frequency response analysis of soil-structure interaction for concrete gravity dams

Frequency response analysis of soil-structure interaction for concrete gravity dams Frequency response analysis of soil-structure interaction for concrete gravity dams Anna De Falco 1, Matteo Mori 2 and Giacomo Sevieri 3 1 Dept. of Energy, Systems, Territory and Construction Engineering,

More information

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface

Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield Surface Transaction A: Civil Engineering Vol. 16, No. 6, pp. 512{519 c Sharif University of Technology, December 2009 Research Note Nonlinear FE Analysis of Reinforced Concrete Structures Using a Tresca-Type Yield

More information

Limit 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 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 information