首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
In this work we studied the performance of different numerical approaches to simulate the large amplifications of long period earthquake ground motion within the Gubbio plain, a closed-shape intra-mountain alluvial basin of extensional tectonic origin in Central Italy, observed during the Umbria-Marche 1997 seismic sequence. Particularly, referring to the Sep 26 1997 Mw6.0 mainshock, we considered the following numerical approximations: (a) 3D model, including a kinematic model of the extended seismic source, a layered crustal structure, and the basin itself with a simplified homogeneous velocity profile; (b) 2D model of a longitudinal and transversal cross-section of the basin, subject to vertical and oblique incidence of plane waves with time dependence at bedrock obtained by the 3D simulations; (c) 1D model. 3D and 2D numerical simulations were carried out using the spectral element code GeoELSE, exploiting in 3D its implementation in parallel computer architectures. 3D numerical simulations were successful to predict the observed large amplification of ground motion at periods beyond about 1 s, due to the prominent onset of surface waves originated at the southern edge of the basin and propagating northwards. More specifically, the difference of 3D vs 2D results is remarkable, since the latter ones fail to approach such large amplification levels, even when an oblique incidence of plane waves is considered.  相似文献   

2.
Azimuth dependent wave amplification in alluvial valleys   总被引:1,自引:0,他引:1  
An extension of the indirect boundary element method (IBEM) to three-dimensional scattering by two-dimensional alluvial valleys is presented. While the valley is two-dimensional, the incident plane waves can arrive outside the 2D plane so the scattering is three-dimensional with coupling of P---SV---SH waves. Such a method makes it possible to take earthquake location into account in the estimation of site effects in alluvial valleys. The method is validated by transparency tests, by comparison with 2D simulations, and by comparison with results of other authors. The advantage of the method is that is combines high accuracy with cost-efficiency in terms of computer-time. It is applied to theoretically estimate site effects across a simplified model of an alluvial valley in the French Alps where azimuth dependence of local amplification has been observed. A parametric study with simulations for a range of azimuths and incidence angles shows that (1) the local amplification depends strongly on both azimuth and incidence of the incoming waves, (2) the global pattern of amplification across the valley is very complex for all azimuths, and (3) it is not possible to predict the 3D response of the valley from 2D modeling. Theoretical spectral ratios are in approximate agreement with observed ones for a station in the center of the valley where the local structure justifies use of a simplified model for the comparison.  相似文献   

3.
基于动力子结构方法的场地地震反应分析方法   总被引:2,自引:1,他引:1  
本文验证了将约束子结构法引入到复杂场地地震反应分析的可行性,并在此基础上提出了两种进一步提高计算效率的简化措施。结合土层有限元模型的自身特点,提出了标准子结构的概念,通过设置标准子结构,可减少相同子结构的重复计算。基于约束子结构本身的性质,假定土层计算区域为局部非线性,可简化等效线性化分析过程,通过数值试验给出了局部非线性区域的取值范围。算例表明,上述简化措施在提高大规模复杂场地地震反应分析计算效率方面具有显著的优势。  相似文献   

4.
In soil‐structure interaction modeling of systems subjected to earthquake motions, it is classically assumed that the incoming wave field, produced by an earthquake, is unidimensional and vertically propagating. This work explores the validity of this assumption by performing earthquake soil‐structure interaction modeling, including explicit modeling of sources, seismic wave propagation, site, and structure. The domain reduction method is used to couple seismic (near‐field) simulations with local soil‐structure interaction response. The response of a generic nuclear power plant model computed using full earthquake soil‐structure interaction simulations is compared with the current state‐of‐the‐art method of deconvolving in depth the (simulated) free‐field motions, recorded at the site of interest, and assuming that the earthquake wave field is spatially unidimensional. Results show that the 1‐D wave‐field assumption does not hold in general. It is shown that the way in which full 3‐D analysis results differ from those which assume a 1‐D wave field is dependent on fault‐to‐site geometry and motion frequency content. It is argued that this is especially important for certain classes of soil‐structure systems of which nuclear power plants subjected to near‐field earthquakes are an example.  相似文献   

