Implementation of a coupled plastic damage distinct lattice spring model for dynamic crack propagation in geomaterials

This paper studies dynamic crack propagation by employing the distinct lattice spring model (DLSM) and 3‐dimensional (3D) printing technique. A damage‐plasticity model was developed and implemented in a 2D DLSM. Applicability of the damage‐plasticity DLSM was verified against analytical elastic solutions and experimental results for crack propagation. As a physical analogy, dynamic fracturing tests were conducted on 3D printed specimens using the split Hopkinson pressure bar. The dynamic stress intensity factors were recorded, and crack paths were captured by a high‐speed camera. A parametric study was conducted to find the influences of the parameters on cracking behaviors, including initial and peak fracture toughness, crack speed, and crack patterns. Finally, selection of parameters for the damage‐plasticity model was determined through the comparison of numerical predictions and the experimentally observed cracking features.     

Dust supply from river floodplains: the case of the lower Huang He (Yellow River) recorded in a loess–palaeosol sequence from the Mangshan Plateau

The Mangshan Plateau is located on the south bank of the Huang He (Yellow River) just west of the city of Zhengzhou, well outside the Loess Plateau in central China. Mixing models of the grain‐size data indicate that the loess deposits are mixtures of three loess components. Comparison of the mixing model with existing models established for a series of loess–palaeosol sequences from the Loess Plateau indicates that the Mangshan loess has been supplied from a proximal dust source, the Huang He floodplain, during major dust outbreaks. The high accumulation rates, the composition of the loess components, and especially the high proportions of a sandy loess component support this. Owing to the exceptionally high accumulation rates, the Mangshan grain size, magnetic susceptibility and carbonate records provide a high‐resolution archive of environmental and climate change. Copyright © 2008 John Wiley & Sons, Ltd.     

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.     

Dynamic analysis of three-dimensional marine cables

Marine cables are widely used in a variety of offshore activities. In this paper, a numerical approach is presented which is capable of predicting the dynamic behaviour of three-dimensional marine cables. The solution of this problem is of practical importance in the estimation of dybamic loading and motion, and thus has direct application to the enhancement of safety and the effectiveness of the offshore activities. The approach is based upon the lump-mass-and-spring model and the finite difference method. Mathematical analysis is provided on the basic characteristics of the model and the numerical stability. A numerical example is given to demonstrate the validity of the approach.     

Optical Flash of GRB 990123: Constraints on the Physical Parameters of the Reverse Shock

The optical flash accompanying GRB 990123 is believed to be powered by the reverse shock of a thin shell. With the best-fit physical parameters for GRB 990123 and the assumption that the parameters in the optical flash are the same as in the afterglow, we show that: 1) the shell is thick rather than thin, and we have provided the light curve for the thick shell case which coincides with the observation; 2) the theoretical peak flux of the optical flash accounts for only 3×10~-4 of the observed. In order to remove this discrepancy, the physical parameters, the electron energy and magnetic ratios, εe and εB, should be 0.61 and 0.39, which are very different from their values for the late afterglow.     

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.     

A numerical method for differential equations of second order

A new algorithm of order five is presented for the solution of the initial value problem where the system of ordinary differential equations is of second order and does not contain the first derivative.     

Runoff responses to afforestation in a watershed of the Loess Plateau,China

Afforestation has been suggested as a means of improving soil and water conservation in north‐western China, especially on the Loess Plateau. Understanding of the hydrological responses to afforestation will help us develop sustainable watershed management strategies. A study was conducted during the period of 1956 to 1980 to evaluate runoff responses to afforestation in a watershed on the Loess Plateau with an area of 1·15 km², using a paired watershed approach. Deciduous trees, including locust (locusta L.), apricot (praecox L.) and elm (ulmus L.), were planted on about 80% of a treated watershed, while a natural grassland watershed remained unchanged. It was estimated that cumulative runoff yield in the treated watershed was reduced by 32% as a result of afforestation. A significant trend was also observed that shows annual runoff reduction increases with the age of the trees planted. Reduction in monthly runoff occurred mainly from June to September, which was ascribed to greater rainfall and utilization by trees during this period. Afforestation also resulted in reduction in the volume and peak flow of storm runoff events in the treated watershed with greater reduction in peak flow. Copyright © 2003 John Wiley & Sons, Ltd.