冻土三轴蠕变非线性数学模型研究

采用MTS-810液压伺服材料试验机, 对人工配制的冻黏土试件进行了三轴蠕变试验, 获得了冻黏土在复杂应力状态下的蠕变曲线. 结果表明: 冻土的蠕变变形具有较强的温度敏感性, 温度越高, 这种温度敏感性越强; 相同温度下, 荷载越大, 变形越大. 运用相关理论, 推导了冻黏土在复杂应力状态下的三轴蠕变非线性数学模型, 采用MATLAB软件的数据拟合功能得到了模型方程参数的数值, 模型参数与温度之间存在密切关系, 建立了二者之间的数学表达式. 冻土三轴蠕变非线性数学模型的曲线与试验曲线拟合精度较高, 建立的数学模型可以精确体现冻土的蠕变规律, 能够为实际冻土工程的变形发展预测提供有效的理论指导.     

降雨过程中边坡临界滑动场

降雨入渗过程中孔隙水压力的升高与基质吸力的降低引起边坡稳定性的下降,是导致边坡滑塌的主要诱导因素。利用饱和-非饱和渗流有限元计算得到的孔隙水压力场,基于Fredlund提出的非饱和土抗剪强度理论,对边坡临界滑动场进行改进,提出可以考虑降雨过程的边坡临界滑动场数值模拟方法,能够方便、快速地计算出边坡局部、整体安全系数和相对应的临界滑动面在降雨过程中的变化历程。将该法用于一个典型均质边坡和一个非均质边坡在降雨过程中的稳定性计算,分析降雨持续时间、降雨强度和非饱和强度参数取值等因素对边坡稳定性的影响,并将计算结果与其他方法进行比较,结果表明临界滑动场方法能搜索任意形状最危险滑面,计算的安全系数合理。     

平面应变下结构性对黄土基坑稳定性的影响

鉴于以往很少考虑黄土结构性对基坑开挖的影响,针对目前普遍存在的基坑稳定性问题,开展室内平面应变条件下的真三轴试验,建立黄土结构性与强度之间关系,通过FISH语言将其引入到基坑开挖数值模拟中,分析结构性对基坑开挖的影响。研究表明,结构性参数与强度指标的变化规律与以往研究结果基本一致,表现为结构性参数与黏聚力呈双曲线关系,但参数略有不同,其受摩擦角的影响较小;结构性影响下基坑侧壁出现明显滑带,且此时连续墙的侧向位移大于不考虑其影响时的位移,坑底隆起量亦如此。研究结果对黄土基坑工程的设计和施工具有一定的指导意义。     

基于三维地质模型的排土场边坡整体稳定性探究

目前所广泛采用的平面应变分析手段,针对排土场边坡稳定性评价时有时会出现某剖面安全系数高于设定允许值,实际中却遭破坏的现象,除了常被提及的力学参数、本构模型等因素,排土场凹形夹持效应和凸形发散作用突出以及滑体端部效应的影响,是分析结果不能准确反映实际的根本原因之一。鉴于此,从平面高程数据、地质剖面、钻孔资料等条件,配合插值效果较为稳定的Kriging法对空间数据进行插值,建立排土场三维地质整体模型。然后,采用自编FISH语言接口程序,将地质模型导入FLAC3D中离散化并计算,较为简单、有效地实现三维地质建模、模型可视化、岩土工程计算分析三位一体功能,有效地克服了平面分析无法解决的上述问题。以攀钢兰尖铁矿尖山第7排土场为实,对其整体稳定性进行应力、应变分析,采用三维强度折减法计算排土场三维整体安全系数,与现场调查及平面计算结果较为一致。基于该模型的计算结果可以有效避免排土场空间效应和滑体端部效应,得到的排土场全域应力-应变规律,对露天矿排土场的生产管理,具有一定的指导意义。     

基于加卸载响应比理论的矿柱动力稳定性分析

利用有限差分程序FLAC建立连续和并排的两类非均质矿柱模型,对矿柱在动力载荷循环作用下的破坏过程进行全时程动力分析。将加卸载响应比(LURR)理论引入矿柱动力稳定性分析中,以施加的动力载荷对矿柱作用的周期性变化作为加载量,以矿柱特征点的位移变化和破坏单元的弹性应变能的演化规律作为响应参数,建立矿柱的加卸载响应比模型,确定矿柱在动力载荷作用期间的加卸载响应比时间序列。计算结果表明,矿柱整体上处于稳定状态和趋于失稳时,其加卸载响应比Y值表现出不同的特征,通过对Y值的分析对矿柱的稳定状态进行评价,预测是否临近峰值承载能力,表明在矿柱动力破裂过程数值模拟分析中引入加卸载响应比是矿柱破坏力学行为研究的一种方便而有效的方法。     

乌东德高拱坝极限抗震能力研究

乌东德水电站是金沙江下游河段上4个水电梯级的最上游梯级,大坝为265 m高的双曲拱坝。采用非线性有限元法,将拱坝坝体-坝基作为一个整体系统,在考虑坝体横缝影响及无限地基辐射阻尼影响基础上,模拟坝肩控制性滑块各滑裂面在地震作用下张开、滑移的非线性力学行为,研究设计与校核地震作用下大坝的抗震性能,探讨大坝-坝基整体系统在强震作用下的抗震潜力。研究结果表明：设计地震条件下,坝体仅在拱端应力集中区域拉应力超过抗拉强度,压应力则均满足强度要求,坝肩块体稳定,大坝-坝基体系的整体安全性可以保证;校核地震条件下,坝体应力值略有增加,能够满足不溃坝的设防要求;以大坝-坝基代表性的位移、变形量变化曲线的拐点作为评价极限抗震能力的标准,可知拱坝-坝基体系的抗震超载安全系数为2.7,相应的基岩水平向峰值加速度0.729g,拱坝-坝基体系抗剪参数强度储备安全系数为2.0,乌东德高拱坝具备较强的极限抗震能力。     

