Seismic risk mitigation technique for art objects: experimental evaluation and numerical modelling of double concave curved surface sliders

The opportunity and the advantages of applying base isolation technologies developed for civil structures to small objects are investigated and some considerations related to the scale effects are discussed. A series of experimental tests carried out at SRMD facility in San Diego are presented with the aim to investigate the dynamic response of a system designed to simulate a sculpture isolated by means of Double Concave Curved Surface Sliders. Such devices are specifically re-designed to fit the peculiar situation of light weight objects. In particular the samples used in the experimentation are designed with geometric and inertia characteristics fit for representing the Michelangelo’s sculptures exposed at the “Galleria dell’Accademia” in Florence, Italy. Finally, the main aspects concerning the numerical simulation of the seismic response of the proposed isolated system are discussed, also to investigate the efficiency of the existing numerical models, which are developed for traditional devices, when they are applied on such small and particular objects. The results of the experimental/numerical campaign show a general efficiency of the isolation system in terms of limiting the transferred action. Further study is needed to understand more clearly some local unexpected phenomena emerged from the experimental tests.     

仿古建筑的隔震加固研究

为对具有传统文化的仿古建筑在不破坏其原有建筑风貌的基础上进行抗震加固,将隔震技术应用于仿古建筑。采用数值模拟的方法研究仿古建筑隔震,建立隔震与非隔震两种有限元计算分析模型,并对加固前、后模型进行动力时程分析。通过输入不同频谱的El Centro Site波、Taft波和RH1TG055波,对比分析两种模型在地震作用下的位移、加速度时程曲线。结果表明：采用隔震技术加固方案后,结构的地震响应得到明显的降低,从而显著地提高了结构的抗震能力。利用该方法加固,节约造价,施工方便,可供同类工程参考。     

Seismic behaviour of shallow foundations: Shaking table experiments vs numerical modelling

The capability of a simplified approach to model the behaviour of shallow foundations during earthquakes is explored by numerical simulation of a series of shaking table tests performed at the Public Works Research Institute, Tsukuba, Japan. After a summary of the experimental work, the numerical model is introduced, where the whole soil–foundation system is represented by a multi‐degrees‐of‐freedom elasto‐plastic macro‐element, supporting a single degree‐of‐freedom superstructure. In spite of its simplicity and of the large intensity of the excitation involving a high degree of nonlinearity in the foundation response, the proposed approach is found to provide very satisfactory results in predicting the rocking behaviour of the system and the seismic actions transmitted to the superstructure. The agreement is further improved by introducing a simple degradation rule of the foundation stiffness parameters, suitable to capture even some minor details of the observed rocking response. On the other hand, the performance of the model is not fully satisfactory in predicting vertical settlements. Copyright © 2007 John Wiley & Sons, Ltd.     

Shake-table tests and numerical simulation of an innovative isolation system for highway bridges

Damage investigation of small to medium-span highway bridges in Wenchuan earthquake revealed that typical damage of these bridges included: sliding between laminated-rubber bearings and bridge girders, concrete shear keys failure, excessive girder displacements and even span collapse. However, the bearing sliding could actually act as a seismic isolation for piers, and hence, damage to piers for these bridges was minor during the earthquake. Based on this concept, an innovative solation system for highway bridges with laminated-rubber bearings is developed. The system is comprised of typical laminated-rubber bearings and steel dampers. Bearing sliding is allowed during an earthquake to limit the seismic forces transmitting to piers, and steel dampers are applied to restrict the bearing displacements through hysteretic energy dissipation. As a major part of this research, a quarter-scale, two-span bridge model was constructed and tested on the shake tables to evaluate the performance of this isolation system. The bridge model was subjected to a Northridge and an artificial ground motion in transverse direction. Moreover, numerical analyses were conducted to investigate the seismic performance of the bridge model. Besides the test bridge model, a benchmark model with the superstructure fixed to the substructure in transverse direction was also included in the numerical analyses. Both the experimental and the numerical results showed high effectiveness of this proposed isolation system in the bridge model. The system was found to effectively control the pier-girder relative displacements, and simultaneously, protect the piers from severe damage. Numerical analyses also validated that the existing finite element methods are adequate to estimate the seismic response of bridges with this isolation system.     

Test and numerical simulation of a new type of hybrid control technique

