scholarly journals Response Analysis of the Free Field under Fault Movements

2018 ◽  
Vol 2018 ◽  
pp. 1-15
Author(s):  
H. L. Qu ◽  
Y. Wu ◽  
B. K. Zhang ◽  
Q. D. Hu ◽  
Z. L. Xiao

A quasistatic simulation of highly nonlinear problems under fault movements was carried out using the EXPLICIT module of ABAQUS. Combined with the secondary development program of the software, the application of the strain softening Mohr–Coulomb model in the simulation was realized. Free field-fault systems were simulated with two types of fault types (normal and reverse faults), four fault dip angles (45°, 60°, 75°, and 90°), and two kinds of soil (sand and clay). Moreover, the rupture laws and sensitivities of the sand and clay were studied with different soil thicknesses and different fault dip angles in the free field. The results show that the width of the ground zone with obvious deformation, which represents the point of the fault outcrop, the critical displacement of the fault, and the rupture characteristics of the overlying soil are closely related to the fault type and soil parameters. The critical displacement of the reverse fault is larger than that of the normal fault. The width of the ground zone with obvious deformation varies from 0.65 to 1.3 and does not exhibit a regular relationship with the type of soil. Compared with a normal fault, the rupture of a reverse fault is not prone to exposure at the surface.

2019 ◽  
Vol 24 (1) ◽  
pp. 311-331 ◽  
Author(s):  
Prashant Kumar ◽  
Carmen Rodrigo ◽  
Francisco J. Gaspar ◽  
Cornelis W. Oosterlee

AbstractWe present a multilevel Monte Carlo (MLMC) method for the uncertainty quantification of variably saturated porous media flow that is modeled using the Richards equation. We propose a stochastic extension for the empirical models that are typically employed to close the Richards equations. This is achieved by treating the soil parameters in these models as spatially correlated random fields with appropriately defined marginal distributions. As some of these parameters can only take values in a specific range, non-Gaussian models are utilized. The randomness in these parameters may result in path-wise highly nonlinear systems, so that a robust solver with respect to the random input is required. For this purpose, a solution method based on a combination of the modified Picard iteration and a cell-centered multigrid method for heterogeneous diffusion coefficients is utilized. Moreover, we propose a non-standard MLMC estimator to solve the resulting high-dimensional stochastic Richards equation. The improved efficiency of this multilevel estimator is achieved by parametric continuation that allows us to incorporate simpler nonlinear problems on coarser levels for variance reduction while the target strongly nonlinear problem is solved only on the finest level. Several numerical experiments are presented showing computational savings obtained by the new estimator compared with the original MC estimator.


2014 ◽  
Vol 543-547 ◽  
pp. 414-419
Author(s):  
Liang Liang Chen ◽  
Hui Long Ren ◽  
Guo Qing Feng ◽  
Fang Duan

Due to the air cushion and splash the slamming into the water of Trimaran connecting bridge are highly nonlinear, this paper discusses the problem of typical Trimaran sections slamming into the water by use of general finite element software MSC.Dytran. Through secondary development program implemented the initialization of water pressure, considered the coupling structural-air-water, a simulation method for connecting bridge slamming into the water has developed. The results show that the method meets the engineering precision and can give reference for Trimaran construction and operation.


Author(s):  
Po Ting Lin ◽  
Wei-Hao Lu ◽  
Shu-Ping Lin

In the past few years, researchers have begun to investigate the existence of arbitrary uncertainties in the design optimization problems. Most traditional reliability-based design optimization (RBDO) methods transform the design space to the standard normal space for reliability analysis but may not work well when the random variables are arbitrarily distributed. It is because that the transformation to the standard normal space cannot be determined or the distribution type is unknown. The methods of Ensemble of Gaussian-based Reliability Analyses (EoGRA) and Ensemble of Gradient-based Transformed Reliability Analyses (EGTRA) have been developed to estimate the joint probability density function using the ensemble of kernel functions. EoGRA performs a series of Gaussian-based kernel reliability analyses and merged them together to compute the reliability of the design point. EGTRA transforms the design space to the single-variate design space toward the constraint gradient, where the kernel reliability analyses become much less costly. In this paper, a series of comprehensive investigations were performed to study the similarities and differences between EoGRA and EGTRA. The results showed that EGTRA performs accurate and effective reliability analyses for both linear and nonlinear problems. When the constraints are highly nonlinear, EGTRA may have little problem but still can be effective in terms of starting from deterministic optimal points. On the other hands, the sensitivity analyses of EoGRA may be ineffective when the random distribution is completely inside the feasible space or infeasible space. However, EoGRA can find acceptable design points when starting from deterministic optimal points. Moreover, EoGRA is capable of delivering estimated failure probability of each constraint during the optimization processes, which may be convenient for some applications.


