Validation Data and Model Development for Fuel Assembly Response to Seismic Loads

DOE Office of Scientific and Technical Information (OSTI.GOV)

Bardet, Philippe; Ricciardi, Guillaume

2016-01-31

Vibrations are inherently present in nuclear reactors, especially in cores and steam generators of pressurized water reactors (PWR). They can have significant effects on local heat transfer and wear and tear in the reactor and often set safety margins. The simulation of these multiphysics phenomena from first principles requires the coupling of several codes, which is one the most challenging tasks in modern computer simulation. Here an ambitious multiphysics multidisciplinary validation campaign is conducted. It relied on an integrated team of experimentalists and code developers to acquire benchmark and validation data for fluid-structure interaction codes. Data are focused on PWRmore » fuel bundle behavior during seismic transients.« less

Self-induced seismicity due to fluid circulation along faults

NASA Astrophysics Data System (ADS)

Aochi, Hideo; Poisson, Blanche; Toussaint, Renaud; Rachez, Xavier; Schmittbuhl, Jean

2014-03-01

In this paper, we develop a system of equations describing fluid migration, fault rheology, fault thickness evolution and shear rupture during a seismic cycle, triggered either by tectonic loading or by fluid injection. Assuming that the phenomena predominantly take place on a single fault described as a finite permeable zone of variable width, we are able to project the equations within the volumetric fault core onto the 2-D fault interface. From the basis of this `fault lubrication approximation', we simulate the evolution of seismicity when fluid is injected at one point along the fault to model-induced seismicity during an injection test in a borehole that intercepts the fault. We perform several parametric studies to understand the basic behaviour of the system. Fluid transmissivity and fault rheology are key elements. The simulated seismicity generally tends to rapidly evolve after triggering, independently of the injection history and end when the stationary path of fluid flow is established at the outer boundary of the model. This self-induced seismicity takes place in the case where shear rupturing on a planar fault becomes dominant over the fluid migration process. On the contrary, if healing processes take place, so that the fluid mass is trapped along the fault, rupturing occurs continuously during the injection period. Seismicity and fluid migration are strongly influenced by the injection rate and the heterogeneity.

Seperating Long-term Hydrological Loading and Tectonic Deformation Observed with Multi-temporal SAR Interferometry and GPS in Qinghai-Tibetan Plateau

NASA Astrophysics Data System (ADS)

LI, G.; Lin, H.

2014-12-01

From 2000 till present, most endorheic lakes in Tibetan plateau experienced quick increasing. Several largest lakes, gathered several meters depth water during one decade. Such massive mass increasing will lead to elastic and visco-elastic deformation of the ground. Qinghai-Tibetan Plateau is one the most active tectonic places in the world; monitoring its ground deformation is essential, when loading effect is a nuisance item. Due to the sparse distribution of GPS sites and most are roving sites, it is hard to distinguish tectonic component from mass loading effect. In this research we took Selin Co Lake located at Nujiang-Bangoin suture zone and evaluated long time ground deformation at hundred kilometers scale by multi-temporal SAR interferometry and simulate the ground deformation by loading history evaluated by multi mission satellite altimetry and optical images observation. At Nujiang-Bangoin suture zone, where GPS presented the maximum ground subsidence in Qinghai-Tibetan Plateau of 3.6mm/a which was found at the shore of Selin Co Lake from 1999 to 2011, when it experienced water level increasing of 0.7m/a. A model of elastic plate lying over Newtonian viscous half-space matches well with the results of multi-temporal SAR interferometry and GPS observations. We concluded that near Selin Co Lake area, ground deformation is composed by both tectonic and hydrological loading part. As SAR image coverage is much smaller than tectonic scale, we contribute the deformation detected by InSAR to loading effect. After evaluating and removing the hydrological loading effect, we founds that Nujiang-Bangoin suture zone did not experience quick subsidence, but only limited to 0.5mm/a. Selin Co Lake's quick volume increasing caused 3mm/a subsidence rate to the nearest GPS site. The Second nearest site showed the 1.4mm/a subsidence totally, which were composed by 1.05mm/a hydrological loading effect and the rest was tectonic. We also found that Young's Modulus is the most essential parameter for loading effect simulation, and our simulation gave the similar Young's Modulus as the previous seismic tomographic INDEPTH III program did. Therefore with accurate seismic tomographic results and loading history detected by remote sensing could accurately simulate ground deformation caused by hydrological loading.

Seismic performance of geosynthetic-soil retaining wall structures

NASA Astrophysics Data System (ADS)

Zarnani, Saman

Vertical inclusions of expanded polystyrene (EPS) placed behind rigid retaining walls were investigated as geofoam seismic buffers to reduce earthquake-induced loads. A numerical model was developed using the program FLAC and the model validated against 1-g shaking table test results of EPS geofoam seismic buffer models. Two constitutive models for the component materials were examined: elastic-perfectly plastic with Mohr-Coulomb (M-C) failure criterion and non-linear hysteresis damping model with equivalent linear method (ELM) approach. It was judged that the M-C model was sufficiently accurate for practical purposes. The mechanical property of interest to attenuate dynamic loads using a seismic buffer was the buffer stiffness defined as K = E/t (E = buffer elastic modulus, t = buffer thickness). For the range of parameters investigated in this study, K ≤50 MN/m3 was observed to be the practical range for the optimal design of these systems. Parametric numerical analyses were performed to generate design charts that can be used for the preliminary design of these systems. A new high capacity shaking table facility was constructed at RMC that can be used to study the seismic performance of earth structures. Reduced-scale models of geosynthetic reinforced soil (GRS) walls were built on this shaking table and then subjected to simulated earthquake loading conditions. In some shaking table tests, combined use of EPS geofoam and horizontal geosynthetic reinforcement layers was investigated. Numerical models were developed using program FLAC together with ELM and M-C constitutive models. Physical and numerical results were compared against predicted values using analysis methods found in the journal literature and in current North American design guidelines. The comparison shows that current Mononobe-Okabe (M-O) based analysis methods could not consistently satisfactorily predict measured reinforcement connection load distributions at all elevations under both static and dynamic loading conditions. The results from GRS model wall tests with combined EPS geofoam and geosynthetic reinforcement layers show that the inclusion of a EPS geofoam layer behind the GRS wall face can reduce earth loads acting on the wall facing to values well below those recorded for conventional GRS wall model configurations.

Staging of the Acoustic Response at Laboratory Modelling of Tidal Influence upon Seismicity

NASA Astrophysics Data System (ADS)

Saltykov, Vadim; Patonin, Andrey; Kugaenko, Yulia

2010-05-01

INTRODUCTION The seismic radiation is varied through the wide range of seismic energy from seismic emission (high-frequency seismic noise, HFSN) to earthquakes. Some features of external influence response on the different scales allow to consider the medium as a single whole seismoactive object. Earth tide is a bright example of external excited field. Tidal topic has long history in seismology. Results obtained by different scientists are contradictory and ambiguous often. We denoted instability of tidal effect manifestation as possible reason of this situation. In view of the aforesaid it is significant, that tidal effects in weak seismicity and HFSN prove more strongly in the stage of large earthquake preparation [Rykunov et al., 1998, Saltykov et al., 2004, 2007]. It is presumed that the metastable medium has more high tidal sensitivity. For example, sources of prepared earthquakes and extensive near-surface zones of micro-fissuring and dilatancy, which appear during source formation and stretch far enough. [Alekseev et all., 2001, Goldin, 2004, 2005]. Common features of observed effects allow to suggest existence of tidal modulation mechanism, which is similar (may be single) for different seismic scales. Modelling of these processes can improve our understanding of tidal effect nature. LABORATORY EXPERIMENT Results of rock sample destruction experiments under controlling are presented. Acoustic emission (AE) pulses act as analogue of seismic events. Tides are simulated by weak long-period variations added to quasi-stationary subcritical loading. The results of tidal modeling confirmed AE intensity synchronization with external periodic influence with large (5-10%) variations of loading are known [Lockner, Beeler, 1999, Ponomarev et al., 2007]. But real (in nature) tidal strain&stress variations are much less and equal to splits of percent. Therefore, investigation of weak modulation influence upon deformed rock is one of main proposed purposes. Used software-programmable electro-hydraulic system INOVA [Patonin, 2006], can provide various procedures of experiment, among them programmable modulatory action. Axial deformation with stable strain rate and additional action of meander with specified period and amplitude was chosen as mode of operation. The relation between background and periodic strains reaches three orders, which corresponds to real relation between maximal tectonic and tidal strains. RESULTS For detection of periodic loading modulation of AE we used procedure based on Rayleigh criteria of uniformity and considered uniformity of AE impulses distribution on time interval, multiple to period of loading. Moreover, the predominant phase of periodical loading, corresponding to maximal AE activity, was calculated in sliding time window. In all experiments we observed instability of modulation effects. So the following stages were distinguished: - synchronization of AE and periodic loading at the initial part of test; - absence of synchronization at the elastic stage; - resumption of synchronization during plastic deformation. Stability of phase corresponding to maximal AE activity was discovered within the initial part and plastic deformation stage. Absolute values of phase for initial loading and during plastic deformation are different. CONCLUSION Now we regard revealed staging of AE response to weak periodical loading as our main result of these experiments. Different stages of AE response are connected with different state of rock samples during loading and destruction. Observed effects of synchronization can be considered as analogue of tidal modulation of HFSN and appearance of "tidal" seismicity in source zone of prepared large earthquake. This investigation was supported by RFBR, grant 08-05-00692.

Analytical simulation of nonlinear response to seismic test excitations of HDR-VKL (Heissdampfreaktor-Versuchskreislauf) piping system

DOE Office of Scientific and Technical Information (OSTI.GOV)

Srinivasan, M.G.; Kot, C.A.; Mojtahed, M.

The paper describes the analytical modeling, calculations, and results of the posttest nonlinear simulation of high-level seismic testing of the VKL piping system at the HDR Test Facility in Germany. One of the objectives of the tests was to evaluate analytical methods for calculating the nonlinear response of realistic piping systems subjected to high-level seismic excitation that would induce significant plastic deformation. Two out of the six different pipe-support configurations, (ranging from a stiff system with struts and snubbers to a very flexible system with practically no seismic supports), subjected to simulated earthquakes, were tested at very high levels. Themore » posttest nonlinear calculations cover the KWU configuration, a reasonably compliant system with only rigid struts. Responses for 800% safe-shutdown-earthquake loading were calculated using the NONPIPE code. The responses calculated with NONPIPE were found generally to have the same time trends as the measurements but contained under-, over-, and correct estimates of peak values, almost in equal proportions. The only exceptions were the peak strut forces, which were underestimated as a group. The scatter in the peak value estimate of displacements and strut forces was smaller than that for the strains. The possible reasons for the differences and the effort on further analysis are discussed.« less

Safety Identifying of Integral Abutment Bridges under Seismic and Thermal Loads

PubMed Central

Easazadeh Far, Narges; Barghian, Majid

2014-01-01

Integral abutment bridges (IABs) have many advantages over conventional bridges in terms of strength and maintenance cost. Due to the integrity of these structures uniform thermal and seismic loads are known important ones on the structure performance. Although all bridge design codes consider temperature and earthquake loads separately in their load combinations for conventional bridges, the thermal load is an “always on” load and, during the occurrence of an earthquake, these two important loads act on bridge simultaneously. Evaluating the safety level of IABs under combination of these loads becomes important. In this paper, the safety of IABs—designed by AASHTO LRFD bridge design code—under combination of thermal and seismic loads is studied. To fulfill this aim, first the target reliability indexes under seismic load have been calculated. Then, these analyses for the same bridge under combination of thermal and seismic loads have been repeated and the obtained reliability indexes are compared with target indexes. It is shown that, for an IAB designed by AASHTO LRFD, the indexes have been reduced under combined effects. So, the target level of safety during its design life is not provided and the code's load combination should be changed. PMID:25405232

Behavior of Steel-Sheathed Shear Walls Subjected to Seismic and Fire Loads.

PubMed

Hoehler, Matthew S; Smith, Christopher M; Hutchinson, Tara C; Wang, Xiang; Meacham, Brian J; Kamath, Praveen

2017-07-01

A series of tests was conducted on six 2.7 m × 3.7 m shear wall specimens consisting of cold-formed steel framing sheathed on one side with sheet steel adhered to gypsum board and on the opposite side with plain gypsum board. The specimens were subjected to various sequences of simulated seismic shear deformation and fire exposure to study the influence of multi-hazard interactions on the lateral load resistance of the walls. The test program was designed to complement a parallel effort at the University of California, San Diego to investigate a six-story building subjected to earthquakes and fires. The test results reported here indicate that the fire exposure caused a shift in the failure mode of the walls from local buckling of the sheet steel in cases without fire exposure, to global buckling of the sheet steel with an accompanying 35 % reduction in lateral load capacity after the wall had been exposed to fire. This behavior appears to be predictable, which is encouraging from the standpoint of residual lateral load capacity under these severe multi-hazard actions.

Modelling strong seismic ground motion: three-dimensional loading path versus wavefield polarization

NASA Astrophysics Data System (ADS)

Santisi d'Avila, Maria Paola; Lenti, Luca; Semblat, Jean-François

2012-09-01

Seismic waves due to strong earthquakes propagating in surficial soil layers may both reduce soil stiffness and increase the energy dissipation into the soil. To investigate seismic wave amplification in such cases, past studies have been devoted to one-directional shear wave propagation in a soil column (1D-propagation) considering one motion component only (1C-polarization). Three independent purely 1C computations may be performed ('1D-1C' approach) and directly superimposed in the case of weak motions (linear behaviour). This research aims at studying local site effects by considering seismic wave propagation in a 1-D soil profile accounting for the influence of the 3-D loading path and non-linear hysteretic behaviour of the soil. In the proposed '1D-3C' approach, the three components (3C-polarization) of the incident wave are simultaneously propagated into a horizontal multilayered soil. A 3-D non-linear constitutive relation for the soil is implemented in the framework of the Finite Element Method in the time domain. The complex rheology of soils is modelled by mean of a multisurface cyclic plasticity model of the Masing-Prandtl-Ishlinskii-Iwan type. The great advantage of this choice is that the only data needed to describe the model is the modulus reduction curve. A parametric study is carried out to characterize the changes in the seismic motion of the surficial layers due to both incident wavefield properties and soil non-linearities. The numerical simulations show a seismic response depending on several parameters such as polarization of seismic waves, material elastic and dynamic properties, as well as on the impedance contrast between layers and frequency content and oscillatory character of the input motion. The 3-D loading path due to the 3C-polarization leads to multi-axial stress interaction that reduces soil strength and increases non-linear effects. The non-linear behaviour of the soil may have beneficial or detrimental effects on the seismic response at the free surface, depending on the energy dissipation rate. Free surface time histories, stress-strain hysteresis loops and in-depth profiles of octahedral stress and strain are estimated for each soil column. The combination of three separate 1D-1C non-linear analyses is compared to the proposed 1D-3C approach, evidencing the influence of the 3C-polarization and the 3-D loading path on strong seismic motions.

Seismic Performance of Composite Shear Walls Constructed Using Recycled Aggregate Concrete and Different Expandable Polystyrene Configurations.

PubMed

Liu, Wenchao; Cao, Wanlin; Zhang, Jianwei; Qiao, Qiyun; Ma, Heng

2016-03-02

The seismic performance of recycled aggregate concrete (RAC) composite shear walls with different expandable polystyrene (EPS) configurations was investigated. Six concrete shear walls were designed and tested under cyclic loading to evaluate the effect of fine RAC in designing earthquake-resistant structures. Three of the six specimens were used to construct mid-rise walls with a shear-span ratio of 1.5, and the other three specimens were used to construct low-rise walls with a shear-span ratio of 0.8. The mid-rise and low-rise shear walls consisted of an ordinary recycled concrete shear wall, a composite wall with fine aggregate concrete (FAC) protective layer (EPS modules as the external insulation layer), and a composite wall with sandwiched EPS modules as the insulation layer. Several parameters obtained from the experimental results were compared and analyzed, including the load-bearing capacity, stiffness, ductility, energy dissipation, and failure characteristics of the specimens. The calculation formula of load-bearing capacity was obtained by considering the effect of FAC on composite shear walls as the protective layer. The damage process of the specimen was simulated using the ABAQUS Software, and the results agreed quite well with those obtained from the experiments. The results show that the seismic resistance behavior of the EPS module composite for shear walls performed better than ordinary recycled concrete for shear walls. Shear walls with sandwiched EPS modules had a better seismic performance than those with EPS modules lying outside. Although the FAC protective layer slightly improved the seismic performance of the structure, it undoubtedly slowed down the speed of crack formation and the stiffness degradation of the walls.

Seismic Performance of Composite Shear Walls Constructed Using Recycled Aggregate Concrete and Different Expandable Polystyrene Configurations

PubMed Central

Liu, Wenchao; Cao, Wanlin; Zhang, Jianwei; Qiao, Qiyun; Ma, Heng

2016-01-01

The seismic performance of recycled aggregate concrete (RAC) composite shear walls with different expandable polystyrene (EPS) configurations was investigated. Six concrete shear walls were designed and tested under cyclic loading to evaluate the effect of fine RAC in designing earthquake-resistant structures. Three of the six specimens were used to construct mid-rise walls with a shear-span ratio of 1.5, and the other three specimens were used to construct low-rise walls with a shear-span ratio of 0.8. The mid-rise and low-rise shear walls consisted of an ordinary recycled concrete shear wall, a composite wall with fine aggregate concrete (FAC) protective layer (EPS modules as the external insulation layer), and a composite wall with sandwiched EPS modules as the insulation layer. Several parameters obtained from the experimental results were compared and analyzed, including the load-bearing capacity, stiffness, ductility, energy dissipation, and failure characteristics of the specimens. The calculation formula of load-bearing capacity was obtained by considering the effect of FAC on composite shear walls as the protective layer. The damage process of the specimen was simulated using the ABAQUS Software, and the results agreed quite well with those obtained from the experiments. The results show that the seismic resistance behavior of the EPS module composite for shear walls performed better than ordinary recycled concrete for shear walls. Shear walls with sandwiched EPS modules had a better seismic performance than those with EPS modules lying outside. Although the FAC protective layer slightly improved the seismic performance of the structure, it undoubtedly slowed down the speed of crack formation and the stiffness degradation of the walls. PMID:28773274

The Slip Behavior and Source Parameters for Spontaneous Slip Events on Rough Faults Subjected to Slow Tectonic Loading

NASA Astrophysics Data System (ADS)

Tal, Yuval; Hager, Bradford H.

2018-02-01

We study the response to slow tectonic loading of rough faults governed by velocity weakening rate and state friction, using a 2-D plane strain model. Our numerical approach accounts for all stages in the seismic cycle, and in each simulation we model a sequence of two earthquakes or more. We focus on the global behavior of the faults and find that as the roughness amplitude, br, increases and the minimum wavelength of roughness decreases, there is a transition from seismic slip to aseismic slip, in which the load on the fault is released by more slip events but with lower slip rate, lower seismic moment per unit length, M0,1d, and lower average static stress drop on the fault, Δτt. Even larger decreases with roughness are observed when these source parameters are estimated only for the dynamic stage of the rupture. For br ≤ 0.002, the source parameters M0,1d and Δτt decrease mutually and the relationship between Δτt and the average fault strain is similar to that of a smooth fault. For faults with larger values of br that are completely ruptured during the slip events, the average fault strain generally decreases more rapidly with roughness than Δτt.

Dynamic Sliding Analysis of a Gravity Dam with Fluid-Structure-Foundation Interaction Using Finite Elements and Newmark's Sliding Block Analysis

NASA Astrophysics Data System (ADS)

Goldgruber, Markus; Shahriari, Shervin; Zenz, Gerald

2015-11-01

To reduce the natural hazard risks—due to, e.g., earthquake excitation—seismic safety assessments are carried out. Especially under severe loading, due to maximum credible or the so-called safety evaluation earthquake, critical infrastructure, as these are high dams, must not fail. However, under high loading local failure might be allowed as long as the entire structure does not collapse. Hence, for a dam, the loss of sliding stability during a short time period might be acceptable if the cumulative displacements after an event are below an acceptable value. This performance is not only valid for gravity dams but also for rock blocks as sliding is even more imminent in zones with higher seismic activity. Sliding modes cannot only occur in the dam-foundation contact, but also in sliding planes formed due to geological conditions. This work compares the qualitative possible and critical displacements for two methods, the well-known Newmark's sliding block analysis and a Fluid-Foundation-Structure Interaction simulation with the finite elements method. The results comparison of the maximum displacements at the end of the seismic event of the two methods depicts that for high friction angles, they are fairly close. For low friction angles, the results are differing more. The conclusion is that the commonly used Newmark's sliding block analysis and the finite elements simulation are only comparable for high friction angles, where this factor dominates the behaviour of the structure. Worth to mention is that the proposed simulation methods are also applicable to dynamic rock wedge problems and not only to dams.

Hydrologically-driven crustal stresses and seismicity in the New Madrid Seismic Zone.

PubMed

Craig, Timothy J; Chanard, Kristel; Calais, Eric

2017-12-15

The degree to which short-term non-tectonic processes, either natural and anthropogenic, influence the occurrence of earthquakes in active tectonic settings or 'stable' plate interiors, remains a subject of debate. Recent work in plate-boundary regions demonstrates the capacity for long-wavelength changes in continental water storage to produce observable surface deformation, induce crustal stresses and modulate seismicity rates. Here we show that a significant variation in the rate of microearthquakes in the intraplate New Madrid Seismic Zone at annual and multi-annual timescales coincides with hydrological loading in the upper Mississippi embayment. We demonstrate that this loading, which results in geodetically observed surface deformation, induces stresses within the lithosphere that, although of small amplitude, modulate the ongoing seismicity of the New Madrid region. Correspondence between surface deformation, hydrological loading and seismicity rates at both annual and multi-annual timescales indicates that seismicity variations are the direct result of elastic stresses induced by the water load.

Instant Variations in Velocity and Attenuation of Seismic Waves in a Friable Medium Under a Vibrational Dynamic Loading

NASA Astrophysics Data System (ADS)

Geza, N.; Yushin, V.

2007-12-01

Instant variations of the velocities and attenuation of seismic waves in a friable medium subjected to dynamic loading have been studied by new experimental techniques using a powerful seismic vibrator. The half-space below the operating vibrator baseplate was scanned by high-frequency elastic waves, and the recorded fluctuations were exposed to a stroboscopic analysis. It was found that the variations of seismic velocities and attenuation are synchronous with the external vibrational load but have phase shift from it. Instant variations of the seismic waves parameters depend on the magnitude and absolute value of deformation, which generally result in decreasing of the elastic-wave velocities. New experimental techniques have a high sensitivity to the dynamic disturbance in the medium and allow one to detect a weak seismic boundaries. The relaxation process after dynamic vibrational loading were investigated and the results of research are presented.

Analysis and seismic tests of composite shear walls with CFST columns and steel plate deep beams

NASA Astrophysics Data System (ADS)

Dong, Hongying; Cao, Wanlin; Wu, Haipeng; Zhang, Jianwei; Xu, Fangfang

2013-12-01

A composite shear wall concept based on concrete filled steel tube (CFST) columns and steel plate (SP) deep beams is proposed and examined in this study. The new wall is composed of three different energy dissipation elements: CFST columns; SP deep beams; and reinforced concrete (RC) strips. The RC strips are intended to allow the core structural elements — the CFST columns and SP deep beams — to work as a single structure to consume energy. Six specimens of different configurations were tested under cyclic loading. The resulting data are analyzed herein. In addition, numerical simulations of the stress and damage processes for each specimen were carried out, and simulations were completed for a range of location and span-height ratio variations for the SP beams. The simulations show good agreement with the test results. The core structure exhibits a ductile yielding mechanism characteristic of strong column-weak beam structures, hysteretic curves are plump and the composite shear wall exhibits several seismic defense lines. The deformation of the shear wall specimens with encased CFST column and SP deep beam design appears to be closer to that of entire shear walls. Establishing optimal design parameters for the configuration of SP deep beams is pivotal to the best seismic behavior of the wall. The new composite shear wall is therefore suitable for use in the seismic design of building structures.

Combined seismic plus live-load analysis of highway bridges.

DOT National Transportation Integrated Search

2011-10-01

"The combination of seismic and vehicle live loadings on bridges is an important design consideration. There are well-established design : provisions for how the individual loadings affect bridge response: structural components that carry vertical li...

Numerical Simulations of Mechanical Erosion from below by Creep on Rate-State Faults

NASA Astrophysics Data System (ADS)

Werner, M. J.; Rubin, A. M.

2012-04-01

The aim of this study is to increase our understanding of how earthquakes nucleate on frictionally-locked fault patches that are loaded by the growing stress concentrations at their boundaries due to aseismic creep. Such mechanical erosion from below of locked patches has previously been invoked by Gillard et al. (1996) to explain accelerating seismicity and increases in maximum earthquake magnitude on a strike-slip streak (a narrow ribbon of tightly clustered seismicity) in Kilauea's East rift, and it might also play a role in the loading of major locked strike-slip faults by creep from below the seismogenic zone. Gillard et al. (1996) provided simple analytical estimates of the size of and moment release within the eroding edge of the locked zone that matched the observed seismicity in Kilauea's East rift. However, an obvious, similar signal has not consistently been found before major strike-slip earthquakes. Here, we use simulations to determine to what extent the simple estimates by Gillard et al. survive a wider range of geometric configurations and slip histories. The boundary between the locked and creeping sections at the base of the seismogenic zone is modeled as a gradual, continuous transition between steady-state velocity-strengthening at greater depth to velocity-weakening surroundings at shallow depth, qualitatively consistent with laboratory estimates of the temperature dependence of (a-b). The goal is to expand the range of possible outcomes to broaden our range of expectations for the behavior of the eroding edge of the locked zones.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Kot, C.A.; Srinivasan, M.G.; Hsieh, B.J.

As part of the Phase II testing at the HDR Test Facility in Kahl/Main, FRG, two series of high-level seismic/vibrational experiments were performed. In the first of these (SHAG) a coast-down shaker, mounted on the reactor operating floor and capable of generating 1000 tonnes of force, was used to investigate full-scale structural response, soil-structure interaction (SSI), and piping/equipment response at load levels equivalent to those of a design basis earthquake. The HDR soil/structure system was tested to incipient failure exhibiting highly nonlinear response. In the load transmission from structure to piping/equipment significant response amplifications and shifts to higher frequencies occurred.more » The performance of various pipe support configurations was evaluated. This latter effort was continued in the second series of tests (SHAM), in which an in-plant piping system was investigated at simulated seismic loads (generated by two servo-hydraulic actuators each capable of generating 40 tonnes of force), that exceeded design levels manifold and resulted in considerable pipe plastification and failure of some supports (snubbers). The evaluation of six different support configurations demonstrated that proper system design (for a given spectrum) rather than number of supports or system stiffness is essential to limiting pipe stresses. Pipe strains at loads exceeding the design level eightfold were still tolerable, indicating that pipe failure even under extreme seismic loads is unlikely inspite of multiple support failures. Conservatively, an excess capacity (margin) of at least four was estimated for the piping system, and the pipe damping was found to be 4%. Comparisons of linear and nonlinear computational results with measurements showed that analytical predictions have wide scatter and do not necessarily yield conservative responses, underpredicting, in particular, peak support forces.« less

Structural evaluation of the 2736Z Building for seismic loads

DOE Office of Scientific and Technical Information (OSTI.GOV)

Giller, R.A.

The 2736Z building structure is evaluated for high-hazard loads. The 2736Z building is analyzed herein for normal and seismic loads and is found to successfully meet the guidelines of UCRL-15910 along with the related codes requirements.

Seasonal variations in shallow Alaska seismicity and stress modulation from GRACE derived hydrological loading

NASA Astrophysics Data System (ADS)

Johnson, C. W.; Fu, Y.; Burgmann, R.

2017-12-01

Shallow (≤50 km), low magnitude (M≥2.0) seismicity in southern Alaska is examined for seasonal variations during the annual hydrological cycle. The seismicity is declustered with a spatio-temporal epidemic type aftershock sequence (ETAS) model. The removal of aftershock sequences allows detailed investigation of seismicity rate changes, as water and ice loads modulate crustal stresses throughout the year. The GRACE surface loads are obtained from the JPL mass concentration blocks (mascons) global land and ocean solutions. The data product is smoothed with a 9˚ Gaussian filter and interpolated on a 25 km grid. To inform the surface loading model, the global solutions are limited to the region from -160˚ to -120˚ and 50˚ to 70˚. The stress changes are calculated using a 1D spherical layered earth model at depth intervals of 10 km from 10 - 50 km in the study region. To evaluate the induced seasonal stresses, we use >30 years of earthquake focal mechanisms to constrain the background stress field orientation and assess the stress change with respect to the principal stress orientation. The background stress field is assumed to control the preferred orientation of faulting, and stress field perturbations are expected to increase or decrease seismicity. The number of excess earthquakes is calculated with respect to the background seismicity rates. Here, we present preliminary results for the shallow seismicity variations and quantify the seasonal stresses associated with changes in hydrological loading.

Earthquake triggering by transient and static deformations

USGS Publications Warehouse

Gomberg, J.; Beeler, N.M.; Blanpied, M.L.; Bodin, P.

1998-01-01

Observational evidence for both static and transient near-field and far-field triggered seismicity are explained in terms of a frictional instability model, based on a single degree of freedom spring-slider system and rate- and state-dependent frictional constitutive equations. In this study a triggered earthquake is one whose failure time has been advanced by ??t (clock advance) due to a stress perturbation. Triggering stress perturbations considered include square-wave transients and step functions, analogous to seismic waves and coseismic static stress changes, respectively. Perturbations are superimposed on a constant background stressing rate which represents the tectonic stressing rate. The normal stress is assumed to be constant. Approximate, closed-form solutions of the rate-and-state equations are derived for these triggering and background loads, building on the work of Dieterich [1992, 1994]. These solutions can be used to simulate the effects of static and transient stresses as a function of amplitude, onset time t0, and in the case of square waves, duration. The accuracies of the approximate closed-form solutions are also evaluated with respect to the full numerical solution and t0. The approximate solutions underpredict the full solutions, although the difference decreases as t0, approaches the end of the earthquake cycle. The relationship between ??t and t0 differs for transient and static loads: a static stress step imposed late in the cycle causes less clock advance than an equal step imposed earlier, whereas a later applied transient causes greater clock advance than an equal one imposed earlier. For equal ??t, transient amplitudes must be greater than static loads by factors of several tens to hundreds depending on t0. We show that the rate-and-state model requires that the total slip at failure is a constant, regardless of the loading history. Thus a static load applied early in the cycle, or a transient applied at any time, reduces the stress at the initiation of failure, whereas static loads that are applied sufficiently late raise it. Rate-and-state friction predictions differ markedly from those based on Coulomb failure stress changes (??CFS) in which ??t equals the amplitude of the static stress change divided by the background stressing rate. The ??CFS model assumes a stress failure threshold, while the rate-and-state equations require a slip failure threshold. The complete rale-and-state equations predict larger ??t than the ??CFS model does for static stress steps at small t0, and smaller ??t than the ??CFS model for stress steps at large t0. The ??CFS model predicts nonzero ??t only for transient loads that raise the stress to failure stress levels during the transient. In contrast, the rate-and-state model predicts nonzero ??t for smaller loads, and triggered failure may occur well after the transient is finished. We consider heuristically the effects of triggering on a population of faults, as these effects might be evident in seismicity data. Triggering is manifest as an initial increase in seismicity rate that may be followed by a quiescence or by a return to the background rate. Available seismicity data are insufficient to discriminate whether triggered earthquakes are "new" or clock advanced. However, if triggering indeed results from advancing the failure time of inevitable earthquakes, then our modeling suggests that a quiescence always follows transient triggering and that the duration of increased seismicity also cannot exceed the duration of a triggering transient load. Quiescence follows static triggering only if the population of available faults is finite.

Masonry Vaults Subjected To Horizontal Loads: Experimental and Numerical Investigations to Evaluate the Effectiveness of A GFRM Reinforcement

NASA Astrophysics Data System (ADS)

Gattesco, Natalino; Boem, Ingrid

2017-10-01

The paper investigates the effectiveness of a modern reinforcement technique based on a Glass Fiber-Reinforced Mortar (GFRM) for the enhancement of the performances of existing masonry vaults subjected to horizontal seismic actions. In fact, the authors recently evidenced, through numerical simulations, that the typical simplified loading patterns generally adopted in the literature for the experimental tests, based on concentrated vertical loads at 1/4 of the span, are not reliable for such a purpose, due to an unrealistic stress distribution. Thus, experimental quasi-static cyclic tests on full-scale masonry vaults based on a specific setup, designed to apply a horizontal load pattern proportional to the mass, were performed. Three samples were tested: an unreinforced vault, a vault reinforced at the extrados and a vault reinforced at the intrados. The experimental results demonstrated the technique effectiveness in both strength and ductility. Moreover, numerical simulations were performed by adopting a simplified FE, smear-crack model, evidencing the good reliability of the prediction by comparison with the experimental results.

Seismic Performance Evaluation of Reinforced Concrete Frames Subjected to Seismic Loads

NASA Astrophysics Data System (ADS)

Zameeruddin, Mohd.; Sangle, Keshav K.

2017-06-01

Ten storied-3 bays reinforced concrete bare frame designed for gravity loads following the guidelines of IS 456 and IS 13920 for ductility is subjected to seismic loads. The seismic demands on this building were calculated by following IS 1893 for response spectra of 5% damping (for hard soil type). Plastic hinges were assigned to the beam and column at both ends to represent the failure mode, when member yields. Non-linear static (pushover) analysis was performed to evaluate the performance of the building in reference to first (ATC 40), second (FEMA 356) and next-generation (FEMA 440) performance based seismic design procedures. Base shear against top displacement curve of structure, known as pushover curve was obtained for two actions of plastic hinge behavior, force-controlled (brittle) and deformation-controlled (ductile) actions. Lateral deformation corresponding to performance point proves the building capability to sustain a certain level of seismic loads. The failure is represented by a sequence of formation of plastic hinges. Deformation-controlled action of hinges showed that building behaves like strong-column-weak-beam mechanism, whereas force-controlled action showed formation of hinges in the column. The study aims to understand the first, second and next generation performance based design procedure in prediction of actual building responses and their conservatism into the acceptance criteria.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Ghosh, A.; Hsiung, S.M.; Chowdhury, A.H.

Long-term stability of emplacement drifts and potential near-field fluid flow resulting from coupled effects are among the concerns for safe disposal of high-level nuclear waste (HLW). A number of factors can induce drift instability or change the near-field flow patterns. Repetitive seismic loads from earthquakes and thermal loads generated by the decay of emplaced waste are two significant factors. One of two key technical uncertainties (KTU) that can potentially pose a high risk of noncompliance with the performance objectives of 10 CFR Part 60 is the prediction of thermal-mechanical (including repetitive seismic load) effects on stability of emplacement drifts andmore » the engineered barrier system. The second KTU of concern is the prediction of thermal-mechanical-hydrological (including repetitive seismic load) effects on the host rock surrounding the engineered barrier system. The Rock Mechanics research project being conducted at the Center for Nuclear Waste Regulatory Analyses (CNWRA) is intended to address certain specific technical issues associated with these two KTUs. This research project has two major components: (i) seismic response of rock joints and a jointed rock mass and (ii) coupled thermal-mechanical-hydrological (TMH) response of a jointed rock mass surrounding the engineered barrier system (EBS). This final report summarizes the research activities concerned with the repetitive seismic load aspect of both these KTUs.« less

Mechanical Erosion of the Seismogenic Zone by Creep from below on Rate-State Faults

NASA Astrophysics Data System (ADS)

Werner, M. J.; Rubin, A. M.

2012-12-01

The aim of this study is to increase our understanding of how earthquakes nucleate on frictionally-locked fault patches that are loaded by the growing stress concentrations at their boundaries due to aseismic creep. Such mechanical erosion from below of locked patches has previously been invoked by Gillard et al. (1996) to explain accelerating seismicity and increases in maximum earthquake magnitude on a strike-slip streak (a narrow ribbon of tightly clustered seismicity) in Kilauea's East rift, and it might also play a role in the loading of major locked strike-slip faults by creep from below the seismogenic zone. Gillard et al. (1996) provided simple analytical estimates of the size of and of the time-dependence of the moment release within the eroding edge of the locked zone that matched observed seismicity in Kilauea's East rift. However, an obvious, similar signal has not consistently been found before major strike-slip earthquakes on presumably analogous faults. Here, we use simulations to determine to what extent the simple estimates by Gillard et al. survive a wider range of geometric configurations, friction laws and slip histories. We model the boundary between the locked and creeping sections at the base of the seismogenic zone as a transition between steady-state velocity-strengthening behavior at greater depth to velocity-weakening at shallow depth, qualitatively consistent with laboratory estimates of the temperature dependence of (a-b). As we increase the ratio of the size H of the velocity-weakening region over the nucleation length scale L∞ [Rubin & Ampuero, 2005] from tens to more appropriate hundreds and thousands in our 1D model, we observe an increasing number of creep fronts that march into the seismogenic zone from the weakening/strengthening transition between surface-rupturing seismic events. For moderate to large ratios of H/L∞ , these fronts begin to appear while afterslip from large events is still propagating down-dip below the transition and reloading the seismogenic zone near the transition. Some creep fronts slow down and fail, especially during very early deep afterslip, while later creep fronts tend to accelerate to seismic speeds (modeled using radiation-damping) after traveling a distance above the transition that appears to be controlled by L∞ . The ensuing earthquakes propagate updip and reach only a fraction of the seismogenic height, until eventually a surface-rupturing event occurs. As we increase a/b towards 1 in the velocity-weakening zone, we observe some large events nucleating well above the transition at stress concentrations left over from earlier seismic events that did not reach the surface. In our simulations thus far, total moment release rates over the earthquake cycle tend to decay with afterslip, occasionally punctuated by the smaller earthquakes, which decrease in frequency over the cycle and seem mostly driven by deep afterslip. It thus appears possible that, at least on homogeneous (aging) rate-state faults, post-seismic afterslip is masking (some of) the seismicity patterns expected from Gillard et al.'s analysis, although a wider range of simulations is warranted.

Liquid Rocket Booster (LRB) for the Space Transportation System (STS) systems study. Appendix E: Pressure-fed booster test bed for the liquid rocket booster study

NASA Technical Reports Server (NTRS)

1989-01-01

The stress analysis/structural design of the Pressure-Fed Booster Engine Test Bed using the existing F-1 Test Facility Test Stand at Huntsville, Alabama is described. The analysis has been coded and set up for solution on NASTRAN. A separate stress program was established to take the NASTRAN output and perform stress checks on the members. Joint checks and other necessary additional checks were performed by hand. The notes include a brief description of other programs which assist in reproducing and reviewing the NASTRAN results. The redesign of the test stand members and the stress analysis was performed per the A.I.S.C. Code. Loads on the stand consist of the loaded run tanks; wind loads; seismic loads; live loads consisting of snow and ice: live and dead loads of steel; and loaded pressurant bottle. In combining loads, wind loads and seismic loads were each combined with full live loads. Wind and seismic loads were not combined. No one third increase in allowables was taken for the environmental loads except at decks 147 and 214, where the increase was used when considering the stay rods, brackets and stay beams. Wind and seismic loads were considered from each of the four coordinate directions (i.e. N,S,E,W) to give eight basic conditions. The analysis was run with the pressurant tank mounted at level 125. One seismic condition was also run with the tank mounted at levels 169 and 214. No failures were noted with mounting at level 169, but extensive deck failure with mounting at level 214 (the loadsets used are included on the tape, but no detailed results are included in the package). Decking support beams at levels 147 and 214 are not included in the model. The stress program thus does not reduce strut lengths to the length between support beams (the struts are attached to the beams at intersection points) and gives stress ratios larger than one for some of the struts. The affected members were therefore checked by hand.

Cyclic performance of concrete-filled steel batten built-up columns

NASA Astrophysics Data System (ADS)

Razzaghi, M. S.; Khalkhaliha, M.; Aziminejad, A.

2016-03-01

Steel built-up batten columns are common types of columns in Iran and some other parts of the world. They are economic and have acceptable performance due to gravity loads. Although several researches have been conducted on the behavior of the batten columns under axial loads, there are few available articles about their seismic performance. Experience of the past earthquakes, particularly the 2003 Bam earthquake in Iran, revealed that these structural members are seismically vulnerable. Thus, investigation on seismic performance of steel batten columns due to seismic loads and providing a method for retrofitting them are important task in seismic-prone areas. This study aims to investigate the behavior of concrete-filled batten columns due to combined axial and lateral loads. To this end, nonlinear static analyses were performed using ANSYS software. Herein, the behaviors of the steel batten columns with and without concrete core were compared. The results of this study showed that concrete-filled steel batten columns, particularly those filled with high-strength concrete, may cause significant increases in energy absorption and capacity of the columns. Furthermore, concrete core may improve post-buckling behavior of steel batten columns.

Centrifuge modeling of rocking-isolated inelastic RC bridge piers

PubMed Central

Loli, Marianna; Knappett, Jonathan A; Brown, Michael J; Anastasopoulos, Ioannis; Gazetas, George

2014-01-01

Experimental proof is provided of an unconventional seismic design concept, which is based on deliberately underdesigning shallow foundations to promote intense rocking oscillations and thereby to dramatically improve the seismic resilience of structures. Termed rocking isolation, this new seismic design philosophy is investigated through a series of dynamic centrifuge experiments on properly scaled models of a modern reinforced concrete (RC) bridge pier. The experimental method reproduces the nonlinear and inelastic response of both the soil-footing interface and the structure. To this end, a novel scale model RC (1:50 scale) that simulates reasonably well the elastic response and the failure of prototype RC elements is utilized, along with realistic representation of the soil behavior in a geotechnical centrifuge. A variety of seismic ground motions are considered as excitations. They result in consistent demonstrably beneficial performance of the rocking-isolated pier in comparison with the one designed conventionally. Seismic demand is reduced in terms of both inertial load and deck drift. Furthermore, foundation uplifting has a self-centering potential, whereas soil yielding is shown to provide a particularly effective energy dissipation mechanism, exhibiting significant resistance to cumulative damage. Thanks to such mechanisms, the rocking pier survived, with no signs of structural distress, a deleterious sequence of seismic motions that caused collapse of the conventionally designed pier. © 2014 The Authors Earthquake Engineering & Structural Dynamics Published by John Wiley & Sons Ltd. PMID:26300573

Centrifuge modeling of rocking-isolated inelastic RC bridge piers.

PubMed

Loli, Marianna; Knappett, Jonathan A; Brown, Michael J; Anastasopoulos, Ioannis; Gazetas, George

2014-12-01

Experimental proof is provided of an unconventional seismic design concept, which is based on deliberately underdesigning shallow foundations to promote intense rocking oscillations and thereby to dramatically improve the seismic resilience of structures. Termed rocking isolation , this new seismic design philosophy is investigated through a series of dynamic centrifuge experiments on properly scaled models of a modern reinforced concrete (RC) bridge pier. The experimental method reproduces the nonlinear and inelastic response of both the soil-footing interface and the structure. To this end, a novel scale model RC (1:50 scale) that simulates reasonably well the elastic response and the failure of prototype RC elements is utilized, along with realistic representation of the soil behavior in a geotechnical centrifuge. A variety of seismic ground motions are considered as excitations. They result in consistent demonstrably beneficial performance of the rocking-isolated pier in comparison with the one designed conventionally. Seismic demand is reduced in terms of both inertial load and deck drift. Furthermore, foundation uplifting has a self-centering potential, whereas soil yielding is shown to provide a particularly effective energy dissipation mechanism, exhibiting significant resistance to cumulative damage. Thanks to such mechanisms, the rocking pier survived, with no signs of structural distress, a deleterious sequence of seismic motions that caused collapse of the conventionally designed pier. © 2014 The Authors Published by John Wiley & Sons Ltd.

Comparison among different retrofitting strategies for the vulnerability reduction of masonry bell towers

NASA Astrophysics Data System (ADS)

Milani, Gabriele; Shehu, Rafael; Valente, Marco

2017-11-01

This paper investigates the effectiveness of reducing the seismic vulnerability of masonry towers by means of innovative and traditional strengthening techniques. The followed strategy for providing the optimal retrofitting for masonry towers subjected to seismic risk relies on preventing active failure mechanisms. These vulnerable mechanisms are pre-assigned failure patterns based on the crack patterns experienced during the past seismic events. An upper bound limit analysis strategy is found suitable to be applied for simplified tower models in their present state and the proposed retrofitted ones. Taking into consideration the variability of geometrical features and the uncertainty of the strengthening techniques, Monte Carlo simulations are implemented into the limit analysis. In this framework a wide range of idealized cases are covered by the conducted analyses. The retrofitting strategies aim to increase the shear strength and the overturning load carrying capacity in order to reduce vulnerability. This methodology gives the possibility to use different materials which can fulfill the structural implementability requirements.

Experimental investigation of damage behavior of RC frame members including non-seismically designed columns

NASA Astrophysics Data System (ADS)

Chen, Linzhi; Lu, Xilin; Jiang, Huanjun; Zheng, Jianbo

2009-06-01

Reinforced concrete (RC) frame structures are one of the mostly common used structural systems, and their seismic performance is largely determined by the performance of columns and beams. This paper describes horizontal cyclic loading tests of ten column and three beam specimens, some of which were designed according to the current seismic design code and others were designed according to the early non-seismic Chinese design code, aiming at reporting the behavior of the damaged or collapsed RC frame strctures observed during the Wenchuan earthquake. The effects of axial load ratio, shear span ratio, and transverse and longitudinal reinforcement ratio on hysteresis behavior, ductility and damage progress were incorporated in the experimental study. Test results indicate that the non-seismically designed columns show premature shear failure, and yield larger maximum residual crack widths and more concrete spalling than the seismically designed columns. In addition, longitudinal steel reinforcement rebars were severely buckled. The axial load ratio and shear span ratio proved to be the most important factors affecting the ductility, crack opening width and closing ability, while the longitudinal reinforcement ratio had only a minor effect on column ductility, but exhibited more influence on beam ductility. Finally, the transverse reinforcement ratio did not influence the maximum residual crack width and closing ability of the seismically designed columns.

Earthquake stress drop and laboratory-inferred interseismic strength recovery

USGS Publications Warehouse

Beeler, N.M.; Hickman, S.H.; Wong, T.-F.

2001-01-01

We determine the scaling relationships between earthquake stress drop and recurrence interval tr that are implied by laboratory-measured fault strength. We assume that repeating earthquakes can be simulated by stick-slip sliding using a spring and slider block model. Simulations with static/kinetic strength, time-dependent strength, and rate- and state-variable-dependent strength indicate that the relationship between loading velocity and recurrence interval can be adequately described by the power law VL ??? trn, where n=-1. Deviations from n=-1 arise from second order effects on strength, with n>-1 corresponding to apparent time-dependent strengthening and n<-1 corresponding to weakening. Simulations with rate and state-variable equations show that dynamic shear stress drop ????d scales with recurrence as d????d/dlntr ??? ??e(b-a), where ??e is the effective normal stress, ??=??/??e, and (a-b)=d??ss/dlnV is the steady-state slip rate dependence of strength. In addition, accounting for seismic energy radiation, we suggest that the static shear stress drop ????s scales as d????s/dlntr ??? ??e(1+??)(b-a), where ?? is the fractional overshoot. The variation of ????s with lntr for earthquake stress drops is somewhat larger than implied by room temperature laboratory values of ?? and b-a. However, the uncertainty associated with the seismic data is large and the discrepancy between the seismic observations and the rate of strengthening predicted by room temperature experiments is less than an order of magnitude. Copyright 2001 by the American Geophysical Union.

Transient triggering of near and distant earthquakes

USGS Publications Warehouse

Gomberg, J.; Blanpied, M.L.; Beeler, N.M.

1997-01-01

We demonstrate qualitatively that frictional instability theory provides a context for understanding how earthquakes may be triggered by transient loads associated with seismic waves from near and distance earthquakes. We assume that earthquake triggering is a stick-slip process and test two hypotheses about the effect of transients on the timing of instabilities using a simple spring-slider model and a rate- and state-dependent friction constitutive law. A critical triggering threshold is implicit in such a model formulation. Our first hypothesis is that transient loads lead to clock advances; i.e., transients hasten the time of earthquakes that would have happened eventually due to constant background loading alone. Modeling results demonstrate that transient loads do lead to clock advances and that the triggered instabilities may occur after the transient has ceased (i.e., triggering may be delayed). These simple "clock-advance" models predict complex relationships between the triggering delay, the clock advance, and the transient characteristics. The triggering delay and the degree of clock advance both depend nonlinearly on when in the earthquake cycle the transient load is applied. This implies that the stress required to bring about failure does not depend linearly on loading time, even when the fault is loaded at a constant rate. The timing of instability also depends nonlinearly on the transient loading rate, faster rates more rapidly hastening instability. This implies that higher-frequency and/or longer-duration seismic waves should increase the amount of clock advance. These modeling results and simple calculations suggest that near (tens of kilometers) small/moderate earthquakes and remote (thousands of kilometers) earthquakes with magnitudes 2 to 3 units larger may be equally effective at triggering seismicity. Our second hypothesis is that some triggered seismicity represents earthquakes that would not have happened without the transient load (i.e., accumulated strain energy would have been relieved via other mechanisms). We test this using two "new-seismicity" models that (1) are inherently unstable but slide at steady-state conditions under the background load and (2) are conditionally stable such that instability occurs only for sufficiently large perturbations. For the new-seismicity models, very small-amplitude transients trigger instability relative to the clock-advance models. The unstable steady-state models predict that the triggering delay depends inversely and nonlinearly on the transient amplitude (as in the clock-advance models). We were unable to generate delayed triggering with conditionally stable models. For both new-seismicity models, the potential for triggering is independent of when the transient load is applied or, equivalently, of the prestress (unlike in the clock-advance models). In these models, a critical triggering threshold appears to be inversely proportional to frequency. Further advancement of our understanding will require more sophisticated, quantitative models and observations that distinguish between our qualitative, yet distinctly different, model predictions.

Effect of URM infills on seismic vulnerability of Indian code designed RC frame buildings

NASA Astrophysics Data System (ADS)

Haldar, Putul; Singh, Yogendra; Paul, D. K.

2012-03-01

Unreinforced Masonry (URM) is the most common partitioning material in framed buildings in India and many other countries. Although it is well-known that under lateral loading the behavior and modes of failure of the frame buildings change significantly due to infill-frame interaction, the general design practice is to treat infills as nonstructural elements and their stiffness, strength and interaction with the frame is often ignored, primarily because of difficulties in simulation and lack of modeling guidelines in design codes. The Indian Standard, like many other national codes, does not provide explicit insight into the anticipated performance and associated vulnerability of infilled frames. This paper presents an analytical study on the seismic performance and fragility analysis of Indian code-designed RC frame buildings with and without URM infills. Infills are modeled as diagonal struts as per ASCE 41 guidelines and various modes of failure are considered. HAZUS methodology along with nonlinear static analysis is used to compare the seismic vulnerability of bare and infilled frames. The comparative study suggests that URM infills result in a significant increase in the seismic vulnerability of RC frames and their effect needs to be properly incorporated in design codes.

Reliability assessment of an OVH HV power line truss transmission tower subjected to seismic loading

NASA Astrophysics Data System (ADS)

Winkelmann, Karol; Jakubowska, Patrycja; Soltysik, Barbara

2017-03-01

The study focuses on the reliability of a transmission tower OS24 ON150 + 10, an element of an OVH HV power line, under seismic loading. In order to describe the seismic force, the real-life recording of the horizontal component of the El Centro earthquake was adopted. The amplitude and the period of this excitation are assumed random, their variation is described by Weibull distribution. The possible space state of the phenomenon is given in the form of a structural response surface (RSM methodology), approximated by an ANOVA table with directional sampling (DS) points. Four design limit states are considered: stress limit criterion for a natural load combination, criterion for an accidental combination (one-sided cable snap), vertical and horizontal translation criteria. According to these cases the HLRF reliability index β is used for structural safety assessment. The RSM approach is well suited for the analysis - it is numerically efficient, not excessively time consuming, indicating a high confidence level. Given the problem conditions, the seismic excitation is shown the sufficient trigger to the loss of load-bearing capacity or stability of the tower.

Experimental Study on the Seismic Performance of Recycled Concrete Brick Walls Embedded with Vertical Reinforcement.

PubMed

Cao, Wanlin; Zhang, Yongbo; Dong, Hongying; Zhou, Zhongyi; Qiao, Qiyun

2014-08-19

Recycled concrete brick (RCB) is manufactured by recycled aggregate processed from discarded concrete blocks arising from the demolishing of existing buildings. This paper presents research on the seismic performance of RCB masonry walls to assess the applicability of RCB for use in rural low-rise constructions. The seismic performance of a masonry wall is closely related to the vertical load applied to the wall. Thus, the compressive performance of RCB masonry was investigated firstly by constructing and testing eighteen RCB masonry compressive specimens with different mortar strengths. The load-bearing capacity, deformation and failure characteristic were analyzed, as well. Then, a quasi-static test was carried out to study the seismic behavior of RCB walls by eight RCB masonry walls subjected to an axial compressive load and a reversed cyclic lateral load. Based on the test results, equations for predicting the compressive strength of RCB masonry and the lateral ultimate strength of an RCB masonry wall were proposed. Experimental values were found to be in good agreement with the predicted values. Meanwhile, finite element analysis (FEA) and parametric analysis of the RCB walls were carried out using ABAQUS software. The elastic-plastic deformation characteristics and the lateral load-displacement relations were studied.

Experimental Study on the Seismic Performance of Recycled Concrete Brick Walls Embedded with Vertical Reinforcement

PubMed Central

Cao, Wanlin; Zhang, Yongbo; Dong, Hongying; Zhou, Zhongyi; Qiao, Qiyun

2014-01-01

Recycled concrete brick (RCB) is manufactured by recycled aggregate processed from discarded concrete blocks arising from the demolishing of existing buildings. This paper presents research on the seismic performance of RCB masonry walls to assess the applicability of RCB for use in rural low-rise constructions. The seismic performance of a masonry wall is closely related to the vertical load applied to the wall. Thus, the compressive performance of RCB masonry was investigated firstly by constructing and testing eighteen RCB masonry compressive specimens with different mortar strengths. The load-bearing capacity, deformation and failure characteristic were analyzed, as well. Then, a quasi-static test was carried out to study the seismic behavior of RCB walls by eight RCB masonry walls subjected to an axial compressive load and a reversed cyclic lateral load. Based on the test results, equations for predicting the compressive strength of RCB masonry and the lateral ultimate strength of an RCB masonry wall were proposed. Experimental values were found to be in good agreement with the predicted values. Meanwhile, finite element analysis (FEA) and parametric analysis of the RCB walls were carried out using ABAQUS software. The elastic-plastic deformation characteristics and the lateral load-displacement relations were studied. PMID:28788170

Medium effect on the characteristics of the coupled seismic and electromagnetic signals.

PubMed

Huang, Qinghua; Ren, Hengxin; Zhang, Dan; Chen, Y John

2015-01-01

Recently developed numerical simulation technique can simulate the coupled seismic and electromagnetic signals for a double couple point source or a finite fault planar source. Besides the source effect, the simulation results showed that both medium structure and medium property could affect the coupled seismic and electromagnetic signals. The waveform of coupled signals for a layered structure is more complicated than that for a simple uniform structure. Different from the seismic signals, the electromagnetic signals are sensitive to the medium properties such as fluid salinity and fluid viscosity. Therefore, the co-seismic electromagnetic signals may be more informative than seismic signals.

Medium effect on the characteristics of the coupled seismic and electromagnetic signals

PubMed Central

HUANG, Qinghua; REN, Hengxin; ZHANG, Dan; CHEN, Y. John

2015-01-01

Recently developed numerical simulation technique can simulate the coupled seismic and electromagnetic signals for a double couple point source or a finite fault planar source. Besides the source effect, the simulation results showed that both medium structure and medium property could affect the coupled seismic and electromagnetic signals. The waveform of coupled signals for a layered structure is more complicated than that for a simple uniform structure. Different from the seismic signals, the electromagnetic signals are sensitive to the medium properties such as fluid salinity and fluid viscosity. Therefore, the co-seismic electromagnetic signals may be more informative than seismic signals. PMID:25743062

Considering the ranges of uncertainties in the New Probabilistic Seismic Hazard Assessment of Germany - Version 2016

NASA Astrophysics Data System (ADS)

Grunthal, Gottfried; Stromeyer, Dietrich; Bosse, Christian; Cotton, Fabrice; Bindi, Dino

2017-04-01

The seismic load parameters for the upcoming National Annex to the Eurocode 8 result from the reassessment of the seismic hazard supported by the German Institution for Civil Engineering . This 2016 version of hazard assessment for Germany as target area was based on a comprehensive involvement of all accessible uncertainties in models and parameters into the approach and the provision of a rational framework for facilitating the uncertainties in a transparent way. The developed seismic hazard model represents significant improvements; i.e. it is based on updated and extended databases, comprehensive ranges of models, robust methods and a selection of a set of ground motion prediction equations of their latest generation. The output specifications were designed according to the user oriented needs as suggested by two review teams supervising the entire project. In particular, seismic load parameters were calculated for rock conditions with a vS30 of 800 ms-1 for three hazard levels (10%, 5% and 2% probability of occurrence or exceedance within 50 years) in form of, e.g., uniform hazard spectra (UHS) based on 19 sprectral periods in the range of 0.01 - 3s, seismic hazard maps for spectral response accelerations for different spectral periods or for macroseismic intensities. The developed hazard model consists of a logic tree with 4040 end branches and essential innovations employed to capture epistemic uncertainties and aleatory variabilities. The computation scheme enables the sound calculation of the mean and any quantile of required seismic load parameters. Mean, median and 84th percentiles of load parameters were provided together with the full calculation model to clearly illustrate the uncertainties of such a probabilistic assessment for a region of a low-to-moderate level of seismicity. The regional variations of these uncertainties (e.g. ratios between the mean and median hazard estimations) were analyzed and discussed.

Tectonic state: its significance and characterization in the assessment of seismic effects associated with reservoir impounding

USGS Publications Warehouse

Castle, R.O.; Clark, M.M.; Grantz, A.; Savage, J.C.

1980-01-01

Any analysis of seismicity associated with the filling of large reservoirs requires an evaluation of the natural tectonic state in order to determine whether impoundment is the basic source, a mechanically unrelated companion feature, or a triggering stimulus of the observed seismicity. Several arguments indicate that the associated seismicity is usually a triggered effect. Among the elements of tectonic state considered here (existing fractures, accumulated elastic strain, and deformational style), deformational style is especially critical in forecasting the occurrence of impoundment-induced seismicity. The observational evidence indicates that seismicity associated with impounding generally occurs in areas that combine steeply dipping faults, relatively high strain rates, and either extensional or horizontal-shear strain. Simple physical arguments suggest: (1) that increased fluid pressures resulting from increased reservoir head should enhance the likelihood of seismic activity, whatever the tectonic environment; (2) that stress changes resulting from surface loading may increase the likelihood of crustal failure in areas of normal and transcurrent faulting, whereas they generally inhibit failure in areas of thrust faulting. Comparisons with other earthquake-producing artificial and natural processes (underground explosions, fluid injection, underground mining, fluid extraction, volcanic emissions) indicate that reservoir loading may similarly modify the natural tectonic state. Subsurface loading resulting from fluid extraction may be a particularly close analogue of reservoir loading; "seismotectonic" events associated with fluid extraction have been recognized in both seismically active and otherwise aseismic regions. Because the historic record of seismicity and surface faulting commonly is short in comparison with recurrence intervals of earthquake and fault-slip events, tectonic state is most reliably appraised through combined studies of historic seismicity and faulting, instrumentally measured strain, and the geological record, especially that of the Quaternary. Experience in California and elsewhere demonstrates that the character and activity of recognized faults can be assessed by means of: instrumental earthquake investigations, repeated geodetic measurements, written history, archeological studies, fault topography, and local stratigraphic relations. Where faults are less easily distinguished, appraisals of tectonic state may be based on both the regional seismicity and the regional history of vertical movement as shown by: repeated levelling and sea-level measurements, written history, archeologic investigations, terrace and shoreline deformation, and denudation and sedimentation studies. ?? 1980.

The Second International Piping Integrity Research Group (IPIRG-2) program. Final report, October 1991--April 1996

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hopper, A.; Wilowski, G.; Scott, P.

1997-03-01

The IPIRG-2 program was an international group program managed by the US NRC and funded by organizations from 15 nations. The emphasis of the IPIRG-2 program was the development of data to verify fracture analyses for cracked pipes and fittings subjected to dynamic/cyclic load histories typical of seismic events. The scope included: (1) the study of more complex dynamic/cyclic load histories, i.e., multi-frequency, variable amplitude, simulated seismic excitations, than those considered in the IPIRG-1 program, (2) crack sizes more typical of those considered in Leak-Before-Break (LBB) and in-service flaw evaluations, (3) through-wall-cracked pipe experiments which can be used to validatemore » LBB-type fracture analyses, (4) cracks in and around pipe fittings, such as elbows, and (5) laboratory specimen and separate effect pipe experiments to provide better insight into the effects of dynamic and cyclic load histories. Also undertaken were an uncertainty analysis to identify the issues most important for LBB or in-service flaw evaluations, updating computer codes and databases, the development and conduct of a series of round-robin analyses, and analyst`s group meetings to provide a forum for nuclear piping experts from around the world to exchange information on the subject of pipe fracture technology. 17 refs., 104 figs., 41 tabs.« less

