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.     

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.     

Fracture energy release and size effect in borehole breakout

The paper presents a simple approximate analytical solution of the remote stresses that cause the collapse of a borehole or other circular cylindrical cavity in an infinite elastic space. Regions of parallel equidistant splitting cracks are assumed to form on the sides of the cavity. Their boundary is assumed to be an ellipse of a growing horizontal axis, the other axis remaining equal to the borehole diameter. The slabs of rock between the splitting cracks are assumed to buckle as slender columns, and their post-critical stress is considered as the residual stress in the cracked rock. The buckling of these slab columns is assumed to be resisted not only by their elastic bending stiffness but also shear stresses produced on rough crack faces by relative shear displacements. The energy release from the infinite medium caused by the growth of the elliptical cracking region is evaluated according to Eschelby's theorem. This release is set equal to the energy dissipated by the formation of all the splitting cracks, which is calculated under the assumption of constant fracture energy. This yields the collapse stress as a function of the elastic moduli, fracture energy, ratio of the remote principal stresses, crack shear resistance characteristic and borehole diameter. The collapse stress as a function of crack spacing is found to have a minimum, and the correct crack spacing is determined from this minimum. For small enough diameters, the crack spacing increases as the (4/5)-power of the borehole diameter, while for large enough diameters a constant spacing is approached. In contrast to plastic solutions, the breakout stress exhibits a size effect, such that for small enough diameters the breakout stress decreases as the (? 2/5)-power of the borehole diameter, while for large enough diameters a constant limiting value is approached. Finally, some numerical estimates are given and the validity of various simplifying assumptions made is discussed.     

Energetic solar electrons in the interplanetary medium

ISEE-3 measurements extending down to 2 keV energy have provided a new perspective on energetic solar electrons in the interplanetary medium. Impulsive solar electron events are observed, on average, several times a day near solar maximum, with 40% detected only below 15 keV. The electron energy spectra have a nearly power-law shape extending smoothly down to 2 keV, indicating that the origin of these events is high in the corona. These coronal flare-like events often produced ³He-rich particle events.In large solar flares which accelerate electrons and ions to relativistic energies, the electron spectrum appears to be modified by a second acceleration which results in a double power-law shape above 10 keV with a break near 100 keV and flattening from 10–100 keV. Large flares result in long-lived (many days) streams of outflowing electrons which dominate the interplanetary fluxes at low energies. Even in the absence of solar activity, significant fluxes of low energy electrons flow out from the Sun.Solar type-III radio bursts are produced by the escaping 2–10² keV electrons through a beam-plasma instability. The detailed ISEE-3 measurements show that electron plasma waves are generated by the bump-on-tail distribution created by the faster electrons running ahead of the slower ones. These plasma waves appear to be converted into radio emission by nonlinear wave-wave interactions.     

LZY1—1型光电浊度仪

本文介绍了LZY_(1-1)型光电浊度仪的设计原理,标定方法及使用结果,该仪器用经典法校准后,能快速测量大批量样品的悬浮物含量。     

Paraconchoecia diacanthus n. sp., a new species of planktonic Ostracoda

Ｐａｒａｃｏｎｃｈｏｅｃｉａｄｉａｃａｎｔｈｕｓｎ．ｓｐ．,ａｎｅｗｓｐｅｃｉｅｓｏｆｐｌａｎｋｔｏｎｉｃＯｓｔｒａｃｏｄａ￥ＣｈｅｎＲｕｉｘｉａｎｇａｎｄＬｉｎＪｉｎｇｈｏｎｇ（ＲｅｃｅｉｖｅｄＡｐｒｉｌ１６,１９９３１ａｃｃｅｐｔｅｄＪｕｌｙ２０...     

Interfacial geochemistry of barium in suspension and sediment from Zhujiang Estuary

Chemical forms of Ba are determined in samples of suspension and sediment from the Zhujiang Estuary. Their interfacial geochemical processes are discussed. Total suspended Ba content is between 2. 4 and 40. 4 μg/L, and mostly exists in the crystalline form (43. 5%), secondly in the Fe-Mn oxidative form (23.1%). Percentages of organic, carbonate and exchangeable forms are 14. 8%, 11. l%and 7.4%, respectively. Total content of Ba in the sediment is between 158. 6 and 48. 0 ug/g. Percentages of crystalline form, Fe-Mn oxidative form and carbonate form are 78. 4% ,13. 5%and 8. 2%, respectively. Organic and exchangeable forms are not detected. The study on the mechanism of interfacial movement suggests that the salinity range of 10 is the turning point for the varied distribution of Ba. The subsidence of crystalline form affects the decrease of Ba content in suspension. The decrease mostly takes place in the salinity range lower than 10, which corresponds to the high value of Ba content in the sediment. The     

On the recovery of surface wave by pressure transfer function

Except the commonly selected pressure transfer function derived from the linear wave theory, a previous study on the pressure transfer function for recovering surface wave from underwater pressure transducer suggested that the pressure transfer function is a function of frequency parameter only. With careful analysis, this study showed that the pressure transfer function should include a transducer submergence parameter as that given by the linear theory. It was found that the previously suggested empirical formula should be restricted to measurements with the pressure transducer close to the surface; otherwise overestimation of wave height would result. Field measurements were carried out with an acoustic wave gauge and a synchronized pressure transducer located at various depths with submergence parameter close to 1 (near the sea floor). It was shown that the previous one-parameter empirical formula might overestimate the significant wave height by more than 30%. This study found that with deep-water wave bursts excluded, the transfer function based on the linear wave theory provided a fairly good estimation on the significant wave heights, with an average deviation of 3.6%.     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.