In this paper, a new hybrid control technique, based on a combination of base-isolation and semi-active variable stiffness/damping in a superstructure, is presented. To illustrate the efficiency of the proposed control system, model tests on a mini-electromagnetic shaking table and a numerical simulation were performed. The test and numerical calculation results indicate that this new hybrid control mode with additional damping and smaller additional stiffness can achieve a better control efficiency.     

Active Tectonics and Seismic Zonation of the Urban Area of Florence, Italy

—?The city of Florence possesses a concentration of cultural and artistic treasures which is unique in the world. In this sense it has a particularly high seismic exposure and a potentially high vulnerability. In order to better evaluate its seismic hazard and risk, we analyzed the seismic response of the urban area of Florence by performing a multidisciplinary study on the effects of earthquakes on the city. By a computer aided methodology we re-evaluated the seismic intensity reports of the May 18 and June 6, 1895 earthquakes in different parts of the city and compared these data with recent studies on surface geology, active tectonics and actual fault movements in the Florence basin. We concluded that more detailed studies of soil response are needed to form a basis for public policy.     

Seismic ground motion amplification pattern induced by a subway tunnel: Shaking table testing and numerical simulation

A series of 1 g shaking table tests, followed by the numerical simulations, is performed to investigate the effect of a circular subway tunnel on the ground motion amplification pattern. Effects of various parameters, including shear wave velocity of soil, frequency content of input motion, flexibility ratio and tunnel depth on the amplification pattern is investigated. Experimental study revealed that the tunnel did not affect free field response at dimensionless period greater than 10. Subsequent parametric study demonstrated that the amount of amplifications were mainly controlled by dimensionless period, dimensionless depth and flexibility ratio. Tunnel effect on the amplification pattern is more prominent for dimensionless period between 3 to 10, flexibility ratio greater than 1 and dimensionless depth less than 3. The study revealed that subway tunnel influences the seismic response of low period buildings, constructed above the tunnel.     

Sediment dynamics in the Black Sea: numerical modelling and remote sensing observations

Here, we address the sediment dynamics in the Black Sea based on analysis of remote sensing data from the Medium Resolution Imaging Spectrometer and numerical simulations with Nucleus for European Modelling of the Ocean model. Boundary conditions consist of realistic meteorological forcing, including significant wave height generated by wave prediction model. A number of sensitivity runs was analysed with the aim to find the most suitable parameters governing sediment fluxes. The comparison between numerical simulations and remote sensing data gives credibility to the quality of simulations. The combined effect of wind waves and currents in the bed layer controls the sediment resuspension that appears to be the major basin-wide source of sediment. Sensitivity experiments included or excluded different forcing terms, e.g. sediment flux from rivers enable to determine the spatial extensions of different point sources. It is concluded that wind-wave forcing is manifested in the sediment dynamics through episodic high energy events contributing to the increase of horizontal sediment fluxes over the northwestern shelf. Both satellite images and numerical model simulations demonstrated that the penetration of suspended sediment into the basin interior was governed by the dynamics of coastal and open-ocean eddies. While fine sediment at sea surface could cross the continental slope propagating into the open ocean, coarser fractions follow the bottom and their penetration into the open ocean is limited. The conclusion is thus that the deposition patterns correlate with the specific shape of Black Sea topography, and the largest depositions are observed in the area of continental slope.     

Seismic reliability assessment and the nonergodicity in the modelling parameter uncertainties

Modelling uncertainty can significantly affect the structural seismic reliability assessment. However, the limit state excursion due to this type of uncertainty may not be described by a Poisson process as it lacks renewal properties with the occurrence of each earthquake event. Furthermore, considering uncertainties related to ground motion representation by employing recorded ground motions together with modelling uncertainties is not a trivial task. Robust fragility assessment, proposed previously by the authors, employs the structural response to recorded ground motion as data in order to update prescribed seismic fragility models. Robust fragility can be extremely efficient for considering also the structural modelling uncertainties by creating a dataset of one-to-one assignments of structural model realizations and as-recorded ground motions. This can reduce the computational effort by more than 1 order of magnitude. However, it should be kept in mind that the fragility concept itself is based on the underlying assumption of Poisson-type renewal. Using the concept of updated robust reliability, considering both the uncertainty in ground motion representation based on as-recorded ground motion and non ergodic modelling uncertainties, the error introduced through structural reliability assessment by using the robust fragility is quantified. It is shown through specific application to an existing RC frame that this error is quite small when the product of the time interval and the standard deviation of failure rate is small and is on the conservative side.     

