154 resultados para numerical simulations
Resumo:
Large earthquakes, such as the Chile earthquake in 1960 and the Sumatra-Andaman earthquake on Dec 26, 2004 in Indonesia, have generated the Earth’s free oscillations. The eigenfrequencies of the Earth’s free oscillations are closely related to the Earth’s internal structures. The conventional methods, which mainly focus on calculating the eigenfrequecies by analytical ways, and the analysis on observations can not easily study the whole processes from earthquake occurrence to the Earth’s free oscillation inspired. Therefore, we try to use numerical method incorporated with large-scale parallel computing to study on the Earth’s free oscillations excited by giant earthquakes. We first give a review of researches and developments of the Earth’s free oscillation, and basical theories under spherical coordinate system. We then give a review of the numerical simulation of seismic wave propagation and basical theories of spectral element method to simulate global seismic wave propagation. As a first step to study the Earth’s free oscillations, we use a finite element method to simulate the propagation of elastic waves and the generation of oscillations of the chime bell of Marquis Yi of Zeng, by striking different parts of the bell, which possesses the oval crosssection. The bronze chime bells of Marquis Yi of Zeng are precious cultural relics of China. The bells have a two-tone acoustic characteristic, i.e., striking different parts of the bell generates different tones. By analysis of the vibration in the bell and the spectrum analysis, we further help the understanding of the mechanism of two-tone acoustic characteristics of the chime bell of Marquis Yi of Zeng. The preliminary calculations have clearly shown that two different modes of oscillation can be generated by striking different parts of the bell, and indicate that finite element numerical simulation of the processes of wave propagation and two-tone generation of the chime bell of Marquis Yi of Zeng is feasible. These analyses provide a new quantitative and visual way to explain the mystery of the two-tone acoustic characteristics. The method suggested by this study can be applied to simulate free oscillations excited by great earthquakes with complex Earth structure. Taking into account of such large-scale structure of the Earth, small-scale low-precision numerical simulation can not simply meet the requirement. The increasing capacity in high-performance parallel computing and progress on fully numerical solutions for seismic wave fields in realistic three-dimensional spherical models, Spectral element method and high-performance parallel computing were incorporated to simulate the seismic wave propagation processes in the Earth’s interior, without the effects of the Earth’s gravitational potential. The numerical simulation shows that, the results of the toroidal modes of our calculation agree well with the theoretical values, although the accuracy of our results is much limited, the calculated peaks are little distorted due to three-dimensional effects. There exist much great differences between our calculated values of spheroidal modes and theoretical values, because we don’t consider the effect the Earth’ gravitation in numerical model, which leads our values are smaller than the theoretical values. When , is much smaller, the effect of the Earth’s gravitation make the periods of spheroidal modes become shorter. However, we now can not consider effects of the Earth’s gravitational potential into the numerical model to simulate the spheroidal oscillations, but those results still demonstrate that, the numerical simulation of the Earth’s free oscillation is very feasible. We make the numerical simulation on processes of the Earth’s free oscillations under spherically symmetric Earth model using different special source mechanisms. The results quantitatively show that Earth’s free oscillations excited by different earthquakes are different, and oscillations at different locations are different for free oscillation excited by the same earthquake. We also explore how the Earth’s medium attenuation will take effects on the Earth’s free oscillations, and take comparisons with the observations. The medium attenuation can make influences on the Earth’s free oscillations, though the effects on lower-frequency fundamental oscillations are weak. At last, taking 2008 Wenchuan earthquake for example, we employ spectral element method incorporated with large-scale parallel computing technology to investigate the characteristics of seismic wave propagation excited by Wenchuan earthquake. We calculate synthetic seismograms with one-point source model and three-point source model respectively. Full 3-D visualization of the numerical results displays the profile of the seismic wave propagation with respect to time. The three-point source, which was proposed by the latest investigations through field observation and reverse estimation, can better demonstrate the spatial and temporal characteristics of the source rupture processes than one-point source. Primary results show that those synthetic signals calculated from three-point source agree well with the observations. This can further reveal that the source rupturing process of Wenchuan earthquake is a multi-rupture process, which is composed by at least three or more stages of rupture processes. In conclusion, the numerical simulation can not only solve some problems concluding the Earth’s ellipticity and anisotropy, which can be easily solved by conventional methods, but also finally solve the problems concluding topography model and lateral heterogeneity. We will try to find a way to fully implement self-gravitation in spectral element method in future, and do our best to continue researching the Earth’s free oscillations using the numerical simulations to see how the Earth’ lateral heterogeneous will affect the Earth’s free oscillations. These will make it possible to bring modal spectral data increasingly to bear on furthering our understanding of the Earth’s three-dimensional structure.
Resumo:
As a marginal subject, dynamic responses of slopes is not only an important problem of engineering geology (Geotechnical problem), but also of other subjects such as seismology, geophysics, seismic engineering and engineering seismic and so on. Owning to the gulf between different subjects, it is arduous to study dynamic responses of slopes and the study is far from ripeness. Studying on the dynamic responses of slopes is very important in theories as well as practices. Supported by hundreds of bibliographies, this paper systemically details the development process of this subject, introduces main means to analyze this subject, and then gives brief remarks to each means respectively. Engineering geology qualitative analysis is the base of slopes dynamic responses study. Because of complexity of geological conditions, engineering geology qualitative analysis is very important in slopes stability study, especially to rock slopes with complex engineering geology conditions. Based on research fruits of forerunners, this paper summarizes factors influencing slopes dynamic stability into five aspects as geology background, stratums, rock mass structure, and topography as well as hydrogeology condition. Based on rock mass structure controlling theory, engineering geology model of the slope is grouped into two typical classes, one is model with obvious controlling discontinuities, which includes horizontal bedded slope, bedding slope, anti-dip slope, slide as well as slope with base rock and weathered crust; the other is model without obvious controlling discontinuities, which includes homogeneous soil slope, joint rock mass slope. Study on slope failure mechanism under dynamic force, the paper concludes that there are two effects will appear in slope during strong earthquake, one is earthquake inertia force, the other is ultra pore pressure buildup. The two effects lead to failure of the slope. To different types of slope failure, the intensity of two effects acting on the slope is different too. To plastic flow failure, pore pressure buildup is dominant; to falling rock failure and toppling failure, earthquake inertia force is dominant in general. This paper briefly introduces the principle of Lagrangian element method. Through a lot of numerical simulations with FLAC3D, the paper comprehensively studies dynamic responses of slopes, and finds that: if the slope is low, displacement, velocity and acceleration are linear enlarging with elevation increasing in vertical direction; if the slope is high enough, displacement, velocity and acceleration are not linear with elevation any more, on the other hand, they fluctuate with certain rhythm. At the same time, the rhythm appears in the horizontal direction in the certain area near surface of the slope. The distribution form of isoline of displacement, velocity and acceleration in the section of the slope is remarkably affected by the slope angle. In the certain area near the slope surface, isoline of displacement, velocity and acceleration is parallel to the surface of the slope, in the mean time, the strike direction of the extreraum area is parallel to the surface of the slope too. Beyond this area, the isoline direction and the strike direction of the extremum area turn to horizontal with invariable distance. But the rhythm appearing or not has nothing to with the slope angle. The paper defines the high slope effect and the low slope effect of slopes dynamic responses, discusses the threshold height H^t of the dynamic high slope effect, and finds that AW is proportional to square root of the dynamic elastic moduli El P , at the same time, it is proportional to period Tof the dynamic input. Thus, the discriminant of H^t is achieved. The discriminant can tell us that to a slope, if its height is larger than one fifth of the wavelength, its response regular will be the dynamic high slope effect; on the other hand, its response regular will be the dynamic low slope effect. Based on these, the discriminant of different slopes taking on same response under the same dynamic input is put forward in this paper. At the same time, the paper studies distribution law of the rhythm extremum point of displacement, velocity and acceleration, and finds that there exists relationship of N = int among the slope height H, the number of the rhythm extremum
VHlhro)
point N and ffthre- Furthermore, the paper points out that if N^l, the response of the slope will be dynamic high slope effect; \fN
Resumo:
In order to realize fast development of the national economy in a healthy way and coordinate progress with whole society, the country has implemented the strategy of development of the western region. An important action of finishing this strategic task is to accelerate the highway construction in the western region, join the western region and places along the coast, the river, the border with goods and materials, technology, and personnel interchanges, and then drive development of the local economy.The western region was influenced by the Himalaya Tectonization in Cenozoic, and the crust rose and became the plateau. In the course of rising, rivers cut down sharply to form a lot of high mountains and gorges.Because of topography and geomorphology, bridges in the traffic construction in the alpine gorge area are needed. Rivers have characteristics of large flow, fast velocity and high and steep river valley, so building a pier in the river is not only very difficult, but also making the cost increase. At the same time, the impact that the pier is corroded and the bridge base that is drawn to be empty by flow are apt to cause destruction of the pier. For those reasons, suspending bridge and cable-stay bridge are usually adopted with the single and large span. For the large span bridge, the pier foundation could receive ten thousand and more vertical strength, bending moment and near kiloton horizontal thrust.Because bank slope in the alpine gorge district is cut deeply and unsettled big, natural stability is worse under endogenic and exogenic force. When bank slope bears heavy vertical strength, bending moment and horizontal thrust facing the river, it will inevitably make the balance state of rock and soil mass change, bridge bank slope deform, and even destroyed. So the key problem at the time of the large span's bridge construction in the alpine gorge area is how to make it stable.So based on the spot investigation, the Engineering Geology Analysis Method is very important to grasp the bank slope stability. It can provide the bank slope stability macroscopic ally and qualitatively, and reference to the indoor calculation. The Engineering Geology Analysis Method is that by way of analyzing and investigating terms of bank slope instability, stability development trend, the ancient rock slide and devolution in the site, stability comprehensive evaluation primarily, current and future stability of bank slope is gotten, realizing the intention to serving the concrete engineering.After the Engineering Geology Analysis Method is applied to project instances of BeiPan River Bridge and BaLin River Bridge, results are accord with bank slope actual conditions, which proves sites are suited to building bridges from site stability.we often meet bank slope stability issues in the traffic construction in the alpine gorge areao Before the evaluation of the bank slope stability, the engineering geological condition is investigated first. After that, the next exploration target and geology measures are decided. So, the Engineering Geology Analysis Method that the investigation of the engineering geological condition is the main content is quite important in practice. The other evaluations of the bank slope stability are based on it. Because foundation receives very heavy load, for the big span's bridge in the alpine gorge area, a long pile of the large diameter (D^0.8m) is usually selected. In order to reflect rock mass's deformation properties under rock-socketed pile function, the author has used the FLAG30 software for rock and soil mass and done many numerical simulations. By them, the author launches the further investigation on deformation properties of bank slope under different slope angle, pile length, diameter, elastic modulus, load, bank slope's structure, etc. Some conclusion meaningful to the design and produce are obtained.
Resumo:
Earlier studies on the distribution of geological environmental indicators in China revealed drastic changes from a zonal climate pattern (planetary-wave-dominant pattern) in the Paleogene to a monsoon-dominant one in the Neogene, which suggested an inception of the initial East-Asian summer monsoon. However, there are different views about the time and causes of the changes.Here, we attempt to compile a series of paleoenvironmental maps based on newly collected climate indicators from the literatures and chronologically constrained evidence of geological maps in order to re-examine the temporal and spatial evolution of climate belts in China during the Cenozoic with special emphasis on the changes of the arid belt. These indicators include mammalian fauna, coal, carbonate concretions, jarosite, salt, gypsum deposits and pollen assemblages etc, with chronological controls that we believe reliable. Pollen assemblages and mammalian fauna have been classified into three categories (arid, semi-arid/sub-humid, humid) to reflect the intensity of aridity/humidity. Salt, jarosite and gypsum deposits are classified as the arid indicators. Carbonate concretions and coal are classified into the semi-arid/sub-humid and humid one respectively. Paleoenvironmental maps at 8 time slices have been reconstructed. They are the Paleocene, Eocene, Oligocene, Miocene, Early Miocene, Middle Miocene, Late Miocene and Pliocene.And furthermore, we attempt to use IAP^AGCM to simulate the evolution of climate belts in emphasizing on the changes of the rain band, and compare the results with the paleoenvironmental maps in order to examine the causes of the drastic paleoenvironmental changes near the Oligocene/Miocene boundary. 36 sensitive numerical experiments are carried out using the IAP__AGCM to analyze the impacts of the uplift of the Himalayan-Tibetan complex, shrinkage of the Paratethys Sea, expansion of the South China Sea and the development of the polar ice sheets on rain band in China.The main conclusions are as follows:The obtained results essentially confirm the earlier conclusions about a zonal climate pattern in the Paleogene and a different pattern in the Neogene, and illustrate that a monsoon-dominant environmental pattern with inland aridity formed by the Early Miocene, which is temporally consist with the onset of eolian deposits in China.Cenozoic cooling and the formation of polar ice sheets are unlikely the main causes to the changes of environmental patterns mentioned above in China. But northern hemispheric cooling and the ice-sheets can intensify the Siberian High Pressure, and strengthen the winter monsoon circulations and enhance the aridity in the west part of China. These results support the earlier studies.Shrinkage of the Paratethys Sea and uplift of the Himalayan-Tibetan complex played important roles in strengthening the East Asian monsoon and induceing the above changes of environmental pattern, which is consistent with the earlier studies. Furthermore, "the monsoon-dominant pattern" appears when the Himalayan-Tibetan complex reaches to about 1000-2000 meters high and the Paratethys Sea retreats to the Turan Plate.4) Expansion of the South China Sea is another significant factor that drives the evolution of environmental patterns. We believe that the above three factors co-act and drive the change of the environmental patterns from a planetary-wave-dominant one to a monsoon-dominant one. However, the impacts of each factor vary by regions. The uplift mainly increases the humidity in Southwestern China and the aridity in northwestern country. The shrinkage mainly increases the humidity in Northern China and also enhances the aridity in the northwestern country. The expansion greatly increases the humidity in the south part of China.
