A mathematical model for the waste region of a long-wall mineworking

The mechanical behaviour of the broken rock material occupying the waste region of a long-wall mineworking is investigated. A mathematical model is proposed, developed from physical considerations, which is based upon the continuum, plane strain, rigid-plastic theory for double-shearing flows of compressible granular materials. The model comprises the stress equilibrium equations, the Coulomb yield criterion, kinematic equations due to Mehrabadi and Cowin (1978) and the continuity equation, together with the waste region geometry and the boundary conditions imposed upon the field variables. Statical indeterminancy is a property of the formulation and this results in a connection between the stress and velocity fields through the boundary conditions. A strategy is presented for the iterative construction of the stress and velocity fields and approximations to solutions of the equations of the model are obtained numerically.     

A matrix formulation for the dynamic analysis of planar mechanisms using point coordinates and velocity transformation

This paper presents a matrix formulation for the dynamic analysis of planar mechanisms consisting of interconnected rigid bodies. The formulation initially uses the rectangular Cartesian coordinates of an equivalent constrained system of particles to define the configuration of the mechanical system. This results in a simple and straightforward procedure for generating the equations of motion. The equations of motion are then derived in terms of relative joint coordinates through the use of a velocity transformation matrix. The velocity transformation matrix relates the relative joint velocities to the Cartesian velocities. For the open-loop case, this process automatically eliminates all of the non-working constraint forces and leads to an efficient integration of the equations of motion. For the closed-loop case, suitable joints should be cut and few cut-joints constraint equations should be included for each closed loop. Two examples are used to demonstrate the generality and efficiency of the proposed method.     

Plane deformations of ideal granular materials3

The flow rule based on a double-shearing model developed by A. J. M. Spencer, and M. M. Mehrabadi and S. C. Cowin is intended for use in analyzing fully developed plane flows of rate-independent, isotropic, rigidperfectly plastic frictional and dilatant granular materials. With the view of extending the range of applicability of this physically based flow rule to the prepeak (prelocalization) regime, we systematically introduce isotropic hardening (softening) via the usual consideration of scalar hardening variables in the yield function. However, for this new model to be in accord with some recent experimental observations of I. Vardoulakis on shear-band localizations, it is found that we have to abandon the classical assumption of the coincidence of stress and velocity characteristics.     

Dynamic modelling of the double wishbone motor-vehicle suspension system

In this paper the dynamic analysis of the double wishbone motor-vehicle suspension system using the point-joint coordinates formulation is presented. The mechanical system is replaced by an equivalent constrained system of particles and then the laws of particle dynamics are used to derive the equations of motion. Due to the presence of large number of geometric and kinematic constraints the velocity transformation approach is used to eliminate some constraints. The equations of motion in terms of the Cartesian coordinates of the particles are transformed to a reduced set in terms of relative joint variables by defining differential-algebraic equations in terms of the joint variables are equal to the number of degrees of freedom of the whole system plus the number of cut-joint constraints corresponding to cut of kinematical closed loops. Use of both the Cartesian and relative joint variables produces an efficient set of equations without loss of generality. The chosen suspension includes open and closed loops with quarter-car model.     

Three-dimensional constitutive relations for granular materials based on the dilatant double shearing mechanism and the concept of fabric

