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.     

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.     

广西中三叠统浊流流向及坡向

通过对广西中三叠统浊积岩系中溢堤流沉积的识别,结合同一鲍马层序Tc段斜层理倾向与其底部槽模走向不一致的分析,以及横向冲刷模、U形障碍冲刷痕等冲蚀构造、大型滑动模等压刻构造的认识,筛选出宽阔水道充填沉积、扇状体及分支水道沉积中底模构造的指向数据,来指示浊流主体流向,并区分出溢堤流流向。从而确定了本区中三叠世时的浊流主体流向为NW向,平均为300°。根据滑塌沉积中滑塌包卷层轴面、滑塌面的测定,确定凤山县金牙处的古坡向平均为322°,与浊流主体流向基本一致。还认为中三叠世时右江盆地沉积物源来自桂东和桂东南蚀源区,海水由其南的越南境内进入。     

Can Brans-Dicke universes in five dimensions be interpreted as 4-dimensional general relativity cosmologies?

The 5-dimensional Jordan-Brans-Dicke cosmologies in vacuum are found for the Bianchi type I metric, their relation with general relativity cosmologies is studied. Two solutions are possible, both produce effective pressure and energy density in the 4-dimensional G.R.-universes. One is a power-law relation, with two cases, the first one is forp _eff=_eff and the other forp _eff=_eff(– 1 < < 1) has a behaviour as the open flat universe. The second solution is an exponential only valid forp _eff=–_eff. In all cases the three-space expansion reaches infinity ast and the fifth dimension can be made to decrease approaching zero. The scalar field can increase or decrease with time.     

A possible mechanism of formation of faculae

We propose a new heating mechanism of faculae. We think that the formation of faculae is a result of the Joule dissipation of the Hall current generated by the interaction of the convection field of granules in an active region and the inter-granular magnetic field. For a region to generate effectively Hall current, its characteristic length must be such that the magnetic Reynolds number is less than 1. The equation of energy balance in the facula region is

$\text{16σT}{}^3{}_{\mathrm{p}}\text{(T}{}_{\mathrm{l}}\text{? T}{}_{\mathrm{p}}\text{)}{}^{\mathrm{n}}\text{HP}{}_{\mathrm{s}}{}^{\mathrm{a}}\text{H}{}^{\mathrm{?}}\text{=}\text{Qn}{}_{\mathrm{s}}\text{m}{}_{\mathrm{i}}\text{u}{}_{\mathrm{x}}{}^2\text{(τ}{}^2{}_{\mathrm{in}}\text{ω}{}_{\mathrm{i}}\text{)}$

.For five observational models of faculae, we calculated the corresponding velocity fields, and the results are in basic agreement with the observed fields. The present mechanism explains the dependence of the facula brightness on the magnetic and velocity fields, the apparent distribution of the faculae on the solar disk and suggest a possible interpretation of the five structures of faculae.     

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.     

Comparative analyses of East Texas forest cover maps generated from Landsat and AVHRR data

Comparing satellite data derived map products are affected by differences in data characteristics, image acquisition dates, processing techniques, and classification schemes used for assigning pixels to a thematic class. By comparing two forest maps generated from Landsat Enhanced Thematic Mapper Plus (ETM+) and Advanced Very High Resolution Radiometer (AVHRR) images acquired on the same day, and processed using identical classification scheme and methods these differences were minimized. The ETM+ derived map had higher classification accuracy values and more precise area estimates than the AVHRR derived map. In the ETM+ derived map, 87 of the 599 verification data were misclassified, whereas in the AVHRR derived map, 155 of the 469 verification data were misclassified. Detailed error analyses by land cover class revealed that a land use based definition of forest accounted for 74% (64 out of 87) and 57% (89 out of 155) of the classification errors in ETM+ and AVHRR derived maps, respectively.