地震作用下立式储液罐罐壁“象足”变形仿真分析

基于ANSYS软件建立了考虑液体晃动和罐底提离立式储液罐有限元模型,分别进行了水平地震和竖向地震作用下罐壁“象足”变形分析。分析表明:立式储液罐罐壁“象足”变形主要是由罐壁纵向压应力超过临界应力而产生的局部屈曲破坏,并非强度破坏。因罐底提离导致的罐底与基础反复撞击加大了作用在罐壁上的应力,使罐壁底部“象足”变形不断发展,最终导致罐壁撕裂。在完全相同地震加速度作用下,水平地震作用比竖向地震作用罐壁更早更容易进入屈曲状态,产生“象足”变形。     

地震作用下立式储液罐罐壁"象足"变形仿真分析

基于ANSYS软件建立了考虑液体晃动和罐底提离立式储液罐有限元模型,分别进行了水平地震和竖向地震作用下罐壁"象足"变形分析.分析表明:立式储液罐罐壁"象足"变形主要是由罐壁纵向压应力超过临界应力而产生的局部屈曲破坏,并非强度破坏.因罐底提离导致的罐底与基础反复撞击加大了作用在罐壁上的应力,使罐壁底部"象足"变形不断发展,最终导致罐壁撕裂.在完全相同地震加速度作用下,水平地震作用比竖向地震作用罐壁更早更容易进入屈曲状态,产生"象足"变形.     

位移谐波激励下拱顶锚固储罐的动力屈曲研究

地震下储罐罐壁发生屈曲是其主要的震害形式,研究储罐的动力屈曲特性尤为重要.以位移简谐波为水平激励,运用ADINA有限元软件对高径比分别为0.82、1.24和1.85的3种拱顶锚固储罐进行动力屈曲分析.根据Budiansky-Roth动力屈曲准则,通过谐波位移幅值对应的加速度峰值PGA和储罐位移响应,绘制特征节点的"伪平衡路径"轨迹,以确定临界加速度峰值,判断储罐动力屈曲.3种储罐模型在位移谐波激励下均为弹性屈曲.通过动力屈曲分析得出结论:储罐屈服的位置随高径比的不同有所不同,屈服多发生在径向变形最大处,罐壁内外面的塑性总应变值和屈服应力值大小不同等.最后以储罐模型B为例,将以"伪平衡路径"轨迹确定的临界应力与规范公式的计算结果进行了比较.     

竖向地震作用对储液罐地震响应的影响分析

地震作用下大型储液罐的安全问题日益引起重视。基于ANSYS软件建立储罐液体耦合有限元模型,考虑罐底非线性接触效应,以El-Centro南北向和竖直向记录地震波为输入,研究水平激励以及水平和竖向同时激励两种工况下储罐的动力响应。研究结果表明,两种工况下靠近罐底1.2m处均发生了"象足"变形,竖向激励下水平相对位移增加了14%。竖向激励使得罐壁环向应力和轴向压应力均有不同程度的增加。竖向地震激励对液面的竖向晃动影响较小。储液罐底板在地震作用下发生了竖向提离和永久滑移,竖向激励时增长幅度均在10%左右。同时罐体基底剪力在竖向地震作用下也有所增大。储罐抗震设计时应考虑竖向地震分量的影响,研究结论可为立式储罐的抗震设计提供一定的参考和依据。     

位移谐波激励下拱顶锚固储罐的动力响应分析

地震中拱顶锚固储罐的破坏将造成严重的后果,研究其地震作用下的动力响应特性具有一定工程意义。以位移简谐波为水平激励,运用AD INA有限元软件对一座3 000m3的拱顶锚固储罐模型进行动力数值模拟,系统分析了储罐内液面晃动波高、动水压力极值、罐壁应力、拱顶应力、拱顶位移等动力响应结果,将以位移谐波激励和以加速度谐波激励的储罐动力响应结果进行比较。分析表明,位移谐波激励下拱顶锚固储罐的液面晃动波高受场地影响较大;储罐应力分布具有一定规律性;采用位移谐波激励时其应力、位移和动水压力等较采用加速度谐波激励时大,液体晃动波高较小。     

基于概率估计方法的立式储罐地震易损性研究

结合储罐的震害资料和震害特点,提出了基于概率估计方法的立式储罐地震易损性研究.考虑了结构的随机性及荷载作用的随机性,给出了各随机变量的概率分布模型及给定烈度下地震发生的概率.从可靠度指标的几何意义入手,提出了储罐的地震易损性计算方法.通过2×103m3和5×104m32种高矮罐的地震易损性研究,结果表明:通过概率方法和绘制地震易损性曲线的方式,更能直观的给出储罐各失效模式的地震易损性.储罐的地震易损性与地震烈度、高径比的关系很大,地震易损性随着地震烈度的增大逐渐增大;与2×103m3储罐相比,5×104m3储罐各失效模式地震易损性明显较小;从多种失效模式进行储罐地震易损性研究更具有科学性.     

