A simple procedure to simulate a smooth elastic-plastic transition in Cam-Clay models

A simple procedure is proposed to simulate a smooth transition from elastic to elastoplastic behaviour in Cam-Clay models. The procedure consists of the definition of an external constitutive surface where full yield is assumed, and an internal one that allows the definition of the position in which plastic strains start to appear before the external yield surface is reached. The comparison of the model results with different laboratory tests shows the validity of the procedure. The method considers one additional parameter with regard to a “standard” critical state model, and it can easily be implemented in existing integration modules.     

不同强度El Nino的衰减过程.Ⅰ,强El Nino的衰减过程

将1951～2004年期间的12次El Nino事件分为强、中等、较弱和弱4类,结果发现,强和较弱El Nino衰减进入La Nina,但是中等和弱El Nino衰减进入平常态.因此,El Nino的衰减结果与自身强度之间存在密切的非线性关系.进一步的研究表明,负异常信号自西太平洋向中东太平洋的东传主导了强El Nino事件向La Nina的转变过程,其具体动力过程类似于西太平洋振子理论.热带西北太平洋（WNP）大气异常反气旋在强El Nino位相转变中起核心作用,它的维持和缓慢东移是赤道东风异常维持和发展的原因,而后者通过激发Kelvin波导致了ENSO从El Nino向La Nina的转变.     

中地壳的水和水岩相互作用实验及其地球物理涵义

本文重点报道了高温高压下流体与流体-岩石相互作用实验结果,提供了中地壳条件下流体性质和水岩反应速率数据.这些数据有助于理解中地壳的一些地球物理现象.作者进行了25℃~435℃和22~39 MPa条件下水-岩相互作用反应动力学实验.同时,研究水在近临界区至超临界区的性质.一般地说,中地壳大致位于10(15)至25 km的深度范围.各地的地壳厚度不同,但是中地壳高导-低速层的深度范围十分相似.中地壳的顶界温度处于300℃,底界大致为450℃范围,压力高达200 MPa以上.流体-岩石相互作用实验表明:硅酸盐矿物和岩石的硅最大溶解速率出现在300℃~400℃.此时,硅酸盐矿物格架解体.通常,地壳里普遍存在水、流体.地壳构造活动导致断裂空隙、减压、流体流动.这时,有可能导致中地壳处于300℃~450℃流体的压力减低,由超临界区进入临界态、亚临界态.这会引发强烈流动的水与岩石相互作用.溶解反应导致岩层的硅淋失,硅的强烈淋失又会导致硅酸盐矿物格架解体,岩石崩塌.同时,进一步促进流体的流动.实验表明300℃~400℃下的强烈水岩相互作用促进了岩石破坏,并有可能影响岩层的地球物理性质,如高导层出现.另外,实验和理论研究表明处于300℃~400℃流体具有高电导率性质.这些水岩相互作用会使中地壳出现高导-低速层.     

A mathematical model for steady‐state regolith production at constant erosion rate

It has been hypothesized that many soil profiles reach a steady‐state thickness. In this work, such profiles were simulated using a one‐dimensional model of reaction with advective and diffusive solute transport. A model ‘rock’ is considered, consisting of albite that weathers to kaolinite in the presence of chemically inert quartz. The model yields three different steady‐state regimes of weathering. At the lowest erosion rates, a local‐equilibrium regime is established where albite is completely depleted in the weathering zone. This regime is equivalent to the transport‐limited regime described in the literature. With an increase in erosion rate, transition and kinetic regimes are established. In the transition regime, both albite and kaolinite are present in the weathering zone, but albite does not persist to the soil–air interface. In the weathering‐limited regime, here called the kinetic regime, albite persists to the soil–air interface. The steady‐state thickness of regolith decreases with increasing erosion rate in the local equilibrium and transition regimes, but in the kinetic regime, this thickness is independent of erosion rate. Analytical expressions derived from the model are used to show that regolith production rates decrease exponentially with regolith thickness. The steady‐state regolith thickness increases with the Darcy velocity of the pore fluid, and in the local equilibrium regime may vary markedly with small variations in this velocity and erosion rate. In the weathering‐limited regime, the temperature dependences for chemical weathering rates are related to the activation energy for the rate constant for mineral reaction and to the ΔH of dissolution, while for local equilibrium regimes they are related to the ΔH only. The model illustrates how geochemical and geomorphological observations are related for a simple compositional system. The insights provided will be useful in interpreting natural regolith profiles. Copyright © 2010 John Wiley & Sons, Ltd.     

具有随机状态转移矩阵的Kalman滤波

讨论状态转移矩阵为随机性矩阵时的Kalman滤波方法,并给出滤波的具体算法及分析滤波的收敛特性和参数的取值范围。通过实际算例证实了该滤波方法的有效性。     