Effects of Irregular Bridge Columns and Feasibility of Seismic Regularity

NASA Astrophysics Data System (ADS)

Thomas, Abey E.

2018-05-01

Bridges with unequal column height is one of the main irregularities in bridge design particularly while negotiating steep valleys, making the bridges vulnerable to seismic action. The desirable behaviour of bridge columns towards seismic loading is that, they should perform in a regular fashion, i.e. the capacity of each column should be utilized evenly. But, this type of behaviour is often missing when the column heights are unequal along the length of the bridge, allowing short columns to bear the maximum lateral load. In the present study, the effects of unequal column height on the global seismic performance of bridges are studied using pushover analysis. Codes such as CalTrans (Engineering service center, earthquake engineering branch, 2013) and EC-8 (EN 1998-2: design of structures for earthquake resistance. Part 2: bridges, European Committee for Standardization, Brussels, 2005) suggests seismic regularity criterion for achieving regular seismic performance level at all the bridge columns. The feasibility of adopting these seismic regularity criterions along with those mentioned in literatures will be assessed for bridges designed as per the Indian Standards in the present study.

3D Numerical Simulation on the Sloshing Waves Excited by the Seismic Shacking

NASA Astrophysics Data System (ADS)

Zhang, Lin; Wu, Tso-Ren

