Application of micropolar plasticity to post failure analysis in geomechanics

A micropolar elastoplastic model for soils is formulated and a series of finite element analyses are employed to demonstrate the use of a micropolar continuum in overcoming the numerical difficulties encountered in application of finite element method in standard Cauchy–Boltzmann continuum. Three examples of failure analysis involving a deep excavation, shallow foundation, and a retaining wall are presented. In all these cases, it is observed that the length scale introduced in the polar continuum regularizes the incremental boundary value problem and allows the numerical simulation to be continued until a clear collapse mechanism is achieved. The issue of grain size effect is also discussed. Copyright © 2004 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.     

PROPAGATION OF 30—60 DAY LOW FREQUENCY OSCILLATIONS AND THEIR INFLUENCE UPON THE SUBTROPICAL WESTERLIES JET STREAM DURING NORTHERN HEMISPHERE WINTER

Based on daily ECMWF gridpoint data of two winters during 1981—1983 including an ENSOyear,propagation of low frequency oscillations(LFO)during Northern Hemisphere winters andtheir influences upon 30—60 day oscillations of the subtropical jet stream are studied with the sta-tistical methods as complex empirical orthogonal function(CEOF)and so on.Results show that inthe winter of a normal year(1981—1982),30—60 day oscillations in the subtropical zone aremainly in the northern and southern flanks of exit region of jet stream.In the ENSO year(1982—1983),they are mainly in the vicinity of entrance and exit regions of jet stream.Intraseasonalchanges of subtropical jet stream manifested themselves as latitudinal fluctuation or longitudinalprogression or regression of about 40 day period.There are marked differences between propagat-ing passages of low frequency modes responsible for changes of subtropical jet stream in the normalyear(1981—1982)and in the ENSO year(1982—1983).Changes of oscillation amplitude showobvious phases.In general,the one in late winter is stronger than that in early winter,strongestone occurs in February.     

A note on turbulent perturbations of salinity in a partially mixed estuary

Measurements of salinity perturbations in a partially mixed estuary have been used to evaluate the usefulness of an inductive salinometer and to determine some of the characteristics of the salinity perturbations. The salinometer performed satisfactorily under most conditions. Although internal wave like effects were present, the turbulence fluctuations were dominant. The salinity fluctuations and the turbulent fluxes $\text{sw}$ and $\text{su}$ were found to behave in a manner similar to the density fluctuations in a thermally stratified atmospheric boundary layer and a laboratory open channel flow. A quadrant analysis suggested that the contribution of each quadrant to the turbulent flux changed with Ri. The turbulence parameters ν and c_γ were found to decrease and increase respectively as Ri increases.     

图象的模型及其统计特性

本文对图象的统计特性及其适用模型作了分析。讨论了以前常用的平稳图象模型的缺陷,对实际图象的统计特性作了分析,认为它是非平稳的,不满足各态历经性,同时在空间上是高度相关的;讨论了相应的随机参数统计模型和描述性统计模型。指出作用于整幅图象的“全局”图象复原算法比每次只独立计算单个象元的“点”图象复原算法更为优越。     

Discrimination of sedimentary facies by association analysis

The ability of association analysis to discriminate sedimentary facies was tested on Purdy's modal analyses of modern sediments of the Great Bahama Bank. Purdy's data set has served in the past as a standard reference for evaluating various multivariate classification algorithms. In order to adapt Purdy's data to association analysis, the percent abundance of the 12 constituents was converted to binary form by dichotomizing each variable on its mean value. The results obtained by association analysis are virtually identical to those obtained by Purdy and other authors. The same four main sedimentary facies were discriminated; 86% of the samples were identically classified (97% when misclassified borderline cases are counted as matches); the total partition variance of the classification is only negligibly greater (4%); and the grouping of the variables yielded the same four groups. The rank order of the three division-attributes responsible for the sample classification is fines, oolites, and corals. Association analysis has been employed by other authors to differentiate meaningful facies groups in studies of ancient reef carbonates, modern reef sediments, and heavy minerals in stream sediments. In all these studies, the results were found to be compatible with those obtained by using the continuous quantitative measurements, indicating that qualitative binary data may often be sufficient for the purpose of facies discrimination in many branches of geology and that association analysis is an effective method for this purpose.     

The Statistical Analysis of Migration of Strong Earthquakes-Taking the North China Region as an Example

The migration of strong earthquakes is an important research topic because the migration phenomena reflect partly the seismic mechanism and involve the prediction of tendency of seismic activity. Research on migration of strong earthquakes has mostly focused on finding the phenomena. Some attempts on getting regularity were comparatively subjective. This paper suggests that there should be indices of migration in earthquake dataset and the indexes should have statistical meaning if there is regularity in the migration of strong earthquakes. In this study, three derivative attributes of migration, i.e., migration orientation, migration distance and migration time interval, were statistically analyzed. Results in the North China region show that the migration of strong earthquakes has statistical meaning. There is a dominant migration orientation (W by S to E by N), a dominant distance (≤100km and on the confines of 300～700km), and a dominant time interval (≤1a and on the confines of 3～4a). The results also show that the migration will differ slightly with different magnitude range or earthquake activity phase.