بررسی اثرات توپوگرافی در ساختگاه سد پاکویما با روش اجزای مرزی سه‌بعدی در حوزه زمان

نوع مقاله : Articles

نویسندگان

1 دانشکده زمین‌شناسی، پردیس علوم، دانشگاه تهران، تهران، ایران

2 دانشکده مهندسی عمران، دانشگاه تبریز، تبریز، ایران

3 دانشکده زمین شناسی، پردیس علوم، دانشگاه تهران، تهران، ایران

چکیده

توپوگرافی ناهموار یکی از عوامل ایجاد تفاوت جنبش در نقاط تکیه‌گاهی سازه‌های بزرگ نظیر سدها می‌باشد. در این مقاله از روش عددی اجزای مرزی سه‌بعدی جهت حل پاسخ لرزه‌ای دره محل سد پاکویما استفاده شده و نتایج حاصله در کنار تحلیل‌های ارائه شده، با مقادیر ثبت شده طی زمین‌لرزه سال 2001 سد پاکویما مقایسه شده است. از مؤلفه‌های افقی زمین‌لرزه واقعی ثبت شده در کف دره پس از پردازش‌های لازم، به‌عنوان موج ورودی جهت تحلیل استفاده گردیده و تاریخچه تغییر مکان و بزرگنمایی نقاط مختلف دره محل سد ارائه و بررسی شده است. نتایج نشانگر آن است که دامنه بیشینه تغییر مکان با افزایش ارتفاع نقاط در دو سوی دره کاهش می‌یابد؛ با این حال در حوزه فرکانس، بزرگنمایی در فرکانس‌های بالا با افزایش ارتفاع نقاط دره افزایش می‌یابد. نسبت دامنه تغییر مکان‌ها در طول دره محل سد تا حدود 3 نیز می‌رسد که این مسئله لزوم توجه به جنبش‌های غیریکنواخت در سازه‌های بزرگ را نشان می‌دهد.

کلیدواژه‌ها


عنوان مقاله [English]

Topography Effects in Pacoima Dam Site Using Time-Domain Three-Dimensional BEM

نویسندگان [English]

  • Abdollah Sohrabi-Bidar 1
  • Mohsen Isari 2
  • Reza Tarinejad 2
  • Shahram Maghami 3
1 School of Geology, College of Science, University of Tehran, Tehran, Iran
2 Faculty of Civil Engineering, University of Tabriz, Tabriz, Iran
3 School of Geology, College of Science, University of Tehran, Tehran, Iran
چکیده [English]