5.
Two hilly sites were selected to study seismic site response due to topography effects. The sites were selected in a manner to be as much as possible homogenous and free of the soft soil layers effects. The hills were instrumented by nine velocimetric stations to record microtremors and the obtained data were analyzed using horizontal to vertical spectral ratios. Some standard spectral ratio tests were performed on noise as well. Then the instrumented hills were modeled (both 2D and 3D) assuming a linear elastic constitutive behavior subjected to vertically propagating SV and P Ricker wavelets. All calculations were performed in time domain using direct boundary element method. Different transfer function components, amplification patterns and spectral ratios were calculated in frequency domain. The frequency of vibration, obtained by experimental studies, is between 4 and 5 Hz for both of the hills. The spectral ratios derived by numerical simulations were compared with the observed spectral ratios. They show relatively good similarities between the results of these two methods. The frequencies of vibration derived from different methods seem to be nearly identical. The agreement in term of resonance frequency between microtremors and numerical modeling suggests that noise measurements could represent a simple, even if preliminary, tool in order to identify possible topographic amplification.  相似文献   

6.
In the present work the seismic site response of Narni ridge (Central Italy) is evaluated by comparing experimental results and numerical simulations. The inhabited village of Narni is located in central Apennines at the top of a steep massive limestone ridge. From March to September 2009 the site was instrumented with 10 weak-motion stations, 3 of which located at the base of the ridge and 7 at the top. The velocimetric network recorded 642 events of ML up to 5.3 and hypocentral distance up to about 100 km. The great amount of data are related to the April 2009 L’Aquila sequence. The site response was analyzed using both reference (standard spectral ratio, SSR) and non reference spectral techniques (horizontal to vertical spectral ratio, HVSR). Moreover directional analyses were performed in order to evaluate the influence of the ridge orientation with respect to the selected source-site paths. In general the experimental results show amplification factors for frequencies between 4 and 5 Hz for almost all stations installed along the crest. The SSR technique provides amplification factors up to 4.5 in a direction perpendicular to the main elongation of the ridge. The results obtained from the data analyses were used as a target for bidimensional and tridimensional numerical simulations, performed using a hybrid finite-boundary element method and a boundary element method for 2D and 3D modelling, respectively. In general, the results obtained through numerical simulation fit well the experimental data in terms of range of amplified frequencies, but they underestimate by a factor of about 2 the observed amplifications.  相似文献   

7.
Observations from many recent strong motion events have shown the importance of local soil conditions and non-linear soil behaviour on the seismic ground response (site effects). As demonstrated by previous seismic microzoning studies (Lebrun et al.) [1]), as well as by at least three historical major earthquakes, Pointe-à-Pitre is prone to strong site effects, due to the particular geology of the area. In this paper, we present a comparison between the strong-motion data available from the stations operating on the swampy site of Pointe-à-Pitre airport and the ground motions derived from 1D non-linear finite element simulations.Results show that, for moderate to strong ground motions, 1D simulations reproduce the main characteristics of site response in terms of duration, energy distribution, amplitude and frequency content. It also shows the importance of very superficial soft layers as peat or saturated mud in low frequency site effects simulations. This point is important for further engineering studies since such very soft formations overlain by stiffer landfills are commonly expected in the Antilles context. Our work also shows that Anderson's criteria, used to quantify the goodness-of-fit of simulated ground motions to the observed ones, appear to be an interesting diagnostic tool for testing the quality of numerical simulations from an engineering point of view.  相似文献   