One decade of multi-objective calibration approaches in hydrological modelling: a review

Abstract

One decade after the first publications on multi-objective calibration of hydrological models, we summarize the experience gained so far by underlining the key perspectives offered by such approaches to improve parameter identification. After reviewing the fundamentals of vector optimization theory and the algorithmic issues, we link the multi-criteria calibration approach with the concepts of uncertainty and equifinality. Specifically, the multi-criteria framework enables recognition and handling of errors and uncertainties, and detection of prominent behavioural solutions with acceptable trade-offs. Particularly in models of complex parameterization, a multi-objective approach becomes essential for improving the identifiability of parameters and augmenting the information contained in calibration by means of both multi-response measurements and empirical metrics (“soft” data), which account for the hydrological expertise. Based on the literature review, we also provide alternative techniques for dealing with conflicting and non-commeasurable criteria, and hybrid strategies to utilize the information gained towards identifying promising compromise solutions that ensure consistent and reliable calibrations.

Citation Efstratiadis, A. & Koutsoyiannis, D. (2010) One decade of multi-objective calibration approaches in hydrological modelling: a review. Hydrol. Sci. J. 55(1), 58–78.     

周至地电阻率观测环境干扰数值模拟

周至地震台NS、NE和NW向测道地电阻率观测值2012年初开始出现显著趋势减小变化,经多次现场调查落实,认为变化与近年来农民在测区大规模搭建猕猴桃架有关.针对周至地震台地电阻率连续观测存在的环境干扰问题,建立环境干扰物理模型,模拟分析环境干扰对周至地电阻率观测的影响,并分析2012年以来趋势减小的可能干扰源.结果显示：2012年新搭建的猕猴桃架会造成周至地电阻率3测道观测值的趋势减小.当干扰电阻R=10 Ω时,模拟结果与实际观测结果符合良好.     

Three-dimensional numerical thermal and rheological modelling in the central Fennoscandian Shield

A three-dimensional model for the central Fennoscandian Shield was constructed for analysing the thermal, the rheological and the structural conditions in the lithosphere. The mesh covers a rectangular area in the southern Finland with horizontal dimensions of 500 km × 400 km and a depth extent of 100 km. Structural boundaries are derived from the several deep seismic soundings carried out in the area. Constructed model is first used in the calculation of the thermal and the rheological models and secondly in analysing the stress and the deformational conditions with the obtained rheology. Thermal and structural models are solved with the finite element method. The calculated surface HFD is between 40 and 48 mW m⁻² in the Proterozoic southern part and below 40 mW m⁻² in the older and northern Archaean part of the model. The calculated rheological strength shows a layered structure with two individual rheologically weak layers in the crust and strong layer in the upper part of the lower crust. The minimum brittle–ductile transition (BDT) depth is around 10 km in the southern part of the model while in the north and north-eastern parts the BDT depth is around 45–50 km. Comparison with the focal depth data shows that as most of the earthquakes occur no deeper than the depth of 10 km are they located in the brittle regime. Resulting stress conditions and possible regions of deformation after the model is subjected to pressure of 50 MPa reveals that the stress field is quite uniformly distributed in different crustal layers and that the elastic parameters control more the state of the stress than the applied rheological structure. In the upper crust, the stress intensity has values between 42 and 45 MPa whereas in the middle crust the values are around 50 MPa. Comparison of the 3-D model with earlier 2-D models shows that some differences in the results are to be expected.     

Three-dimensional numerical thermal and rheological modelling in the central Fennoscandian Shield

