Dynamic Response of Wind Turbines taking into account the Rotational Component of Ground Motions
|
|
- Willis Doyle
- 6 years ago
- Views:
Transcription
1 Dynamic Response of Wind Turbines taking into account the Rotational Component of Ground Motions L. Hermanns, L.E. Quiros, J. Gaspar-Escribano & B. Benito Universidad Politecnica de Madrid,Spain M.A. Santoyo Universidad Complutense de Madrid, Spain J. Vega CEMIM-F2I2, Madrid, Spain SUMMARY: This work studies the dynamic response of a wind turbine structure under seismic loading, taking into account the rotational component of ground shaking. The study covers both the definition of synthetic input ground motion, and the analysis of the dynamic response. Based on previous studies, carried out in the framework of a research project, and the common seismological characteristics of shallow earthquakes, the synthetic seismograms were generated using an appropriate set of source properties (e.g. focal mechanism, source time function, etc...) and an earth subsurface velocity structure for the region, taking into account the tectonic environment of the modelling site. Synthetic time histories have been computed using the Discrete Wavenumber Method, assuming a point source and a horizontal layered earth structure. Results will be analyzed in order to assess the factors that lead to scenarios in which the contribution of the rotational components to the overall response cannot be neglected. Keywords: Rotational component, Wind Turbine, Discrete Wavenunber Method 1. INTRODUCTION Rotational ground motion analysis is a novel discipline of seismology that still presents significant unexplored fields. Accordingly, the development of instruments for recording rotational motions is limited. Thus, strong motion networks do not routinely measure these parameters. In this paper a set of complete synthetic ground motion time histories (i. e., including both translational and rotational components) was used in order to evaluate the seismic response of a wind turbine. This structure is interesting because it presents a large inertial mass at the top and a rotating component around a horizontal axis. Different methods to calculate the dynamic response of the structure are examined. At this stage of the analysis, a Finite Element model seems a good choice for performing transient dynamic calculations. The structural response of the wind turbine is characterized by the maximum value of Von Mises stress at different levels of the tower. The target turbine site corresponds to a hypothetic site located at Madrid Region (Spain), where moderate events could take place. Values of the corresponding seismological parameters are selected accordingly. However, the upper soil layer of the Earth velocity model is very stiff. This should always be remembered when interpreting the results. 2. SEISMICITY AND SEISMIC HAZARD IN CENTRAL SPAIN Central Spain is characterized by a relatively low seismic activity. The largest magnitude ever recorded corresponds to the 2007, m blg 5.1 Pedro Muñoz earthquake and the largest intensity ever reported is the V and corresponds to different events occurred in the province of Cuenca during the first part of the XX century (Garaballa 1919, Motilla del Palancar 1929). Accordingly, seismic hazard is low and current regulations do not enforce the development of any special plan for civil protection
2 against seismic risk. However, there are tectonic structures with neotectonic activity in the area (Giner et al. 1996) with potential of generating a moderate event, as the faults bounding the Altomira Range (Figure 1). Figure 1. Geologic map of the study area in Central Spain. A limitation for ground motion characterization exists due to the limited availability of data on earthquake sources (e. g., focal mechanisms and fault slip rates) and of strong motion records. Therefore, synthetic ground-motion time histories are developed to model the ground shaking (both translation and rotation motions) at the wind turbine site. 3. SYNTHETIC GROUND MOTIONS 3.1. Ground Motion Definition The seismic inputs used in this work (translational and rotational time histories) were theoretically generated for an earth layered structure embedding a double couple point source at shallow depths. The model takes into account the geometric and the inelastic attenuation of the medium (Q) for both P and S waves. Seismic events were characterized by a specific focal mechanism, moment magnitude and a Source-Time Function (STF) equivalent to an earthquake with circular finite rupture. The synthetic seismograms (accelerograms) were computed using the discrete wave-number method described by Bouchon and Aki (1977) and Bouchon (1979) for point dislocations. The velocity model assumed corresponds to the crustal properties of the south-eastern region of the Madrid Region, Spain (Corchete and Chourak, 2010; Table K). Table 3.1. Crustal Structure used in this study Layer Z (m) Vp (m/s) Vs (m/s) ρ (kg/m 3 ) Qp Qs The set of synthetic seismograms used here were generated for four different simulation groups (A, B, C and D) each of them with 30 different STF. Group A was simulated also for 30 different focal mechanisms randomly generated. Groups B, C and D were simulated using the mechanisms shown in table 3.2.
3 Table 3.2. Focal mechanisms for simulations B, C and D. Group Strike Dip Rale B C D The epicentre of the earthquake was located 23 km northeast of the wind turbine site. A Mw=5.5 earthquake at 6.5 km depth was assumed, which is the mean depth of seismicity in this region. The set of STF were randomly generated taking into account the mean properties of shallow crustal earthquakes of this magnitude (Tanioka and Ruff, 1997; Houston, 2001). STF were constructed assuming a lognormal shape (Liao y Huang, 2005), with skewness 0.15<S k <0.35 (Houston, 2001). They were added with 10% and 20% of white noise, low pass filtered at a f c =25 hz and tapered at its extremes to avoid possible alias effects. The maximum total duration of these functions was τ=2.0s, with a mean duration of τ=1.5s at the 90% of the total moment release. Focal mechanisms were generated taking into account the tectonic setting of the region. In Figure 1 an example of two of the STF used in the study is shown. Figure 2. Source Time Functions for two of the simulated earthquakes. Seismograms were synthesized in accelerations and displacements for each of the four groups and each STF simulation at four triaxial receivers located at the wind turbine placement. Seismograms were computed up to a frequency of 75.0 Hz and low pass filtered with a 4 pole Butterworth filter with a cut-off frequency of 40 Hz. From here the highest usable frequency is about 35.0 Hz. Once seismograms were generated, the time histories of rotations and angular accelerations were obtained by a first-order finite difference approach for the computation of the spatial derivatives of the ground motion. Here, the complete tensor of synthetic displacement gradients G i,j =Δu i /Δx j, where Δu i,i =x,y,z, are the differences in the synthetic displacements in the three spatial directions and Δx j, j=1,2,3 are the space increment in x, y and z directions respectively. The synthetic displacement time histories were computed again by the discrete wave-number method (Bouchon 1979). The spacing among computing points was taken as Δx j = λ min /100=V s /100 f N, where Δx j =Δx 1 =Δx 2 =Δx 3, V s = minimum S wave velocity of the considered velocity structure, and f N =maximum analyzed frequency. The time increment for this analysis Δt=1/f s =1/2f N. Once obtained the components of G i,j, these were used to derive the rigid body rotations by: 1 u u i j ω = i, j (5.1) 2 x j xi and the angular accelerations by:
4 2 ωi, j 2 ω i, j = (5.2) t where i=1,2,3; j=1,2,3; u 1 =u, u 2 =v, u 3 =w; u, v and w are the displacements in the x, y and z directions at a given time. At the surface, due to the stress free boundary conditions, three components of G i,j are not independent: u 1 / x 3 = - u 3 / x 1, u 2 / x 3 = - u 3 / x 2 and u 3 / x 3 = η ( u 2 / x 2 + u 1 / x 1 ) and η=-λ/ (λ+2µ); λ and µ are the Lamé parameters. From here, ε 1,3 =ε 3,1 = ε 2,3 = ε 3,2 =0. Cotton and Coutant (1997) tested this methodology with respect to the solution using analytic spatial derivatives, obtaining an excellent agreement with almost undistinguishable results between both techniques Sensitivity Analysis A number of model parameters characterizing the fault, such as type of focal mechanism, source-time function and expected magnitudes are varied in order to study the sensitivity of the model to these input parameters. All cases tested show that the predominant rotational component is around the Y axis, i. e., the axis connecting source and site locations (Figure 3). Figure 3. Rotational ground-motion time histories for a magnitude 6.0 event (source-to-site distance and depth are 20 and 7 km, respectively. The rotational effect varies with magnitude and source-to-site distance. Rotation around Y axis for a magnitude 6.0 event is about one order of magnitude higher than for a magnitude 5.0 event. For relatively short distances (about 10 km), permanent rotation might be reached for moderate-high magnitude earthquakes (Figure 4). Figure 4. Rotational ground-motion time histories as a function of epicentral distance, for a magnitude 6.0 event and a depth of 7 km. The impact of different focal mechanisms on modeled ground motions turned out to be significant and made it difficult to identify any variation pattern. Subsequent tests were made fixing the focal mechanism and varying source time functions, ensuring that the final STF was consistent with the
5 chosen magnitude of FINITE ELEMENT MODEL OF THE WIND TURBINE Of course the simplest dynamical model has only 1 degree of freedom (dof). However, this type of model is not suited to represent adequately the dynamic response of a wind turbine during a seismic event. Much more sophisticated models that include several soil layers represent the opposite extreme leading to prohibitive calculation times. In the latter case 3d continuum elements are used for the soil, shell elements for the tower as well as the blades and concentrated masses to represent heavy equipment located in the nacelle. Fig. 5 presents some views of a wireframe model and a photograph. A good compromise between accuracy and calculation times can be achieved by choosing a Finite Element (FE) model that represents a medium level of sophistication. To this end the soil s stiffness and damping properties are approximated by means of spring and damper elements that are located at the bottom of the footing. Apart from the spring and damper elements additional masses can also be included in order to account for the surrounding soil that accompanies the motion of the footing. The tower and the blades are modelled using beam elements. If the forces in the blades are not of special interest, rotor, nacelle and the heavy equipment may be modelled by means of concentrated masses at the top of the tower. The resulting FE model permits to study the contribution of the different bending modes to the overall response and allows studying the effects of changes in footing dimensions or soil properties. Therefore this type of model has been used for the simulations. Figure 5. Wireframe model and photograph of the wind turbine MADE AE-46/1 As already mentioned in previous sections, the upper soil layer is very, very stiff. This should always be remembered when interpreting results presented in the present paper. The highest excitation frequency is limited to approximately 30 Hz. With this in mind, a convergence study has been performed in order to end up with a suitable FE mesh. The varying diameter and thickness of the tower lead to continuously changing mass and stiffness properties between top and bottom. The used beam elements have a constant cross-section which increases the required number of elements. Finally, a mesh with 45 elements of equal length has been used for the simulations. The characteristics of the fist eight modes of vibration are presented in table 4.1. Table 4.2. Characteristics of the modes of vibration of the FE model. Number Type Frequency (Hz) Period (s) 1 Bending x direction Bending z direction Bending x direction Bending z direction Axial transl. y dirct Bending x direction Bending z direction Torsion rot y
6 Rayleigh damping equivalent to 1 % of critical damping has been used using the frequencies of the first and third modes. 5. RESULTS 5.1. Results considering the first set of ground motions A total of 60 simulations have been carried out, two for each of the 30 sets of time history accelerations: once considering that the ground motion produces only motion on the translational dof at the base, and a second one with additional motion on the rotational dof. For a given seismic input, the results with or without rotational seismic input may be very similar. Comparisons may be established in terms of the maximum value of Von Mises stresses at a given section of the tower. Fig.6 (a) shows this magnitude as a function of the height of the tower considering both kinds of seismic input for a given set of time history accelerations. In this case differences are smaller than 10%. Fig. 6 (b) displays the same plot for a different set. This time differences are much more important, indicating that for this set the base rotations do have a big influence on the response of the tower. Figure 6. Comparison of the distribution of maximum Von Mises stress along the tower's height with or without rotational input for records: (a) no.12 and (b) no.18. The observed differences between responses caused by both records can be easily explained by the means of response spectra. They are displayed in figure 7 for horizontal ground motions and rotations around the horizontal axis. The differences in the radiation pattern between both cases as well as resonance effects between vibration frequencies of the top soil layer and main structural modes may explain the discrepancies. Figure 7. Comparison of the distribution of maximum Von Mises stress along the tower's height with or without rotational input for records: (a) no.12 and (b) no.18.
7 In order to assess the differences obtained between simulations with or without the rotational seismic input, results from all simulations have been analyzed globally. For each section (height) the maximum value of the Von Mises stress resulting from each simulation has been considered as a realization of a random variable X. It has been assumed that X distributes as a lognormal random variable, and its two parameters have been estimated from the 30 samples available. Fig. 8 shows plots of the probability density functions for both cases at three different heights. At the base (0m) and around 30m both distributions are quite different. However this is not the case at 15m. The reason is that that modes number 3 and 4, which are the ones that are excited significantly by base rotations, have reduced moments at this height Figure 8. Comparison of probability density functions of maximum Von Mises stress at sections placed at different heights of the tower. The obtained values indicate that there is a large scatter among results. Some records display important differences while others do not fit this point. It is not clear whether important differences are associated to particular focal mechanism or to particular SF. This question is addressed next Results considering sets B, C and D For each record of the new sets, again two simulations have been done: one considering the rotational seismic components, and another without them. The importance of the rotational components will be assessed using the following ratio max max ( σ VM, disp rot ) ( σ VM, disp ) +. (5.3) High values of the ratio indicate that rotations are important, while values close to unity indicate the opposite. Figure 9 displays the cumulative distribution function (CDF) of the ratios for the three sets. Set C has an almost vertical CDF indicating that all the ratios are close to 1. Discrepancies between results with or without rotational components are noticeable for the records in set B, and remarkably important for those in set D.
8 Figure 9. Comparison of cumulative distribution functions or the ratio between the maximum Von Mises stresses obtained with and without rotational seismic components. 6. CONCLUSIONS The analysis of synthetic time histories of translational and rotational motions indicates the significant effect of the rotational components of motion, especially for short distances and moderate magnitudes. The results obtained in this study indicate the importance of the rotational components of ground motions on wind turbine response patterns. The importance of the rotational components is directly related to its magnitude, which is highly dependent on the focal mechanisms. Source time functions have proven to also have an effect on it, but it only arises with suitable focal mechanisms. AKCNOWLEDGEMENT The work presented in this paper is partly funded by the Universidad Politécnica de Madrid and by the Comunidad Autónoma de Madrid. REFERENCES Bouchon, M. and K. Aki (1977). Discrete wave-number representation of seismic-source wave fields. Bulletin of the Seismological Society of America 67:2, Bouchon, M. (1979), discrete wavenumber representation of elastic wave fields in three-space dimensions, Journal of Geophysical Research 84, Cotton F., Coutant O. (1997) Dynamic stress variations due to shear faults in a plane-layered medium. Geophys. J. Int.,128, Corchete V. and M. Chourak (2011). Shear-wave velocity structure of the south-eastern part of the Iberian Peninsula from Rayleigh wave analysis. International Journal of Earth Sciences, 100, Giner, J. L., G. De Vicente, A. C. Pérez González, J. G. Sánchez Cabañero, and L. Pinilla (1996). Crisis tectónicas cuaternarias en la Cuenca de Madrid. Geogaceta 20:4, Houston, H. (2001). Influence of depth, focal mechanism, and tectonic setting on the shape and duration of earthquake source time functions. Journal of Geophysical Research 106:b6, 11,137-11,150. Liao B. Y. and Huang H.C. (2005). Estimation of the source time function based on blind deconvolution with Gaussian mixtures. Pure Appl. Geophys. PAGEOPH 162, doi /s z; /05/ Tanioka, Y. and L. J. Ruff (1997) Source Time Functions. Seismological Research Letters 68:3,
Outstanding Problems. APOSTOLOS S. PAPAGEORGIOU University of Patras
NEAR-FAULT GROUND MOTIONS: Outstanding Problems APOSTOLOS S. PAPAGEORGIOU University of Patras Outline Characteristics of near-fault ground motions Near-fault strong ground motion database A mathematical
More informationDynamic Loads CE 543. Examples. Harmonic Loads
CE 543 Structural Dynamics Introduction Dynamic Loads Dynamic loads are time-varying loads. (But time-varying loads may not require dynamic analysis.) Dynamics loads can be grouped in one of the following
More informationFOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES
FOCAL MECHANISM DETERMINATION USING WAVEFORM DATA FROM A BROADBAND STATION IN THE PHILIPPINES Vilma Castillejos Hernandez Supervisor: Tatsuhiko Hara MEE10508 ABSTRACT We performed time domain moment tensor
More informationDamping Scaling of Response Spectra for Shallow CCCCCCCCCrustalstallPaper Crustal Earthquakes in Active Tectonic Title Line Regions 1 e 2
Damping Scaling of Response Spectra for Shallow CCCCCCCCCrustalstallPaper Crustal Earthquakes in Active Tectonic Title Line Regions 1 e 2 S. Rezaeian U.S. Geological Survey, Golden, CO, USA Y. Bozorgnia
More informationPart 2 - Engineering Characterization of Earthquakes and Seismic Hazard. Earthquake Environment
Part 2 - Engineering Characterization of Earthquakes and Seismic Hazard Ultimately what we want is a seismic intensity measure that will allow us to quantify effect of an earthquake on a structure. S a
More informationSTUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING
STUDYING THE IMPORTANT PARAMETERS IN EARTHQUAKE SIMULATION BASED ON STOCHASTIC FINITE FAULT MODELING H. Moghaddam 1, N. Fanaie 2* and H. Hamzehloo 1 Professor, Dept. of civil Engineering, Sharif University
More informationSEISMIC RESPONSE OF SINGLE DEGREE OF FREEDOM STRUCTURAL FUSE SYSTEMS
3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 377 SEISMIC RESPONSE OF SINGLE DEGREE OF FREEDOM STRUCTURAL FUSE SYSTEMS Ramiro VARGAS and Michel BRUNEAU
More informationEXAMPLE OF PILED FOUNDATIONS
EXAMPLE OF PILED FOUNDATIONS The example developed below is intended to illustrate the various steps involved in the determination of the seismic forces developed in piles during earthquake shaking. The
More informationVALIDATION AGAINST NGA EMPIRICAL MODEL OF SIMULATED MOTIONS FOR M7.8 RUPTURE OF SAN ANDREAS FAULT
VALIDATION AGAINST NGA EMPIRICAL MODEL OF SIMULATED MOTIONS FOR M7.8 RUPTURE OF SAN ANDREAS FAULT L.M. Star 1, J. P. Stewart 1, R.W. Graves 2 and K.W. Hudnut 3 1 Department of Civil and Environmental Engineering,
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 informationIZMIT BAY BRIDGE SOUTH APPROACH VIADUCT: SEISMIC DESIGN NEXT TO THE NORTH ANATOLIAN FAULT
Istanbul Bridge Conference August 11-13, 2014 Istanbul, Turkey IZMIT BAY BRIDGE SOUTH APPROACH VIADUCT: SEISMIC DESIGN NEXT TO THE NORTH ANATOLIAN FAULT A. Giannakou 1, J. Chacko 2 and W. Chen 3 ABSTRACT
More informationEARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE MECHANISMS SHOW MOTION
6-1 6: EARTHQUAKE FOCAL MECHANISMS AND PLATE MOTIONS Hebgen Lake, Montana 1959 Ms 7.5 1 Stein & Wysession, 2003 Owens Valley, California 1872 Mw ~7.5 EARTHQUAKE LOCATIONS INDICATE PLATE BOUNDARIES EARTHQUAKE
More informationSynthetic Seismicity Models of Multiple Interacting Faults
Synthetic Seismicity Models of Multiple Interacting Faults Russell Robinson and Rafael Benites Institute of Geological & Nuclear Sciences, Box 30368, Lower Hutt, New Zealand (email: r.robinson@gns.cri.nz).
More informationDynamic behavior of turbine foundation considering full interaction among facility, structure and soil
Dynamic behavior of turbine foundation considering full interaction among facility, structure and soil Fang Ming Scholl of Civil Engineering, Harbin Institute of Technology, China Wang Tao Institute of
More informationSpatial Correlation of Ground Motions in Seismic Hazard Assessment
Spatial Correlation of Ground Motions in Seismic Hazard Assessment Taojun Liu tliu82@uwo.ca Department of Civil & Environmental Engineering University of Western Ontario London, Ontario, Canada 1 Outline
More informationCodal 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 informationON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL
1 Best Practices in Physics-based Fault Rupture Models for Seismic Hazard Assessment of Nuclear ON NEAR-FIELD GROUND MOTIONS OF NORMAL AND REVERSE FAULTS FROM VIEWPOINT OF DYNAMIC RUPTURE MODEL Hideo AOCHI
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lecture 20, 30 Nov. 2017 www.geosc.psu.edu/courses/geosc508 Seismic Spectra & Earthquake Scaling laws. Seismic Spectra & Earthquake Scaling laws. Aki, Scaling law
More informationEmpirical Green s Function Analysis of the Wells, Nevada, Earthquake Source
Nevada Bureau of Mines and Geology Special Publication 36 Empirical Green s Function Analysis of the Wells, Nevada, Earthquake Source by Mendoza, C. 1 and Hartzell S. 2 1 Centro de Geociencias, Universidad
More informationCE6701 STRUCTURAL DYNAMICS AND EARTHQUAKE ENGINEERING QUESTION BANK UNIT I THEORY OF VIBRATIONS PART A
CE6701 STRUCTURAL DYNAMICS AND EARTHQUAKE ENGINEERING QUESTION BANK UNIT I THEORY OF VIBRATIONS PART A 1. What is mean by Frequency? 2. Write a short note on Amplitude. 3. What are the effects of vibration?
More informationDynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement
Emily Morton Geop 523: Theoretical Seismology 29 April 2011 Introduction Dynamic Earthquake Triggering Due to Stress from Surface Wave Particle Displacement Earthquakes can be triggered by other earthquakes.
More informationTime-Varying Dynamic Properties of Offshore Wind Turbines Evaluated by Modal Testing
Time-Varying Dynamic Properties of Offshore Wind Turbines Evaluated by Modal Testing M. Damgaard *, L.V. Andersen Ɨ, L.B. Ibsen Ɨ * Technology and Engineering Solutions, Vestas Wind Systems A/S, Denmark
More informationDIRECT HAZARD ANALYSIS OF INELASTIC RESPONSE SPECTRA
DIRECT HAZARD ANALYSIS OF INELASTIC RESPONSE SPECTRA ABSTRACT Y. Bozorgnia, M. Hachem, and K.W. Campbell Associate Director, PEER, University of California, Berkeley, California, USA Senior Associate,
More informationOn the Horizontal-to-Vertical Spectral Ratio in Sedimentary Basins
Bulletin of the Seismological Society of America, 90, 4, pp. 1101 1106, August 2000 On the Horizontal-to-Vertical Spectral Ratio in Sedimentary Basins by Zakaria Al Yuncha and Francisco Luzón Abstract
More informationRelevance of Fault-Normal/Parallel and Maximum Direction Rotated Ground Motions on Nonlinear Behavior of Multi-Story Buildings
Relevance of Fault-Normal/Parallel and Maximum Direction Rotated Ground Motions on Nonlinear Behavior of Multi-Story Buildings J.C. Reyes Universidad de los Andes, Bogotá, Colombia E. Kalkan United States
More information7 SEISMIC LOADS. 7.1 Estimation of Seismic Loads. 7.2 Calculation of Seismic Loads
1 7 SEISMIC LOADS 7.1 Estimation of Seismic Loads 7.1.1 Seismic load and design earthquake motion (1) For ordinary buildings, seismic load is evaluated using the acceleration response spectrum (see Sec.7.2)
More informationEMEA. Liudmila Feoktistova Engineer Atomenergoproekt
EMEA Liudmila Feoktistova Engineer Atomenergoproekt Dr. Adolf Birbraer, Atomenergoproekt Dr. Alexandr Roleder, Atomenergoproekt Dmitry Mikhaluk, St. Petersburg State Polytechnical University Finite element
More informationDr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum
STRUCTURAL DYNAMICS Dr.Vinod Hosur, Professor, Civil Engg.Dept., Gogte Institute of Technology, Belgaum Overview of Structural Dynamics Structure Members, joints, strength, stiffness, ductility Structure
More informationRESPONSE SPECTRUM METHOD FOR ESTIMATION OF PEAK FLOOR ACCELERATION DEMAND
RESPONSE SPECTRUM METHOD FOR ESTIMATION OF PEAK FLOOR ACCELERATION DEMAND Shahram Taghavi 1 and Eduardo Miranda 2 1 Senior catastrophe risk modeler, Risk Management Solutions, CA, USA 2 Associate Professor,
More informationStructural Dynamics Lecture 4. Outline of Lecture 4. Multi-Degree-of-Freedom Systems. Formulation of Equations of Motions. Undamped Eigenvibrations.
Outline of Multi-Degree-of-Freedom Systems Formulation of Equations of Motions. Newton s 2 nd Law Applied to Free Masses. D Alembert s Principle. Basic Equations of Motion for Forced Vibrations of Linear
More informationFrequency 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 informationPHASE ANGLE PROPERTIES OF EARTHQUAKE STRONG MOTIONS: A CRITICAL LOOK
565 PHASE ANGLE PROPERTIES OF EARTHQUAKE STRONG MOTIONS: A CRITICAL LOOK B TILIOUINE 1, M HAMMOUTENE And P Y BARD 3 SUMMARY This paper summarises the preliminary results of an investigation aimed at identifying
More informationRepresentative ground-motion ensembles for several major earthquake scenarios in New Zealand
Representative ground-motion ensembles for several major earthquake scenarios in New Zealand K. Tarbali & B.A. Bradley Department of Civil and Natural Resources Engineering, University of Canterbury, Christchurch.
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 33 Probabilistic methods in earthquake engineering-2 So, we have
More informationResponse Spectrum Analysis Shock and Seismic. FEMAP & NX Nastran
Response Spectrum Analysis Shock and Seismic FEMAP & NX Nastran Table of Contents 1. INTRODUCTION... 3 2. THE ACCELEROGRAM... 4 3. CREATING A RESPONSE SPECTRUM... 5 4. NX NASTRAN METHOD... 8 5. RESPONSE
More informationDesign Procedures For Dynamically Loaded Foundations
Design Procedures For Dynamically Loaded Foundations 1/11 Choice of parameters for equivalent lumped systems Lumped mass : the mass of the foundation and supported machinery Damping : 1 Geometrical(or
More informationSPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES
2010/2 PAGES 1 8 RECEIVED 21. 9. 2009 ACCEPTED 20. 1. 2010 Y. KOLEKOVÁ, M. PETRONIJEVIĆ, G. SCHMID SPECIAL DYNAMIC SOIL- STRUCTURE ANALYSIS PROCEDURES DEMONSTATED FOR TWO TOWER-LIKE STRUCTURES ABSTRACT
More informationCHARACTERISTICS OF SOURCE SPECTRA OF SMALL AND LARGE INTERMEDIATE DEPTH EARTHQUAKES AROUND HOKKAIDO, JAPAN
13 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1861 CHARACTERISTICS OF SOURCE SPECTRA OF SMALL AND LARGE INTERMEDIATE DEPTH EARTHQUAKES AROUND HOKKAIDO,
More informationKeywords: Drab-e-Astane earthquake, Empirical green's function, Kostrov slip function, Simulation.
Synthesizing the 2006 Silakhor (Darb-e-Astane) earthquake (Iran), using Empirical Green s function method A. Nicknam, R. Abbasnia, Y. Eslamian, M. Bozorgnasab, fariborz Nicknam, M.A.Sadredini 1. Assistant
More informationEXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION
EXAMINATION ON CONSECUTIVE RUPTURING OF TWO CLOSE FAULTS BY DYNAMIC SIMULATION M. Muto 1, K. Dan 1, H. Torita 1, Y. Ohashi 1, and Y. Kase 2 1 Ohsaki Research Institute, Inc., Tokyo, Japan 2 National Institute
More informationVertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes
Vertical to Horizontal (V/H) Ratios for Large Megathrust Subduction Zone Earthquakes N.J. Gregor Consultant, Oakland, California, USA N.A. Abrahamson University of California, Berkeley, USA K.O. Addo BC
More informationOVERVIEW INTRODUCTION 3 WHAT'S MISSING? 4 OBJECTIVES 5
OVERVIEW INTRODUCTION 3 WHAT'S MISSING? 4 OBJECTIVES 5 DISTORTION OF SEISMIC SOURCE SPECTRUM 6 PRINCIPLE 7 SEISMIC SOURCE SPECTRUM 8 EFFECT OF RECORDING INSTRUMENTS 9 SEISMOMETERS 9 CORRECTION FOR FREQUENCY
More information7 Ground Motion Models
7 Ground Motion Models 7.1 Introduction Ground motion equations are often called attenution relations but they describe much more than just the attenutation of the ground motion; they describe the probability
More informationSURFACE WAVE MODELLING USING SEISMIC GROUND RESPONSE ANALYSIS
43 SURFACE WAVE MODELLING USING SEISMIC GROUND RESPONSE ANALYSIS E John MARSH And Tam J LARKIN SUMMARY This paper presents a study of surface wave characteristics using a two dimensional nonlinear seismic
More informationElastic rebound theory
Elastic rebound theory Focus epicenter - wave propagation Dip-Slip Fault - Normal Normal Fault vertical motion due to tensional stress Hanging wall moves down, relative to the footwall Opal Mountain, Mojave
More informationStochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore
Stochastic Structural Dynamics Prof. Dr. C. S. Manohar Department of Civil Engineering Indian Institute of Science, Bangalore Lecture No. # 32 Probabilistic Methods in Earthquake Engineering-1 (Refer Slide
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 informationModelling Strong Ground Motions for Subduction Events in the Wellington Region, New Zealand
Proceedings of the Ninth Pacific Conference on Earthquake Engineering Building an Earthquake-Resilient Society 14-16 April, 2011, Auckland, New Zealand Modelling Strong Ground Motions for Subduction Events
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 informationANALYSIS OF THE CORRELATION BETWEEN INSTRUMENTAL INTENSITIES OF STRONG EARTHQUAKE GROUND MOTION
ANALYSIS OF THE CORRELATION BETWEEN INSTRUMENTAL INTENSITIES OF STRONG EARTHQUAKE GROUND MOTION J.Enrique Martinez-Rueda 1, Evdokia Tsantali 1 1 Civil Engineering & Geology Division School of Environment
More informationWidespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake
Widespread Ground Motion Distribution Caused by Rupture Directivity during the 2015 Gorkha, Nepal Earthquake Kazuki Koketsu 1, Hiroe Miyake 2, Srinagesh Davuluri 3 and Soma Nath Sapkota 4 1. Corresponding
More informationGround-Motion Attenuation Relationships for Subduction- Zone Earthquakes in Northern Taiwan
Ground-Motion Attenuation Relationships for Subduction- Zone Earthquakes in Northern Taiwan Lin, P.S., Lee, C.T. Bulletin of the Seismology Society of America (2008) Presenter: Yang Pei-Xin Adviser: Lee
More informationSynthetic Near-Field Rock Motions in the New Madrid Seismic Zone
Synthetic Near-Field Rock Motions in the New Madrid Seismic Zone Genda Chen*, Ph.D., P.E., and Mostafa El-Engebawy Engebawy,, Ph.D. *Associate Professor of Civil Engineering Department of Civil, Architecture
More informationEARTHQUAKE CLUSTERS, SMALL EARTHQUAKES
EARTHQUAKE CLUSTERS, SMALL EARTHQUAKES AND THEIR TREATMENT FOR HAZARD ESTIMATION Gary Gibson and Amy Brown RMIT University, Melbourne Seismology Research Centre, Bundoora AUTHORS Gary Gibson wrote his
More informationSimulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the Pseudo Point-Source Model
6 th International Conference on Earthquake Geotechnical Engineering -4 November 25 Christchurch, New Zealand Simulation of Strong Ground Motions for a Shallow Crustal Earthquake in Japan Based on the
More informationEarthquake stress drop estimates: What are they telling us?
Earthquake stress drop estimates: What are they telling us? Peter Shearer IGPP/SIO/U.C. San Diego October 27, 2014 SCEC Community Stress Model Workshop Lots of data for big earthquakes (rupture dimensions,
More informationSOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION
Synopses of Master Papers Bulletin of IISEE, 47, 19-24, 2013 SOURCE PROCESS OF THE 2003 PUERTO PLATA EARTHQUAKE USING TELESEISMIC DATA AND STRONG GROUND MOTION SIMULATION Fabricio Moquete Everth* Supervisor:
More informationDepartment of Civil Engineering, Serbia
FACTA UNIVERSITATIS Series: Architecture and Civil Engineering Vol. 10, N o 2, 2012, pp. 131-154 DOI: 10.2298/FUACE1202131B TOWARDS PREPARATION OF DESIGN SPECTRA FOR SERBIAN NATIONAL ANNEX TO EUROCODE
More informationGrandstand Terraces. Experimental and Computational Modal Analysis. John N Karadelis
Grandstand Terraces. Experimental and Computational Modal Analysis. John N Karadelis INTRODUCTION Structural vibrations caused by human activities are not known to be particularly damaging or catastrophic.
More informationEarthquake Stochastic Response and Soil Type Variable Effects
Australian Journal of Basic and Applied Sciences, 5(8): 581-586, 211 ISSN 1991-8178 Earthquake Stochastic Response and Soil Type Variable Effects Zaniar Tokmechi Department of Civil Engineering, Science
More informationStructural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian
Structural Dynamics Lecture Eleven: Dynamic Response of MDOF Systems: (Chapter 11) By: H. Ahmadian ahmadian@iust.ac.ir Dynamic Response of MDOF Systems: Mode-Superposition Method Mode-Superposition Method:
More informationEstimation of Peak Ground Acceleration for Delhi Region using Finsim, a Finite Fault Simulation Technique
215 Estimation of Peak Ground Acceleration for Delhi Region using Finsim, a Finite Fault Simulation Technique NEELIMA SATYAM. D* and K. S. RAO** * Earthquake Engineering Research Centre, International
More informationEARTHQUAKE SOURCE PARAMETERS OF MODERATELY EARTHQUAKE IN THE SOUTH EASTERN IRAN BASED ON TELESEISMIC AND REGIONAL DISTANCES
EARTHQUAKE SOURCE PARAMETERS OF MODERATELY EARTHQUAKE IN THE SOUTH EASTERN IRAN BASED ON TELESEISMIC AND REGIONAL DISTANCES M. Mostafazadeh 1 and S.Ashkpour Motlagh 2 1 Assistant Professor,Seismology Research
More informationFINITE FAULT MODELING OF FUTURE LARGE EARTHQUAKE FROM NORTH TEHRAN FAULT IN KARAJ, IRAN
Proceedings of the 31 st Conference on Earthquake Engineering, JSCE, 2011 FINITE FAULT MODELING OF FUTURE LARGE EARTHQUAKE FROM NORTH TEHRAN FAULT IN KARAJ, IRAN Meghdad Samaei 1, Masakatsu Miyajima 2,
More informationIntroduction to Geotechnical Earthquake Engineering
Module 1 Introduction to Geotechnical Earthquake Engineering by Dr. Deepankar Choudhury Professor Department of Civil Engineering IIT Bombay, Powai, Mumbai 400 076, India. Email: dc@civil.iitb.ac.in URL:
More information2 Approaches To Developing Design Ground Motions
2 Approaches To Developing Design Ground Motions There are two basic approaches to developing design ground motions that are commonly used in practice: deterministic and probabilistic. While both approaches
More informationFoundation models for the dynamic response of offshore wind turbines
Marine Renewable Energy Conference (MAREC), Newcastle, UK, September 00. Foundation models for the dynamic response of offshore wind turbines. M.B. Zaaijer, MSc Delft University of Technology, The Netherlands
More informationTHE USE OF INPUT ENERGY FOR SEISMIC HAZARD ASSESSMENT WITH DIFFERENT DUCTILITY LEVEL
th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, Paper No. 8 THE USE OF INPUT ENERGY FOR SEISMIC HAZARD ASSESSMENT WITH DIFFERENT DUCTILITY LEVEL Mao-Sheng GONG And Li-Li
More informationModeling Of Earthquake Ground Motion for Tbilisi Region
P a g e 2 Vol. 10 Issue 3 (Ver 1.0) July 2010 Global Journal of Researches in Engineering Modeling Of Earthquake Ground Motion for Tbilisi Region P. Rekvava 1 k. Mdivani 2 GJRE Classification E (FOR) 090702,
More informationIntroduction to Displacement Modeling
Introduction to Displacement Modeling Introduction Deformation on the Earth surface informs us about processes and material properties below surface Observation tools: GPS (static-dynamic) InSAR (static)
More informationME 475 Modal Analysis of a Tapered Beam
ME 475 Modal Analysis of a Tapered Beam Objectives: 1. To find the natural frequencies and mode shapes of a tapered beam using FEA.. To compare the FE solution to analytical solutions of the vibratory
More informationThe 2016, October 26, Central Italy Earthquake Origin Time 17:10:36 UTC, M L(ISNet) =5.3; M W(ISNet) =5.6
The 2016, October 26, Central Italy Earthquake Origin Time 17:10:36 UTC, M L(ISNet) =5.3; M W(ISNet) =5.6 RISSC-Lab: Laboratorio di RIcerca in Sismologia Sperimentale e Computazionale The event as seen
More informationProgram System for Machine Dynamics. Abstract. Version 5.0 November 2017
Program System for Machine Dynamics Abstract Version 5.0 November 2017 Ingenieur-Büro Klement Lerchenweg 2 D 65428 Rüsselsheim Phone +49/6142/55951 hd.klement@t-online.de What is MADYN? The program system
More informationWind tunnel sectional tests for the identification of flutter derivatives and vortex shedding in long span bridges
Fluid Structure Interaction VII 51 Wind tunnel sectional tests for the identification of flutter derivatives and vortex shedding in long span bridges J. Á. Jurado, R. Sánchez & S. Hernández School of Civil
More informationEarthquake patterns in the Flinders Ranges - Temporary network , preliminary results
Earthquake patterns in the Flinders Ranges - Temporary network 2003-2006, preliminary results Objectives David Love 1, Phil Cummins 2, Natalie Balfour 3 1 Primary Industries and Resources South Australia
More informationMicrotremor survey methods in the Tamar Valley, Launceston, Tasmania: Evidence of 2D resonance from microtremor observations.
Microtremor survey methods in the Tamar Valley, Launceston, Tasmania: Evidence of 2D resonance from microtremor observations. Maxime Claprood 1, Michael W. Asten 2 1. Corresponding author. PhD candidate,
More informationEffects of Surface Geology on Seismic Motion
4 th IASPEI / IAEE International Symposium: Effects of Surface Geology on Seismic Motion August 23 26, 2011 University of California Santa Barbara TUNING THE DEEP VELOCITY STRUCTURE MODEL OF THE TOKYO
More informationStrong Ground Motion Prediction of Future Large Earthquake from Niavaran Fault in Tehran, Iran by Finite Fault Method
Strong Ground Motion Prediction of Future Large Earthquake from Niavaran Fault in Tehran, Iran by Finite Fault Method M. Samaei & M. Miyajima Kanazawa University, Japan M. Tsurugi Geo-Research Institute,
More informationVIBRATION PROBLEMS IN ENGINEERING
VIBRATION PROBLEMS IN ENGINEERING FIFTH EDITION W. WEAVER, JR. Professor Emeritus of Structural Engineering The Late S. P. TIMOSHENKO Professor Emeritus of Engineering Mechanics The Late D. H. YOUNG Professor
More informationof other regional earthquakes (e.g. Zoback and Zoback, 1980). I also want to find out
4. Focal Mechanism Solutions A way to investigate source properties of the 2001 sequence is to attempt finding well-constrained focal mechanism solutions to determine if they are consistent with those
More informationNONLINEAR STRUCTURAL DYNAMICS USING FE METHODS
NONLINEAR STRUCTURAL DYNAMICS USING FE METHODS Nonlinear Structural Dynamics Using FE Methods emphasizes fundamental mechanics principles and outlines a modern approach to understanding structural dynamics.
More informationSeismic Source Mechanism
Seismic Source Mechanism Yuji Yagi (University of Tsukuba) Earthquake Earthquake is a term used to describe both failure process along a fault zone, and the resulting ground shaking and radiated seismic
More informationTopography on Earthquake Motions in Sedimentary Basins
TRANSPORTATION RESEARCH RECORD 1411 Effects of Three-Dimensional Bedrock Topography on Earthquake Motions in Sedimentary Basins ARTHUR FRANKEL Work being done at the U.S. Geological Survey on 3-D simulations
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 informationAverage response spectra from some Australian earthquakes
Average response spectra from some Australian earthquakes T.I. Allen 1, G. Gibson 1,2 & J.P. Cull 1 1 School of Geosciences, Monash University, Melbourne, Australia 2 Seismology Research Centre, Environmental
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 informationMain Menu. Summary. Introduction
Kyosuke Okamoto *, JSPS Research Fellow, Kyoto University; Ru-shan Wu, University of California, Santa Cruz; Hitoshi Mikada, Tada-nori Goto, Junichi Takekawa, Kyoto University Summary Coda-Q is a stochastic
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 informationKinematic Waveform Inversion Study of Regional Earthquakes in Southwest Iberia
Kinematic Waveform Inversion Study of Regional Earthquakes in Southwest Iberia Ana Domingues Under supervision of Prof. João Fonseca and Dr. Susana Custódio Dep. Physics, IST, Lisbon, Portugal November
More informationEarthquakes. Forces Within Eartth. Faults form when the forces acting on rock exceed the rock s strength.
Earthquakes Vocabulary: Stress Strain Elastic Deformation Plastic Deformation Fault Seismic Wave Primary Wave Secondary Wave Focus Epicenter Define stress and strain as they apply to rocks. Distinguish
More informationPseudo-dynamic tests in centrifugal field for structure-foundation-soil systems
Pseudo-dynamic tests in centrifugal field for structure-foundation-soil systems Yoshikazu Takahashi Kyoto University, Kyoto, Japan Masako Kodera Tokyo Electric Power Company, Kawasaki, Japan SUMMARY: The
More informationSIMULATION OF NEAR FIELD STRONG GROUND MOTIONS AT TOMBAK SITE IN IRAN USING HYBRID METHOD
SIMULATION OF NEAR FIELD STRONG GROUND MOTIONS AT TOMBAK SITE IN IRAN USING HYBRID METHOD ABSTRACT : B. Ebrahimian 1, M.S. Sahraeian 2 and As. Noorzad 3 1 Ph.D. Candidate, School of Civil Engineering,
More informationStochastic Dynamics of SDOF Systems (cont.).
Outline of Stochastic Dynamics of SDOF Systems (cont.). Weakly Stationary Response Processes. Equivalent White Noise Approximations. Gaussian Response Processes as Conditional Normal Distributions. Stochastic
More informationPreliminary Examination - Dynamics
Name: University of California, Berkeley Fall Semester, 2018 Problem 1 (30% weight) Preliminary Examination - Dynamics An undamped SDOF system with mass m and stiffness k is initially at rest and is then
More informationDistribution Restriction Statement Approved for public release; distribution is unlimited.
CECW-ET Engineer Manual 1110-2-6050 Department of the Army U.S. Army Corps of Engineers Washington, DC 20314-1000 EM 1110-2-6050 30 June 1999 Engineering and Design RESPONSE SPECTRA AND SEISMIC ANALYSIS
More informationEffective stress analysis of pile foundations in liquefiable soil
Effective stress analysis of pile foundations in liquefiable soil H. J. Bowen, M. Cubrinovski University of Canterbury, Christchurch, New Zealand. M. E. Jacka Tonkin and Taylor Ltd., Christchurch, New
More informationACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT
ACCOUNTING FOR SITE EFFECTS IN PROBABILISTIC SEISMIC HAZARD ANALYSIS: OVERVIEW OF THE SCEC PHASE III REPORT Edward H FIELD 1 And SCEC PHASE III WORKING GROUP 2 SUMMARY Probabilistic seismic hazard analysis
More informationInclusion of a Sacrificial Fuse to Limit Peak Base-Shear Forces During Extreme Seismic Events in Structures with Viscous Damping
Inclusion of a Sacrificial Fuse to Limit Peak Base-Shear Forces During Extreme Seismic Events in Structures with Viscous Damping V. Simon, C. Labise, G.W. Rodgers, J.G. Chase & G.A. MacRae Dept. of Civil
More informationlog 4 0.7m log m Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Module 1 Seismology Exercise Problems :
Seismic Analysis of Structures by TK Dutta, Civil Department, IIT Delhi, New Delhi. Module Seismology Exercise Problems :.4. Estimate the probabilities of surface rupture length, rupture area and maximum
More informationFault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction
Fault Specific, Dynamic Rupture Scenarios for Strong Ground Motion Prediction H. Sekiguchi Disaster Prevention Research Institute, Kyoto University, Japan Blank Line 9 pt Y. Kase Active Fault and Earthquake
More information