岩石的临界孔隙度反演及横波预测(英文)

临界孔隙度模型是利用岩石的临界孔隙度来计算岩石骨架的弹性模量,岩石的临界孔隙度值受到很多因素的影响,而实际应用中通常无法获得准确的临界孔隙度值,只能选取经验临界孔隙度值,就会给岩石物理建模带来误差。本文提出了一种利用纵波速度反演岩石的临界孔隙度的方法,并且把它应用于横波预测中。实验室和测井数据应用结果表明本文提出的方法可以降低以往选取经验值带来的不确定性,并且能够为横波预测提供准确的临界孔隙度值,提高了横波预测的精度。     

Analytical modelling of a large‐scale dynamic testing facility

An analytical model, which aims at reproducing the response of a large‐scale dynamic testing facility, that is a system composed of the specimen/shaking table/reaction‐mass/airbags/dampers/soil is developed. The Lagrangian of the system is derived, under the assumption of large displacements and rotations. A set of four nonlinear differential equations is obtained and solved with numerical methods. Preliminary verifications of the derived model are carried out by reproducing both well‐known results in the literature as well as those of a lumped model employed in the design of an existing dynamic testing facility. The case‐study for validating the nonlinear equations of motion is the shaking table of the EUCENTRE Laboratory. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of irregular wave overtopping against a smooth sea dike

Based on the filtered Navier-Stokes equations and Smagorinsky turbulence model,a numerical wave flume is developed to investigate the overtopping process of irregular waves over smooth sea dikes.Simulations of fully nonlinear standing wave and regular wave’s run-up on a sea dike are carried out to validate the implementation of the numerical wave flume with wave generation and absorbing modules.To model stationary ergodic stochastic processes,several cases with different random seeds are computed for each specified irregular wave spectrum.It turns out that the statistical mean overtopping discharge shows good agreement with empirical formulas,other numerical results and experimental data.     

Wave height transformation and set-up between a submerged permeable breakwater and a seawall

In this study,we investigated wave transformation and wave set-up between a submerged permeable breakwater and a seawall.Modified time-dependent mild-slope equations,which involve parameters of the porous medium,were used to calculate the wave height transformation and the mean water level change around a submerged breakwater.The numerical solution is verified with experimental data.The simulated results show that modulations of the wave profile and wave set-up are clearly observed between the submerged breakwater and the seawall.In contrast to cases without a seawall,the node or pseudo-node of wave height evolution can be found between the submerged breakwater and the seawall.Higher wave set-up occurs if the nodal or pseudo-nodal point appears near the submerged breakwater.We also examined the influence of the porosity and friction factor of the submerged permeable breakwater on wave transformation and set-up.     

Evaluation of gene expression programming approaches for estimating daily evaporation through spatial and temporal data scanning

Evaporation of water from free water surfaces or from land surfaces is one of the main components of the hydrological cycle, and its occurrence is governed by various meteorological and physical factors. There is a multitude of models developed for estimating daily evaporation values by using weather data. This paper evaluates a Gene Expression Programming (GEP) model for estimating evaporation through spatial and temporal data scanning techniques. It is by using ‘leave‐one‐out’ procedures, a complete scan of the possible train and test set configurations is carried out according to temporal and spatial criteria. Comparison of the GEP model with empirical‐physical models shows that daily evaporation values computed by the GEP model are more accurate. Further, local calibration of the GEP model may not be needed, if enough climatic data are available at other stations. The performance of the GEP model fluctuates throughout the period of study and between stations. A suitable assessment of the model should consider a complete temporal and/or spatial scan of the data set used. Copyright © 2012 John Wiley & Sons, Ltd.     

基于内变量和张量函数表示定理的本构方程