2018 ◽  
Vol 15 (6) ◽  
pp. 661-677 ◽  
Author(s):  
Toufiq Ouzandja ◽  
Mohamed Hadid

Purpose This paper aims to present the investigation of the linear and nonlinear seismic site response of a saturated inhomogeneous poroviscoelastic soil profile for different soil properties, such as pore-water saturation, non-cohesive fines content FC, permeability k, porosity n and coefficient of uniformity Cu. Design/methodology/approach The inhomogeneous soil profile is idealized as a multi-layered saturated poroviscoelastic medium and is characterized by the Biot’s theory, with a shear modulus varying continuously with depth according to the Wichtmann’s model. Seismic response analysis has been evaluated through a computational model, which is based on the exact stiffness matrix method formulated in the frequency domain assuming that the incoming seismic waves consist of inclined P-SV waves. Findings Unlike the horizontal seismic response, the results indicate that the vertical one is strongly affected by the pore water saturation. Moreover, in the case of fully saturated soil profile, the same vertical response spectra are found for the two cases of soil behavior, linear and nonlinear. Originality/value This research is a detailed study of the geotechnical soil properties effect on the bi-directional seismic response of saturated inhomogeneous poroviscoelastic soil profile, which has not been treated before; the results are presented in terms of the peak acceleration ratio, as well as the free-field response spectra and the spectral ratio (V/H).


2000 ◽  
Vol 37 (6) ◽  
pp. 1368-1382 ◽  
Author(s):  
Kevin J Bentley ◽  
M Hesham El Naggar

Recent destructive earthquakes have highlighted the need for increased research into the revamping of design codes and building regulations to prevent further catastrophic losses in terms of human life and economic assets. The present study investigated the response of single piles to kinematic seismic loading using the three-dimensional finite element program ANSYS. The objectives of this study were (i) to develop a finite element model that can accurately model the kinematic soil–structure interaction of piles, accounting for the nonlinear behaviour of the soil, discontinuity conditions at the pile–soil interface, energy dissipation, and wave propagation; and (ii) to use the developed model to evaluate the kinematic interaction effects on the pile response with respect to the input ground motion. The static performance of the model was verified against exact available solutions for benchmark problems including piles in elastic and elastoplastic soils. The geostatic stresses were accounted for and radiating boundaries were provided to replicate actual field conditions. Earthquake excitation with a low predominant frequency was applied as an acceleration–time history at the base bedrock of the finite element mesh. To evaluate the effects of the kinematic loading, the responses of both the free-field soil (with no piles) and the pile head were compared. It was found that the effect of the response of piles in elastic soil was slightly amplified in terms of accelerations and Fourier amplitudes. However, for elastoplastic soil with separation allowed, the pile head response closely resembled the free-field response to the low-frequency seismic excitation and the range of pile and soil parameters considered in this study.Key words: numerical modelling, dynamic, lateral, piles, kinematic, seismic.


2013 ◽  
Vol 405-408 ◽  
pp. 1939-1944
Author(s):  
Gui Lan Tao ◽  
Li Zhang

Spring-damper units were set on the boundaries to absorb incident waves and reflected scattering waves to realize viscoelastic artificial boundary (VAB). The equivalent node load input method was used to simulate the VAB and viscoelastic boundary element wave input. Programming is based on APDL secondary development language with ANSYS finite element software. Considering the interaction between chamber structure and the surrounding soil, docking chamber structure dynamic model is established based on the VAB. The linear elastic model was used for concrete structure. The D-P nonlinear model was used for the back soil calculation. Docking chamber structure dynamic analysis under conditions of fixed boundaries and viscoelastic boundaries were conducted. The result indicated that under the viscoelastic boundary conditions, dynamic acceleration response is significant on the top of the lock wall, which is approximately 2.5 times of the value on the bottom of the lock wall. The maximum response stress appears near the cross point of the lock wall and the bottom floor with value of approximately 5620 kPa;.The chamber bottom floor is subjected to tension and maximum stress with the value of approximately 6180 kPa. Usually, the structure response under the fixed boundary conditions is higher than the structure response under the viscoelastic boundary conditions.


