Weak Zone Related Seismic Cycles in Progressive Failure Leading to Collapse in Brittle Crust

Until quite recently our understanding of the basic mechanical process responsible for earthquakes and faulting was not well known. It can be argued that this was partly a consequence of the complex nature of fracture in crust and in part because evidence of brittle phenomena in the natural laboratory of the earth is often obliterated or obscured by other geological processes. While it is well understood that the spatial and temporal complexity of earthquakes and the fault structures emerge from geometrical and material built-in heterogeneities, one important open question is how the shearing becomes localized into a band of intense fractures. Here the authors address these questions through a numerical approach of a tectonic plate by considering rockmass heterogeneity both in microscopic scale and in mesoscopic scale. Numerical simulations of the progressive failure leading to collapse under long-range slow driving forces in the far-field show earthquake-like rupture behavior. $En Echelon$ crack-arrays are reproduced in the numerical simulation. It is demonstrated that the underlying fracturing induced acoustic emissions (or seismic events) display self-organized criticality------from disorder to order. The seismic cycles and the geometric structures of the fracture faces, which are found greatly depending on the material heterogeneity (especially on the macroscopic scale), agree with that observed experimentally in real brittle materials. It is concluded that in order to predict a main shock, one must have extremely detailed knowledge on very minor features of the earth's crust far from the place where the earthquake originated. If correct, the model proposed here seemingly provides an explanation as to why earthquakes to date are not predicted so successfully. The reason is not that the authors do not understand earthquake mechanisms very well but that they still know little about our earth's crust.

Instability from steady and axisymmetric to steady and asymmetric floating half zone convection in a fat liquid bridge of larger Prandtl number

The linear instability analysis of the present paper shows that the thermocapillary convection in a half floating zone of larger Prandtl number has a steady instability mode w(i) = 0 and m = 1 for a fat liquid bridge V = 1.2 with small geometrical aspect ratio A = 0.6. This conclusion is different from the usual idea of hydrothermal instability, and implies that the instability of the system may excite a steady and axial asymmetric state before the onset of oscillation in the ease of large Prandtl number.

Dynamic Behaviour of Nanoscale Electrostatic Actuators

The dynamic behaviour for nanoscale electrostatic actuators is studied. A two Parameter mass-spring model is shown to exhibit a bifurcation from the case excluding an equilibrium point to the case including two equilibrium points as the geometrical dimensions of the device are altered. Stability analysis shows that one is a stable Hopf bifurcation point and the other is an unstable saddle point. In addition, we plot the diagram phases, which have periodic orbits around the Hopf point and a homoclinic orbit passing though the unstable saddle point.

On the Formation of Periodic Segmentation Cracks in a Coating Plated on a Substrate with Periodic Subsurface Inclusions

The mechanism of the formation of periodic segmentation cracks of a coating plated on a substrate with periodic subsurface inclusions (PSI) is investigated. The internal stress in coating and subsequently the strain energy release rate (SERR) of the segmentation cracks are computed with finite element method (FEM). And the effect of the geometrical parameters of the PSI is studied. The results indicate that the ratio of the width of the inclusion to the period of the repeated structure has an optimum value, at which the maximum internal tensile stress and SERR arise. On the other hand, the ratio of the max-thickness of the inclusion to the thickness of the coating has a threshold value, above which the further increase of this ratio should seldom influence the internal stress or the SERR.

Influence of Liquid Bridge Volume on the Floating Zone Convection

The liquid bridge volume is a critical geometrical parameter in addition to the aspect ratio for onset of oscillation in the floating zone convection. The oscillatory features are generally divided into two characteristic regions: slender liquid bridge region and fat liquid bridge region. The oscillatory modes in two regions are discussed in the present paper.

Molecular Dynamics Study Of The Mechanical Behavior Of Nickel Nanowire: Strain Rate Effects

We present the analysis of uniaxial deformation of nickel nanowires using molecular dynamics simulations, and address the strain rate effects on mechanical responses and deformation behavior. The applied strain rate is ranging from 1 x 10(8) s(-1) to 1.4 x 10(11) s(-1). The results show that two critical strain rates, i.e., 5 x 10(9) s(-1) and 8 x 10(10) s(-1), are observed to play a pivotal role in switching between plastic deformation modes. At strain rate below 5 x 10(9) s(-1), Ni nanowire maintains its crystalline structure with neck occurring at the end of loading, and the plastic deformation is characterized by {111} slippages associated with Shockley partial dislocations and rearrangements of atoms close to necking region. At strain rate above 8x10(10) s(-1), Ni nanowire transforms from a fcc crystal into a completely amorphous state once beyond the yield point, and hereafter it deforms uniformly without obvious necking until the end of simulation. For strain rate between 5 x 10(9) s(-1) and 8 x 10(10) s(-1), only part of the nanowire exhibits amorphous state after yielding while the other part remains crystalline state. Both the {111} slippages in ordered region and homogenous deformation in amorphous region contribute to the plastic deformation. (C) 2007 Published by Elsevier B.V.

Plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant materials

The plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant material are investigated in this paper. Applying the stress-strain relation for the pressure-sensitive dilatant material, we have obtained an exact asymptotic solution for the plane strain tip fields for two types of cracks, one of which lies in the pressure-sensitive dilatant material and the other in the elastic material and their tips touch both the bimaterial interface. In cases, numerical results show that the singularity and the angular variations of the fields obtained depend on the material hardening exponent n, the pressure sensitivity parameter mu and geometrical parameter lambda.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

Dynamic equations for curved submerged floating tunnel

In virtue of reference Cartesian coordinates, geometrical relations of spatial curved structure are presented in orthogonal curvilinear coordinates. Dynamic equations for helical girder are derived by Hamilton principle. These equations indicate that four generalized displacements are coupled with each other. When spatial structure degenerates into planar curvilinear structure, two generalized displacements in two perpendicular planes are coupled with each other. Dynamic equations for arbitrary curvilinear structure may be obtained by the method used in this paper.

Waverider configurations derived from general conical flowfields

A method based on the computational fluid dynamics (CFD) is presented for a flexible waverider's design. The generating bodies of this method could be any cones. In addition, either the leading edge or the profile of the scramjet's inlet is used as the waverider's definition curve, parameterized by the quadric function, the sigmoid function or the B-spline function. Furthermore, several numerical examples are carried out to validate the method and the relevant codes. The CFD results of the configurations show that all the designs are successful. Moreover, primary suggestions are proposed for practical design by comparing the geometrical and aerodynamic performances of the cone-derived waveriders at Mach 6.

Investigation on a simulation model of floating half zone convection .2. Experiment

A simulation model of floating half zone with non-uniform temperature distribution at the upper rod and uniform temperature distribution at lower rod was discussed by numerical investigation in a previous paper. In the present paper, the experimental investigation of the simulation model is given generally. The results of the present model show that the temperature profile is quite different and the critical applied temperature difference is lower than the one of usual model with same geometrical parameters in most cases. The features of critical Marangoni number depending on the liquid bridge volume are also different from the ones of usual model.

Elastodynamic stress-intensity factors for a semiinfinite crack under 3-d combined mode loading

The dynamic stress intensity factor histories for a half plane crack in an otherwise unbounded elastic body are analyzed. The crack is subjected to a traction distribution consisting of two pairs of suddenly-applied shear point loads, at a distance L away from the crack tip. The exact expression for the combined mode stress intensity factors as the function of time and position along the crack edge is obtained. The method of solution is based on the direct application of integral transforms together with the Wiener-Hopf technique and the Cagniard-de Hoop method, which were previously believed to be inappropriate. Some features of solutions are discussed and the results are displayed in several figures.

The inhomogeneity of plastic-deformation in ductile single-crystals

A discrete slip model which characterizes the inhomogeneity of material properties in ductile single crystals is proposed in this paper. Based on this model rate-dependent finite element investigations are carried out which consider the finite deformation, finite rotation, latent hardening effect and elastic anisotropy. The calculation clearly exhibits the process from microscopic inhomogeneous and localized deformation to necking and the formation of LSBS and reveals several important features of shear localization. For example, the inhomogeneous deformation is influenced by the imperfections and initial non-uniformities of material properties. The inhomogeneous deformation may either induce necking which results in the lattice rotation and leads to geometrical softening, which in turn promotes the formation of CSBS, or induces heavily localized deformation. The microscopic localized deformation eventually develops into the LSBS and results in a failure. These results are in close agreement with experiment. Our calculations also find that the slip lines on the specimen's surface at necking become curved and also find that if the necking occurs before the formation of LSBS, this band must be misoriented from the operative slip systems. In this case, the formation of LSBS must involve non-crystallographic effects. These can also be indirectly confirmed by experiment. All these suggest that our present discrete slip model offers a correct description of the inhomogeneous deformation characterization in ductile crystals.

Elastodynamic stress intensity factor history for a semi-infinite crack under three-dimensional transient loading

The dynamic stress intensity factor history for a semi-infinite crack in an otherwise unbounded elastic body is analyzed. The crack is subjected to a pair of suddenly-applied point loadings on its faces at a distance L away from the crack tip. The exact expression for the mode I stress intensity factor as a function of time is obtained. The method of solution is based on the direct application of integral transforms, the Wiener-Hopf technique and the Cagniard-de Hoop method. Due to the existence of the characteristic length in loading this problem was long believed a knotty problem. Some features of the solutions are discussed and graphical result for numerical computation is presented.

3-dimensional transient wave response in a cracked elastic solid

The three-dimensional transient wave response problem is presented for an infinite elastic medium weakened by a plane crack of infinite length and finite width. Tractions are applied suddenly to the crack, which simulates the case of impact loading. The integral transforms are utilized to reduce the problem to a standard Fredholm integral equation in the Laplace transform variable and sequentially invert the Laplace transforms of the stress components by numerical inversion method. The dynamic mode I stress intensity factors at the crack tip are obtained and some numerical results are presented in graphical form.