The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Three-Dimensional Network Model And Its Parameter Calibration

A material model, whose framework is parallel spring-bundles oriented in 3-D space, is proposed. Based on a discussion of the discrete schemes and optimum discretization of the solid angles, a 3-D network cell consisted of one-dimensional components is developed with its geometrical and physical parameters calibrated. It is proved that the 3-D network model is able to exactly simulate materials with arbitrary Poisson ratio from 0 to 1/2, breaking through the limit that the previous models in the literature are only suitable for materials with Poisson ratio from 0 to 1/3. A simplified model is also proposed to realize high computation accuracy within low computation cost. Examples demonstrate that the 3-D network model has particular superiority in the simulation of short-fiber reinforced composites.

General Solution to Two-Dimensional Nonslipping JKR Model with a Pulling Force in an Arbitrary Direction

In this paper, a generalized JKR model is investigated, in which an elastic cylinder adhesively contacts with an elastic half space and the contact region is assumed to be perfect bonding. An external pulling force is acted on the cylinder in an arbitrary direction. The contact area changes during the pull-off process, which can be predicted using the dynamic Griffith energy balance criterion as the contact edge shifts. Full coupled solution with an oscillatory singularity is obtained and analyzed by numerical calculations. The effect of Dundurs' parameter on the pull-off process is analyzed, which shows that a nonoscillatory solution can approximate the general one under some conditions, i.e., larger pulling angle (pi/2 is the maximum value), smaller a/R or larger nondimensional parameter value of Delta gamma/E*R. Relations among the contact half width, the external pulling force and the pulling angle are used to determine the pull-off force and pull-off contact half width explicitly. All the results in the present paper as basic solutions are helpful and applicable for experimenters and engineers.

Embedded cell model of three-dimensional two-phase composites

An embedded cell model is presented to obtain the effective elastic moduli and the elastic-plastic stress-strain relations of three-dimensional two-phase particulate composites. Each cell consists of an ellipsoidal inclusion surrounded by a finite ellipsoidal matrix that embedded in an infinite matrix. When both matrix and particle are elastic, the effective elastic moduli are derived which is an exact analytic formula without any simplified approximation that can be expressed in an explicit form. Further, the elastic-plastic stress-strain relations are obtained for spherical cells and oblate spheroid cells, in which the matrix is elastic and the particle is elastic-plastic. In addition, the macroscopic elastic-plastic constitutive relation of particle reinforced composites (PRC) is investigated by a systematic approach [1] in which the matrix is elastic-plastic and the particle is elastic.

Application of Three-Dimensional Discrete Element Face-to-Face Contact Model with Fissure Water Pressure to Stability Analysis of Landslide in Panluo Iron Mine

Three-dimensional discrete element face-to-face contact model with fissure water pressure is established in this paper and the model is used to simulate three-stage process of landslide under fissure water pressure in the opencast mine, according to the actual state of landslide in Panluo iron mine where landslide happened in 1990 and was fathered in 1999. The calculation results show that fissure water pressure on the sliding surface is the main reason causing landslide and the local soft interlayer weakens the stability of slope. If the discrete element method adopts the same assumption as the limit equilibrium method, the results of two methods are in good agreement; while if the assumption is not adopted in the discrete element method, the critical phi numerically calculated is less than the one calculated by use of the limit equilibrium method for the same C. Thus, from an engineering point of view, the result from the discrete element model simulation is safer and has more widely application since the discrete element model takes into account the effect of rock mass structures.

Two-Dimensional Kinematic Wave Model of Overland-Flow

A two-dimensional kinematic wave model was developed for simulating runoff generation and flow concentration on an experimental infiltrating hillslope receiving artificial rainfall. Experimental observations on runoff generation and flow concentration on irregular hillslopes showed that the topography of the slope surface controlled the direction and flow lines of overland flow. The model-simulated results satisfactorily compared with experimental observations. The erosive ability of the concentrated flow was found to mainly depend on the ratio of the width and depth of confluent grooves.

Two-dimensional fluid model simulation of bell jar top inductively coupled plasma

In the present paper, argon (Ar) plasmas in a bell jar inductively coupled plasma (ICP) source are systematically studied over pressures from 5 to 20 mtorr and power inputs from 0.2 to 0.5 kW. In this study, both a two-dimensional (2-D) fluid model simulation and global model calculation are compared, The 2-D fluid model simulation with a self-consistent power deposition is developed to describe the Ar plasma behavior as well as predict the plasma parameter distributions, Finally, a quantitative comparison between the global model and the fluid model is made to test their validity.

The Displacement and Strain Field of Three-Dimensional Rheologic Model of Earthquake Preparation