The reliable estimation of seismic loads on a structure is required in order to earthquakes resistant design of the structure. The difference in seismic loading in different support points of the structure is important in large and long structures. In general, the lack of access to the reliable time histories in different support points of the structure is the main problem of performing non-uniform excitation analyses. Numerical analyses and calculation of ground motion at different points of the foundation of the structure is one of the ways to achieve the non-uniform support excitation. This paper aims to evaluate the seismic response of Pacoima dam site by performing three-dimensional boundary element analyses in the time domain. The pattern of displacement and amplification due to seismic waves scattering in the dam site are evaluated, and calculated results are compared to the recorded ground motions.
The numerical modeling has been executed using the time-domain boundary element that is based on the boundary integral equation of the wave motion. To transform the governing integral equation into the ideal form, it has been discretized in both time and spatial domains. Finally, the obtained equations have been expressed in the matrix form and have been implemented in a computer code named as BEMSA. Earlier, several different examples of wave scattering have been solved in order to show the accuracy and efficiency of the implemented BE algorithm in carrying out the site response analysis of topographic structures.
Pacoima dam is a concrete arch dam located in the San Gabriel Mountains in Los Angeles County. The height and the length of the crest of the dam are 113 m and 180 m, respectively. The dam is instrumented by use of 17 accelerometers at different elevations on the dam body and its abutments. For site response analyses, the dam site has been subjected to vertically propagating recorded motions of the Pacoima dam 2001 earthquake with a magnitude of 4.3, depth of about 9 and epicentral distance of about 6 km south of the dam. The medium assumed to be homogeneous linear elastic with density of 2.64 ton/m3, shear wave velocity of 2000 m/s and Poisson's ratio of 0.25. The 3D topographic model has been generated up to a radius of 5000 meters, using 1218 eight-node quadrilateral isoparametric elements with the average effective element size of 25 m in the center part of the model.
In order to investigate the seismic response of the canyon, a couple of points at four levels have been considered on both sides of the canyon and the results analyzed in time and the frequency domains. Despite the actual record earthquake motions, which includes the effects of the interaction between the foundation and the dam structure as well as the lake behind the dam, the calculated motions include only the wave scattering by the topography of the canyon. Therefore, although the exact matching of the recorded and calculated motions are not expected, comparison of the motions show that the patterns of the displacements are close together. This phenomenon indicates the importance of valley shape and its important contribution to the dynamic response of the dam site. Assessment of the displacement time histories in various points at both sides of the canyon indicates that the amplitude of the motions decreases when the height of the point increase. Besides, the comparison between the motions of the left and right sides points show have a higher frequency content and a higher shear-wave velocity. 4) In all ten soil groups, the shear wave velocity
that due to the non-symmetricity of the canyon, displacement amplitudes in the left side are larger than the right side. Based on the calculated displacements on the various points, the maximum amplitude along the canyon would be changed up to three times.
In the frequency domain, different points of the canyon surface have generally the similar amplification patterns. There are two main peaks of amplification in the frequency range of 3-5 Hz and the frequency range of 6-8 Hz. In both sides of the canyon by increasing the height of the points amplification is increased, especially in the frequency range of 6-8 Hz. Moreover, at the same elevation points, the amplification value in the left side is higher than the right side. Comparison of amplification curves of recorded and calculated motions, show the appearance of new peaks of amplification in higher frequency, which could be related to the real conditions of the dam site.
Finally, although the motion amplitude in time domain decreases by height increasing on both sides of the canyon, as expected, the amplification in the frequency domain, especially in high frequencies, increases. This insists that the amplification characteristics of a site should be considered and interpreted as a frequency dependence phenomenon. Moreover, the results indicate the spatial variation of the motion due to the topography effect along the canyon, in which the amplitude of peak ground displacements along the canyon has been changed up to three times.