8.
This study assesses the 3D amplification effects in shallow basins and quantifies the effects of site‐city interaction (SCI) on high‐rise buildings. A regional‐scale 3D spectral element simulation is conducted on the Tuen Mun‐Yuen Long basin, which contains multiple subbasins with heterogeneous and nonlinear soil profiles, while 3D city models with various building layouts are fully integrated into the basin model for our SCI study. We found a good correlation between spectral amplification factors and soil depths. Site response is significantly amplified at basin edges and centers due to surface waves generated at basin edges and the focusing effects stemming from 3D basin geometry. Transfer functions of 3D basins can be up to fourfold at fundamental frequencies as compared to 1D response, and further amplifications occur at high frequencies due to surface waves. In the SCI simulations, we observe wave trapping in the open space amid buildings resulting in energy concentration and up to twofold PGA amplifications. The wave trapping effect diminishes as the space between buildings increase beyond their range of influence (~100 m). The SCI analyses show that destructive kinetic energy in superstructures increases 28% in one horizontal direction but decreases 22% in the other. Our study concluded that, 1D site response analysis can significantly underestimate the seismic demand in shallow basins. Site‐city interaction of high‐rise buildings increases the short‐period spectra of ground motions, leading to an increase in their story accelerations by up to 50% and to a substantial decrease in the seismic safety of short structures in their vicinity.  相似文献   

9.
This paper presents a 2D model of the geological structure of Thessaloniki city and results of empirical and theoretical approaches for the evaluation of site response due to complex site effects. The construction of the 2D model is based on the available geophysical and geotechnical data in terms of the most important parameters needed to model site response. The well-known subsoil structure, despite the existence of some local uncertainties, gave the possibility to correlate the geometry and the dynamic properties of the 2D model with the results of site response determined from the analysis of one event in frequency and time domains and 1D numerical modelling. The study of site response shows the effect of the lateral variations on ground motion and suggests that the contribution of locally generated surface waves to the resonant peak may be important. In this case history, the limitations of the 1D approximation to simulate ground motion under complex soil conditions in both frequency and time domains are also shown. This paper lays the ground for a companion article dealing with 2D site effects.  相似文献   

10.
The transient response of large embedded foundation elements of length-to-diameter aspect ratio D/B=2–6 is characterized by a complex stress distribution at the pier–soil interface that cannot be adequately represented by means of existing models for shallow foundations or flexible piles. On the other hand, while three-dimensional (3D) numerical solutions are feasible, they are infrequently employed in practice due to their associated cost and effort. Prompted by the scarcity of simplified models for design in current practice, we here develop an analytical model that accounts for the multitude of soil resistance mechanisms mobilized at their base and circumference, while retaining the advantages of simplified methodologies for the design of non-critical facilities. The characteristics of soil resistance mechanisms and corresponding complex spring functions are developed on the basis of finite element simulations, by equating the stiffness matrix terms and/or overall numerically computed response to the analytical expressions derived by means of the proposed Winkler model. Sensitivity analyses are performed for the optimization of the truncated numerical domain size, the optimal finite element size and the far-field dynamic boundary conditions to avoid spurious wave reflections. Numerical simulations of the transient system response to vertically propagating shear waves are next successfully compared to the analytically predicted response. Finally, the applicability of the method is assessed for soil profiles with depth-varying properties. The formulation of frequency-dependent complex spring functions including material damping is also described, while extension of the methodology to account for nonlinear soil behavior and soil–foundation interface separation is described in the conclusion and is being currently investigated.  相似文献   

11.
12.
This paper presents results of numerical modelling of site response for Thessaloniki, obtained with two different 2D methods; a finite difference and a finite element method. Ground motion across a 2D model of the subsoil of the city has been simulated for vertically incident SH waves. The predominance of locally generated surface waves is very clear in the synthetic seismograms of a weak event and of stronger ones. These results are then compared with the observations in time domain and frequency domain. The role of the soil formations with high attenuation in the lateral propagation and the effect of the differential motion close to the lateral variations are also pinpointed. The stronger events were finally used to compute strong ground motion in order to reveal and to discuss practical engineering aspects such as peak ground acceleration value, the most familiar indicator in seismic norms, the soil to rock spectral coefficients for the period bandwidth of interest, and the aggravation factor in terms of 2D to 1D response spectra as a useful ruler to account for complex site effects.  相似文献   

