岩石单轴压缩下能量与损伤演化规律研究

能量变化是材料各参数变化的本质特征。通过颗粒流模拟岩石单轴压缩试验,从细观力学角度分析能量转化、裂纹扩展及损伤演化规律。颗粒流利用细观颗粒的运动克服了宏观力学理论不易实现试件多裂纹破坏形态的缺点,以此研究能量变化更加合理。应力-应变各阶段能量与微裂纹及损伤存在着相互对应的发展关系。通过3种工况论证了宏-细观力学参数对应关系,采用最小二乘法拟合得出微裂纹与轴向应变呈幂次函数关系。采用割线模量定义损伤变量,取模量加速下降处（对应裂纹加速扩展）为损伤门槛,对应损伤阈值为0.158。3种工况下微裂纹数量与损伤呈线性发展关系,为损伤演化的深入研究提供参考。     

地球物理资料群体智能反演(英文)

复杂地球物理资料的反演问题往往是一个求解多参数非线性多极值的最优解问题。而鸟和蚂蚁等群体觅食的过程,正好与寻找地球物理反演最优解的过程相似。基于自然界群体协调寻优的思想,本文提出了交叉学科的群体智能地球物理资料反演方法,并给出了其对应的数学模型。用一个有无限多个局部最优解的已知模型对该类方法进行了试验。然后,将它们应用到了不同的复杂地球物理反演问题中:(1)对噪声敏感的线性问题;(2)非线性和线性同步反演问题;(3)非线性问题。反演结果表明,群体智能反演是可行的。与常规遗传算法和模拟退火法相比,该类方法有收敛速度相对快、收敛精度相对高等优点;与拟牛顿法和列文伯格一马夸特法相比,该类方法有能跳出局部最优解等优点。     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

Dynamics of a particle in a gravitational field of a homogeneous annulus disk

We formulate the problem of an infinitesimal particle moving in the space under the influence of the gravitational force induced by a homogeneous annulus disk fixed on a plane. We compute using different coordinates and in terms of an elliptic integral the potential associated to this problem. Also we study the symmetries of the associated potential. After that, we look at the dynamics in some particular cases, namely: in the line perpendicular to the plane that contains the annulus disk and passes through the center of the ring and on the plane that contains the massive annulus.      

Development dynamics of intense solar proton flares

We consider the relationship of electromagnetic radiation in the three most intense flares of solar cycle 23, more specifically, those of October 28, 2003, January 20, 2005, and September 7, 2005, to the acceleration and release of protons into interplanetary space. The impulsive phase of these flares lasted ～ 20 min and consisted of at least three energy release episodes, which differed by their manifestation in the soft (1–8 Å, GOES) and hard (>150 keV, INTEGRAL) X-ray ranges as well as at radio frequencies of 245 MHz and 8.8 GHz. The protons and electrons were accelerated in each episode, but with a different efficiency; the relativistic protons were accelerated only after 5–6min of impulsive-phase development after the onset of a coronal mass ejection. It is at this time that maximum hard X-ray fluxes were observed in the September 7, 2005 event, which exceeded severalfold those for the other two flares considered. We associate the record fluxes of protons with energies > 200MeV observed in the heliosphere in the September 7, 2005 event with the dynamics of the impulsive phase. The extreme intensities of the microwave emission in the October 28, 2003 and January 20, 2005 events were probably attributable to the high-energy electron trapping conditions and did not reflect the acceleration process.     

Composite suspended sediment particles in river systems: their incidence,dynamics and physical characteristics

Most of the existing data on the effective particle size characteristics of fluvial suspended sediment derive from instantaneous sampling methods that may not be representative of the overall suspended sediment loads. This presents difficulties when there is a need to incorporate effective particle size data into numerical models of floodplain sedimentation and sediment‐associated contaminant transfer. We have used a field‐based water elutriation apparatus (WEA) to assemble a large (36 flood) database on the time‐integrated nature of the effective and absolute particle size characteristics of suspended sediment in four subcatchments of the River Exe basin of southwest England. These catchments encompass a wide range of terrains and fluvial environments that are broadly representative of much of the UK and temperate, low relief northwest Europe. The WEA provides important data on the physical characteristics of composite particles that are not attainable using other methods. This dataset has allowed, for the first time, detailed interbasin comparisons of the time‐integrated particle size characteristics of suspended sediment and reliable estimates of the contribution of five effective size classes to the mean annual suspended sediment load of the study catchments. The suspended sediment load of each river is dominated by composite rather than primary particles, with, for example, almost 60% (by mass) of the sediment load of the River Exe at Thorverton transported as composite particles > 16 µm in size. All the effective size classes contain significant clay components. A key outcome of this study is the recognition that each catchment has a distinctive time‐integrated effective particle size signature. In addition, the time‐integrated effective particle size characteristics of the suspended loads in each of the catchments display much greater spatial variability than the equivalent absolute particle size distributions. This indicates that the processes producing composite particles vary significantly between these catchments, and this has important implications for our understanding of the dynamics of suspended sediment properties. Copyright © 2007 John Wiley & Sons, Ltd.     

