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猜想等价。     

POS与DEM辅助机载SAR多普勒参数估计

由于机载SAR多普勒参数受载机不稳定飞行和复杂地形影响明显,目前仍缺少有效的矢量估计法。本文对POS和DEM数据及其误差对机载多普勒参数的影响进行了分析,以距离-共面方程构建的SAR几何关系为基础,在POS和DEM数据支持下,建立了一种多普勒参数矢量估计法,并通过仿真对提出的方法进行了验证和分析。理论和仿真表明,本文方法在高山地形、不稳飞行及大斜视角等复杂条件下,亦能够为机载SAR瞬时多普勒中心提供高精度估计。     

侧视雷达影像距离-共面模型严密定位

本文对现有侧视雷达遥感影像的几何构像模型在摄影测量应用中的不足进行了理论分析,以外方位元素作为定向参数,依次建立了像点坐标显函数形式的距离-共面(R-Cp)方程、地心直角坐标系中星载侧视雷达影像的距离-共面方程、以及星载侧视雷达遥感影像的严密定位模型,并进行了稀少控制点条件的星载TerraSAR-X影像定位试验。理论和试验表明,本文所建立的R-Cp严密定位模型,使得雷达影像像点量测坐标作为摄影测量观测值的属性得到体现,容易吸收和利用光学遥感影像基于共线方程模型发展起来的成熟的数字摄影测量数据处理方法,有利于提高雷达遥感影像摄影测量数据处理的严密性、稳定性、可靠性及数据处理效率。     

缺少控制点的卫星遥感对地目标定位

从单线阵推扫式传感器的成像机理出发,利用6个卫星轨道开普勒参数和3个传感器姿态角建立了推扫式卫星遥感影像坐标与其地面点在地心坐标系下的坐标关系式,即构像方程。按照所建立的构像方程,对某地区一景SPOT-5影像进行对地目标定位,获得了实地上83.392m的平面精度;利用单个地面控制点对卫星轨道开普勒参数和传感器姿态实施调整后,目标定位精度提高到14.217m。试验证实,所建立的构像方程是正确的,在卫星遥感对地目标定位中有较好的应用前景。     

共线条件方程式教学中的几个问题

指出了共线条件方程式教学中应注意的一些问题：共线条件方程式是联立的两个平面方程式，存在双主距（fx,fy)时的几何概念，以及它的变换式与直接线性变换关系式的异同点。     

旋转矩阵表达方法对相机检校的影响

在大倾角像片的相机检校中,经典共线方程是采用欧拉角（妒,纠,Ⅳ）描述旋转矩阵的,该方法受限于无法获取位置与姿态的初始值,结果可能导致迭代不收敛。共线方程中的旋转矩阵还可以直接利用方向余弦或单位四元数来描述,相机检校时可以无需依赖位置与姿态的初始值。在相同实验数据、初始值和收敛条件的实验中,直接利用方向余弦描述旋转矩阵的方法明显优于单位四元数方法,主要体现在收敛情况和计算结果接近欧拉角方法两个方面。建议在非量测相机检校时,最好选用直接利用方向余弦描述旋转矩阵的共线方程方法。