Multi-scale variation characteristics of polar vortex at 10 hPa in the Southern Hemisphere

By use of 1948-2007 NCEP/NCAR reanalysis monthly geopotential data, a set of circulation indices are defined to characterize the polar vortex at 10 hPa in the Southern Hemisphere, including area-(S), intensity-(P) and centre position-(λc , φc) indices. Sea-sonal variation, interannual anomalies and their possible causes of 10 hPa polar vortex in the Southern Hemisphere are analyzed by using these indices, the relationship between 10 hPa polar vortex strength and the Antarctic Oscillation are analyzed as well. The results show that: (1) the polar region at 10 hPa in the Southern Hemisphere is controlled by anticyclone (cyclone) from Dec. to Jan. (from Mar. to Oct.), Feb. and Nov. are circulation transition seasons. (2) Intensity index (P) and area index (S) of anticy-clone (cyclone) in Jan. (Jul.) show a significant spike in the late 1970s, the anticyclone (cyclone) enhances (weakens) from ex-tremely weak (strong) oscillation to near the climatic mean before a spike, anticyclone tends to the mean state from very strong oscillation and cyclone oscillates in the weaker state after the spike. (3) There is significant interdecadal change for the anticyclone center in Jan., while markedly interannual variation for cyclone center in July. (4) The ozone anomalies can cause the interannual anomaly of the polar anticyclone at 10 hPa in the Southern Hemisphere in Jan. (positive correlation between them), but it is not related to the polar cyclone anomalies. (5) There is notable negative correlation between the polar vortex intensity index P and the Antarctic Oscillation index (AAOI), thus AAOI can be represented by P.     

循环荷载下砂土的剪切硬化边界面本构模型

基于临界状态土力学框架,建立了一个适用于往返循环荷载作用的砂土边界面本构模型。采用无纯弹性域假设,认为受到反向荷载的瞬时土体就产生塑性变形,砂土的弹性区域退化为一个点。屈服面为倒子弹头型,由于砂土孔隙比与压力之间不存在惟一对应的关系,使得屈服面大小无法与体积应变直接耦合,故采用塑性偏应变而不是剑桥模型那种塑性体应变作为硬化参数。流动法则采用加入状态参数的修正的Rowe应力剪胀关系,体现了依赖状态的剪胀思想。屈服面大小的比值 反映了塑性模量的演化,并推导了 的表达式。只用1套参数,该模型就能合理地模拟砂土在不同密度和固结压力下循环荷载的应力-应变关系曲线。     

山西国际旅游市场营销目标选择研究

在梳理旅游市场营销目标选择方法的基础上,对原有竞争态模型进行扩展,分析山西省主体客源国的亲景度和市场发展指数,并基于两个指数的组合,建立旅游市场相对竞争态模型,将山西的主体客源国市场分为兴旺型、崛起型、成熟型和停滞型,并结合旅游市场绝对竞争格局,依据各类型市场的特征,确定山西省国际旅游市场的营销目标。     

主应力轴循环旋转下高密实粉土稳定性影响研究

为研究主应力轴循环旋转对粉土性状影响,对长江入海口高密实饱和粉土进行了主应力轴180°转幅的循环旋转试验。试验发现,等向固结的高密实粉土在恒定动剪应力比的主应力轴循环旋转下能发生孔压达到初始围压水平的液化破坏,但液化前孔压开展会因动应力水平的不同而呈现以渐变相态点和激变相态点区分的3阶段变化形式,或仅以激变相态点区分的两阶段变化形式。试样应变开展则始终呈现由崩塌点区分的两阶段变化形式,且孔压激变相态点与应变崩塌点同步反映了土体结构崩塌状态,崩塌状态时主应变差双向幅值限制在0.2%～0.4%这一较小范围内。试样进行主应力轴旋转达到崩塌时的应力状态与频率无显著关系,可由 - 空间中一条准不稳定状态线归一化。而以崩塌状态为终点进行归一化的主应力轴循环旋转下孔压振次曲线,可用修正的Seed公式进行较好拟合     

三维应力状态下砂的临界状态本构模型分叉分析

就服从非关联流动法则软化型砂的临界状态本构模型,在不同应力路径下的分叉进行了理论和数值分析。理论分析表明：分叉现象强烈依赖于应力路径,即当平均应力p一定时,应力路径在洛德角在-25o～15o之间时,均会在应力-应变关系硬化阶段出现分叉现象,而其他应力路径下不会产生分叉。并且分叉时对应的应力比、主应变及剪切带倾角都随洛德角的增加而变化,其值是先增加而后减小。采用回映应力更新算法,编写了该本构模型材料子程序,借助有限元软件ABAQUS及材料子程序,通过数值计算方法预测到了分叉点所对应的应力状态,表明了分叉现象在数值计算中的存在性,并通过数值预测值和理论解的对比,两者结果基本一致。