2003 ◽  
Vol 125 (3) ◽  
pp. 533-539 ◽  
Author(s):  
Zekai Ceylan ◽  
Mohamed B. Trabia

Welded cylindrical containers are susceptible to stress corrosion cracking (SCC) in the closure-weld area. An induction coil heating technique may be used to relieve the residual stresses in the closure-weld. This technique involves localized heating of the material by the surrounding coils. The material is then cooled to room temperature by quenching. A two-dimensional axisymmetric finite element model is developed to study the effects of induction coil heating and subsequent quenching. The finite element results are validated through an experimental test. The container design is tuned to maximize the compressive stress from the outer surface to a depth that is equal to the long-term general corrosion rate of the container material multiplied by the desired container lifetime. The problem is subject to several geometrical and stress constraints. Two different solution methods are implemented for this purpose. First, an off-the-shelf optimization software is used. The results however were unsatisfactory because of the highly nonlinear nature of the problem. The paper proposes a novel alternative: the Successive Heuristic Quadratic Approximation (SHQA) technique. This algorithm combines successive quadratic approximation with an adaptive random search within varying search space. SHQA promises to be a suitable search method for computationally intensive, highly nonlinear problems.


1994 ◽  
Vol 37 (6) ◽  
Author(s):  
G. C. P. King ◽  
R. M. Wood

The character of the hydrological changes that follow major earthquakes has been investigated and found to be critically dependent on the style of fault displacement. In areas where fracture-flow in the crystalline crust communicates uninterrupted with the surface the most significant response is found to accompany major normal fault earthquakes. Increases in spring and river discharges peak a few days after the earthquake and typically excess flow is sustained for a period of 4 12 months. Rainfall equivalent discharges, have been found to ceed 100 mm close to the fault and remain above 10 mm at distances greater than 50 km. The total volume of water released in two M 7 normal fault earthquakes in the Western U.S.A. was 0.3-0.5 km3. In contrast, hydroIogical changes accompanying reverse fault earthquakes are either undetected or else involve falls in well-levels and spring-flows. The magnitude and distribution of the water-discharge for these events is compared with deformation models calibrated from seismic and geodetic information, and found to correlate with the crustal volume strain down to a depth of at least 5 km. Such relatively rapid drainage is only possible if the fluid was formerly contained in high aspect ratio fissures interconnected throughout much of the seismogenic upper crust. The rise and decay times of the discharge are shown to be critically dependent on crack widths, for which the «characteristic» or dominant cracks cannot be wider than 0.03 mm. These results suggest that fluid-filled cracks are ubiquitous throughout the brittle continental crust, and that these cracks open and close through the earthquake cycle. Seismohydraulic fluid flows have major implications for our understanding of the mechanical and chemical behaviour of crustal rocks, of the tectonic controls of fluid flow associated with petroleum migration, hydrothermal mineralisation and a significant hazard for underground waste disposal.


Author(s):  
Miftahul Jannah ◽  
Adi Suryadi ◽  
Muchtar Zafir ◽  
Randi Saputra ◽  
Ihsanul Hakim ◽  
...  

On the study area there are three types of structure, those are fault, fold and joint. Types of fault were found  in the study area, reverse fault with the strike/dip is N215oE/75o, normal fault has a fault directions N22oE and N200oE with pitch 35o, and dextral fault with pitch 10o and strike N219oE. Fold and joint structures used to determine the direction of the main stress on the study area. Further, an analysis used stereonet for data folds and joints. So that from the data got three directions of main stress, those are Northeast – Southwest (T1), North – South (T2) and Southeast – Northwest (T3). On the Northeast – Southwest (T1) stress there are four geological structures, anticline fold at ST.3 , syncline folds at ST. 13a, ST. 13b, ST. 13c and ST. 33, chevron fold at ST. 44 and joint at ST. 2. On the North – South (T2) stress there are three geological structures, those are syncline fold at ST. 35, anticline fold at ST. 54 and joints at ST. 41, ST. 46 and ST. 47. On the Southeast – Northwest (T3) stress were also three geological structures, those are chevron fold at ST 42a, overturned fold at ST. 42b, syncline fold at ST. 42c and joints at ST. 5 and ST. 34.


Sign in / Sign up

Export Citation Format

Share Document