考虑桩土-隔震层-罐体相互作用的隔震储罐地震响应分析

为了研究桩土、隔震层与罐体相互作用对15×104m3隔震储罐地震响应的影响,桩与隔震层采用弹簧-阻尼单元,罐体采用壳单元,流体采用势流体单元,基于ADINA建立数值仿真分析模型,在8度E1 Centro地震动激励下,应用Newmark数值积分方法进行地震响应分析.结果表明:隔震储罐基底剪力和基底弯矩相对于非隔震储罐明显降低,储罐罐壁加速度沿罐壁高度方向近乎平动,但对波高的控制效果不明显;除罐壁顶部外,动液压力明显降低;隔震后,罐壁有效应力和环向应力减震效应明显,对轴向应力和径向应力影响较小,建议从安全角度考虑,储罐非隔震设计时不考虑桩土作用,隔震设计时考虑桩土对罐体的影响.     

立式浮顶储罐并联基础隔震模拟地震振动台试验研究

立式浮顶储罐在强震作用下会发生罐壁屈曲、卡顶、浮顶下沉等震害。为降低储罐地震响应,以1 000 m3储罐为原型,按1:5的几何相似比设计缩尺模型,并以隔震周期为0.5 s进行并联隔震装置设计,给出其力学本构关系。采用模拟地震振动台试验方法,选取四种不同频谱特性的地震动对模型罐进行地震响应分析,结果表明:隔震层顶部加速度卓越频率有所降低,周期延长;隔震后储罐加速度、层间位移及动应力地震响应明显降低,但晃动波高有所放大;三向地震激励与单向地震激励相比,加速度地震响应明显放大,对晃动波高影响较小。建议储罐抗震减震设计时应考虑三向地震作用。     

考虑SSI效应的15×10~4 m~3储罐基础隔震数值仿真分析

为考虑土与结构相互作用(SSI)对15×104m3大型立式储罐基础隔震效应的影响,采用弹簧-阻尼系统模拟地基土和隔震层,罐壁及底板采用壳单元,流体采用势流体单元,利用ADINA建立15×104m3储罐有限元模型,在峰值加速度0.2g El Centro波地震激励下,应用Newmark数值积分方法进行地震响应分析,结果表明:考虑SSI效应时,非隔震储罐的地震响应有所减小,而基础隔震时地震响应有放大效应。储罐抗震减震设计时,中软地基土上储罐从结构设计安全角度来说需要考虑土与结构的相互作用。     

立式储罐与地基相互作用地震反应分析

在考虑地基与储罐相互作用的情况下,采用有限元法对储罐住水平地震荷载作用下的反应进行了计算。推导了系统刚度矩阵及质量矩阵的计算过程。通过对四种场地条件下的地震响应分析,得出了储罐各种参数对地震反应的影响。分析表明：不同的地震激励对储罐系统的反应是不同的;地基刚度对系统的动力响应的影响也很大：即地基越“柔”,基底弯矩、罐壁轴向应力越小,但环向应力和底板变形却在增加。     

A Monte Carlo study of rainfall forecasting with a stochastic model

A procedure for short-term rainfall forecasting in real-time is developed and a study of the role of sampling on forecast ability is conducted. Ground level rainfall fields are forecasted using a stochastic space-time rainfall model in state-space form. Updating of the rainfall field in real-time is accomplished using a distributed parameter Kalman filter to optimally combine measurement information and forecast model estimates. The influence of sampling density on forecast accuracy is evaluated using a series of a simulated rainfall events generated with the same stochastic rainfall model. Sampling was conducted at five different network spatial densities. The results quantify the influence of sampling network density on real-time rainfall field forecasting. Statistical analyses of the rainfall field residuals illustrate improvement in one hour lead time forecasts at higher measurement densities.     

Earthquake Engineering and Engineering Vibration Aims And Scope

正This journal is established by the Institute of Engineering Mechanics(IEM),China Earthquake Administration,to promote scientific exchange between Chinese and foreign scientists and engineers so as to improve the theory and practice of earthquake hazards mitigation,preparedness,and recovery.To accomplish this purpose,the journal aims to attract a balanced number of papers between Chinese and     

Solute concentrations in the lower niger river and the source rock contribution

