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.     

Fluid infiltration during contact metamorphism of interbedded marble and calc-silicate hornfels, Twin Lakes area, central Sierra Nevada, California

Abstract In the Twin Lakes area, central Sierra Nevada, California, most contact metamorphosed marbles contain calcite + dolomite + forsterite ± diopside ± phlogopite ± tremolite, and most calc-silicate hornfelses contain calcite + diopside + wollastonite + quartz ± anorthite ± K-feldspar ± grossular ± titanite. Mineral-fluid equilibria involving calcite + dolomite + tremolite + diopside + forsterite in two marble samples and wollastonite + anorthite + quartz + grossular in three hornfels samples record P± 3 kbar and T± 630° C. Various isobaric univariant assemblages record CO₂-H₂O fluid compositions of χCO₂= 0.61–0.74 in the marbles and χCO₂= 0.11 in the hornfelses. Assuming a siliceous dolomitic limestone protolith consisting of dolomite + quartz ° Calcite ± K-feldspar ± muscovite ± rutile, all plausible prograde reaction pathways were deduced for marble and hornfels on isobaric T-XCO₂ diagrams in the model system K₂O-CaO-MgO-Al₂O₃-SiO₂-H₂O-CO₂. Progress of the prograde reactions was estimated from measured modes and mass-balance calculations. Time-integrated fluxes of reactive fluid which infiltrated samples were computed for a temperature gradient of 150 °C/km along the fluid flow path, calculated fluid compositions, and estimated reaction progress using the mass-continuity equation. Marbles and hornfelses record values in the range 0.1–3.6 × 10⁴ cm³/cm² and 4.8–12.9 × 10⁴ cm³/cm², respectively. For an estimated duration of metamorphism of 10⁵ years, average in situ metamorphic rock permeabilities, calculated from Darcy's Law, are 0.1–8 × 10^?6 D in the marbles and 10–27 × 10^?6 D in the hornfelses. Reactive metamorphic fluids flowed up-temperature, and were preferentially channellized in hornfelses relative to the marbles. These results appear to give a general characterization of hydrothermal activity during contact metamorphism of small pendants and screens (dimensions ± 1 km or less) associated with emplacement of the Sierra Nevada batholith.     

Hα Photometry of Two Contact Binaries

We carried out optical and Hα photometry of two contact binaries (V861 Herculis, EQ Tauri). The light curve modeling revealed stellar spots in both contact systems and strong Hα excess in the position of the observed stellar spots. A correlation was found between the V−R and R−H_α colour indices of V861 Her.     

On the Period Variation of the W UMa-type Contact Binary V502 Ophiuchi

The variation in the orbital period of the W UMa type contact binary V502 Oph is analyzed. The orbital period exhibits a wavelike variation with a periodicity of 23.0 years and an amplitude of △P = 1.24×10~(-6) days superimposed on secular decrease of dP/dt = 1.68×10-7 day per year. The long-term decrease may be accompanied by the contraction of the secondary at a rate of 83 m per year and a mass transfer rate from the primary to the secondary of 4.28×10~8 M per year. The short-term oscillation may be explained by the presence of a third component. Orbital elements of the third body and its possible mass are presented.     

海水冷却技术

城市用水中约 5 0 %是工业冷却水。海水代替淡水作为工业冷却用水——海水冷却 ,是解决我国沿海城市和地区淡水资源危机问题的重要途径之一。海水冷却技术包括海水直流冷却技术和海水循环冷却技术 ,目前我国基本具备了具有自主知识产权的海水(直流、循环 )冷却关键技术。寻求政府在资金、政策等方面的大力支持 ,加强有关行业、院所、企业的“强强联合”,适时借鉴国外先进的海水冷却塔技术 ,加快产业化技术示范 ,是推进海水循环冷却技术快速发展的必由之路     

Random wave-current interactions in water of varying depth

Refraction of incoherent random gravity waves with currents and bottom topography results in spatial variations in the spectral characteristics of the free surface. Prediction of such variations based on the radiation transfer equation is in a simple analytic form for the case of one dimensional inhomogeneities in currents and topography. This analytic form is examined in terms of two-dimensional wave number- and polar frequency-direction spectra along the associated dynamic and kinematic constraints relevant to wave breaking and reflection. Results are specialized to the simplest case of horizontal shear currents in deep and shallow water with explicit examples to illustrate the relative and combined effects of currents and topography on free surface spectra.     

基于坐标归一化和奇异值分解的直接线性变换解法

直接线性变换方法是数字摄影测量和计算机视觉领域最常用的解析处理方法之一,特别适用于未检校数码相机的三维测量,但直接线性参数间的相关性、物方控制的约束和设计矩阵元素数量级的较大差异,均可导致法矩阵严重病态,从而影响解的稳定性。本文借鉴改进的八点基本矩阵估计算法,采用基于坐标归一化和奇异值分解的解法,即首先将像点和物点坐标进行相似变换得到归一化坐标后组成法矩阵,其次利用矩阵奇异值分解方法代替常规的最小二乘方法,模拟和真实数据表明,此方法可以有效提高解算精度和稳定性。     

PRELIMINARY NUMERICAL STUDY OF TOPOGRAPHIC EFFECTS OF THE TIBETAN PLATEAU ON SURFACE DIRECT RADIATION*

To improve the simulating ability of a model,this paper presents a scheme of calculating direct radiation at land surface with topography in the model.A numerical study is conducted for the topographic effects of the Tibetan Plateau on the direct radiation using NCEP terrain data.Results show that,after taking account into the topographic radiation effect,the regional average of the radiation over the Plateau obviously increases in the local early morning and late afternoon,but changes less around noon.The effect is stronger in winter than that in summer.And heterogeneous topography has also affected the distribution of the radiation in this area.A simple numerical experiment shows that considering the effect will lead ground temperature to increase on the slope having more sunshine,and vice versa.