2016-04-01

In the event of 2015 Nepal earthquake, a video clip broadcasted worldwide showed a violent water spilling in a hotel swimming pool. This sloshing phenomenon indicates a potential water loss in the sensitive facilities, e.g. the spent fuel pools in nuclear power plant, has to be taken into account carefully under the consideration of seismic-induced ground acceleration. In the previous studies, the simulation of sloshing mainly focused on the pressure force on the structure by using a simplified Spring-Mass Method developed in the field of solid mechanics. However, restricted by the assumptions of plane water surface and limited wave height, significant error will be made in evaluating the amount of water loss in the tank. In this paper, the computational fluid dynamical model, Splash3D, was adopted for studying the sloshing problem accurately. Splash3D solved 3D Navier-Stokes Equation directly with Large-Eddy Simulation (LES) turbulent closure. The Volume-of-fluid (VOF) method with piecewise linear interface calculation (PLIC) was used to track the complex breaking water surface. The time series acceleration of a design seismic was loaded to excite the water. With few restrictions from the assumptions, the accuracy of the simulation results were improved dramatically. A series model validations were conducted by compared to a 2D theoretical solution, and a 3D experimental data. Good comparisons can be seen. After the validation, we performed the simulation for considering a sloshing case in a rectangular water tank with a dimension of 12 m long, 8 m wide, 8 m deep, which contained water with 7 m in depth. The seismic movement was imported by considering time-series acceleration in three dimensions, which were about 0.5 g to 1.2 g in the horizontal directions, and 0.3 g to 1 g in the vertical direction. We focused the discussions on the kinematics of the water surface, wave breaking, velocity field, pressure field, water force on the side walls, and, most importantly, the amount of water loosed in the event. The simulated water movement excited by the seismic acceleration was visually similar to the video clip mentioned before. From the simulation results, we observed that the water was mainly leaked at the corner of the water tank with a nonlinear curve of the free-surface. This phenomenon can't be found in the conventional studies with acceleration in a sole direction. We also studied the effect from a porous body placed on the lower part of the tank. Detailed results and discussion will be presented in the full paper. Keywords Sloshing, Splash3D, LES, Breaking waves, VOF, spent fuel pool, Nuclear power plant