Based on the three-dimensional elastic inclusion model proposed by Dobrovolskii, we developed a rheological inclusion model to study earthquake preparation processes. By using the Corresponding Principle in the theory of rheologic mechanics, we derived the analytic expressions of viscoelastic displacement U(r, t) , V(r, t) and W(r, t), normal strains epsilon(xx) (r, t), epsilon(yy) (r, t) and epsilon(zz) (r, t) and the bulk strain theta (r, t) at an arbitrary point (x, y, z) in three directions of X axis, Y axis and Z axis produced by a three-dimensional inclusion in the semi-infinite rheologic medium defined by the standard linear rheologic model. Subsequent to the spatial-temporal variation of bulk strain being computed on the ground produced by such a spherical rheologic inclusion, interesting results are obtained, suggesting that the bulk strain produced by a hard inclusion change with time according to three stages (alpha, beta, gamma) with different characteristics, similar to that of geodetic deformation observations, but different with the results of a soft inclusion. These theoretical results can be used to explain the characteristics of spatial-temporal evolution, patterns, quadrant-distribution of earthquake precursors, the changeability, spontaneity and complexity of short-term and imminent-term precursors. It offers a theoretical base to build physical models for earthquake precursors and to predict the earthquakes.

Two-Dimensional Hillslope Scale Soil Erosion Model

On a hillslope, overland flow first generates sheet erosion and then, with increasing flux, it causes rill erosion. Sheet erosion (interrill erosion) and rill erosion are commonly observed to coexist on hillslopes. Great differences exist between both the intensities and incidences of rill and interrill erosion. In this paper, a two-dimensional rill and interrill erosion model is developed to simulate the details of the soil erosion process on hillslopes. The hillslope is treated as a combination of a two-dimensional interrill area and a one-dimensional rill. The rill process, the interrill process, and the joint occurrence of rill and interrill areas are modeled, respectively. Thus, the process of sheet flow replenishing rill flow with water and sediment can be simulated in detail, which may possibly render more truthful results for rill erosion. The model was verified with two sets of data and the results seem good. Using this model, the characteristics of soil erosion on hillslopes are investigated. Study results indicate that (1) the proposed model is capable of describing the complex process of interrill and rill erosion on hillslopes; (2) the spatial distribution of erosion is simulated on a simplified two-dimensional hillslope, which shows that the distribution of interrill erosion may contribute to rill development; and (3) the quantity of soil eroded increases rapidly with the slope gradient, then declines, and a critical slope gradient exists, which is about 15-20 degrees for the accumulated erosion amount.

BRST invariant theory of a generalized 1+1 Dimensional nonlinear sigma model with topological term

We give a generalized Lagrangian density of 1 + 1 Dimensional O( 3) nonlinear sigma model with subsidiary constraints, different Lagrange multiplier fields and topological term, find a lost intrinsic constraint condition, convert the subsidiary constraints into inner constraints in the nonlinear sigma model, give the example of not introducing the lost constraint. N = 0, by comparing the example with the case of introducing the lost constraint, we obtain that when not introducing the lost constraint, one has to obtain a lot of various non-intrinsic constraints. We further deduce the gauge generator, give general BRST transformation of the model under the general conditions. It is discovered that there exists a gauge parameter beta originating from the freedom degree of BRST transformation in a general O( 3) nonlinear sigma model, and we gain the general commutation relations of ghost field.

Friction phenomena and phase transition in the underdamped two-dimensional Frenkel-Kontorova model

Locked-to-sliding phase transition has been studied in the driven two-dimensional Frenkel-Kontorova model with the square symmetric substrate potential. It is found that as the driving force increases, the system transfers from the locked state to the sliding state where the motion of particles is in the direction different from that of driving force. With the further increase in driving force, at some critical value, the particles start to move in the direction of driving force. These two critical forces, the static friction or depinning force, and the kinetic friction force for which particles move in the direction of driving force have been analyzed for different system parameters. Different scenarios of phase transitions have been examined and dynamical phases are classified. In the case of zero misfit angle, the analytical expressions for static and kinetic friction force have been obtained.

Seasonal and interannual variability of primary and export production in the South China Sea: a three-dimensional physical-biogeochemical model study

To investigate the seasonal and interannual variations in biological productivity in the South China Sea (SCS), a Pacific basin-wide physical - biogeochemical model has been developed and used to estimate the biological productivity and export flux in the SCS. The Pacific circulation model, based on the Regional Ocean Model Systems (ROMS), is forced with daily air-sea fluxes derived from the NCEP (National Centers for Environmental Prediction) reanalysis between 1990 and 2004. The biogeochemical processes are simulated with a carbon, Si(OH)(4), and nitrogen ecosystem (CoSiNE) model consisting of silicate, nitrate, ammonium, two phytoplankton groups (small phytoplankton and large phytoplankton), two zooplankton grazers (small micrograzers and large mesozooplankton), and two detritus pools. The ROMS-CoSiNE model favourably reproduces many of the observed features, such as ChI a, nutrients, and primary production (PP) in the SCS. The modelled depth-integrated PP over the euphotic zone (0-125 m) varies seasonally, with the highest value of 386 mg C m (-2) d (-1) during winter and the lowest value of 156 mg C m (-2) d (-1) during early summer. The annual mean value is 196 mg C m (-2) d (-1). The model-integrated annual mean new production (uptake of nitrate), in carbon units, is 64.4 mg C m (-2) d (-1) which yields an f-ratio of 0.33 for the entire SCS. The modelled export ratio (e-ratio: the ratio of export to PP) is 0.24 for the basin-wide SCS. The year-to-year variation of biological productivity in the SCS is weaker than the seasonal variation. The large phytoplankton group tends to dominate over the smaller phytoplankton group, and likely plays an important role in determining the interannual variability of primary and new production.