A three-dimensional model for the central Fennoscandian Shield was constructed for analysing the thermal, the rheological and the structural conditions in the lithosphere. The mesh covers a rectangular area in the southern Finland with horizontal dimensions of 500 km × 400 km and a depth extent of 100 km. Structural boundaries are derived from the several deep seismic soundings carried out in the area. Constructed model is first used in the calculation of the thermal and the rheological models and secondly in analysing the stress and the deformational conditions with the obtained rheology. Thermal and structural models are solved with the finite element method. The calculated surface HFD is between 40 and 48 mW m⁻² in the Proterozoic southern part and below 40 mW m⁻² in the older and northern Archaean part of the model. The calculated rheological strength shows a layered structure with two individual rheologically weak layers in the crust and strong layer in the upper part of the lower crust. The minimum brittle–ductile transition (BDT) depth is around 10 km in the southern part of the model while in the north and north-eastern parts the BDT depth is around 45–50 km. Comparison with the focal depth data shows that as most of the earthquakes occur no deeper than the depth of 10 km are they located in the brittle regime. Resulting stress conditions and possible regions of deformation after the model is subjected to pressure of 50 MPa reveals that the stress field is quite uniformly distributed in different crustal layers and that the elastic parameters control more the state of the stress than the applied rheological structure. In the upper crust, the stress intensity has values between 42 and 45 MPa whereas in the middle crust the values are around 50 MPa. Comparison of the 3-D model with earlier 2-D models shows that some differences in the results are to be expected.     

基于数值模拟和地震活动性统计的混合地震预测: 在中国地震科学实验场的应用

地震数值预测是近年来最为前沿的科学研究问题之一,但是大部分地震数值预测的探索都受困于数值模拟步长的问题.本文引入混合预测方法,探索了较短预测时长的地震数值预测方法,并在中国地震科学实验场区中进行应用.该混合预测方法主要建立在活动块体理论和构造力学分析的基础上,利用数值模拟推测未来长时间的应力变化量,并引入上一次7级地震之后数十年间,中等强度地震同震应力的影响,最后通过重分类和加权,在应力变化量上叠加年尺度小地震活动异常的统计,从而实现区域地震危险程度的混合评估.基于该方法,我们将地震数值预测的预测时间尺度缩短到年尺度,以年为单位滚动预测了中国地震科学实验场自2021—2030年或更长时间段内可能出现中强地震的区域.2022年芦山6.1级地震、2021年漾濞6.4级地震、2021年泸州6.0级地震和2022年泸定6.8级地震均发生在危险性较高的区域,而2022年马尔康6.0级地震则不在预测区域内.本文展示的混合预测尝试为地震数值预测研究提供了一个可供参考的例证,也为解决数值预测步长控制难的问题提供了一个可能的探索途径.

     

Estimation of topographical effects at Narni ridge (Central Italy): comparisons between experimental results and numerical modelling

In the present work the seismic site response of Narni ridge (Central Italy) is evaluated by comparing experimental results and numerical simulations. The inhabited village of Narni is located in central Apennines at the top of a steep massive limestone ridge. From March to September 2009 the site was instrumented with 10 weak-motion stations, 3 of which located at the base of the ridge and 7 at the top. The velocimetric network recorded 642 events of M_L up to 5.3 and hypocentral distance up to about 100 km. The great amount of data are related to the April 2009 L’Aquila sequence. The site response was analyzed using both reference (standard spectral ratio, SSR) and non reference spectral techniques (horizontal to vertical spectral ratio, HVSR). Moreover directional analyses were performed in order to evaluate the influence of the ridge orientation with respect to the selected source-site paths. In general the experimental results show amplification factors for frequencies between 4 and 5 Hz for almost all stations installed along the crest. The SSR technique provides amplification factors up to 4.5 in a direction perpendicular to the main elongation of the ridge. The results obtained from the data analyses were used as a target for bidimensional and tridimensional numerical simulations, performed using a hybrid finite-boundary element method and a boundary element method for 2D and 3D modelling, respectively. In general, the results obtained through numerical simulation fit well the experimental data in terms of range of amplified frequencies, but they underestimate by a factor of about 2 the observed amplifications.     

Dynamics of the Oder river plume in the Southern Baltic Sea: satellite data and numerical modelling

