共查询到20条相似文献，搜索用时 93 毫秒
1.
This paper presents a new analysis method, called macro–micro analysis method (MMAM) for numerical simulation of wave propagation in stochastic media, which could be used to predict distribution of earthquake strong motion with high accuracy and spatial resolution. This MMAM takes advantage of the bounding medium theory (BMT) and the singular perturbation expansion (SPE). BMT can resolve uncertainty of soil and crust structures by obtaining optimistic and pessimistic estimates of expected strong motion distribution. SPE leads to efficient multi‐scale analysis for reducing a huge amount of computation. The MMAM solution is given as the sum of waves of low resolution covering a whole city and waves of high resolution for each city portion. This paper presents BMT and SPE along with the formulation of MMAM for wave propagation in three‐dimensional elastic media. Application examples are presented to verify the validity of the MMAM and demonstrate potential usefulness of this approach. In a companion paper (Earthquake Engng. Struct. Dyn., this issue) application examples of earthquake strong motion prediction are also presented. Copyright © 2005 John Wiley & Sons, Ltd. 相似文献
2.
《地震工程与结构动力学》2018,47(5):1309-1328
A direct finite element method for nonlinear earthquake analysis of 2‐dimensional dam–water–foundation rock systems has recently been presented. The analysis procedure uses standard viscous‐damper absorbing boundaries to model the semi‐unbounded foundation‐rock and fluid domains and specifies the seismic input as effective earthquake forces at these boundaries. Presented in this paper is a generalization of the direct finite element method with viscous‐damper boundaries to 3‐dimensional dam–water–foundation rock systems. Step‐by‐step procedures for determining the effective earthquake forces starting from a ground motion specified at a control point on the foundation‐rock surface is developed, and several numerical examples are computed and compared with independent benchmark solutions to demonstrate the effectiveness of the analysis procedure for modeling 3‐dimensional systems. 相似文献
3.
T. Ichimura M. Hori P. E. Quinay M. L. L. Wijerathne T. Suzuki S. Noguchi 《地震工程与结构动力学》2012,41(4):795-811
The seismic structural response is affected by temporal and spatial variations in strong ground motion. It can be evaluated through the fault‐structure system: the fault mechanism, wave propagation through the crust, amplification near the surface, and soil‐structure interaction. To analyze this system at high resolution and accuracy, we previously proposed a new multiscale analysis method and numerically verified its validity. However, the problem of the extremely large computation cost of constructing a three‐dimensional numerical model and solving the discretized governing equations still remains. Here, we introduce a new method to resolve these difficulties. By combining this new method with our multiscale analysis, we developed a tool for fault‐structure system analysis. The accuracy of this tool is verified by comparing it to a Green's function solution. Finally, we demonstrate the potential utility of the method by estimating the seismic response of a large and complex underground highway junction in a given earthquake scenario. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
4.
Current state of integrated earthquake simulation for earthquake hazard and disaster 总被引：1，自引：0，他引：1
This paper presents the current state of integrated simulation for earthquake hazard and disaster. This simulation takes advantage
of the macro–micro analysis method; this method estimates an earthquake’s strong motion with high spatial and temporal resolution,
using the bounding medium theory to obtain optimistic and pessimistic estimates of expected strong motion distribution and
the singular perturbation expansion that results in an efficient multi-scale analysis. Integrated earthquake simulation calculates
seismic responses for all structures in a target area, inputting simulated strong ground motion to a structure analysis method
that is plugged into the system by means of a wrapper; a suitable method, linear or nonlinear, is chosen depending on the
type of the structure. The results of all simulations are visualized so that residences and government officials can share
a common recognition of earthquake hazard and disaster. Two examples of this integrated earthquake simulations are presented;
one is made by plugging nonlinear structure analysis methods into the system, and the other is made for an actual city, the
computer model of which is constructed with the help of available geographical information systems. 相似文献
5.
In numerical simulation of wave scattering under oblique incident body waves using the finite element method, the free field motion at the incident lateral boundary induced by the background layered half-space complicates the computational area. In order to replace the complex frequency domain method, a time-domain method to calculate the free field motion of a layered half-space subjected to oblique incident body waves is developed in this paper. The new method decouples the equations of motion used in the finite element method and offers an interpolation formula of the free field motion. This formula is based on the fact that the apparent horizontal velocity of the free field motion is constant and can be calculated exactly. Both the theoretical analysis and numerical results demonstrate that the proposed method offers a high degree of accuracy. 相似文献
6.
The paper presents a numerical model for the dynamic analysis of pile groups with inclined piles in horizontally layered soil deposits. Piles are modelled with Euler–Bernoulli beams, while the soil is supposed to be constituted by independent infinite viscoelastic horizontal layers. The pile–soil–pile interaction as well as the hysteretic and geometric damping is taken into account by means of two‐dimensional elastodynamic Green's functions. Piles cap is considered by introducing a rigid constraint; the condensation of the problem permits a consistent derivation of both the dynamic impedance matrix of the soil–foundation system and the foundation input motion. These quantities are those used to perform inertial soil–structure interaction analyses in the framework of the substructure approach. Furthermore, the model allows evaluating the kinematic stress resultants in piles resulting from waves propagating in the soil deposit, taking into account the pile–soil–pile interactions. The model validation is carried out by performing accuracy analyses and comparing results in terms of dynamic impedance functions, kinematic response parameters and pile stress resultants, with those furnished by 3D refined finite element models. To this purpose, classical elastodynamic solutions are adopted to define the soil–pile interaction problem. The model results in low computational demands without significant loss of precision, compared with more rigorous approaches or refined finite element models. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
7.
区域性地震波强地面运动的量化数值模拟分析结果可以用来指导城市规划建设,并在防震减灾中发挥重要作用. 本文采用有限元数值模拟方法来模拟由地震激发的区域地震波强地面运动过程,并得到地表峰值速度和加速度的分布特征.考虑到福州是东南沿海的重要省会城市,其明显的盆地结构特征使它很容易遭受强地震灾害的影响. 因此本文以福州盆地为主要研究对象,假定盆地的主要断层某一位置在未来可能发生灾害性地震,则该地震会激发地震波的强地面运动,并由于地形、沉积层等主要影响,在盆地内部发生放大.通过对地震波传播过程的数值模拟和后处理分析,给出了该盆地的地表峰值位移、峰值速度和峰值加速度分布图谱.该图谱可为未来的福州城市规划和抗震减灾提供定量科学参考依据. 相似文献
8.
Characterizing the pore space of rock samples using three‐dimensional (3D) X‐ray computed tomography images is a crucial step in digital rock physics. Indeed, the quality of the pore network extracted has a high impact on the prediction of rock properties such as porosity, permeability and elastic moduli. In carbonate rocks, it is usually very difficult to find a single image resolution which fully captures the sample pore network because of the heterogeneities existing at different scales. Hence, to overcome this limitation a multiscale analysis of the pore space may be needed. In this paper, we present a method to estimate porosity and elastic properties of clean carbonate (without clay content) samples from 3D X‐ray microtomography images at multiple resolutions. We perform a three‐phase segmentation to separate grains, pores and unresolved porous phase using 19 μm resolution images of each core plug. Then, we use images with higher resolution (between 0.3 and 2 μm) of microplugs extracted from the core plug samples. These subsets of images are assumed to be representative of the unresolved phase. We estimate the porosity and elastic properties of each sample by extrapolating the microplug properties to the whole unresolved phase. In addition, we compute the absolute permeability using the lattice Boltzmann method on the microplug images due to the low resolution of the core plug images. In order to validate the results of the numerical simulations, we compare our results with available laboratory measurements at the core plug scale. Porosity average simulations for the eight samples agree within 13%. Permeability numerical predictions provide realistic values in the range of experimental data but with a higher relative error. Finally, elastic moduli show the highest disagreements, with simulation error values exceeding 150% for three samples. 相似文献
9.
A numerical method has been developed for the dynamic analysis of a tall building structure with viscous dampers. Viscous dampers are installed between the top of an inverted V‐shaped brace and the upper beam on each storey to reduce vibrations during strong disturbances like earthquakes. Analytically, it is modelled as a multi‐degree‐of freedom (MDOF) system with the Maxwell models. First, the computational method is formulated in the time domain by introducing a finite element of the Maxwell model into the equation of motion in the discrete‐time system, which is based on the direct numerical integration. Next, analyses for numerical stability and accuracy of the proposed method are discussed. The results show its numerical stability. Finally, the proposed method is applied to the numerical analysis of a realistic building structure to demonstrate its practical validity. 相似文献
10.
利用模拟长周期地震动的显式有限元技术和并行计算技术、场地一维线性反演和等效线性化方法以及有限断层模型,针对一个具体的城市地下结构模型,本文给出了一套近断层地震动预测的方法和相应的预测结果,并与近断层地震动的特征对比,分析了预测结果的合理性。 相似文献
11.
Modified image algorithm to simulate seismic channel waves in 3D tunnel model with rugged free surfaces 总被引：1，自引：0，他引：1 
下载免费PDF全文
下载免费PDF全文 The existence of rugged free‐surface three‐dimensional tunnel conditions in the coal seams, caused either by geological or mining processes, will inevitably influence wave propagation characteristics when the seismic waves go through the coal mines. Thus, a modified image algorithm has been developed to account for seismic channel waves propagating through this complicated topography with irregular free surfaces. Moreover, the seismic channel waves commonly exhibit damped and dispersive signatures, which is not only because of their own unique sandwich geometry of rock–coal–rock but also because of the viscoelastic behavior of coal. Considering the complexity of programming in three‐dimensional tunnel models with rugged free surfaces, an optimized vacuum grid search algorithm, enabling to model highly irregular topography and to compute efficiently, is also proposed when using high‐order staggered finite‐difference scheme to simulate seismic channel wave propagations in viscoelastic media. The numerical simulations are implemented to investigate the accuracy and stability of the method and the impact of coal's viscoelastic behavior on seismic channel wave propagation characteristics. The results indicate that the automatic vacuum grid search algorithm can be easily merged into high‐order staggered finite‐difference scheme, which can efficiently be applied to calculate three‐dimensional tunnel models with rugged free surfaces in the viscoelastic media. The simulation also suggests that the occurrence of a three‐dimensional tunnel with free surfaces has a remarkable influence on the seismic channel wave propagation characteristics and elastic energy distribution. 相似文献
12.
Intensity measures for probabilistic assessment of non‐structural components acceleration demand 下载免费PDF全文
下载免费PDF全文 Marco De Biasio Stephane Grange Frederic Dufour Frederic Allain Ilie Petre‐Lazar 《地震工程与结构动力学》2015,44(13):2261-2280
A fundamental issue in the framework of seismic probabilistic risk analysis is the choice of ground motion intensity measures (IMs). Based on the floor response spectrum method, the present contribution focuses on the ability of IMs to predict non‐structural components (NSCs) horizontal acceleration demand. A large panel of IMs is examined and a new IM, namely equipment relative average spectral acceleration (E‐ASAR), is proposed for the purpose of NSCs acceleration demand prediction. The IMs efficiency and sufficiency comparisons are based on (i) the use of a large dataset of recorded earthquake ground motions; (ii) numerical analyses performed on three‐dimensional numerical models, representing actual structural wall and frame buildings; and (iii) systematic statistical analysis of the results. From the comparative study, the herein introduced E‐ASAR shows high efficiency with respect to the estimation of maximum floor response spectra ordinates. Such efficiency is particularly remarkable in the case of structural wall buildings. Besides, the sufficiency and the simple formulation allowing the use of existing ground motion prediction models make the E‐ASAR a promising IMs for seismic probabilistic risk assessment. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
13.
Ground motion prediction is important for earthquake early warning systems, because the region’s peak ground motion indicates the potential disaster. In order to predict the peak ground motion quickly and precisely with limited station wave records, we propose a real-time numerical shake prediction and updating method. Our method first predicts the ground motion based on the ground motion prediction equation after P waves detection of several stations, denoted as the initial prediction. In order to correct the prediction error of the initial prediction, an updating scheme based on real-time simulation of wave propagation is designed. Data assimilation technique is incorporated to predict the distribution of seismic wave energy precisely. Radiative transfer theory and Monte Carlo simulation are used for modeling wave propagation in 2-D space, and the peak ground motion is calculated as quickly as possible. Our method has potential to predict shakemap, making the potential disaster be predicted before the real disaster happens. 2008 M S8.0 Wenchuan earthquake is studied as an example to show the validity of the proposed method. 相似文献
14.
The impact of the incident angle of earthquake motion on the seismic response of the long lined tunnels is studied. Based on the time‐domain finite element method with the viscous‐spring artificial boundary condition, the earthquake motion of oblique incidence is transformed into the equivalent nodal forces acting on the truncated boundary of finite element model. In the present work, the formulas of equivalent nodal forces for the plane P wave with arbitrary incident angle are deduced and implemented into the commercial software abaqus 1 . The effectiveness of the formulas and its implementation are demonstrated by two numerical examples with the reference solutions. The proposed method is applied to investigate the seismic responses of the long lined tunnels under the obliquely incident P waves. The numerical results indicate that the seismic responses of the long lined tunnels are highly affected by the incident angles of P waves. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
15.
2008年11月10日在青海柴达木盆地北缘发生了大柴旦M_W6.3地震,为了研究该地震的区域地震波传播与地面运动特征,本文利用地质资料和地壳速度结构研究成果,构建了柴达木盆地及周边区域三维传播介质模型,采用有限差分方法模拟了大柴旦地震波场传播过程以及地面运动分布特征.结果表明,柴达木盆地对波场传播有明显影响,表现为地震波传入盆地后在边界产生次生面波,盆地沉积物对地震波具有围陷作用,地震地面运动在盆地内振幅增大、持时延长.模拟结果给出的地震地面运动峰值速度分布以及理论地震图均和观测结果符合较好,反映数值模拟较好地给出了观测地面运动的主要特征以及传播介质模型的合理性. 相似文献
16.
大地震等诸多激励均能激发全球自由振荡现象,通常表现为驻波形式的全球整体振荡.现有的地震波数值模拟方法多为非保结构方法,无法压制长时程计算中的积累误差.本文采用优化的三阶辛格式谱元法,对地球自由振荡及全球尺度的地震波传播进行了长时程模拟.通过与传统的基于Newmark算法的谱元方法结果对比分析,明确验证了本文所得优化的三阶辛格式谱元法在模拟地球自由振荡等大规模长时程问题上的优越性和准确性.上述进展在方法论层面为今后探测、刻画全球尺度地球非均匀结构的驻波数值方法奠定了部分基础,并为相关研究领域提供了新的选择. 相似文献
17.
Scattering of incident plane harmonic pseudo P‐, SH‐, and SV‐waves by a two‐dimensional basin of arbitrary shape is investigated by using an indirect boundary integral equation approach. The basin and surrounding half‐space are assumed to be generally anisotropic, homogeneous, linearly elastic solids. No material symmetries are assumed. The unknown scattered waves are expressed as linear combinations of full‐space time‐harmonic two‐dimensional Green functions. Using the Radon transform, the Green functions are obtained in the form of finite integrals over a unit circle. An algorithm for the accurate and efficient numerical evaluation of the Green functions is discussed. A detailed convergence and parametric analysis of the problem is presented. Excellent agreement is obtained with isotropic results available in the literature. Steady‐state surface ground motion is presented for semi‐circular basins with generally anisotropic material properties. The results show that surface motion strongly depends upon the material properties of the basin as well as the angle of incidence and frequency of the incident wave. Significant mode conversion can be observed for general triclinic materials which are not present in isotropic models. Comparison with an isotropic basin response demonstrates that anisotropy is very important for assessing the nature of surface motion atop basins. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
18.
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. 相似文献
19.
In order to account for the site-response part of the seismic ground motion coherency for sites of interest, an analytical stochastic methodology is proposed in this paper. By combining the pseudo-excitation method with wave motion finite element simulation techniques, a numerical approach for the computation of the coherency function between observation points is developed firstly. Then the orthogonal expansion method is introduced into this approach to study the effect of the uncertainty in soil properties on the coherency function. Finally some numerical examples are given to show the applicability of the methodology. The computational results demonstrate that the lagged coherency values tend to decrease in the vicinity of the resonant frequencies of the site. 相似文献
20.
近场强地震动除受场地条件的影响外,还受到震源破裂面上子源的空间分布特点、子源破裂先后顺序的强烈控制,基于数值格林函数法的近场强地震动数值模拟方法可以综合考虑震源、传播途径及局部场地条件的影响,对计算过程进行合理简化,分2步完成地震动模拟:第1步,在介质均匀区采用矩张量的解析解计算所有子源在盖层底面的位移,形成下一步有限元计算的输入场;第2步,在盖层介质不均匀区,结合局部人工透射边界技术,采用时、空解耦的波动显式有限元方法计算地表强地震动。在有限断层模型中,采用具有9个力偶的等效地震矩张量表达断层产状、滑动方向等的影响,采用Brune模型定义各子源的滑动时间函数,描述滑动的时、空不均匀分布特征,从而细化震源模型。通过对Northridge地震中4个基岩台站地表地震动的模拟结果和强震记录,验证了此简化计算方法的可行性 相似文献