Bond rolling resistance and its effect on yielding of bonded granulates by DEM analyses

A discrete element modelling of bonded granulates and investigation on the bond effect on their behaviour are very important to geomechanics. This paper presents a two‐dimensional (2‐D) discrete element theory for bonded granulates with bond rolling resistance and provides a numerical investigation into the effect of bond rolling resistance on the yielding of bonded granulates. The model consists of mechanical contact models and equations governing the motion of bonded particles. The key point of the theory is that the assumption in the original bond contact model previously proposed by the authors (55th CSCE‐ASCE Conference, Hamilton, Ont., Canada, 2002; 313–320; J. Eng. Mech. (ASCE) 2005; 131 (11):1209–1213) that bonded particles are in contact at discrete points, is here replaced by a more reliable assumption that bonded particles are in contact over a width. By making the idealization that the bond contact width is continuously distributed with the normal/tangential basic elements (BE) (each BE is composed of spring, dashpot, bond, slider or divider), we establish a bond rolling contact model together with bond normal/tangential contact models, and also relate the governing equations to local equilibrium. Only one physical parameter β needs to be introduced in the theory in comparison to the original bond discrete element model. The model has been implemented into a 2‐D distinct element method code, NS2D. Using the NS2D, a total of 86 1‐D, constant stress ratio, and biaxial compressions tests have been carried out on the bonded granular samples of different densities, bonding strengths and rolling resistances. The numerical results show that: (i) the new theory predicts a larger internal friction angle, a larger yielding stress, more brittle behaviour and larger final broken contact ratio than the original bond model; (ii) the yielding stress increases nonlinearly with the increasing value of β, and (iii) the first‐yield curve (initiation of bond breakage), which define a zone of none bond breakage and which shape and size are affected by the material density, is amplified by the bond rolling resistance in analogous to that predicted by the original bond model. Copyright © 2006 John Wiley & Sons, Ltd.     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Force-free magnetic fields on curvilinear coordinate surfaces

To seek nonlinear solutions of force-free magnetic fields, some symmetries or approximations are usually invoked. We consider magnetic fields lying on coordinate surfaces of an orthogonal curvilinear coordinate system. We conclude that only fields on parallel planes or spherical shells can be expressed in the form provided by Low in 1980s. These force-free fields are stable against small perturbations with rigid boundaries. Fields on cylindrical shells are also considered.     

On the microwave millisecond spike emission and its associated phenomena during the impulsive phase of large solar flare

We propose a model of two acceleration regions, which can explain, on the basis of microwave maser caused by a “hollow-beam” distribution of electrons, the presence of millisecond spikes in the event of 1981 May 16 and their absence in the event of 1981 October 12, and the enhanced continuous emission in the latter. We have also uncovered relations among the features, e.g. the Type III_G, Type IV_DCIM and hard x-ray bursts, that accompany the microwave millisecond spikes during the impulsive phase of a large flare.     

The global joule heat production rate and the AE index

The accuracy of the AE index as a measure of the joule heat production rate is examined for a typical substorm event on 18 March 1978 by estimating the global joule heat production rate, as a function of time, using data obtained from the IMS six meridian chains. In spite of the fact that the AE index had had an initial slow growth which was followed by a rapid growth, the joule heat production rate attained a high level during the slow growth and thus their variations were considerably different from those of the AE index. Therefore, although the AE index is statistically linearly related to the global joule heat production rate, one should be cautious in assuming that details of time variations of the AE index during individual events represent those of the joule heat production rate.     

Triassic Nb-enriched basalts, magnesian andesites, and adakites of the Qiangtang terrane (Central Tibet): evidence for metasomatism by slab-derived melts in the mantle wedge