A Study on Seismic Hazard Evaluation at the Nagaoka CO2 Storage Site, Japan

NASA Astrophysics Data System (ADS)

Horikawa, S.

2015-12-01

RITE carried out the first Japanese pilot-scale CO2 sequestration project from July, 2003 to January, 2005 in Nagaoka City.Supercritical CO2 was injected into an onshore saline aquifer at a depth of 1,100m. CO2 was injected at a rate of 10,400 tonnes. 'Mid Niigata Prefecture Earthquake in 2004' (Mw6.6) and 'The Niigataken Chuetsu-oki Earthquake in 2007' (Mw6.6) occurred during the CO2 injection-test and after the completion of injection-test. Japan is one of the world's major countries with frequent earthquakes.This paper presents a result of seismic response analysis, and reports of seismic hazard evaluation of a reservoir and a caprock. In advance of dynamic response analysis, the earthquake motion recorded on the earth surface assumed the horizontally layer model, and set up the input wave from a basement layer by SHAKE ( = One-Dimensional Seismic Response Analysis). This wave was inputted into the analysis model and the equation of motion was solved using the direct integral calculus by Newmark Beta Method. In Seismic Response Analysis, authors have used Multiple Yield Model (MYM, Iwata, et al., 2013), which can respond also to complicated geological structure. The intensity deformation property of the foundation added the offloading characteristic to the composition rule of Duncan-Chang model in consideration of confining stress dependency, and used for and carried out the nonlinear repetition model. And the deformation characteristic which made it depend on confining stress with the cyclic loadings and un-loadings, and combined Mohr-Coulomb's law as a strength characteristic.The maximum dynamic shearing strain of caprock was generated about 1.1E-04 after the end of an earthquake. Although the dynamic safety factor was 1.925 on the beginning, after the end of an earthquake fell 0.05 point. The dynamic safety factor of reservoir fell to 1.20 from 1.29. As a result of CO2 migration monitoring by the seismic cross-hole tomography, CO2 has stopped in the reservoir through two earthquakes till the present after injection, and the leak is not accepted till the present. By the result of seismic response simulation, it turned out that the stability of the foundation is not spoiled after the earthquake.