针对各向同性材料,基于张量函数表示定理,建立了本构关系的张量不变性表示,其中,3个不可约基张量取决于应力的0～2次幂,且相互正交,3个系数由塑性应变增量和应力的不变量表示。基于塑性应变增量的不变量定义内变量,本构关系归结为确定内变量的演化。使用张量函数表示定理,给出了内变量演化方程的一般表达式,它取决于应力不变量的增量,因而与主轴旋转无关。讨论了如何根据试验资料和引入适当的假定,确定具体的演化方程。通过与塑性势理论和多重屈服面理论进行比较,表明所建模型是这些理论的最一般表示,且简捷直观、使用方便。     

饱和地基上条形弹性基础的摇摆振动

基于Biot动力方程,研究了饱和均质弹性半空间上弹性条形基础的摇摆振动问题。通过Fourier积分变换求解了饱和土的动力控制方程,然后结合基础底部为混合边界的条件得到了弹性条形基础的摇摆振动对偶积分方程,利用正交多项式将对偶积分方程转化为求解一组线性代数方程组,同时利用复合Simpson法则,得到了动力柔度系数的表达式,通过算例得出了不同参数时地基动力柔度系数随无量纲频率的关系曲线。     

非连续变形分析（DDA）线性方程组的高效求解算法

非连续变形分析(DDA)方法对大规模工程问题的数值模拟耗时太长,其中线性方程组求解耗时可占总计算时间的70%以上,因此,高效的线性方程组解法是重要研究课题。首先,阐述了适用于DDA方法的基于块的行压缩法和基于试验-误差迭代格式的非0位置记录;然后,针对DDA的子矩阵技术,将块雅可比迭代法 (BJ)、预处理的块共轭梯度法 (PCG,包括Jacobi-PCG、SSOR-PCG) 引入DDA方法,重点研究了线性方程组求解过程中的关键运算;最后,通过两个洞室开挖算例,分析了各线性方程组求解算法在DDA中的计算效率。研究表明：与迭代法相比,直解法无法满足大规模工程计算需要;BJ迭代法与块超松弛迭代法(BSOR)的效率差别不大,但明显不如PCG迭代法。因此,建议采用PCG迭代法求解DDA线性方程组,特别是SSOR-PCG值得推广;如果开展并行计算研究,Jacobi-PCG是较好的选择,当刚度矩阵惯性优势明显时,BJ迭代法同样有效。     

墙下条形基础与层状横观各向同性地基共同作用

运用对偶积分方程来求解层状横观各向同性地基与墙下条形基础的共同作用问题。从直角坐标平面应变问题控制方程出发,通过傅里叶（Fourier）变换和层间连续性条件,可以得到层状横观各向同性地基的传递矩阵解。基于该传递矩阵解,并利用条形基础与地基接触的混合边值条件,推导出一组关于基础挠度和地基反力的对偶积分方程。考虑墙下条形基础受到竖向集中荷载的情况,利用弹性薄板理论先求解出条形基础挠度;随后应用雅可比（Jacobi）正交多项式和级数展开的方法,将对偶积分方程转化为线性代数方程组进行求解。编制了相应的计算程序,其计算结果与有限元软件ABAQUS的结果基本吻合,从而验证了所提理论的正确性。算例分析表明,板土相对刚度与地基成层性对地基反力、地表沉降和沿z轴竖向正应力有很大的影响。     

隧道围岩松动圈的现场测试与理论分析

围岩松动圈是制定隧道支护措施和整体稳定性评价的重要参数和依据,如何快速准确地确定松动圈范围已备受关注。以张-石（张家口至石家庄）高速公路上的多条隧道工程为背景,采用声波探测和理论计算相结合的方法,对不同埋深、围岩等级条件下的隧道松动圈进行了研究。基于D-P准则推导了围岩应力松动圈的弹塑性理论和损伤理论计算方法,其计算结果均小于单孔声波法测试结果,但两者变化规律类似。由于考虑了岩体的峰后软化特性,损伤理论分析法计算结果比弹塑性理论更加接近于实测值,运用损伤理论对围岩（特别是低级别围岩）松动圈进行预测是可行的,研究结果对其他类似工程具有借鉴作用。