An extension of a three-dimensional model proposed by Anand and Gu (2000) for amorphous granular materials to include the effects of initial and induced anisotropy is presented in this paper. The proposed model can also be considered as a three-dimensional generalization of a model recently developed by Zhu et al. (2005) for the planar deformation of granular materials. The main ingredients of the model include the dilatant double shearing mechanism (Spencer, 1964, Mehrabadi and Cowin, 1978), the concept of fabric (Oda, 1972), and an extension of the Mohr–Coulomb yield criterion (Shield, 1955, Spencer, 1982) to three dimensions.The constitutive equations are implemented in the finite element program ABAQUS/Explicit (ABAQUS, 2001) by developing a user-material subroutine to conduct numerical triaxial compression tests for samples of granular materials with different initial anisotropy. The numerical results agree with the observed behavior and show that the extended constitutive model is capable of capturing the strength anisotropy of granular materials. Employing the anisotropic model developed here, we have also repeated the numerical simulation of the stress state in a static conical sand pile conducted earlier by Anand and Gu (2000). We find that fabric has little or no influence on the vertical stress distribution except at the base of the sand pile where the peak value of this stress is slightly higher than that predicted by the model of Anand and Gu (2000) which does not include the effects of fabric. We also find that the direction of the principal compressive stress changes from vertical at points away from the center of the pile to almost horizontal at points close to the center of the pile. This result provides a possible explanation for the observed dip in the vertical stress distribution in sand piles.     

Simulating turbulent Dean flow in Cartesian coordinates

A simplified approach to simulate turbulent flows in curved channels is proposed. A set of governing equations of motion in Cartesian coordinates is derived from the full Navier–Stokes equations in cylindrical coordinates. Terms to first order in the dimensionless curvature parameter are retained, whereas higher‐order terms are neglected. The curvature terms are implemented in a conventional Navier–Stokes code using Cartesian coordinates. Direct numerical simulations (DNS) of turbulent flow in weakly curved channels are performed. The pronounced asymmetries in the mean flow and the turbulence statistics observed in earlier DNS studies are faithfully reproduced by the present simplified Navier–Stokes model. It is particularly rewarding that also distinct pairs of counter‐rotating streamwise‐oriented vortices are embedded in the simulated flow field. Copyright © 2008 John Wiley & Sons, Ltd.     

A new method for obtaining the exact solutions of the Navier-Stokes (NS) equations and its applications

In this paper we propose a new method for obtaining the exact solutions of the Mavier-Stokes (NS) equations for incompressible viscous fluid in the light of the theory of simplified Navier-Stokes (SNS) equations developed by the first author^[1,2], Using the present method we can find some new exact solutions as well as the well-known exact solutions of the NS equations. In illustration of its applications, we give a variety of exact solutions of incompressible viscous fluid flows for which NS equations of fluid motion are written in Cartesian coordinates, or in cylindrical polar coordinates, or in spherical coordinates. The project supported by National Natural Science Foundation of China.     

Dynamic equations of a prestressed magnetoelectroelastic medium

The constitutive relations of nonlinear mechanics of a magnetoelectroelastic medium subjected to initial mechanical stresses are linearized in the framework of material (Lagrangian) coordinates. The final expressions are constructed independently of the choice of curvilinear coordinates and are represented in a form convenient for theoretical and applied studies. The constitutive relations for the motion of a prestressed magnetoelectroelastic medium are given in rectangular Cartesian coordinates. The influence of the initial mechanical stresses on piezomagnetoelectric materials of the class 6mm is studied.     

Dynamic Analysis of Spatial Linkages: A Recursive Approach

In this paper, recursive equations of motion of spatial linkages are presented. The method uses the concepts of linear and angular momentums to generate the rigid body equations of motion in terms of the Cartesian coordinates of a dynamically equivalent constrained system of particles, without introducing any rotational coordinates and the corresponding rotational transformation matrix. For the open-chain system, the equations of motion are generated recursively along the serial chains. Closed-chain system is transformed to open-chain by cutting suitable kinematic joints and introducing cut–joint constraints. An example is chosen to demonstrate the generality and simplicity of the developed formulation.     

On the use of linear graph theory in multibody system dynamics