کلیدواژه‌ها [English]

  • Topography Effects
  • Boundary Elements
  • Time-domain
  • Numerical modelling
  • Pacoima Dam
  1. Hariri-Ardebili, M.A. and Mirzabozorg, H. (2010) Nonlinear Seismic Performance Evaluation of Concrete Arch Dams Using Endurance Time Method. M.Sc. Thesis, Civil Engineering, K.N. Toosi University of Technology, Tehran, Iran.
  2. Tarinejad, R., Fatehi, R. and Harichandarn, R.S. (2013) Response of an arch dam to non-uniform excitation generated by a seismic wave scattering model. Soil Dynamics and Earthquake Engineering, 52, 40-54.
  3. Friedman, M.B. and Shaw, R.P. (1962) Diffraction of pulses by cylindrical obstacles of arbitrary cross section. Journal of Applied Mechanics, 29, 40-46.
  4. Niwa, Y., Fukui, T., Kato, S., Fujiki, K. (1980) An application of the integral equation method to two-dimensional elastodynamics. Theory of Applied Mechanics, 28, 281-290.
  5. Mansur, W.J. (1983) A Time-Stepping Technique to Solve Wave Propagation Problems Using the Boundary Element Method. Ph.D. Dissertation, Southampton University.
  6. Antes, H. (1985) A boundary element procedure for transient wave propagation in twodimensional isotropic elastic media. Finite Elemens in Analysis and Design, 1, 313-322.
  7. Mansur, W.J. and Brebbia, C.A. (1985) Transient elastodynamics. Topics in Boundary Element Research. CA Brebbia, ed., Vol. 2: Time-dependent and Vibration Problems, Chap 5, pp. 124-155.
  8. Karabalis, D.L. and Beskos, D.E. (1984) Dynamic response of 3-D rigid surface foundations by time domain boundary element method. Earth. Eng. and Struc. Dyn., 12, 73-93.
  9. Manolis, G.D., Ahmad, S., Banerjee, P.K. (1985) ‘Boundary element method implementation for three-dimensional elasto-dynamics’. In: Developments in Boundary Element Methods: IV. P.K. Banerjee, J.O. Watson, eds. Elsevier Applied Science Publishers: London, 29-63.
  10. Zhao, C., Valliappan, S. and Wang, Y.C. (1992) A numerical model for wave scattering problems in infinite media due to P-and SV-Wave incidences. International Journal for Numerical methods in Engineering, 33, 1661-1682.
  11. Huang, H.C. and Chiu, H.C. (1995) The effect of canyon topography on strong ground motion at Feitsui Damsite: Quantitive Results. Earth. Eng. and Struc. Dyn., 24, 977-990.
  12. Paolucci, R. (2002) Amplification of earthquake ground motion by steep topographic irregularities. Earth. Eng. and Struc. Dyn., 31, 1831-1853.
  13. Álvarez-Rubio, S., Jose Benito, J Sanchez-Sesma, F.J., and Alarcon, E. (2005) The use of direct boundary element method for gaining insight into complex seismic site response. Computers and Structures, 83, 821-835.
  14. Kamalian, M., Jafari, M.K., Sohrabi Bidar, A., Razmkhah, A. and Gatmiri, B. (2006) time-domain two-dimensional site response analysis of non-homogeneous topographic structures by a hybrid FE/BE method. Soil Dynamic and Earthquake Engineering, 26, 753-765.
  15. Kamalian, M., Gatmiri, B., Sohrabi-Bidar, A. and Khalaj, A. (2007) Ampliï‌cation pattern of 2D semi-sine-shaped valleys subjected to vertically propagating incident waves. Communications in Numerical Methods in Engineering, 23, 871-887.
  16. Tarinejad, R., Ahmadi, M.T., and Khaji, N. (2007) Analysis of Topographic Amplification Effects on Canyon Sites using 3D Boundary Element Method. Journal of Seismology and Earthquake Engineering, 9, 25-37.
  17. Gatmiri, B., Arsonb, C., Nguyen, K.V. (2008) Seismic site effects by an optimized 2D BE/FE method I. Theory, numerical optimization and application to topographical irregularities. Soil Dynamics and Earthquake Engineering, 28, 632-645.
  18. Gatmiri, B., Maghoul, P., Arson, C. (2009) Site-speciï‌c spectral response of seismic movement due to geometrical and geotechnical characteristics of sites. Soil Dynamics and Earthquake Engineering, 29, 51-70.
  19. Sohrabi-Bidar, A., Kamalian, M. and Jafari, M.K. (2010) Seismic response of 3-D Gaussian-shaped valleys to vertically propagating incident waves. Geophys. J. Int., 183, 1429-1442.
  20. Sohrabi-Bidar, A. and Kamalian, M. (2013) Effects of three-dimensionality on seismic response of Gaussian-shaped hills for simple incident pulses. Soil Dynamics and Earthquake Engineering, 52, 1-12.
  21. Panji, M., Kamalian, M., Asgari Marnani, J. and Jafari, M.K. (2013) Transient analysis of wave propagation problems by half-plane BEM. Geophys. J. Int., 194, 1849-1865.
  22. Panji, M., Kamalian, M., Asgari Marnani, J., and Jafari, M.K. (2014) Analyzing seismic convex topographies by a half-plane time-domain BEM. Geophys. J. Int., 197, 591-607.
  23. Sohrabi-Bidar, A. (2008) Seismic Behavior Assessment of Surface Topographies Using Time Domain 3D Boundary Elements Method. Ph.D. Dissertation Geophysics-Seismology, International Institute of Earthquake Engineering and Seismology.
  24. Brebbia, C.A., Dominguez, J. (1989) Boundary Elements, an Introductory Course. Computational Mechanics Publications: Boston.
  25. Alves, S.W. (2004) Nonlinear Analysis of Pacoima Dam with Spatially Non-Uniform Ground Motion. In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy. California Institute of Technology Pasadena, California.