Finite element analysis of fault bend influence on stick-slip instability along an intra-plate fault

Earthquakes have been recognized as resulting from stick-slip frictional instabilities along the faults between deformable rocks. A three-dimensional finite-element code for modeling the nonlinear frictional contact behaviors between deformable bodies with the node-to-point contact element strategy has been developed and applied here to investigate the fault geometry influence on the nucleation and development process of the stick-slip instability along an intra-plate fault through a typical fault bend model, which has a pre-cut fault that is artificially bent by an angle of 5.6degrees at the fault center. The numerical results demonstrate that the geometry of the fault significantly affects nucleation, termination and restart of the stick-slip instability along the intra-plate fault, and all these instability phenomena can be well simulated using the current finite-element algorithm.

A unified friction description and its application to the simulation of frictional instability using the finite element method

This article first summarizes some available experimental results on the frictional behaviour of contact interfaces, and briefly recalls typical frictional experiments and relationships, which are applicable for rock mechanics, and then a unified description is obtained to describe the entire frictional behaviour. It is formulated based on the experimental results and applied with a stick and slip decomposition algorithm to describe the stick-slip instability phenomena, which can describe the effects observed in rock experiments without using the so-called state variable, thus avoiding related numerical difficulties. This has been implemented to our finite element code, which uses the node-to-point contact element strategy proposed by the authors to handle the frictional contact between multiple finite-deformation bodies with stick and finite frictional slip, and applied here to simulate the frictional behaviour of rocks to show its usefulness and efficiency.

Parameters affecting the performance of wetting and drying in a two-dimensional finite element long wave hydrodynamic model

The marsh porosity method, a type of thin slot wetting and drying algorithm in a two-dimensional finite element long wave hydrodynamic model, is discussed and analyzed to assess model performance. Tests, including comparisons to simple examples and theoretical calculations, examine the effects of varying the marsh porosity parameters. The findings demonstrate that the wetting and drying concept of marsh porosity, often used in finite element hydrodynamic modeling, can behave in a more complex manner than initially expected.

Prediction of the growth of wear-type rail corrugation

The growth behaviour of the vibrational wear phenomenon known as rail corrugation is investigated analytically and numerically using mathematical models. A simplified feedback model for wear-type rail corrugation that includes a wheel pass time delay is developed with an aim to analytically distil the most critical interaction occurring between the wheel/rail structural dynamics, rolling contact mechanics and rail wear. To this end, a stability analysis on the complete system is performed to determine the growth of wear-type rail corrugations over multiple wheelset passages. This analysis indicates that although the dynamical behaviour of the system is stable for each wheel passage, over multiple wheelset passages, the growth of wear-type corrugations is shown to be the result of instability due to feedback interaction between the three primary components of the model. The corrugations are shown analytically to grow for all realistic railway parameters. From this analysis an analytical expression for the exponential growth rate of corrugations in terms of known parameters is developed. This convenient expression is used to perform a sensitivity analysis to identify critical parameters that most affect corrugation growth. The analytical predictions are shown to compare well with results from a benchmarked time-domain finite element model. (C) 2004 Elsevier B.V. All rights reserved.

Distortional buckling of I-beams by finite element method

Distortional buckling, unlike the usual lateral-torsional buckling in which the cross-section remains rigid in its own plane, involves distortion of web in the cross-section. This type of buckling typically occurs in beams with slender web and stocky flanges. Most of the published studies assume the web to deform with a cubic shape function. As this assumption may limit the accuracy of the results, a fifth order polynomial is chosen here for the web displacements. The general line-type finite element model used here has two nodes and a maximum of twelve degrees of freedom per node. The model not only can predict the correct coupled mode but also is capable of handling the local buckling of the web.