Multibody dynamics involves the generation and solution of the equations of motion for a system of connected material bodies. The subject of this paper is the use of graph-theoretical methods to represent multibody system topologies and to formulate the desired set of motion equations; a discussion of the methods available for solving these differential-algebraic equations is beyond the scope of this work. After a brief introduction to the topic, a review of linear graphs and their associated topological arrays is presented, followed in turn by the use of these matrices in generating various graph-theoretic equations. The appearance of linear graph theory in a number of existing multibody formulations is then discussed, distinguishing between approaches that use absolute (Cartesian) coordinates and those that employ relative (joint) coordinates. These formulations are then contrasted with formal graph-theoretic approaches, in which both the kinematic and dynamic equations are automatically generated from a single linear graph representation of the system. The paper concludes with a summary of results and suggestions for further research on the graph-theoretical modelling of mechanical systems.     

Numerical study of the noninertial systems: application to train coupler systems

Car coupler forces have a significant effect on the longitudinal train dynamics and stability. Because the coupler inertia is relatively small in comparison with the car inertia; the high stiffness associated with the coupler components can lead to high frequencies that adversely impact the computational efficiency of train models. The objective of this investigation is to study the effect of the coupler inertia on the train dynamics and on the computational efficiency as measured by the simulation time. To this end, two different models are developed for the car couplers; one model, called the inertial coupler model, includes the effect of the coupler inertia, while in the other model, called the noninertial model, the effect of the coupler inertia is neglected. Both inertial and noninertial coupler models used in this investigation are assumed to have the same coupler kinematic degrees of freedom that capture geometric nonlinearities and allow for the relative translation of the draft gears and end of car cushioning (EOC) devices as well as the relative rotation of the coupler shank. In both models, the coupler kinematic equations are expressed in terms of the car body and coupler coordinates. Both the inertial and noninertial models used in this study lead to a system of differential and algebraic equations that are solved simultaneously in order to determine the coordinates of the cars and couplers. In the case of the inertial model, the coupler kinematics is described using the absolute Cartesian coordinates, and the algebraic equations describe the kinematic constraints imposed on the motion of the system. In this case of the inertial model, the constraint equations are satisfied at the position, velocity, and acceleration levels. In the case of the noninertial model, the equations of motion are developed using the relative joint coordinates, thereby eliminating systematically the algebraic equations that represent the kinematic constraints. A quasistatic force analysis is used to determine a set of coupler nonlinear force algebraic equations for a given car configuration. These nonlinear force algebraic equations are solved iteratively to determine the coupler noninertial coordinates which enter into the formulation of the equations of motion of the train cars. The results obtained in this study showed that the neglect of the coupler inertia eliminates high frequency oscillations that can negatively impact the computational efficiency. The effect of these high frequencies that are attributed to the coupler inertia on the simulation time is examined using frequency and eigenvalue analyses. While the neglect of the coupler inertia leads, as demonstrated in this investigation, to a much more efficient model, the results obtained using the inertial and noninertial coupler models show good agreement, demonstrating that the coupler inertia can be neglected without having an adverse effect on the accuracy of the solution.     

New mixed plate-bending elements with shear strain interpolations

This paper is concerned with two mixed plate-bending elements with shear strain interpolations, a quadrilateral element and an 8-node serendipity-type element based on discussions on the element proposed in Ref.[1]. The shear strains and inner-forces in the natural coordinates are interpolated in an element and then transformed into Cartesian coordinates in accordance with covariant and contravariant tensor transformations, respectively. The quadrilateral element coincides with the element in Ref.[1] when it is rectangular. Numerical examples show that the two new elements are free from shear locking and spurious kinematic modes under regular and irregular meshes and have the advantages of being insensitive to element distortion and able to give fairly accurate results.The Project supported by National Natural Science Foundation of China.     

Implementation of CLEAR algorithm on non‐orthogonal curvilinear co‐ordinates for solution of incompressible flow and heat transfer

