Numerical prediction of blast‐induced stress wave from large‐scale underground explosion

This paper presents a numerical model for predicting the dynamic response of rock mass subjected to large‐scale underground explosion. The model is calibrated against data obtained from large‐scale field tests. The Hugoniot equation of state for rock mass is adopted to calculate the pressure as a function of mass density. A piecewise linear Drucker–Prager strength criterion including the strain rate effect is employed to model the rock mass behaviour subjected to blast loading. A double scalar damage model accounting for both the compression and tension damage is introduced to simulate the damage zone around the charge chamber caused by blast loading. The model is incorporated into Autodyn3D through its user subroutines. The numerical model is then used to predict the dynamic response of rock mass, in terms of the peak particle velocity (PPV) and peak particle acceleration (PPA) attenuation laws, the damage zone, the particle velocity time histories and their frequency contents for large‐scale underground explosion tests. The computed results are found in good agreement with the field measured data; hence, the proposed model is proven to be adequate for simulating the dynamic response of rock mass subjected to large‐scale underground explosion. Extended numerical analyses indicate that, apart from the charge loading density, the stress wave intensity is also affected, but to a lesser extent, by the charge weight and the charge chamber geometry for large‐scale underground explosions. Copyright © 2004 John Wiley & Sons, Ltd.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

结构内设备的摩擦摆隔震响应分析

基于系统耦和振动微分方程，分析了结构内摩擦摆隔震设备的响应规律和隔震效果。计算结果表明，经过合理设计的摩擦摆系统能够显著控制设备的地震绝对加速度响应从而有效提高设备的抗震可靠度。     

川西亚高山三个森林群落的湿林冠蒸发速率

于3个假说和林冠上方2m处的气象变量，采用Penman-Monteith组合模型估算了一个生长季节内川西亚高山林区分别以云杉(SF)、冷杉(FF)和白桦(BF)为优势树种的3个林分的湿林冠蒸发速率(Er)。研究结果表明，SF、FF和BF的湿林冠蒸发量(E)分别为44．51mm、88．51mm和57．8mm，分别占总降雨量的9．2％、16．6％和10．2％。与SF和BF相比，FF具有最高的月平均Er和蒸发比例。SF、FF和BF的平均Er分别为0．097mm／h(变化范围：0．028-0．487mm／h)、0．242mm／h(变化范围：0．068～0．711mm／h)和0．149mm／h(0．060～0．576mm／h)。最高和最低的月平均Er分别在6月(SF、FF和BF分别为0．120mm／h、0．317mm／h和0．169mm／h)和10月(SF、FF和BF分别为0．083mm／h、0．187mm／h和0．101mm／h)。8：00至16：00期间的平均点Er显著高于0：00至8：00以及16：00至0：00期间的平均Er。Er显著的日变化和月变化主要归因于林冠上方的太阳辐射、空气温度和相对湿度的变化。     

Fermat猜想的两个等价命题

对每个奇素数p,我们提出了关于p的两个命题:(I)不存在正整数m和n,使代数曲线f(y)=y~3-(2m+1)y-2n的三个零点同时为p方有理整数;(I)不存在正整数m,使不定方程x~2+3y~2=4(2m+1)具有三组正整数解x_i,y_i(i=1,2,3),满足x_1=x_2+x_3,x_1~2+x_2~2+x_3~2=6(2m+1)。本文旨在证明(I)及(I)都与关于素数p的Fermat猜想等价。     

变降深抽水附近地下水非稳定运动的积分方程解

根据完整井变流量抽水附近地下水运动的解析解,利用积分方程的方法,推导出了水井水位降深以任意一种函数关系随时间变化时含水层中地下水非稳定运动积分方程形式的解,并提出了该积分解相应的近似计算方法.最后通过一个实例验证了该方法的实用性.     

整体最小二乘的迭代解法

在引入整体最小二乘平差准则的基础上,推导了整体最小二乘的迭代解法;同时,引入多元函数隐函数求导的方法以确定未知参数对观测数据的线性信息,解决了整体最小二乘下的精度评定问题。给出了运用新的解法在拟合函数确定以及坐标转换参数确定等方面的应用实例,验证了新算法的可行性。     

关于动力Alfven波的色散与朗道阻尼

根据完全动力论理论,证明了可由简化漂移动力论方程求得动力Ａｌｆｖｅｎ波色散方程,即横向用磁流体力学方程,纵向用Ｖｌａｓｏｖ方程,在初级近似下由此推导出适合于日冕和太阳风的色散关系和Ｌａｎｄａｕ阻尼,得到在性质上与ＭＨＤ　Ａｌｆｖｅｎ波完全不同的动力Ａｌｆｖｅｎ波,太阳风中Ａｌｆｖｅｎ湍流很容易由动力Ａｌｆｖｅｎ波演化而来。提出由动力Ａｌｆｒｅｎ波构筑太阳风高速流模型将更符合观测结果。     

考虑土-结构相互作用的振动台实时子结构试验仿真分析

本文基于状态空间方程进行了实时子结构试验的初步探索,提出了一种新的实时子结构试验方法。通过simulink仿真发现,这种方法能很好地再现整体分析的结构反应。最后,对考虑土-结构相互作用的振动台实时子结构试验进行了仿真分析。     

Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.