Observation and Simulation of Microseisms Offshore Ireland

NASA Astrophysics Data System (ADS)

Le Pape, Florian; Bean, Chris; Craig, David; Jousset, Philippe; Donne, Sarah; Möllhoff, Martin

2017-04-01

Although more and more used in seismic imagery, ocean induced ambient seismic noise is still not so well understood, particularly how the signal propagates from ocean to land. Between January and September 2016, 10 broadband Ocean Bottom Seismometers (OBSs) stations, including acoustic sensors (hydrophone), were deployed across the shelf offshore Donegal and out into the Rockall Trough. The preliminary results show spatial and temporal variability in the ocean generated seismic noise which holds information about changes in the generation source process, including meteorological information, but also in the geological structure. In addition to the collected OBS data, numerical simulations of acoustic/seismic wave propagation are also considered in order to study the spatio-temporal variation of the broadband acoustic wavefield and its connection with the measured seismic wavefield in the region. Combination of observations and simulations appears significant to better understand what control the acoustic/seismic coupling at the sea floor as well as the effect of the water column and sediments thickness on signal propagation. Ocean generated seismic ambient noise recorded at the seafloor appears to behave differently in deep and shallow water and 3D simulations of acoustic/seismic wave propagation look particularly promising for reconciling deep ocean, shelf and land seismic observations.

Analysis and Simulation of Far-Field Seismic Data from the Source Physics Experiment

DTIC Science & Technology

2012-09-01

ANALYSIS AND SIMULATION OF FAR-FIELD SEISMIC DATA FROM THE SOURCE PHYSICS EXPERIMENT Arben Pitarka, Robert J. Mellors, Arthur J. Rodgers, Sean...Security Site (NNSS) provides new data for investigating the excitation and propagation of seismic waves generated by buried explosions. A particular... seismic model. The 3D seismic model includes surface topography. It is based on regional geological data, with material properties constrained by shallow

Seismic wave propagation in granular media

NASA Astrophysics Data System (ADS)

Tancredi, Gonzalo; López, Francisco; Gallot, Thomas; Ginares, Alejandro; Ortega, Henry; Sanchís, Johnny; Agriela, Adrián; Weatherley, Dion

2016-10-01

Asteroids and small bodies of the Solar System are thought to be agglomerates of irregular boulders, therefore cataloged as granular media. It is a consensus that many asteroids might be considered as rubble or gravel piles.Impacts on their surface could produce seismic waves which propagate in the interior of these bodies, thus causing modifications in the internal distribution of rocks and ejections of particles and dust, resulting in a cometary-type comma.We present experimental and numerical results on the study of propagation of impact-induced seismic waves in granular media, with special focus on behavior changes by increasing compression.For the experiment, we use an acrylic box filled with granular materials such as sand, gravel and glass spheres. Pressure inside the box is controlled by a movable side wall and measured with sensors. Impacts are created on the upper face of the box through a hole, ranging from free-falling spheres to gunshots. We put high-speed cameras outside the box to record the impact as well as piezoelectic sensors and accelerometers placed at several depths in the granular material to detect the seismic wave.Numerical simulations are performed with ESyS-Particle, a software that implements the Discrete Element Method. The experimental setting is reproduced in the numerical simulations using both individual spherical particles and agglomerates of spherical particles shaped as irregular boulders, according to rock models obtained with a 3D scanner. The numerical experiments also reproduces the force loading on one of the wall to vary the pressure inside the box.We are interested in the velocity, attenuation and energy transmission of the waves. These quantities are measured in the experiments and in the simulations. We study the dependance of these three parameters with characteristics like: impact speed, properties of the target material and the pressure in the media.These results are relevant to understand the outcomes of impacts in rubble/gravel pile asteroids.