New chronological, geochemical, and isotopic data are reported for Triassic (219–236 Ma) adakite-magnesian andesite-Nb-enriched basaltic rock associations from the Tuotuohe area, central Qiangtang terrane. The adakites and magnesian andesites are characterized by high Sr/Y (25–45), La/Yb (14–42) and Na₂O/K₂O (12–49) ratios, high Al₂O₃ (15.34–18.28 wt%) and moderate to high Sr concentrations (220–498 ppm) and ε_ND (t) (+0.86 to +1.21) values. Low enrichments of Th, Rb relative to Nb, and subequal normalized Nb and La contents, and enrichments of light rare earth elements combine to distinguish a group of Nb-enriched basaltic rocks (NEBs). They have positive ε_ND (t) (+2.57 to +5.16) values. Positive correlations between Th, La and Nb and an absence of negative Nb anomalies on mantle normalized plots indicate the NEBs are products of a mantle source metasomatized by a slab melt rather than by hydrous fluids. A continuous compositional variation between adakites and magnesian andesites confirms slab melt interaction with mantle peridotite. The spatial association of the NEBs with adakites and magnesian andesites define an “adakitic metasomatic volcanic series” recognized in many demonstrably subduction-related environments (e.g., Mindanao arc, Philippines; Kamchatka arc, Russia; and southern Baja California arc, Mexico). The age of the Touhuohe suite, and its correlation with Triassic NEB to the north indicates that volcanism derived from subduction-modified mantle was abundant prior to 220 Ma in the central Qiangtang terrane.     

海洋走航式单边连续观测地震折射波Δt解释方法--1种新的浅层折射单支时距曲线解释方法

在浅层工程折射地震测量中 ,通常采用相遇观测系统采集数据和相应的解释方法进行解释 ;对于单边观测采集方法由于与其相应的解释方法较少而少被采用。本文提出 1种方法简单、精度可靠、可直接利用同一点相邻 2排列 2次观测的初至折射时间并适合于海洋走航式单边连续观测系统的折射初至资料解释方法。其精度基本与同观测段的采用相遇观测的差异时距方法解释的精度十分接近。该方法即适合于人工解释又适合于计算机自动解释。     

雪宝顶碱性花岗岩岩石地球化学与成矿控制

对雪宝顶碱性花岗岩具有高丰度的W、Sn、Be含量和W、Sn矿化花岗岩的微量元素组合特征进行了分析，结果显示：稀土元素具有∑REE含量与LREE／HREE比值偏低和负Eu异常的特点；元素比值K／Rb，U／Th，Nb／Ta的规律性变化以及矿脉产状特征等表明W-(Sn)-Be脉状矿床形成于岩浆演化晚期，与岩浆期后的热液活动有关。     

青岛文昌鱼的生活史──年龄、生长和死亡研究

于１９８９年８月－１９９０年７月对青岛文昌鱼的生活史进行研究。结果表明：（１）青岛文昌鱼种群由４个齿组组成，它至少生产４ａ以上，甚至可达６ａ左右。（２）在Ⅰ、Ⅱ齿时生长较快，平均年生长１３－１５ｍｍ，Ⅲ龄后生长速度减慢，即使同一年龄组的文晶鱼，在不同季节其生长速度也不同。（３）在冬春采集的文昌鱼中Ⅲ龄组的峰值较小。（４）繁殖季节过后，老龄文昌鱼因衰老而死亡。     

一种腐屑食性鱼类对生态系统氮循环的影响

研究表明 :腐屑食性鱼类——梭鱼对环境中有机氮的摄入能力随体重的增加而增加 ,随温度的升高 ,体重系数显著减少。摄入氮 (CN,mg/ d)与体重 (W,g)和温度 (T,℃ )的数值关系可用如下模型表示 :CN=0 .6 370 W1 .0 642 ln(3.70 0 8T)。体重为 1 0 0 g的梭鱼每日可消耗环境中 332 mg氮有机质。梭鱼氮的吸收效率主要受温度的影响 ,而氮的转化效率与体重呈显著的正相关关系。随体重的增加 ,排泄氮所占比例降低