Optimal precursor of the transition from Kuroshio large meander to straight path

We used the conditional nonlinear optimal perturbation(CNOP) method to explore the optimal precursor of the transition from Kuroshio large meander(LM) to straight path within a barotropic inflowoutflow model,and found that large amplitudes of the optimal precursor are mainly located in the east of Kyushu,which implies that perturbations in the region are important for the transition from LM to straight path.Furthermore,we investigated the transition processes caused by the optimal precursor,and found that these processes could be divided into three stages.In the first stage,a cyclonic eddy is advected to the formation region of the Kuroshio large meander,which enhances the LM path and causes a cyclonic eddy to shed from the Kuroshio mainstream.This process causes the LM path to change into a small meander path.Subsequently,the small meander is maintained for a period because the vorticity advection is balanced by the beta effect in the second stage.In the third stage,the small meander weakens and the straight path ultimately forms.The positive vorticity advecting downstream is responsible for this process.The exploration of the optimal precursor will conduce to improve the prediction of the transition processes from LM path to straight path,and its spatial structure can be used to guide Kuroshio targeted observation studies.     

Optimization of slope angle and its seismic stability: A case study for the proposed open pit coalmine in Phulbari, NW Bangladesh

The present study reflects upon the results of substantial program of two-dimensional Finite Element Method （FEM） numerical analyses of the open pit that links to slope angle optimization associated with the safety factor of the pit slope of a coal mine in Bangladesh. In the present analyses, two types of models have been presented. The first model estimates safety factor without seismic effect on the overall pit slope of the model; the second model incorporates safety factor with seismic stability of the model. The calculated optimum slope angle of the first model is 31% with a rational safety factor of 1.51, prior to the seismic effect. However, the value is reduced to 0.93, 0.82, and 0.72, after we applies the seismic effect in the second model with M6, M6.5, and M7, respectively. Finally, our modeling results emphasize that for the case of the proposed Phulbari coalmine, there is extremely high prospect for causing massive slope failure along the optimum pit slope angle with 31% if the mine area felt seismic shaking, like the Sikkim （in northern India） earthquake with M6.9 on September 18, 2011.     

Analysis of implicit HHT‐α integration algorithm for real‐time hybrid simulation

Real‐time hybrid simulation is a viable experiment technique to evaluate the performance of structures equipped with rate‐dependent seismic devices when subject to dynamic loading. The integration algorithm used to solve the equations of motion has to be stable and accurate to achieve a successful real‐time hybrid simulation. The implicit HHT α‐algorithm is a popular integration algorithm for conducting structural dynamic time history analysis because of its desirable properties of unconditional stability for linear elastic structures and controllable numerical damping for high frequencies. The implicit form of the algorithm, however, requires iterations for nonlinear structures, which is undesirable for real‐time hybrid simulation. Consequently, the HHT α‐algorithm has been implemented for real‐time hybrid simulation using a fixed number of substep iterations. The resulting HHT α‐algorithm with a fixed number of substep iterations is believed to be unconditionally stable for linear elastic structures, but research on its stability and accuracy for nonlinear structures is quite limited. In this paper, a discrete transfer function approach is utilized to analyze the HHT α‐algorithm with a fixed number of substep iterations. The algorithm is shown to be unconditionally stable for linear elastic structures, but only conditionally stable for nonlinear softening or hardening structures. The equivalent damping of the algorithm is shown to be almost the same as that of the original HHT α‐algorithm, while the period elongation varies depending on the structural nonlinearity and the size of the integration time‐step. A modified form of the algorithm is proposed to improve its stability for use in nonlinear structures. The stability of the modified algorithm is demonstrated to be enhanced and have an accuracy that is comparable to that of the existing HHT α‐algorithm with a fixed number of substep iterations. Both numerical and real‐time hybrid simulations are conducted to verify the modified algorithm. The experimental results demonstrate the effectiveness of the modified algorithm for real‐time testing. Copyright © 2011 John Wiley & Sons, Ltd.     

Nonlinear earthquake analysis of high arch dam–water–foundation rock systems

A nonlinear finite element model for earthquake response analysis of arch dam–water–foundation rock systems is proposed in this paper. The model includes dynamic dam–water and dam–foundation rock interactions, the opening of contraction joints, the radiation damping of semi‐unbounded foundation rock, the compressibility of impounded water, and the upstream energy propagating along the semi‐unbounded reservoir. Meanwhile, a new equivalent force scheme is suggested to achieve free‐field input in the model. The effects of the earthquake input mechanism, joint opening, water compressibility, and radiation damping on the earthquake response of the Ertan arch dam (240 m high) in China are investigated using the proposed model. The results show that these factors significantly affect the earthquake response of the Ertan arch dam. Such factors should therefore be considered in the earthquake response analysis and earthquake safety evaluation of high arch dams. Copyright © 2011 John Wiley & Sons, Ltd.