Stress modulation of earthquakes: A study of long and short period stress perturbations and the crustal response

NASA Astrophysics Data System (ADS)

Johnson, Christopher W.

Decomposing fault mechanical processes advances our understanding of active fault systems and properties of the lithosphere, thereby increasing the effectiveness of seismic hazard assessment and preventative measures implemented in urban centers. Along plate boundaries earthquakes are inevitable as tectonic forces reshape the Earth's surface. Earthquakes, faulting, and surface displacements are related systems that require multidisciplinary approaches to characterize deformation in the lithosphere. Modern geodetic instrumentation can resolve displacements to millimeter precision and provide valuable insight into secular deformation in near real-time. The expansion of permanent seismic networks as well as temporary deployments allow unprecedented detection of microseismic events that image fault interfaces and fracture networks in the crust. The research presented in this dissertation is at the intersection of seismology and geodesy to study the Earth's response to transient deformation and explores research questions focusing on earthquake triggering, induced seismicity, and seasonal loading while utilizing seismic data, geodetic data, and modeling tools. The focus is to quantify stress changes in the crust, explore seismicity rate variations and migration patterns, and model crustal deformation in order to characterize the evolving state of stress on faults and the migration of fluids in the crust. The collection of problems investigated all investigate the question: Why do earthquakes nucleate following a low magnitude stress perturbation? Answers to this question are fundamental to understanding the time dependent failure processes of the lithosphere. Dynamic triggering is the interaction of faults and triggering of earthquakes represents stress transferring from one system to another, at both local and remote distances [Freed, 2005]. The passage of teleseismic surface waves from the largest earthquakes produce dynamic stress fields and provides a natural laboratory to explore the causal relationship between low-amplitude stress changes and dynamically triggered events. Interestingly, observations of dynamically triggered M≥5.5 earthquakes are absent in the seismic records [Johnson et al., 2015; Parsons and Velasco, 2011], which invokes questions regarding whether or not large magnitude events can be dynamically triggered. Emerging results in the literature indicate undocumented M≥5.5 events at near to intermediate distances are dynamically triggered during the passage of surface waves but are undetected by automated networks [Fan and Shearer, 2016]. This raises new questions about the amplitude and duration of dynamic stressing for large magnitude events. I used 35-years of global seismicity and find that large event rate increases only occur following a delay from the transient load, suggesting aseismic processes are associated with large magnitude triggered events. To extend this finding I investigated three cases of large magnitude delayed dynamic triggering following the M8.6 2012 Indian Ocean earthquake [Pollitz et al., 2012] by producing microseismicity catalogs and modeling the transient stresses. The results indicate immediate triggering of microseismic events that hours later culminate into a large magnitude event and support the notion that large magnitude events are triggerable by transient loading, but seismic and aseismic processes (e.g. induced creep or fluid mobilization) are contributing to the nucleation process. Open questions remain concerning the source of a nucleation delay period following a stress perturbation that require both geodetic and seismic observations to constrain the source of delayed dynamic triggering and possibly provide insight into a precursory nucleation phase. Induced seismicity has gained much attention in the past 5 years as earthquake rates in regions of low tectonic strain accumulation accelerate to unprecedented levels [Ellsworth, 2013]. The source of the seismicity is attributed to shallow fluid injection associated with energy production. As hydrocarbon extraction continues to increase in the U.S. the deformation and induced seismicity from wastewater injection is providing new avenues to explore crustal properties. The large magnitude events associated with regions of high rate injection support the notion that the crust is critically stressed. Seismic data in these areas provides the opportunity to delineate fault structures in the crust using precise earthquake locations. To augment the studies of transient loading cycles I investigated induced seismicity at The Geysers geothermal field in northern California. Using high-resolution hypocenter data I implement an epidemic type aftershock sequence (ETAS) model to develop seismicity rate time series in the active geothermal field and characterize the migration of fluids from high volume water injection. Subtle stress changes induced by thermo- and poroelastic strains trigger seismicity for 5 months after peak injection at depths 3 km below the main injection interval. This suggests vertical migration paths are maintained in the geothermal field that allows fluid propagation on annual time scales. Fully describing the migration pattern of fluids in the crust and the associated stresses are applicable to tectonic related faulting and triggered seismic activity. Seasonal hydrological loading is a source of annual periodic transient deformation that is ideal for investigating the modulation of seismicity. The initial step in exploring the modulation of seismicity is to validate that a significant annual period does exist in California earthquake records. The periodicity results [Dutilleul et al., 2015] motivate continued investigation of seismically active regions that experience significant seasonal mass loading, i.e. high precipitation and snowfall rates, to quantify the magnitude of seasonal stress changes and possible correlation with seismicity modulation. The implication of this research addresses questions concerning the strength and state of stress on faults. High-resolution water storage time series throughout California are developed using continuous GPS records. The results allow an estimation of the stress changes induced by hydrological loading, which is combined with a detailed focal mechanism analysis to characterize the modulation of seismicity. The hydrologic loading is augmented with the contribution of additional deformation sources (e.g. tidal, atmosphere, and temperature) and find that annual stress changes of 5 kPa are modulating seismicity, most notably on dip-slip structures. These observations suggest that mechanical differences exist between the vertically dipping strike-slip faults and the shallowly dipping oblique structures in California. When comparing all the annual loading cycles it is evident that future studies incorporate all the sources of solid Earth deformation to fully describe the stresses realized on fault systems that respond to seasonal loads.

3-D Voxel FEM Simulation of Seismic Wave Propagation in a Land-Sea Structure with Topography

NASA Astrophysics Data System (ADS)

Ikegami, Y.; Koketsu, K.

2003-12-01

We have already developed the voxel FEM (finite element method) code to simulate seismic wave propagation in a land structure with surface topography (Koketsu, Fujiwara and Ikegami, 2003). Although the conventional FEM often requires much larger memory, longer computation time and farther complicated mesh generation than the Finite Difference Method (FDM), this code consumes a similar amount of memory to FDM and spends only 1.4 times longer computation time thanks to the simplicity of voxels (hexahedron elements). The voxel FEM was successfully applied to inland earthquakes, but most earthquakes in a subduction zone occur beneath a sea, so that a simulation in a land-sea structure should be essential for waveform modeling and strong motion prediction there. We now introduce a domain of fluid elements into the model and formulate displacements in the elements using the Lagrange method. Sea-bottom motions are simulated for the simple land-sea models of Okamoto and Takenaka (1999). The simulation results agree well with their reflectivity and FDM seismograms. In order to enhance numerical stability, not only a variable mesh but also an adaptive time step is introduced. We can now choose the optimal time steps everywhere in the model based the Courant condition. This doubly variable formulation may result in inefficient parallel computing. The wave velocity in a shallow part is lower than that in a deeper part. Therefore, if the model is divided into horizontal slices and they are assigned to CPUs, a shallow slice will consist of only small elements. This can cause unbalanced loads on the CPUs. Accordingly, the model is divided into vertical slices in this study. They also reduce inter-processor communication, because a vertical cross section is usually smaller than a horizontal one. In addition, we will consider higher-order FEM formulation compatible to the fourth-order FDM. We will also present numerical examples to demonstrate the effects of a sea and surface topography on seismic waves and ground motions.

The Indirect Boundary Element Method (IBEM) for Seismic Response of Topographical Irregularities in Layered Media

NASA Astrophysics Data System (ADS)

Contreras Zazueta, M. A.; Perton, M.; Sanchez-Sesma, F. J.; Sánchez-Alvaro, E.

2013-12-01

The seismic hazard assessment of extended developments, such as a dam, a bridge or a pipeline, needs the strong ground motion simulation taking into account the effects of surface geology. In many cases the incoming wave field can be obtained from attenuation relations or simulations for layered media using Discrete Wave Number (DWN). Sometimes there is a need to include in simulations the seismic source as well. A number of methods to solve these problems have been developed. Among them the Finite Element and Finite Difference Methods (FEM and FDM) are generally preferred because of the facility of use. Nevertheless, the analysis of realistic dynamic loading induced by earthquakes requires a thinner mesh of the entire domain to consider high frequencies. Consequently this may imply a high computational cost. The Indirect Boundary Element Method (IBEM) can also be employed. Here it is used to study the response of a site to historical seismic activity. This method is particularly suited to model wave propagation through wide areas as it requires only the meshing of boundaries. Moreover, it is well suited to represent finely the diffraction that can occur on a fault. However, the IBEM has been applied mainly to simple geometrical configurations. In this communication significant refinements of the formulation are presented. Using IBEM we can simulate wave propagation in complex geometrical configurations such as a stratified medium crossed by thin faults or having a complex topography. Two main developments are here described; one integrates the DWN method inside the IBEM in order to represent the Green's functions of stratified media with relatively low computational cost but assuming unbounded parallel flat layers, and the other is the extension of IBEM to deal with multi-regions in contact which allows more versatility with a higher computational cost compared to the first one but still minor to an equivalent FEM formulation. The two approaches are fully described here and their results compared within the hazard studies of CFE-Las Cruces, Nayarit, Mexico, hydroelectrical project. ACKNOWLEDGEMENTS. This study is partially supported by DGAPA-UNAM under Project IN104712.

Wind/seismic comparisons for upgrading existing structures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Giller, R.A.

1989-10-01

This paper depicts the analysis procedures and methods used to evaluate three existing building structures for extreme wind loads. The three structures involved in this evaluation are located at the US Department of Energy's Hanford Site near Richland, Washington. This site is characterized by open flat grassland with few surrounding obstructions and has extreme winds in lieu of tornados as a design basis accident condition. This group of buildings represents a variety of construction types, including a concrete stack, a concrete load-bearing wall structure, and a rigid steel-frame building. The three structures included in this group have recently been evaluatedmore » for response to the design basis earthquake that included non-linear time history effects. The resulting loads and stresses from the wind analyses were compared to the loads and stresses resulting from seismic analyses. This approach eliminated the need to prepare additional capacity calculations that were already contained in the seismic evaluations. 4 refs., 5 figs., 5 tabs.« less

Source characterization and dynamic fault modeling of induced seismicity

NASA Astrophysics Data System (ADS)

Lui, S. K. Y.; Young, R. P.

2017-12-01

In recent years there are increasing concerns worldwide that industrial activities in the sub-surface can cause or trigger damaging earthquakes. In order to effectively mitigate the damaging effects of induced seismicity, the key is to better understand the source physics of induced earthquakes, which still remain elusive at present. Furthermore, an improved understanding of induced earthquake physics is pivotal to assess large-magnitude earthquake triggering. A better quantification of the possible causes of induced earthquakes can be achieved through numerical simulations. The fault model used in this study is governed by the empirically-derived rate-and-state friction laws, featuring a velocity-weakening (VW) patch embedded into a large velocity-strengthening (VS) region. Outside of that, the fault is slipping at the background loading rate. The model is fully dynamic, with all wave effects resolved, and is able to resolve spontaneous long-term slip history on a fault segment at all stages of seismic cycles. An earlier study using this model has established that aseismic slip plays a major role in the triggering of small repeating earthquakes. This study presents a series of cases with earthquakes occurring on faults with different fault frictional properties and fluid-induced stress perturbations. The effects to both the overall seismicity rate and fault slip behavior are investigated, and the causal relationship between the pre-slip pattern prior to the event and the induced source characteristics is discussed. Based on simulation results, the subsequent step is to select specific cases for laboratory experiments which allow well controlled variables and fault parameters. Ultimately, the aim is to provide better constraints on important parameters for induced earthquakes based on numerical modeling and laboratory data, and hence to contribute to a physics-based induced earthquake hazard assessment.

Global seasonal strain and stress models derived from GRACE loading, and their impact on seismicity

NASA Astrophysics Data System (ADS)

Chanard, K.; Fleitout, L.; Calais, E.; Craig, T. J.; Rebischung, P.; Avouac, J. P.

2017-12-01

Loading by continental water, atmosphere and oceans deforms the Earth at various spatio-temporal scales, inducing crustal and mantelic stress perturbations that may play a role in earthquake triggering.Deformation of the Earth by this surface loading is observed in GNSS position time series. While various models predict well vertical observations, explaining horizontal displacements remains challenging. We model the elastic deformation induced by loading derived from GRACE for coefficients 2 and higher. We estimate the degree-1 deformation field by comparison between predictions of our model and IGS-repro2 solutions at a globally distributed network of 700 GNSS sites, separating the horizontal and vertical components to avoid biases between components. The misfit between model and data is reduced compared to previous studies, particularly on the horizontal component. The associated geocenter motion time series are consistent with results derived from other datasets. We also discuss the impact on our results of systematic errors in GNSS geodetic products, in particular of the draconitic error.We then compute stress tensors time series induced by GRACE loads and discuss the potential link between large scale seasonal mass redistributions and seismicity. Within the crust, we estimate hydrologically induced stresses in the intraplate New Madrid Seismic Zone, where secular stressing rates are unmeasurably low. We show that a significant variation in the rate of micro-earthquakes at annual and multi-annual timescales coincides with stresses induced by hydrological loading in the upper Mississippi embayment, with no significant phase-lag, directly modulating regional seismicity. We also investigate pressure variations in the mantle transition zone and discuss potential correlations between the statistically significant observed seasonality of deep-focus earthquakes, most likely due to mineralogical transformations, and surface hydrological loading.

Seismicity rate changes along the central California coast due to stress changes from the 2003 M 6.5 San Simeon and 2004 M 6.0 Parkfield earthquakes

USGS Publications Warehouse

Aron, A.; Hardebeck, J.L.

2009-01-01

We investigated the relationship between seismicity rate changes and modeled Coulomb static stress changes from the 2003 M 6.5 San Simeon and the 2004 M 6.0 Parkfield earthquakes in central California. Coulomb stress modeling indicates that the San Simeon mainshock loaded parts of the Rinconada, Hosgri, and San Andreas strike-slip faults, along with the reverse faults of the southern Los Osos domain. All of these loaded faults, except for the San Andreas, experienced a seismicity rate increase at the time of the San Simeon mainshock. The Parkfield earthquake occurred 9 months later on the loaded portion of the San Andreas fault. The Parkfield earthquake unloaded the Hosgri fault and the reverse faults of the southern Los Osos domain, which both experienced seismicity rate decreases at the time of the Parkfield event, although the decreases may be related to the decay of San Simeon-triggered seismicity. Coulomb stress unloading from the Parkfield earthquake appears to have altered the aftershock decay rate of the southern cluster of San Simeon after-shocks, which is deficient compared to the expected number of aftershocks from the Omori decay parameters based on the pre-Parkfield aftershocks. Dynamic stress changes cannot explain the deficiency of aftershocks, providing evidence that static stress changes affect earthquake occurrence. However, a burst of seismicity following the Parkfield earthquake at Ragged Point, where the static stress was decreased, provides evidence for dynamic stress triggering. It therefore appears that both Coulomb static stress changes and dynamic stress changes affect the seismicity rate.