Models of high velocity impacts into dust-covered ice: Application to Martian northern lowlands

The detection of fresh impact craters with bright floors and ejecta (arising from fresh clean water ice) in the northern lowlands of Mars (Byrne et al., 2009b, Science 325, 1674), together with observations of polygonal structures and evidence from the Phoenix probe, suggests that there are substantial water ice deposits just below the surface over large areas. Specifically in cases of the larger craters observed, the impacts themselves may have raised the temperature and the pressure of the water ice deposits locally to values which allow phase changes. In this paper, we use smoothed particle hydrodynamics to model hyper-velocity impacts. We estimate peak shock pressures in a solid water ice target covered by a layer of loose material, modeled by pre-damaged dunite. In addition, we account for the possibility of a thin layer of sub-surface water ice by using a three-layer model where the ice is surrounded by dunite. We find that the peak shock pressures reached in the simulated events are high enough to produce several 100-1000 kg of liquid water depending upon the impact parameters and the exact shock pressure needed for the phase change. A difficulty remains however in determining whether liquid is generated or whether a type of fluidized ice is produced (or indeed some combination of the two). We also note that the process can become rather complex as the number of layers increases because of reflections of the shock at sub-surface boundaries—a process which should lead to increased fluidization.     

海底双相随机介质声反射的SVM神经网络分类识别研究

考虑海底沉积介质为双相介质,为了更好地模拟实际海底底质的不均匀性,将随机介质理论引入双相介质理论。首先,通过基于随机-双相介质理论的高阶有限差分数值技术模拟计算海底底质分别为泥质砂、泥、泥质砾时的地震反射波信号。然后利用小波变换分别求取不同底质的一次反射波的包络作为其特征向量,最后利用基于粒子群智能算法优化的支持向量机神经网络对这些反射波信号进行分类识别。为了进一步考察所用方法的抗噪能力,对正演得到的海底底质反射波信号分别加入10%、30%、50%的高斯白噪音之后再进行分类,支持向量机仍然取得了较好的分类预测效果。基于上述正演模拟及分类识别方法的论证,提出了一套行之有效的微机软件模拟海底沉积物分类识别的一般化流程,这将有利于开展海底沉积物反射特征的进一步研究。     

河口悬浮粒子Stern电位的测定方法

首次提出了一种间接测定悬浮粒子Stern电位的方法,即通过测定悬浮粒子在不同浓度NaCl(0.01～0.40mol/L)溶液中的ζ-电位,按照Gouy-Chapman理论一般公式之对数变换式,用作图法计算悬浮粒子的Stern电位.用该方法研究了黄河口悬浮粒子Stern电位的变化规律.结果表明,在淡、咸水混合初期,悬浮粒子的Stern电位迅速减小;当S>3.5时,Stern电位不再随盐度的增大而变化.由此提出,对Ca2+和Ms2+的吸附是影响河口悬浮粒子Stern电位变化的主要因素.     

A PACKING GENERATION SCHEME FOR THE GRANULAR ASSEMBLIES WITH PLANAR ELLIPTICAL PARTICLES

In this paper, a new generator algorithm and a computer program PG2D is introduced for 2D numerical simulation of packing configuration in a granular material composed of elliptical particles of different a/b aspect ratios. Each elliptical particle is approximated by four connected arcs. The centre co-ordinates and radius of each arc and co-ordinates of connecting points can be determined from the formulae derived by entering the major axis length, 2a, and the eccentricity. The domain to be filled with particles can be a polygon of any shape. Given the size of the packing domain, geometrical information and numbers of particles to be generated, the packing location of each particle and the co-ordinates of contact points along with contact normal rose diagram can be generated as outputs. Simulation results show that this new algorithm can provide quite a reasonable packing model in accordance with the initial input required for the analysis of the mechanics of granular material. This generation scheme has the potential to cover packing generation and behaviour analysis of 3D sphere or ellipsoidal shaped granular materials. © 1997 by John Wiley & Sons, Ltd.