The Oder river discharge into the Pomeranian Bight of the Baltic Sea was investigated in a combined study using satellite data, numerical modelling and shipborne measurements. The aim was to understand the dynamical processes forming the freshwater distribution patterns during the prevailing winds. From an analysis of typical distribution patterns of the river discharge in relation to the main wind directions and in comparison to seasonal wind statistics, the two main transport directions were determined. The prevailing westerly winds produce an onshore transport and a downwind coastal jet which transports the river water along the Polish coast, in certain cases over a distance of 300 km to the Gdansk Bay. During a period of stable westerly winds in June 1994, the calculated time scale for a water transport over 250 km corresponded to the observed time of 12 d. In spring, the period of maximum river runoff, easterly winds dominate and transport occurs along the German coast into the Arkona Sea. The river water is guided by upwelling processes in front of the Polish coast. During occasional north-easterly winds stable plumes form in front of the Swine river mouth; this occurred in May 1991 for several days. The numerical model showed that the stability of the plume is caused by an interaction between the alignment of the coast, the large-scale circulation in the north, the buoyancy of the freshwater and the Coriolis effect. The underlying anticyclonic eddy is indicated by warm rings in a high resolution Landsat Thematic Mapper scene. From the different datasets the range of the spatial and temporal scales of a stable plume were determined. The volume varied between 0.14 and 0.9 km³, and the suspended matter and chlorophyll load between 1120 and 7200 t and 2.8 and 18 t, respectively. These values are important for ecological budget calculations in turnover process studies.     

Submesoscale tidal eddies in the wake of coral islands and reefs: satellite data and numerical modelling

Interaction of tidal flow with a complex topography and bathymetry including headlands, islands, coral reefs and shoals create a rich submesoscale field of tidal jets, vortices, unsteady wakes, lee eddies and free shear layers, all of which impact marine ecology. A unique and detailed view of the submesoscale variability in a part of the Great Barrier Reef lagoon, Australia, that includes a number of small islands was obtained by using a “stereo” pair of 2-m-resolution visible-band images that were acquired just 54 s apart by the WorldView-3 satellite. Near-surface current and vorticity were extracted at a 50-m-resolution from those data using a cross-correlation technique and an optical-flow method, each yielding a similar result. The satellite-derived data are used to test the ability of the second-generation Louvain-la-Neuve ice-ocean model (SLIM), an unstructured-mesh, finite element model for geophysical and environmental flows, to reproduce the details of the currents in the region. The model succeeds in simulating the large-scale (> 1 km) current patterns, such as the main current and the width and magnitude of the jets developing in the gaps between the islands. Moreover, the order of magnitude of the vorticity and the occurrence of some vortices downstream of the islands are correctly reproduced. The smaller scales (< 500 m) are resolved by the model, although various discrepancies with the data are observed. The smallest scales (< 50 m) are unresolved by both the model- and image-derived velocity fields. This study shows that high-resolution models are able to a significant degree to simulate accurately the currents close to a rugged coast. Very-high-resolution satellite oceanography stereo images offer a new way to obtain snapshots of currents near a complex topography that has reefs, islands and shoals, and is a potential resource that could be more widely used to assess the predictive ability of coastal circulation models.     

Particle path characteristics at the large gravel‐bed river Danube: results from a tracer study and numerical modelling

A tracer study performed on a 3 km long reach of the Danube River in Austria is presented. Forty artificial stones of three different sizes (intermediate b‐axis: 25 mm, 40 mm, 70 mm) were produced and a coded radio acoustic transmitter was implanted. The measurement system had to be improved to be applicable to large rivers with water depths up to 12 m. The positions of the stones were observed approximately once a week, depending on hydrology, over a period of at least one year by radio‐tracking from a boat, including a 15 year flood event. Transport paths and velocities, as well as the incipient motion of bedload transport, could be monitored for the first time on a large gravel‐bed river. The particle paths were found to be mostly bankline‐parallel, even though the stones passed a 30° river bend. The median of the transverse particle displacement was found to be 4% of the longitudinal displacement. Calculations considering both transverse slope and transverse flow velocities showed transverse transport to be 6·6% of the longitudinal transport indicating that marginal lateral transport is mainly influenced by morphology. A three‐dimensional (3D) numerical model using a stochastic particle tracing approach was validated with the data, indicating that the observed positions are well reproduced by the model. Within the observation period, 74% of all stones passed the reach. With more than 1000 detections, particle transport could be characterized by a mean travel velocity of about 10 m per day (variable for the different grain sizes); single tracer stones were transported up to 1000 m during a single flood event. Size‐selective behaviour could be shown and the incipient motion of the large 70 mm gravel was detected at lower discharges than predicted by commonly used uniform bedload transport formulae. Copyright © 2012 John Wiley & Sons, Ltd.