In this paper, the CLEAR (coupled and linked equations algorithm revised) algorithm is extended to non‐orthogonal curvilinear collocated grids. The CLEAR algorithm does not introduce pressure correction in order to obtain an incompressible flow field which satisfies the mass conservation law. Rather, it improves the intermediate velocity by solving an improved pressure equation to make the algorithm fully implicit since there is no term omitted in the derivation process. In the extension of CLEAR algorithm from a staggered grid system in Cartesian coordinates to collocated grids in non‐orthogonal curvilinear coordinates, three important issues are appropriately treated so that the extended CLEAR can lead to a unique solution without oscillation of pressure field and with high robustness. These three issues are (1) solution independency on the under‐relaxation factor; (2) strong coupling between velocity and pressure; and (3) treatment of the cross pressure gradient terms. The flow and heat transfer problems in a rectangular enclosure with an internal eccentric circle and the flow in a lid‐driven inclined cavity are computed by using the extended CLEAR. The results show that the extended CLEAR can guarantee the solution independency on the under‐relaxation factor, the smoothness of pressure profile even at very small under‐relaxation factor and good robustness which leads to a converged solution for the small inclined angle of 5^° only with 5‐point computational molecule while the extended SIMPLE‐series algorithm usually can get a converged solution for the inclined angle larger than 30^° under the same condition. Copyright © 2006 John Wiley & Sons, Ltd.     

Laminar flow of micropolar fluid in rectangular microchannels

Compared with the classic flow on macroscale, flows in microchannels have some new phenomena such as the friction increase and the flow rate reduction. Papautsky and co-workers explained these phenomena by using a micropolar fluid model where the effects of micro-rotation of fluid molecules were taken into account. But both the curl of velocity vector and the curl of micro-rotation gyration vector were given incorrectly in the Cartesian coordinates and then the micro-rotation gyration vector had only one component in the z-direction. Besides, the gradient term of the divergence of micro-rotation gyration vector was missed improperly in the angular moment equation. In this paper, the governing equations for laminar flows of micropolar fluid in rectangular microchannels are reconstructed. The numerical results of velocity profiles and micro-rotation gyrations are obtained by a procedure based on the Chebyshev collocation method. The micropolar effects on velocity and micro-rotation gyration are discussed in detail.The project was supported by the National Natural Science Foundation of China (10472054). The English text was polished by Boyi Wang.     

Implicit application of non‐reflective boundary conditions for Navier–Stokes equations in generalized coordinates

The non‐reflective boundary conditions (NRBC) for Navier–Stokes equations originally suggested by Poinsot and Lele (J. Comput. Phys. 1992; 101 :104–129) in Cartesian coordinates are extended to generalized coordinates. The characteristic form Navier–Stokes equations in conservative variables are given. In this characteristic‐based method, the NRBC is implicitly coupled with the Navier–Stokes flow solver and are solved simultaneously with the flow solver. The calculations are conducted for a subsonic vortex propagating flow and the steady and unsteady transonic inlet‐diffuser flows. The results indicate that the present method is accurate and robust, and the NRBC are essential for unsteady flow calculations. Copyright © 2005 John Wiley & Sons, Ltd.     

Numerical Study of Three-Dimensional Stability Problems for Laminated and Ribbon-Reinforced Composites

Classes of basic problems (classical problems and new problems derived by partial separation of variables) of the theory of elasticity and three-dimensional linearized stability are formulated in Cartesian coordinates. For each class of problems, the base scheme and the base system of equations are constructed and substantiated. They allow us to write explicitly the global mesh equations corresponding to an arbitrary problem from the class under consideration     

Analysis of Laminar Flow and Forced Convection Heat Transfer in a Porous Medium