Water quality analyses for the Niger River for the 1980/81 hydrological year are presented. The samples were collected from the main river at Lokoja, and from two main tributaries, the Kaduna and the Benue Rivers. Different water types were distinguished by the concentrations of major ions. The type Ca > Na > Mg > K - HCO₃ > SO₄ > Cl was represented at all stations during at least part of the year. Chloride was found to dominate the sulphate ion in the Kaduna and Niger, while the Benue maintained a higher concentration of sulphate relative to chloride all year round. Distinct patterns of seasonal variation in the ion concentrations were observed, particularly for the samples collected at Lokoja. Low ion concentrations were prominent during periods of high discharge, while low flow periods coincided with high dissolved ion concentrations. The contribution of rainwater to the total dissolved solids in the river waters was assessed indirectly using rainwater chemistry data from the Gulf of Guinea. The estimated rainwater contribution to the Lower Niger amounts to 5.15 mg 1^?1. Geochemical weathering calculations involving reactions of the four major minerals of granitic rocks - anorthite, biotite, albite, and K-feldspar - with carbon dioxide and water, can account for the average water composition of the Lower Niger. The proportion of the ionic components was also related to the occurrence of the respective element in the minerals.     

Scenario predictions for potential near-field and far-field earthquakes affecting Hong Kong

Historical earthquakes noted in the written records of the South China region, including Hong Kong, are not well delineated along identified prominent fault sources. Despite the lack of any definitive, localised trend in the spatial distribution of seismic activity in the region, there does appear to be some major disparity in the seismic activity rates (especially for large magnitude earthquakes) between the near-field and the far-field regions of Hong Kong. Despite this observation, previous studies of the regional seismic activity and seismic ground motion hazard (the latter using a probabilistic seismic hazard assessment, PSHA) have considered very broad source zone regions, in which uniform levels of seismic activity have been assumed. The present paper further scrutinises this broad source zone (BSZ) approach by adopting a novel expanding circular disc (ECD) method to determine the rates of earthquake recurrence. Such a method is intended to counter-check previously developed models by determining earthquake scenario events in terms of magnitude–distance (M–R) pairs or combinations, having defined values of average return period. Unlike the BSZ approach, the ECD method specifically accounts for the supposed variations in the seismic activity rates between events in the near-field and the far-field of Hong Kong.The form of the developed method is particularly suited to the determination of design-level earthquake ground motions for bedrock sites, since it assumes a directionally-independent attenuation model as described in the companion paper. It is found that, whilst the BSZ approach may indicate the overall average levels of hazard that are representative of the South China region as a whole, it does not capture the large disparity in seismic activity rates between near-field and far-field events. This important feature is expected to have a significant impact upon engineering assessments of the seismic safety of structures in Hong Kong and elsewhere in the South China region. For example, it is found that for events with M≥6, the seismic activity rate (normalised by time and area) in the very far-field is around 3.5 times larger than in the near-field and medium-field of Hong Kong. The resulting design M-R combinations, covering a range of return periods from 70 to 2500 years, are limited, for very long return periods and for distant events, by the maximum credible earthquake (MCE) magnitude. Intensive research to determine this seismic hazard parameter is recommended, in order to refine further the results of the ECD analysis, which presently conservatively assumes the MCE to range between M=6 in the near-field of Hong Kong to M=8 in the very far-field, at distances greater than 280 km from Hong Kong.     

Depositional ages and characteristics of Middle–Upper Jurassic and Lower Cretaceous lacustrine deposits in southeastern Mongolia

Lower Cretaceous lacustrine oil shales are widely distributed in southeastern Mongolia. Due to the high organic carbon content of oil shale, many geochemical studies and petroleum exploration have been conducted. Although most of the oil shales are considered to be Early Cretaceous in age, a recent study reveals that some were deposited in the Middle Jurassic. The present study aims at establishing depositional ages and characteristics of the Jurassic and Cretaceous lacustrine deposits in Mongolia. The Lower Cretaceous Shinekhudag Formation is about 250 m thick and composed of alternating beds of shale and dolomite. The Middle Jurassic Eedemt Formation is about 150 m thick and composed of alternating beds of shale, dolomitic marl, and siltstone. The alternations of shale and dolomite in both formations were formed by lake level changes, reflecting precipitation changes. Shales were deposited in the center of a deep lake during highstand, while dolomites were formed by primary precipitation during lowstand. Based on the radiometric age dating, the Shinekhudag Formation was deposited between 123.8 ±2.0 Ma and 118.5 ±0.9 Ma of the early Aptian. The Eedemt Formation was deposited at around 165–158 Ma of Callovian–Oxfordian. The calculated sedimentation rate of the Shinekhudag Formation is between 4.7 ±2.6 cm/ky and 10.0 ±7.6 cm/ky. Shales in the Shinekhudag Formation show micrometer‐scale lamination, consisting of algal organic matter and detrital clay mineral couplets. Given the average thickness of micro‐laminae and calculated sedimentation rate, the micro‐lamination is most likely of varve origin. Both Middle–Upper Jurassic and Lower Cretaceous lacustrine oil shales were deposited in intracontinental basins in the paleo‐Asian continent. Tectonic processes and basin evolution basically controlled the deposition of these oil shales. In addition, enhanced precipitation under humid climate during the early Aptian and the Callovian–Oxfordian was another key factor inducing the widespread oil shale deposition in Mongolia.     

