Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

One-dimensional consolidation of soil under multistage load based on continuous drainage boundary

In engineering practice, a rapid loading rate can result in ground failure when the strength of soft soils is relatively low, and a multistage loading scheme is always utilized to deal with this situation. Firstly, under a multistage load and the continuous drainage boundary, an analytical solution of excess pore-water pressure and consolidation degree is obtained by virtue of the superposition formula of excess pore-water pressure, and a more general continuous drainage boundary under arbitrary time-dependent load is developed. Then, a comparison with existing analytical solutions is conducted to verify the present solution. A preliminary attempt on applying the continuous drainage boundary into the finite element model is made, and the feasibility of the numerical model for the one-dimensional consolidation under the continuous drainage boundary is verified by comparing the results calculated by FEM with that from present analytical solution. Finally, the consolidation behavior of soil is investigated in detail for different int erface parameters or loading scheme. The results show that, in land reclamation projects, a horizontal drain should be placed close to the boundary with a smaller interface parameter to improve the consolidation efficiency. The degree of consolidation is also related to the applied time-dependent load and interface parameters.     

Enhancement of backtracking search algorithm for identifying soil parameters

In this paper, an enhanced backtracking search algorithm (so-called MBSA-LS) for parameter identification is proposed with two modifications: (a) modifying the mutation of original backtracking search algorithm (BSA) considering the contribution of current best individual for accelerating convergence speed and (b) novelly incorporating an efficient differential evolution (DE) as local search for improving the quality of population. The proposed MBSA-LS is first validated with better performance than the original BSA and some other typical state-of-the-art optimization algorithms on a benchmark of soil parameter identification in terms of effectiveness, efficiency, and robustness. Then, the efficiency of the MBSA-LS is further illustrated by two representative cases: identifying soil parameters from both laboratory tests and field measurements. All comparisons demonstrate that the proposed MBSA-LS algorithm can give accurate results in a short time. Finally, to conveniently solve the problems of parameter identification, a practical tool ErosOpt for parameter identification is developed by integrating the proposed MBSA-LS and some other efficient algorithms for readers to conduct the parameter identification using optimisation algorithms.     

Sensitivity of soil parameters in unsaturated zone modelling and the relation between effective,laboratory and in situ estimates

Simulation of soil moisture content requires effective soil hydraulic parameters that are valid at the modelling scale. This study investigates how these parameters can be estimated by inverse modelling using soil moisture measurements at 25 locations at three different depths (at the surface, at 30 and 60 cm depth) on an 80 by 20 m hillslope. The study presents two global sensitivity analyses to investigate the sensitivity in simulated soil moisture content of the different hydraulic parameters used in a one‐dimensional unsaturated zone model based on Richards' equation. For estimation of the effective parameters the shuffled complex evolution algorithm is applied. These estimated parameters are compared to their measured laboratory and in situ equivalents. Soil hydraulic functions were estimated in the laboratory on 100 cm³ undisturbed soil cores collected at 115 locations situated in two horizons in three profile pits along the hillslope. Furthermore, in situ field saturated hydraulic conductivity was estimated at 120 locations using single‐ring pressure infiltrometer measurements. The sensitivity analysis of 13 soil physical parameters (saturated hydraulic conductivity (K_s), saturated moisture content (θ_s), residual moisture content (θ_r), inverse of the air‐entry value (α), van Genuchten shape parameter (n), Averjanov shape parameter (N) for both horizons, and depth (d) from surface to B horizon) in a two‐layer single column model showed that the parameter N is the least sensitive parameter. K_s of both horizons, θ_s of the A horizon and d were found to be the most sensitive parameters. Distributions over all locations of the effective parameters and the distributions of the estimated soil physical parameters from the undisturbed soil samples and the single‐ring pressure infiltrometer estimates were found significantly different at a 5% level for all parameters except for α of the A horizon and K_s and θ_s of the B horizon. Different reasons are discussed to explain these large differences. Copyright © 2004 John Wiley & Sons, Ltd.     

PEIV模型参数估计新算法