13.
The paper discusses the seismic response of circular tunnels in dry sand and investigates the efficiency of current seismic analysis methods at extreme lining flexibilities. Initially, a dynamic centrifuge test on a flexible circular model tunnel, embedded in dry sand, is analyzed by means of rigorous full dynamic analysis of the coupled soil–tunnel system, applying various non-linear soil and soil–tunnel interface models. The numerical results are compared to the experimental ones, aiming to better understand the recorded response and calibrate the numerical models. Then a series of numerical analyses are conducted using the validated numerical model, in order to investigate the effect of the tunnel lining rigidity on the dynamic response of the soil–tunnel system. In parallel, the accuracy of currently used simplified analysis methods is evaluated, by comparing their predictions with the results of the a priori more accurate and well validated numerical models. The comparative analyses allow us to highlight and discuss several crucial aspects of the soil-tunnel system seismic response, including (1) the post-earthquake residual values of the lining forces, which are amplified with the increase of the flexibility of the tunnel and (2) the importance of the soil-tunnel interface conditions. It is finally concluded that simplified analysis methods may provide a reasonable framework for the analysis at a preliminary stage, under certain conditions.  相似文献   

14.
In recent years, because of a seismic gap of 136 years in the north of Chile, several studies have been carried out with the purpose of characterizing the dynamic properties of the soils of the most populated cities in this area. The purpose of this work is to estimate, by numerical modeling, the seismic amplification of a densely populated area of the city of Arica. Spatial distribution of main soil dynamic properties has been obtained from a detailed geophysical survey, including surface-wave based methods and gravimetry. To estimate the site effects, we solve the wave propagation equation in full heterogeneous media by the spectral element method. This method allows including the topography, the irregular contact between soils and bedrock and heterogeneities of main materials’ properties along the computational domain. The inelastic behavior of the soil has been considered, using stiffness degradation and damping curves. The computations were done using the high performance open-source numerical code SPEED. Results from the distribution of the peak ground acceleration are compared against standard 1D horizontally layered modeling using equivalent linear and nonlinear viscoelastic approximations of nonlinear soil behavior. Linear viscoelastic results indicate that the north area of the city of Arica is not especially prone to tridimensional site effects induced by soil’s natural variability and the apparent bedrock geometry. Nevertheless, non-conservative differences of about 20% between 1D v/s 3D analyses were found under inelastic material behavior. In general, the two 3D approaches followed to include soils’ nonlinearity provide similar site amplification estimates.  相似文献   

15.
El Infiernillo, a 145 m high rockfill dam in Mexico built in a narrow V-shaped canyon, was subjected to eight major earthquakes since its construction. In this study, the dynamic dam response is analysed using (i) the recorded November 15, 1975 bedrock-crest acceleration histories and (ii) the results of a 1970 full-scale test conducted by UNAM, in which eight upstream–downstream and longitudinal resonant frequencies and configurations were observed and documented. These observed and seismically induced dynamic responses are compared herein to predictions of two different numerical models of El Infiernillo dam; a newly developed simplified three-dimensional (3D) model, and a 3D finite element model. The dynamic dam response characteristics are assessed, and performance of the employed numerical models is evaluated. It is found that (i) higher modes of vibration had participated significantly in the recorded seismic response and (ii) upstream–downstream response is well represented by the two numerical models employed. Using the simplified model, the September 19, 1985 earthquake non-linear response is computed and shown to compare satisfactorily with its recorded counterpart in the upstream-downstream and vertical directions. The largest computed dynamic accelerations, stresses and strains are found to occur within the upper third of the dam body.  相似文献   

16.
场地条件对地震动特性影响显著,在抗震设计反应谱的确定过程中,需根据场地条件对加速度反应谱予以相应的调整。已有场地条件影响调整方案研究成果,均基于数值模拟或局部地区强震动记录统计,多数仅给出了峰值加速度PGA场地条件影响调整系数,对非线性的考虑缺乏观测数据依据。为此在全球强震动记录统计获得的PGA归一化加速度反应谱和日本钻井台阵记录获得的加速度反应谱平台值非线性衰减指数的基础上,结合钻孔模型数值模拟和近期研究成果,建立了考虑场地条件影响非线性的地震动加速度反应谱场地条件影响调整系数方案。  相似文献   