THEORETICAL APPROACH TO THE SOLUTION OF THE INFILTRATION PROBLEM

ABSTRACT

The one-dimensional transient downward entry of water in unsaturated soils is investigated theoretically. The mathematical equation describing the infiltration process is derived by combining Darcy's dynamic equation of motion with the continuity and thermodynamic state equations adjusted for the unsaturated flow conditions. The resulting equation together with the corresponding initial and boundary conditions constitues a mathematical initial boundary value problem requiring the solution of a nonlinear partial differential equation of the parabolic type. The volumetric water content is taken as the dependent variable and the time and the position along the vertical direction are taken as the independent variables. The governing equation is of such nature that a solution exists for t > 0 and is uniquely determined if two relationships are defined, together with the specified state of the system, at the initial time t = 0 and at the two boundaries. The two required relations are those of pressure versus permeability and pressure versus volumetric water content.

Since the partial differential equation has strong non-linear terms, a discrete solution is obtained by approximating the derivatives with finite-differences at discrete mesh points in the solution domain and integrated for the corresponding initial and boundary conditions. The use of an implicit difference scheme is employed in order to generate a system of simultaneous non-linear equations that has to be solved for each time increment. For n mesh points the two boundary conditions provide two equations and the repetition of the recurrence formula provides n—2 equations, the total being n equations for each time increment. The solution of the system is obtained by matrix inversion and particularly with a back-substitution technique. The FORTRAN statements used for obtaining the solution with an electronic digital computer (IBM 704) are presented together with the input data.

Analysis of the errors involved in the numerical solution is made and the stability and convergence of the solution of the approximate difference equation to that of the differential equation is investigated. The method applied is that of making a Fourier series expansion of a whole line of errors and then following the progress of the general term of the series expansion and also the behavior of each constituent harmonic. The errors (forming a continuous function of points in an abstract Banach space) are represented by vectors with the Fourier coefficients constituting a second Banach space. The amplification factor of the difference equation is shown to be always less than unity which guarantees the stability of the employed implicit recurrence scheme.

Experiments conducted on a vertical column packed uniformly with very fine sand, show a satisfactory agreement between the theoretically and experimentally obtained values. Many experimental results are shown in an attempt to explain the infiltration phenomenon with emphasis on the shape and movement of the wet front, and the effects of the degree of compaction, initial water content and deaired water on the infiltration rate.     

Late Permian post‐ophiolitic trondhjemites from Central Iran: a mark of subduction role in growth of Paleozoic continental crust

Late Permian trondhjemites in the Anarak area occur as stocks and dykes, which cross cut the Anarak ophiolite and its overlying metasedimentary rocks, and are exposed along the northern Anarak east–west main faults. These leucocratic intrusive bodies have enclaves of all ophiolitic units and metamorphic rocks. They are composed of amphibole, plagioclase (oligoclase), quartz, zircon and muscovite. Secondary minerals are chlorite (pycnochlorite), epidote, albite, magnetite and calcite. Whole‐rock major‐ and trace‐element analyses reveal that they are characterized by high SiO ₂ (67.8–71.0 wt%), Al ₂ O ₃ (14.9–17.1 wt%) and Na ₂ O (5.3–8.6 wt%), low K ₂ O (0.1–1.5 wt%; average: 0.8 wt%), low Rb/Sr ratio (0.01–0.40; average: 0.09), low Y (3–6 ppm), negative Ti, Nb and Ta anomalies, slightly negative or positive Eu anomaly, LREE enrichment and fractionated HREE. These rocks present 2 to 40 times enrichment in inclined chondrite‐normalized REE patterns. Geochemical characteristics of the Anarak trondhjemites all reflect melting of a mafic protolith at more than 10 kbar. The field evidence and whole‐rock chemistry reveal that these rocks have been crystallized from magmas derived from melting of subducted Anarak oceanic crust. This study reveals that melting of garnet amphibolite was an important element of continent formation in the study area.