Seismic behavior of circular reinforced concrete bridge columns under combined loading including torsion.

DOT National Transportation Integrated Search

2009-12-01

Reinforced concrete (RC) columns of skewed and curved bridges with unequal spans and column heights can be subjected to : combined loading including axial, flexure, shear, and torsion loads during earthquakes. The combination of axial loads, shear : ...

Seismic Hazard analysis of Adjaria Region in Georgia

NASA Astrophysics Data System (ADS)

Jorjiashvili, Nato; Elashvili, Mikheil

2014-05-01

The most commonly used approach to determining seismic-design loads for engineering projects is probabilistic seismic-hazard analysis (PSHA). The primary output from a PSHA is a hazard curve showing the variation of a selected ground-motion parameter, such as peak ground acceleration (PGA) or spectral acceleration (SA), against the annual frequency of exceedance (or its reciprocal, return period). The design value is the ground-motion level that corresponds to a preselected design return period. For many engineering projects, such as standard buildings and typical bridges, the seismic loading is taken from the appropriate seismic-design code, the basis of which is usually a PSHA. For more important engineering projects— where the consequences of failure are more serious, such as dams and chemical plants—it is more usual to obtain the seismic-design loads from a site-specific PSHA, in general, using much longer return periods than those governing code based design. Calculation of Probabilistic Seismic Hazard was performed using Software CRISIS2007 by Ordaz, M., Aguilar, A., and Arboleda, J., Instituto de Ingeniería, UNAM, Mexico. CRISIS implements a classical probabilistic seismic hazard methodology where seismic sources can be modelled as points, lines and areas. In the case of area sources, the software offers an integration procedure that takes advantage of a triangulation algorithm used for seismic source discretization. This solution improves calculation efficiency while maintaining a reliable description of source geometry and seismicity. Additionally, supplementary filters (e.g. fix a sitesource distance that excludes from calculation sources at great distance) allow the program to balance precision and efficiency during hazard calculation. Earthquake temporal occurrence is assumed to follow a Poisson process, and the code facilitates two types of MFDs: a truncated exponential Gutenberg-Richter [1944] magnitude distribution and a characteristic magnitude distribution [Youngs and Coppersmith, 1985]. Notably, the software can deal with uncertainty in the seismicity input parameters such as maximum magnitude value. CRISIS offers a set of built-in GMPEs, as well as the possibility of defining new ones by providing information in a tabular format. Our study shows that in case of Ajaristkali HPP study area, significant contribution to Seismic Hazard comes from local sources with quite low Mmax values, thus these two attenuation lows give us quite different PGA and SA values.

Seismic Design of a Single Bored Tunnel: Longitudinal Deformations and Seismic Joints

NASA Astrophysics Data System (ADS)

Oh, J.; Moon, T.

2018-03-01

The large diameter bored tunnel passing through rock and alluvial deposits subjected to seismic loading is analyzed for estimating longitudinal deformations and member forces on the segmental tunnel liners. The project site has challenges including high hydrostatic pressure, variable ground profile and high seismic loading. To ensure the safety of segmental tunnel liner from the seismic demands, the performance-based two-level design earthquake approach, Functional Evaluation Earthquake and Safety Evaluation Earthquake, has been adopted. The longitudinal tunnel and ground response seismic analyses are performed using a three-dimensional quasi-static linear elastic and nonlinear elastic discrete beam-spring elements to represent segmental liner and ground spring, respectively. Three components (longitudinal, transverse and vertical) of free-field ground displacement-time histories evaluated from site response analyses considering wave passage effects have been applied at the end support of the strain-compatible ground springs. The result of the longitudinal seismic analyses suggests that seismic joint for the mitigation measure requiring the design deflection capacity of 5-7.5 cm is to be furnished at the transition zone between hard and soft ground condition where the maximum member forces on the segmental liner (i.e., axial, shear forces and bending moments) are induced. The paper illustrates how detailed numerical analyses can be practically applied to evaluate the axial and curvature deformations along the tunnel alignment under difficult ground conditions and to provide the seismic joints at proper locations to effectively reduce the seismic demands below the allowable levels.

Light Steel-Timber Frame with Composite and Plaster Bracing Panels

PubMed Central

Scotta, Roberto; Trutalli, Davide; Fiorin, Laura; Pozza, Luca; Marchi, Luca; De Stefani, Lorenzo

2015-01-01

The proposed light-frame structure comprises steel columns for vertical loads and an innovative bracing system to efficiently resist seismic actions. This seismic force resisting system consists of a light timber frame braced with an Oriented Strand Board (OSB) sheet and an external technoprene plaster-infilled slab. Steel brackets are used as foundation and floor connections. Experimental cyclic-loading tests were conduced to study the seismic response of two shear-wall specimens. A numerical model was calibrated on experimental results and the dynamic non-linear behavior of a case-study building was assessed. Numerical results were then used to estimate the proper behavior factor value, according to European seismic codes. Obtained results demonstrate that this innovative system is suitable for the use in seismic-prone areas thanks to the high ductility and dissipative capacity achieved by the bracing system. This favorable behavior is mainly due to the fasteners and materials used and to the correct application of the capacity design approach. PMID:28793642

Light Steel-Timber Frame with Composite and Plaster Bracing Panels.

PubMed

Scotta, Roberto; Trutalli, Davide; Fiorin, Laura; Pozza, Luca; Marchi, Luca; De Stefani, Lorenzo

2015-11-03

The proposed light-frame structure comprises steel columns for vertical loads and an innovative bracing system to efficiently resist seismic actions. This seismic force resisting system consists of a light timber frame braced with an Oriented Strand Board (OSB) sheet and an external technoprene plaster-infilled slab. Steel brackets are used as foundation and floor connections. Experimental cyclic-loading tests were conduced to study the seismic response of two shear-wall specimens. A numerical model was calibrated on experimental results and the dynamic non-linear behavior of a case-study building was assessed. Numerical results were then used to estimate the proper behavior factor value, according to European seismic codes. Obtained results demonstrate that this innovative system is suitable for the use in seismic-prone areas thanks to the high ductility and dissipative capacity achieved by the bracing system. This favorable behavior is mainly due to the fasteners and materials used and to the correct application of the capacity design approach.

Compressive Strength Estimation of Marble Specimens using Acoustic Emission Hits in Time and Natural Time Domains: Preliminary Results

NASA Astrophysics Data System (ADS)

Hloupis, George; Stavrakas, Ilias; Vallianatos, Filippos; Triantis, Dimos

2013-04-01

The current study deals with preliminary results of characteristic patterns derived from acoustic emissions during compressional stress. Two loading cycles were applied to a specimen of 4cm x 4cm x 10 cm Dionysos marble while acoustic emissions (AE) were recorded using one acoustic sensor coupled at the expected direction of the main crack (at the center of the specimen). The produced time series comprised from the number of counts per AE hit under increasing and constant load. Processing took place in two domains: in conventional time domain (t), using multiresolution wavelet analysis for the study of temporal variation of the wavelet-coefficients' standard deviation (SDEV) [1] and in natural time domain (χ), using the variance (κ1) of natural-time transformed time-series [2,3]. Results in both cases, dictate that identification of the region where the increasing stress (σ), exceeds 40% of the ultimate compressional strength (σ*), is possible. More specific, in conventional time domain, the temporal evolution of SDEV presents a sharp change around σ* during first loading cycle and less than σ* during second loading cycle. In natural time domain, the κ1 value clearly oscillate around 0.07 at natural time indexes corresponding to σ* during first loading cycle. Merging both results leads to a preliminary observation that we have an identification of the time when the compressional stress exceeds σ*. References [1] Telesca, L., Hloupis, G., Nikolintaga, I., Vallianatos, F.,."Temporal patterns in southern Aegean seismicity revealed by the multiresolution wavelet analysis", Communications in Nonlinear Science and Numerical Simulation, vol. 12, issue 8, pp 1418-1426, 2007 [2] P. A. Varotsos, N. V. Sarlis, and E. S. Skordas, "Natural Time Analysis: The New View of Time. Precursory Seismic Electric Signals, Earthquakes and other Complex Time-Series", Springer-Verlag, Berlin, Heidelberg, 2011. [3] N. V. Sarlis, P. A. Varotsos, and E. S. Skordas, "Flux Avalances in YBa2Cu307-x films and rice piles: natural time domain analysis", Physical Review B, 73, 054504, 2006. Acknowledgements This work was supported by the THALES Program of the Ministry of Education of Greece and the European Union in the framework of the project entitled "Integrated understanding of Seismicity, using innovative Methodologies of Fracture mechanics along with Earthquake and non extensive statistical physics - Application to the geodynamic system of the Hellenic Arc. SEISMO FEAR HELLARC".

Features of Changing Microwave Radiation from Loaded Rock in Elastic Phase

NASA Astrophysics Data System (ADS)

Wu, Lixin; Mao, Wenfei; Huang, Jianwei; Liu, Shanjun; Xu, Zhongying

2017-04-01

Since the discovery of satellite infrared anomaly occurred before some earthquake by Russian geo-scientists in 1980's, both satellite remote sensing on seismic activities and experimental infrared detection on rock physics in process of rock loading were undertaken in many counties including China, Japan, Europe nations and United States. Infrared imager and spectrum instruments were applied to detect the changed infrared radiation from loaded rock to fracturing, which lead to the development of Remote Sensing Rock Mechanics. However, the change of microwave radiation from loaded rock was not so much studied, even if abnormal changes of microwave brightness temperature (MBT) preceding some large earthquakes were observed by satellite sensors such as AMSR-E on boarded Aqua. To monitor rock hazards, seismic activities, and to make earthquake precautions by via of microwave detection or microwave remote sensing, it is fairly demanded to explore the laws of microwave radiation variation with changed stress and to uncover the rock physics. We developed a large scale rock loading system with capability of 500 tons and 10 tons of load, respectively, at two horizontal loading head, and designed a group of microwave detectors in C, K, and Ka bands. To investigate the changed microwave radiation from loaded granite and sandstone in its elastics deformation phase, the first horizontal stress was circularly applied on rock samples of size 10×30×60cm3 at a constant second horizontal stress, and the changes microwave radiation was detected by the detectors hanged overhead the rock sample. The experiments were conducted outdoor at nighttime to keep off environmental radiation and to simulate the satellite observation conditions in background of cool sky. The first horizontal stress and the microwave radiations were synchronically detected and recorded. After reducing the random noise of detected microwave signals with wavelet method, we found the MBT increase with stress rising and decrease with stress dropping, and the correlation factor (R2) of MBT-stress reached 0.88. The experiments and results revealed an important rock physical phenomenon of rock dielectric property changing with stress, which leads to detectable MBT variation.

Bridge seismic retrofit measures considering subduction zone earthquakes.

DOT National Transportation Integrated Search

2015-07-01

Over the years, earthquakes have exposed the vulnerability of reinforced concrete structures under : seismic loads. The recent occurrence of highly devastating earthquakes near instrumented regions, e.g. 2010 Maule, Chile : and 2011 Tohoku, Japan, ha...

Studies of Fault Interactions and Regional Seismicity Using Numerical Simulations

NASA Astrophysics Data System (ADS)

Yikilmaz, Mehmet Burak

Numerical simulations are routinely used for weather and climate forecasting. It is desirable to simulate regional seismicity for seismic hazard analysis. One such simulation tool is the Virtual California earthquake simulator. We have used Virtual California (VC) to study various aspects of fault interaction and analyzed the statistics of earthquake recurrence times and magnitudes generated synthetically. The first chapter of this dissertation investigates the behavior of seismology simulations using three relatively simple models involving a straight strike-slip fault. We show that a series of historical earthquakes observed along the Nankai Trough in Japan exhibit similar patterns to those obtained in our model II. In the second chapter we utilize Virtual California to study regional seismicity in northern California. We generate synthetic catalogs of seismicity using a composite simulation. We use these catalogs to analyze frequency-magnitude and recurrence interval statistics on both a regional and fault specific level and compare our modeled rates of seismicity and spatial variability with observations. The final chapter explores the jump distance for a propagating rupture over a stepping strike-slip fault. Our study indicates that between 2.5 and 5.5 km of the separation distance, the percentage of events that jump from one fault to the next decreases significantly. We find that these step-over distance values are in good agreement with geologically observed values.

A novel EBSD-based finite-element wave propagation model for investigating seismic anisotropy: Application to Finero Peridotite, Ivrea-Verbano Zone, Northern Italy

NASA Astrophysics Data System (ADS)

Zhong, Xin; Frehner, Marcel; Kunze, Karsten; Zappone, Alba

2014-10-01

A novel electron backscatter diffraction (EBSD) -based finite-element (FE) wave propagation simulation is presented and applied to investigate seismic anisotropy of peridotite samples. The FE model simulates the dynamic propagation of seismic waves along any chosen direction through representative 2D EBSD sections. The numerical model allows separation of the effects of crystallographic preferred orientation (CPO) and shape preferred orientation (SPO). The obtained seismic velocities with respect to specimen orientation are compared with Voigt-Reuss-Hill estimates and with laboratory measurements. The results of these three independent methods testify that CPO is the dominant factor controlling seismic anisotropy. Fracture fillings and minor minerals like hornblende only influence the seismic anisotropy if their volume proportion is sufficiently large (up to 23%). The SPO influence is minor compared to the other factors. The presented FE model is discussed with regard to its potential in simulating seismic wave propagation using EBSD data representing natural rock petrofabrics.

Altimetry data and the elastic stress tensor of subduction zones

NASA Technical Reports Server (NTRS)

Caputo, Michele

1987-01-01