17.
We study site effects using 520 weak motion earthquake records from a vertical array in Aegion, Greece. The array is inside a basin, has four stations in soil, and one in bedrock (178 m depth). The site is marked by high seismicity and complex surface geology. We first use the records to establish the downhole accelerometer orientations and their evolution with time. Then we estimate site effects using empirical spectral ratios with and without a reference site (standard and horizontal-to-vertical spectral ratio). We find significant site amplification which cannot be accounted for by 1D model predictions, along with a significant difference in the amplification level between the two horizontal components. These are indications of 2D effects, namely surface waves generated at the basin edge. The difference in amplification between the horizontal components is maximised when these are rotated with respect to the orientation of the basin edge. The strongest amplification takes place in the direction parallel to the basin edge (SH, or out-of-plane motion), and is up to 2 times higher than in the perpendicular direction (SV, or in-plane motion). This directional effect on the amplification is corroborated by numerical 2D modelling using incident SH and SV waves, with the former possibly generating strong Love waves. In the records, the directionality is clear for windows containing the largest amplitudes of the records (S waves and strong surface waves), while it tends to vanish for coda-wave windows. This directionality is also observed when using response spectral ratios rather than Fourier ratios. We compute soil-to-rock amplification factors for peak ground acceleration (PGA) and find it is significantly higher than what is predicted by current design codes. We attribute this difference to the basin edge amplification, linear soil behaviour, and to the inability of simple scalar values like PGA to describe complex amplification effects. Finally, we analyse the earthquake records at a surface station near the slope crest and do not observe significant topographic amplification.  相似文献   

18.
The problem of obtaining field‐scale surface response to rainfall events is complicated by the spatial variability of infiltration characteristics of the soil and rainfall. In this paper, we develop and test a simplified model for generating surface runoff over fields with spatial variation in both rainfall rate and saturated hydraulic conductivities. The model is able to represent the effects of local variation in infiltration, as well as the run‐on effect that controls infiltration of excess water from saturated upstream areas. The effective rainfall excess is routed to the slope outlet using a simplified solution of the kinematic wave approximation. Model results are compared to averaged hydrographs from numerically‐intensive Monte–Carlo simulations for observed and design rainfall events and soil patterns that are typical of Central Italy. The simplified model is found to yield satisfactory results at a relatively small computational expense. A proposal to include a simple channel routing scheme is also presented as a prelude to extend this conceptualization to watershed scales. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

19.
It is well-known that the response of a site to a seismic solicitation depends on local topographical and geotechnical characteristics. Many aspects of seismic site effect still need to be studied in more detail and they can be incorporated in the seismic norms after quantification. The purpose of this paper is to contribute to establishment of a simple method to include complex site effects in a building code. Horizontal ground movements in various points of two-dimensional (2D) irregular configurations subjected to synthetic SV waves of vertical incidence are calculated. The parametric studies are achieved by means of HYBRID program combining finite elements in the near field and boundary elements in the far field (FEM/BEM). The results are shown in the form of pseudo-acceleration response spectra. For the empty valleys, we can classify the spectral response according to a unique geometric criterion: the “surface/angle” ratio, where surface is the area of the valley opening, and angle denotes the angle between the slope and horizontal line in the above corner. To assess the influence of the 2D effect on the spectral response of filled valleys, the response of alluvial basins are compared with the response of one-dimensional columns of soil. Finally, an offset criterion is proposed to choose a relevant computation method for the spectral acceleration at the surface of alluvial basins.  相似文献   

20.
A detailed numerical simulation of the ground motion and a site response analysis for two towns in the Marche Region (Treia and Cagli) is carried out on the basis of structural models deduced from available geological and geophysical data. In both cases, the reference event is an M = 5.7 earthquake associated with a normal fault located beneath each town. The ground motion is computed using the 2D spectral element method (SPEM 2D). The method solves the propagation of the seismic field through complex geological structures and enables an estimate of the effects of deep crustal structure, superficial geology, and topography on ground motion. Numerical simulations of the seismic field are performed along 2D vertical planes containing the seismic source. Strong ground motion has not been yet recorded in the two towns; therefore, the numerical simulation of ground motion represents a way to overcome the lack of instrumental data. The simulations carried out for Treia show that ground motion is influenced by both source mechanism and effects due to propagation through the geological structure, while ground motion in Cagli features strong local effects, caused by the presence of alluvial deposits under a large area of the town.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号