PEIV(Partial Errors-In-Variables)模型是EIV模型的扩展,它能解决系数矩阵含有非随机元素或存在结构特性的问题。针对常规PEIV模型算法的复杂性,提出了一种PEIV模型参数估计的新算法。该算法将系数矩阵含误差的元素看成是一类观测值,与平差模型原观测值构成两类观测值,将PEIV平差模型表示为类似于传统的最小二乘间接平差模型,再通过非线性最小二乘平差理论,推导出了算法的迭代公式和精度评定公式。算法迭代格式与间接平差类似,通过算例验证了算法的可行性和正确性。     

求解GM(1,1)模型新总体最小二乘算法

在GM(1,1)模型中系数矩阵和观测向量都是由原始序列组成的。系数矩阵中同样是有误差的,与观测向量中的误差一样,亦来源于原始序列,即它们误差同源。不同位置的相同元素应该有相同的改正数,采用传统总体最小二乘求解则不能达到此目的。针对这一缺陷,推导了一种新的总体最小二乘算法;并且通过算例验证了新方法的可行性和有效性。     

IEC的风机选型参数估算方法在中国的适用性研究

国际电工委员会编制的《风力发电机组设计要求》(IEC 61400-1)推荐了针对风电机组安全等级评估的极端风速和湍流强度特征值估算方法,因其简单便捷,在风电领域被广泛采用。利用全国风能资源专业观测网的193座测风塔观测数据,对IEC推荐的极端风速计算方法与我国规范推荐的基于极值I型概率分布方法进行比较,发现两种方法计算的193座塔70 m高度层50年一遇10 min平均最大风速,仅有7座测风塔较为一致,差异在±1%;IEC推荐方法的计算结果多数偏小,其中偏小10%以上的测风塔有121座,偏小30%以上的有44座测风塔,而偏大10%以上的只有9座测风塔;IEC方法计算的极值风速大幅度偏小的测风塔主要分布在台风影响的东南沿海地区,偏差较小的测风塔主要分布在西北和华北地形平缓区域,但同时偏大10%以上的测风塔也多分布在这一地区。以目前行业领域普遍采用的以15 m·s~(-1)风速的平均湍流强度作为风电机组选型指标,与严格按照规范,以15 m·s~(-1)风速段所有样本湍流强度的90%分位数处的值作为指标进行风电机组等级确定作对比,发现193座塔中有46座塔的选型是不安全的,甚至相差两个等级。     

高建筑对周围建筑雷击保护距离的模拟

在已有先导连接参数化方案的基础上,选取近地面层为研究区域,保持方案中其他基本参量不变,通过改变闪电的空间形态,在同一建筑分布的背景下进行多次闪电模拟。研究多个建筑之间的屏蔽作用以及建筑雷击保护距离与建筑相关特征参数之间的关系,结果表明:高建筑对矮建筑具有屏蔽作用,并存在一个临界保护距离,当高、矮建筑高度分别为190 m和165 m,宽度均为20 m时,建筑之间的距离在12 m以内,矮建筑受高建筑完全保护不遭受雷击;建筑之间的距离超过12 m,矮建筑遭雷击次数明显增多。     

东山乌礁湾海滩表层沉积物粒度时空变化分析

为了进一步了解海岛岬湾海岸表层沉积物粒度的时空分布特征和形成机制,作者以福建东山岛乌礁湾为典型研究区,通过秋、春、夏3个不同季节海滩表层沉积物的采集和点位测量,综合应用激光法和筛析法进行粒度测定。结果表明,东山乌礁湾海滩表层沉积物以0.16～0.50 mm之间的中、细砂为主,并含有少量的粗砂和细砾,这与区域砂质沉积背景有关;湾内从南到北剖面沉积物粒径由粗砂到细砂逐渐变细,主要受剖面地形、局地物源、季节性风浪作用和近岸往复水动力的影响。时间变化上,各取样站位表现出沉积物粒径粗、细不同程度的多种变化趋势,以秋季为参考,整体上表现为由南到北粒级的变小、增大和稳定,这与东北、南南西季风影响下的浪、潮作用以及沿岸流系格局变化下的水动力环境有关。