The maximum shear stress (mss) field due to mass anomalies is estimated in the Apennines, the Kermadec-Tonga Trench, and the Rio Grande Rift areas and the results for each area are compared to observed seismicity. A maximum mss of 420 bar was calculated in the Kermadec-Tonga Trench region at a depth of 28 km. Two additional zones with more than 300 bar mss were also observed in the Kermadec-Tonga Trench study. Comparison of the calculated mss field with the observed seismicity in the Kermadec-Tonga showed two zones of well correlated activity. The Rio Grande Rift results showed a maximum mss of 700 bar occurring east of the rift and at a depth of 6 km. Recorded seismicity in the region was primarily constrained to a depth of approximately 5 km, correlating well to the results of the stress calculations. Two areas of high mss are found in the Apennine region: 120 bar at a depth of 55 km, and 149 bar at the surface. Seismic events observed in the Apennine area compare favorably with the mss field calculated, exhibiting two zones of activity. The case of loading by seamounts and icecaps are also simulated. Results for this study show that the mss reaches a maximum of about 1/3 that of the applied surface stress for both cases, and is located at a depth related to the diameter of the surface mass anomaly.

Structural pounding of concrete frame structure with masonry infill wall under seismic loading

NASA Astrophysics Data System (ADS)

Ismail, Rozaina; Hasnan, Mohd Hafizudin; Shamsudin, Nurhanis

2017-10-01

Structural pounding is additional problem than the other harmful damage that may occurs due to the earthquake vibrations. A lot of study has been made by past researcher but most of them did not include the walls. The infill masonry walls are rarely involved analysis of structural systems but it does contribute to earthquake response of the structures. In this research, a comparison between adjacent building of 10-storey and 7-storey concrete frame structure without of masonry infill walls and the same dynamic properties of buildings. The diagonal strut approach is adopted for modeling masonry infill walls. This research also focused on finding critical building separation in order to prevent the adjacent structures from pounding. LUSAS FEA v14.03 software has been used for modeling analyzing the behavior of structures due to seismic loading and the displacement each floor of the building has been taken in order to determine the critical separation distance between the buildings. From the analysis that has been done, it is found that masonry infill walls do affect the structures behavior under seismic load. Structures without masonry infill walls needs more distance between the structures to prevent structural pounding due to higher displacement of the buildings when it sways under seismic load compared to structures with masonry infill walls. This shows that contribution of masonry infill walls to the analysis of structures cannot be neglected.

Experimental and analytical studies on the seismic behavior of conventional and hybrid braced frames

NASA Astrophysics Data System (ADS)

Lai, Jiun-Wei

This dissertation summarizes both experimental and analytical studies on the seismic response of conventional steel concentrically braced frame systems of the type widely used in North America, and preliminary studies of an innovative hybrid braced frame system: the Strong-Back System. The research work is part of NEES small group project entitled "International Hybrid Simulation of Tomorrow's Braced Frames." In the experimental phase, a total of four full-scale, one-bay, two-story conventional braced frame specimens with different bracing member section shapes and gusset plate-to-beam connection details were designed and tested at the NEES Berkeley Laboratory. Three braced frame specimens were tested quasi-statically using the same predefined loading protocol to investigate the inelastic cyclic behavior of code-compliant braced frames at both the global and local level. The last braced frame specimen was nearly identical to one of those tested quasi-statically. However, it was tested using hybrid simulation techniques to examine the sensitivity of inelastic behavior on loading sequence and to relate the behavior observed to different levels of seismic hazard. Computer models of the test specimens were developed using two different computer software programs. In the software framework OpenSees fiber-based line elements were used to simulate global buckling of members and yielding and low-cycle fatigue failure at sections. The LS-DYNA analysis program was also used to model individual struts and the test specimens using shell elements with adaptive meshing and element erosion features. This program provided enhanced ability to simulate section local buckling, strain concentrations and crack development. The numerical results were compared with test results to assess and refine and the ability of the models to predict braced frame behavior. A series of OpenSees numerical cyclic component simulations were then conducted using the validated modeling approach. Two hundred and forty pin-ended struts with square hollow structural section shape were simulated under cyclic loading to examine the effect of width-to-thickness ratios and member slenderness ratios on the deformation capacity and energy dissipation characteristics of brace members. The concept of a hybrid system, consisting of a vertical elastic truss or strong-back, and a braced frame that responds inelastically, is proposed herein to mitigate the tendency of weak-story mechanisms to form in conventional steel braced frames. A simple design strategy about member sizing of the proposed Strong-Back System is provided in this study. To assess the ability of the new Strong-Back System to perform well under seismic loading, a series of inelastic analyses were performed considering three six-story hybrid braced frames having different bracing elements, and three six-story conventional brace frames having different brace configurations. Monotonic and cyclic quasi-static inelastic analyses and inelastic time history analyses were carried out. The braced frame system behavior, bracing member force-displacement hysteresis loops, and system residual drifts were the primary response quantities examined. These indicated that the new hybrid system was able to achieve its design goals. Experimental results show for the same loading history that the braced frame specimen using round hollow structural sections as brace members has the largest deformation capacity among the three types of bracing elements studied. Beams connected to gusset plates at the column formed plastic hinges adjacent to the gusset plate. The gusset plates tend to amplify the rotation demands at these locations and stress concentrations tended to result in early fractures of the plastic hinges that form. To remedy this problem, pinned connection details used in the last two specimens; these proved to prevent failures at these locations under both quasi-static and pseudo-dynamic tests. Failure modes observed near the column to base plate connections in all of the specimens suggest the need for further study. Both OpenSees and LS-DYNA models developed in this study predict the global braced frame behavior with acceptable accuracy. In both models, low-cycle fatigue damage models were needed to achieve an acceptable level of fidelity. Shell element models were able to predict local behavior and the mode of failures with greater but not perfect confidence. OpenSees analysis results show that the proposed hybrid braced frames would perform better than conventional braced frames and that the story deformations are more uniform. Finally, future research targets are briefly discussed at the end of this dissertation.

Lattice Boltzmann Simulation of Seismic Mobilization of Residual Oil in Sandstone

NASA Astrophysics Data System (ADS)

Guo, R.; Jiang, F.; Deng, W.

2017-12-01

Seismic stimulation is a promising technology for enhanced oil recovery. However, current mechanism studies are mainly in the single constricted tubes or idealized porous media, and no study has been conducted in real reservoir porous media. We have developed a numerical simulation which uses the lattice Boltzmann method to directly calculate the characteristics of residual oil clusters to quantify seismic mobilization of residual oil in real Berea sandstone in a scale of 400μm x 400μm x 400μm. The residual oil clusters will be firstly obtained by applying the water flooding scheme to the oil-saturated sandstone. Then, we will apply the seismic stimulation to the sandstone by converting the seismic effect to oscillatory inertial force and add to the pore fluids. This oscillatory inertial force causes the mobilization of residual oil by overcoming the capillary force. The response of water and oil to the seismic stimulation will be observed in our simulations. Two seismic oil mobilization mechanisms will be investigated: (1) the passive response of residual oil clusters to the seismic stimulation, and (2) the resonance of oil clusters subject to low frequency seismic stimulation. We will then discuss which mechanism should be the dominant mechanism for the seismic stimulation oil recovery for practical applications.

Forward and adjoint spectral-element simulations of seismic wave propagation using hardware accelerators

NASA Astrophysics Data System (ADS)

Peter, Daniel; Videau, Brice; Pouget, Kevin; Komatitsch, Dimitri

2015-04-01

Improving the resolution of tomographic images is crucial to answer important questions on the nature of Earth's subsurface structure and internal processes. Seismic tomography is the most prominent approach where seismic signals from ground-motion records are used to infer physical properties of internal structures such as compressional- and shear-wave speeds, anisotropy and attenuation. Recent advances in regional- and global-scale seismic inversions move towards full-waveform inversions which require accurate simulations of seismic wave propagation in complex 3D media, providing access to the full 3D seismic wavefields. However, these numerical simulations are computationally very expensive and need high-performance computing (HPC) facilities for further improving the current state of knowledge. During recent years, many-core architectures such as graphics processing units (GPUs) have been added to available large HPC systems. Such GPU-accelerated computing together with advances in multi-core central processing units (CPUs) can greatly accelerate scientific applications. There are mainly two possible choices of language support for GPU cards, the CUDA programming environment and OpenCL language standard. CUDA software development targets NVIDIA graphic cards while OpenCL was adopted mainly by AMD graphic cards. In order to employ such hardware accelerators for seismic wave propagation simulations, we incorporated a code generation tool BOAST into an existing spectral-element code package SPECFEM3D_GLOBE. This allows us to use meta-programming of computational kernels and generate optimized source code for both CUDA and OpenCL languages, running simulations on either CUDA or OpenCL hardware accelerators. We show here applications of forward and adjoint seismic wave propagation on CUDA/OpenCL GPUs, validating results and comparing performances for different simulations and hardware usages.

Global and local waveform simulations using the VERCE platform

NASA Astrophysics Data System (ADS)

Garth, Thomas; Saleh, Rafiq; Spinuso, Alessandro; Gemund, Andre; Casarotti, Emanuele; Magnoni, Federica; Krischner, Lion; Igel, Heiner; Schlichtweg, Horst; Frank, Anton; Michelini, Alberto; Vilotte, Jean-Pierre; Rietbrock, Andreas

2017-04-01

In recent years the potential to increase resolution of seismic imaging by full waveform inversion has been demonstrated on a range of scales from basin to continental scales. These techniques rely on harnessing the computational power of large supercomputers, and running large parallel codes to simulate the seismic wave field in a three-dimensional geological setting. The VERCE platform is designed to make these full waveform techniques accessible to a far wider spectrum of the seismological community. The platform supports the two widely used spectral element simulation programs SPECFEM3D Cartesian, and SPECFEM3D globe, allowing users to run a wide range of simulations. In the SPECFEM3D Cartesian implementation the user can run waveform simulations on a range of pre-loaded meshes and velocity models for specific areas, or upload their own velocity model and mesh. In the new SPECFEM3D globe implementation, the user will be able to select from a number of continent scale model regions, or perform waveform simulations for the whole earth. Earthquake focal mechanisms can be downloaded within the platform, for example from the GCMT catalogue, or users can upload their own focal mechanism catalogue through the platform. The simulations can be run on a range of European supercomputers in the PRACE network. Once a job has been submitted and run through the platform, the simulated waveforms can be manipulated or downloaded for further analysis. The misfit between the simulated and recorded waveforms can then be calculated through the platform through three interoperable workflows, for raw-data access (FDSN) and caching, pre-processing and finally misfit. The last workflow makes use of the Pyflex analysis software. In addition, the VERCE platform can be used to produce animations of waveform propagation through the velocity model, and synthetic shakemaps. All these data-products are made discoverable and re-usable thanks to the VERCE data and metadata management layer. We demonstrate the functionality of the VERCE platform with two use cases, one using the pre-loaded velocity model and mesh for the Maule area of Chile using the SPECFEM3D Cartesian workflow, and one showing the output of a global simulation using the SPECFEM3D globe workflow. It is envisioned that this tool will allow a much greater range of seismologists to access these full waveform inversion tools, and aid full waveform tomographic and source inversion, synthetic shakemap production and other full waveform applications, in a wide range of tectonic settings.

Dynamic load environment of bridge-mounted sign support structures : research implementation plan.

DOT National Transportation Integrated Search

2005-09-01

Welded aluminum highway sign support trusses must withstand in-service dynamic loads, which largely : constitute the fatigue environment. Sources of these dynamic loads include the natural wind and seismic : environment, the artificial wind environme...

Assessment of seismic loading on structures based on airborne LiDAR data from the Kalochori urban area (N. Greece)

NASA Astrophysics Data System (ADS)

Rovithis, Emmanouil; Kirtas, Emmanouil; Marini, Eleftheria; Bliziotis, Dimitris; Maltezos, Evangelos; Pitilakis, Dimitris; Makra, Konstantia; Savvaidis, Alexandros

2016-08-01

Airborne LiDAR monitoring integrated with field data is employed to assess the fundamental period and the seismic loading of structures composing an urban area under prescribed earthquake scenarios. Α piecewise work-flow is adopted by combining geometrical data of the building stock derived from a LiDAR-based 3D city model, structural data from in-situ inspections on representative city blocks and results of soil response analyses. The procedure is implemented in the residential area of Kalochori, (west of Thessaloniki in Northern Greece). Special attention is paid to the in-situ inspection of the building stock in order to discriminate recordings between actual buildings and man-made constructions that do not conform to seismic design codes and to acquire additional building stock data on structural materials, typologies and number of stories which is not feasible by the LiDAR process. The processed LiDAR and field data are employed to compute the fundamental period of each building by means of code-defined formulas. Knowledge of soil conditions in the Kalochoti area allows for soil response analyses to obtain free-field at ground surface under earthquake scenarios with varying return period. Upon combining the computed vibrational characteristics of the structures with the free-field response spectra, the seismic loading imposed on the structures of the urban area under investigation is derived for each one of the prescribed seismic motions. Results are presented in GIS environment in the form of spatially distributed spectral accelerations with direct implications in seismic vulnerability studies of an urban area.

Bed load transport and boundary roughness changes as competing causes of hysteresis in the relationship between river discharge and seismic amplitude recorded near a steep mountain stream

NASA Astrophysics Data System (ADS)

Roth, Danica L.; Finnegan, Noah J.; Brodsky, Emily E.; Rickenmann, Dieter; Turowski, Jens M.; Badoux, Alexandre; Gimbert, Florent

2017-05-01

Hysteresis in the relationship between bed load transport and river stage is a well-documented phenomenon with multiple known causes. Consequently, numerous studies have interpreted hysteresis in the relationship between seismic ground motion near rivers and some measure of flow strength (i.e., discharge or stage) as the signature of bed load transport. Here we test this hypothesis in the Erlenbach stream (Swiss Prealps) using a metric to quantitatively compare hysteresis in seismic data with hysteresis recorded by geophones attached beneath steel plates within the streambed, a well-calibrated proxy for direct sediment transport measurements. We find that while both the geophones and seismometers demonstrate hysteresis, the magnitude and direction of hysteresis are not significantly correlated between these data, indicating that the seismic signal at this site is primarily reflecting hysteresis in processes other than sediment transport. Seismic hysteresis also does not correlate significantly with the magnitude of sediment transport recorded by the geophones, contrary to previous studies' assumptions. We suggest that hydrologic sources and changes in water turbulence, for instance due to evolving boundary conditions at the bed, rather than changes in sediment transport rates, may sometimes contribute to or even dominate the hysteresis observed in seismic amplitudes near steep mountain rivers.