Abstract
Widespread damage caused by some recent earthquakes has posed many new challenges to the research community to study near fault motion and its effects on the ground surface. Earthquakes in different geological regions show drastic variations in their effects, such as large surface upliftment/displacements of soil deposits and distinctly different near fault ground motions in contrast to the far field ground motions. Numerous researchers studied this phenomenon through experiments and tried to establish relationships between various parameters of soil and resulting surface ground motion characteristics. Carrying out the experiments to understand such behaviour itself was difficult, because of complexities involved in replicating the exact site conditions, considering all the necessary parameters of soil deposit and underlying bedrock. Many experiments of high precision are to be conducted to establish relationship between the ground motion and various parameters such as width of affected zone, soil properties, thickness of deposit and dip angle of fault. Moreover, shapes and locations of the surface rupture and its propagation through the soil deposit before reaching the ground surface, are important in city planning, especially in the design of critical structures that may be located near the zone of faulting. Numerical modelling allows investigation of effect of a number of aspects of the fault rupture propagation, which may be difficult to obtain from earthquake case histories or from experiments conducted on physical models. Additionally, the sparse coverage of recording stations limits the ability to capture ground motions close to fault ruptures. Therefore, numerical simulation of near fault ground motions is necessary to understand the severity of near fault ground motion in future earthquakes.
In this thesis, the numerical method called Applied Element Method is explained along with its application in the fault motion analysis. Since the problem is related to the fault rupture propagation, a method is needed, which can handle the discontinuities. The methods based on continuum modeling face lots of complications when we want to apply it to the material that has discontinuities. The Applied Element Method (AEM) has many advantages with respect to the fault rupture problem. Using AEM, the crack initiation and propagation can be modeled in reasonable time by using the available parallel computing power. This method is used for studying the spatial variation of ground motion due to seismic bedrock displacement at the bedrock level.
In the present study, the Applied Element Method is used to study the dip-slip faults subjected to seismic base fault movement. The method used here has been explained and validated with analytical and experimental studies with good comparable results. The preliminary modelling of dip-slip faults has been done using different slip rate values. Increase in the response on the surface due to decrease in rise time values has been seen. Influence of dip angle has been investigated and the change in the response due to the presence of lower velocity layer on the near fault ground motion has been studied. The results of the numerical simulation have explained many dynamic and geometric aspects of the reverse faulting. In all cases with different fault dip angle, there is greater ground motion on the hanging wall side compared to the ground motion of foot wall side. This effect is intrepreted as being due to two important reasons. First, the points on the hanging wall are closer to the fault plane and secondly, the trapped seismic energy in the wedge shape hanging wall leads to multiple reflections. The results from different dip angles indicate that the near fault ground motion is sensitive to the dip angle. Decrease in the ground motion with increase in the dip angle has been seen. The presence of lower velocity layer leads to the reduction in the ground motion on the surface due to inconsistency in the material properties along the rupture direction and also due to the reduction in the stress drop value, whereas the presence of lower velocity block on the hanging wall leads to the amplification of the ground motion at the place of the lower velocity block.
Then, the fault rupture propagation is explored in presence of rock boulders and its effect on the surface peak ground acceleration. Surface faulting has been examined by keeping the boulder at different positions. Presence of boulders in the footwall and in proximity of the rupture zone increases the PGA value on the surface of hanging wall to a large extent, whereas the presence of boulder in hanging wall away from the rupture zone does not affect much on the surface ground motion. It has been seen that there is a considerable change in the rupture pattern in presence of rock boulder in the rupture zone.
Recent ground-motion observations suggest that there is a considerable difference in surface-rupturing earthquakes