The flow of an incompressible Newtonian fluid confined in a planar geometry with different wall temperatures filled with a homogenous and isotropic porous medium is analyzed in terms of determining the unsteady state and steady state velocities, the temperature and the entropy generation rate as function of the pressure drop, the Darcy number, and the Brinkman number. The one-dimensional approximate equation in the rectangular Cartesian coordinates governing the flow of a Newtonian fluid through porous medium is derived by accounting for the order of magnitude of terms as well as accompanying approximations to the full-blown three-dimensional equations by using scaling arguments. The one-dimensional approximate energy and the entropy equations with the viscous dissipation consisting of the velocity gradient and the square of velocity are derived by following the same procedure used in the derivation of velocity expressions. The one-dimensional approximate equations for the velocity, the temperature, and the entropy generation rate are analytically solved to determine the velocity, the temperature, and the entropy distributions in the saturated porous medium as functions of the effective process parameters. It is found that the pressure drop, the Darcy number, and the Brinkman number affect the temperature distribution in the similar way, and besides the above parameters, the irreversibility distribution ratio also affects the entropy generation rate in the similar way.     

A velocity transformation method for the nonlinear dynamic simulation of railroad vehicle systems

The solution of the constrained multibody system equations of motion using the generalized coordinate partitioning method requires the identification of the dependent and independent coordinates. Using this approach, only the independent accelerations are integrated forward in time in order to determine the independent coordinates and velocities. Dependent coordinates are determined by solving the nonlinear constraint equations at the position level. If the constraint equations are highly nonlinear, numerical difficulties can be encountered or more Newton–Raphson iterations may be required in order to achieve convergence for the dependent variables. In this paper, a velocity transformation method is proposed for railroad vehicle systems in order to deal with the nonlinearity of the constraint equations when the vehicles negotiate curved tracks. In this formulation, two different sets of coordinates are simultaneously used. The first set is the absolute Cartesian coordinates which are widely used in general multibody system computer formulations. These coordinates lead to a simple form of the equations of motion which has a sparse matrix structure. The second set is the trajectory coordinates which are widely used in specialized railroad vehicle system formulations. The trajectory coordinates can be used to obtain simple formulations of the specified motion trajectory constraint equations in the case of railroad vehicle systems. While the equations of motion are formulated in terms of the absolute Cartesian coordinates, the trajectory accelerations are the ones which are integrated forward in time. The problems associated with the higher degree of differentiability required when the trajectory coordinates are used are discussed. Numerical examples are presented in order to examine the performance of the hybrid coordinate formulation proposed in this paper in the analysis of multibody railroad vehicle systems.     

The equations of complete depth-averaged turbulence model in general orthogonal coordinates

ＤｉｎｇＹａｎ（丁剡）;ＺｈｏｕＸｕｅｙｉ（周雪漪）;ＹｕＣｈａｎｇｚｈａｏ（余常昭）;ＬｉａｎｇＤｏｎｇ（梁栋）（ＲｅｃｅｉｖｅｄＪｕｎｅ１４,１９９４;ＣｏｍｍｕｎｉｃａｔｅｄｂｙＢｉａｎＹｉｎｇｕｉ）ＴＨＥＥＱＵＡＴＩＯＮＳＯＦＣＯＭＰＬＥＴＥ...     

A Cartesian grid technique based on one‐dimensional integrated radial basis function networks for natural convection in concentric annuli

This paper reports a radial basis function (RBF)‐based Cartesian grid technique for the simulation of two‐dimensional buoyancy‐driven flow in concentric annuli. The continuity and momentum equations are represented in the equivalent stream function formulation that reduces the number of equations from three to one, but involves higher‐order derivatives. The present technique uses a Cartesian grid to discretize the problem domain. Along a grid line, one‐dimensional integrated RBF networks (1D‐IRBFNs) are employed to represent the field variables. The capability of 1D‐IRBFNs to handle unstructured points with accuracy is exploited to describe non‐rectangular boundaries in a Cartesian grid, while the method's ability to avoid the reduction of convergence rate caused by differentiation is instrumental in improving the quality of the approximation of higher‐order derivatives. The method is applied to simulate thermally driven flows in annuli between two circular cylinders and between an outer square cylinder and an inner circular cylinder. High Rayleigh number solutions are achieved and they are in good agreement with previously published numerical data. Copyright © 2007 John Wiley & Sons, Ltd.