Centrifuge and discrete element modelling of tunnelling effects on pipelines

The paper presents centrifuge test data of the problem of tunnelling effects on buried pipelines and compares them to predictions made using DEM simulations. The paper focuses on the examination of pipeline bending moments, their distribution along the pipe, and their development with tunnel volume loss. Centrifuge results are obtained by PIV analysis and compared to results obtained using the DEM model. The DEM model was built to replicate the centrifuge model as closely as possible and included numerical features formulated specially for this task, such as structural elements to replicate the tunnel and pipeline. Results are extremely encouraging, with deviations between DEM and centrifuge test bending moment results being very small. © 2010 Taylor & Francis Group, London.

Representations and algorithms for finite-state bisimulations of linear discrete-time control systems

While a large amount of research over the past two decades has focused on discrete abstractions of infinite-state dynamical systems, many structural and algorithmic details of these abstractions remain unknown. To clarify the computational resources needed to perform discrete abstractions, this paper examines the algorithmic properties of an existing method for deriving finite-state systems that are bisimilar to linear discrete-time control systems. We explicitly find the structure of the finite-state system, show that it can be enormous compared to the original linear system, and give conditions to guarantee that the finite-state system is reasonably sized and efficiently computable. Though constructing the finite-state system is generally impractical, we see that special cases could be amenable to satisfiability based verification techniques. ©2009 IEEE.

Modeling fatigue crack growth in crystalline solids with discrete dislocation plasticity

Analyses of crack growth under cyclic loading conditions are discussed where plastic flow arises from the motion of large numbers of discrete dislocations and the fracture properties are embedded in a cohesive surface constitutive relation. The formulation is the same as used to analyse crack growth under monotonic loading conditions, differing only in the remote loading being a cyclic function of time. Fatigue, i.e. crack growth in cyclic loading at a driving force for which the crack would have arrested under monotonic loading, emerges in the simulations as a consequence of the evolution of internal stresses associated with the irreversibility of the dislocation motion. A fatigue threshold, Paris law behaviour, striations, the accelerated growth of short cracks and the scaling with material properties are outcomes of the calculations. Results for single crystals and polycrystals will be discussed.

The Effect of Motivation on Movement: A Study of Bradykinesia in Parkinson's Disease

Background: Bradykinesia is a cardinal feature of Parkinson's disease (PD). Despite its disabling impact, the precise cause of this symptom remains elusive. Recent thinking suggests that bradykinesia may be more than simply a manifestation of motor slowness, and may in part reflect a specific deficit in the operation of motivational vigour in the striatum. In this paper we test the hypothesis that movement time in PD can be modulated by the specific nature of the motivational salience of possible action-outcomes. Methodology/Principal Findings: We developed a novel movement time paradigm involving winnable rewards and avoidable painful electrical stimuli. The faster the subjects performed an action the more likely they were to win money (in appetitive blocks) or to avoid a painful shock (in aversive blocks). We compared PD patients when OFF dopaminergic medication with controls. Our key finding is that PD patients OFF dopaminergic medication move faster to avoid aversive outcomes (painful electric shocks) than to reap rewarding outcomes (winning money) and, unlike controls, do not speed up in the current trial having failed to win money in the previous one. We also demonstrate that sensitivity to distracting stimuli is valence specific. Conclusions/Significance: We suggest this pattern of results can be explained in terms of low dopamine levels in the Parkinsonian state leading to an insensitivity to appetitive outcomes, and thus an inability to modulate movement speed in the face of rewards. By comparison, sensitivity to aversive stimuli is relatively spared. Our findings point to a rarely described property of bradykinesia in PD, namely its selective regulation by everyday outcomes. © 2012 Shiner et al.

Control of bimanual rhythmic movements: Trading efficiency for robustness depending on the context

This paper investigates how the efficiency and robustness of a skilled rhythmic task compete against each other in the control of a bimanual movement. Human subjects juggled a puck in 2D through impacts with two metallic arms, requiring rhythmic bimanual actuation. The arms kinematics were only constrained by the position, velocity and time of impacts while the rest of the trajectory did not influence the movement of the puck. In order to expose the task robustness, we manipulated the task context in two distinct manners: the task tempo was assigned at four different values (hence manipulating the time available to plan and execute each impact movement individually); and vision was withdrawn during half of the trials (hence reducing the sensory inflows). We show that when the tempo was fast, the actuation was rhythmic (no pause in the trajectory) while at slow tempo, the actuation was discrete (with pause intervals between individual movements). Moreover, the withdrawal of visual information encouraged the rhythmic behavior at the four tested tempi. The discrete versus rhythmic behavior give different answers to the efficiency/robustness trade-off: discrete movements result in energy efficient movements, while rhythmic movements impact the puck with negative acceleration, a property preserving robustness. Moreover, we report that in all conditions the impact velocity of the arms was negatively correlated with the energy of the puck. This correlation tended to stabilize the task and was influenced by vision, revealing again different control strategies. In conclusion, this task involves different modes of control that balance efficiency and robustness, depending on the context. © 2008 Springer-Verlag.

Continuous dynamical systems that realize discrete optimization on the hypercube

We study the problem of finding a local minimum of a multilinear function E over the discrete set {0,1}n. The search is achieved by a gradient-like system in [0,1]n with cost function E. Under mild restrictions on the metric, the stable attractors of the gradient-like system are shown to produce solutions of the problem, even when they are not in the vicinity of the discrete set {0,1}n. Moreover, the gradient-like system connects with interior point methods for linear programming and with the analog neural network studied by Vidyasagar (IEEE Trans. Automat. Control 40 (8) (1995) 1359), in the same context. © 2004 Elsevier B.V. All rights reserved.

Plane-strain discrete dislocation plasticity with climb-assisted glide motion of dislocations

A small-strain two-dimensional discrete dislocation plasticity (DDP) framework is developed wherein dislocation motion is caused by climb-assisted glide. The climb motion of the dislocations is assumed to be governed by a drag-type relation similar to the glide-only motion of dislocations: such a relation is valid when vacancy kinetics is either diffusion limited or sink limited. The DDP framework is employed to predict the effect of dislocation climb on the uniaxial tensile and pure bending response of single crystals. Under uniaxial tensile loading conditions, the ability of dislocations to bypass obstacles by climb results in a reduced dislocation density over a wide range of specimen sizes in the climb-assisted glide case compared to when dislocation motion is only by glide. A consequence is that, at least in a single slip situation, size effects due to dislocation starvation are reduced. By contrast, under pure bending loading conditions, the dislocation density is unaffected by dislocation climb as geometrically necessary dislocations (GNDs) dominate. However, climb enables the dislocations to arrange themselves into lower energy configurations which significantly reduces the predicted bending size effect as well as the amount of reverse plasticity observed during unloading. The results indicate that the intrinsic plasticity material length scale associated with GNDs is strongly affected by thermally activated processes and will be a function of temperature. © 2013 IOP Publishing Ltd.

Crystallographically driven Au catalyst movement during growth of InAs/GaAs axial nanowire heterostructures

The movement of Au catalysts during growth of InAs on GaAs nanowires has been carefully investigated by transmission electron microscopy. It has been found that Au catalysts preferentially stay on { 112 } B GaAs sidewalls. Since a {112} surface is composed of a {111} facet and a {002} facet and since {111} facets are polar facets for the zinc-blende structure, this crystallographic preference is attributed to the different interface energies caused by the different polar facets. We anticipate that these observations will be useful for the design of nanowire heterostructure based devices. © 2009 American Institute of Physics.

The impact of sand slugs against beams and plates: Coupled discrete particle/finite element simulations

The impact of a slug of dry sand particles against a metallic sandwich beam or circular sandwich plate is analysed in order to aid the design of sandwich panels for shock mitigation. The sand particles interact via a combined linear-spring-and-dashpot law whereas the face sheets and compressible core of the sandwich beam and plate are treated as rate-sensitive, elastic-plastic solids. The majority of the calculations are performed in two dimensions and entail the transverse impact of end-clamped monolithic and sandwich beams, with plane strain conditions imposed. The sand slug is of rectangular shape and comprises a random loose packing of identical, circular cylindrical particles. These calculations reveal that loading due to the sand is primarily inertial in nature with negligible fluid-structure interaction: the momentum transmitted to the beam is approximately equal to that of the incoming sand slug. For a slug of given incoming momentum, the dynamic deflection of the beam increases with decreasing duration of sand-loading until the impulsive limit is attained. Sandwich beams with thick, strong cores significantly outperform monolithic beams of equal areal mass. This performance enhancement is traced to the "sandwich effect" whereby the sandwich beams have a higher bending strength than that of the monolithic beams. Three-dimensional (3D) calculations are also performed such that the sand slug has the shape of a circular cylindrical column of finite height, and contains spherical sand particles. The 3D slug impacts a circular monolithic plate or sandwich plate and we show that sandwich plates with thick strong cores again outperform monolithic plates of equal areal mass. Finally, we demonstrate that impact by sand particles is equivalent to impact by a crushable foam projectile. The calculations on the equivalent projectile are significantly less intensive computationally, yet give predictions to within 5% of the full discrete particle calculations for the monolithic and sandwich beams and plates. These foam projectile calculations suggest that metallic foam projectiles can be used to simulate the loading by sand particles within a laboratory setting. © 2013 Elsevier Ltd.

Study on small-strain behaviours of methane hydrate sandy sediments using discrete element method

Methane hydrate bearing soil has attracted increasing interest as a potential energy resource where methane gas can be extracted from dissociating hydrate-bearing sediments. Seismic testing techniques have been applied extensively and in various ways, to detect the presence of hydrates, due to the fact that hydrates increase the stiffness of hydrate-bearing sediments. With the recognition of the limitations of laboratory and field tests, wave propagation modelling using Discrete Element Method (DEM) was conducted in this study in order to provide some particle-scale insights on the hydrate-bearing sandy sediment models with pore-filling and cementation hydrate distributions. The relationship between shear wave velocity and hydrate saturation was established by both DEM simulations and analytical solutions. Obvious differences were observed in the dependence of wave velocity on hydrate saturation for these two cases. From the shear wave velocity measurement and particle-scale analysis, it was found that the small-strain mechanical properties of hydrate-bearing sandy sediments are governed by both the hydrate distribution patterns and hydrate saturation. © 2013 AIP Publishing LLC.

Multi-scale modelling of granular avalanches

Avalanches, debris flows, and landslides are geophysical hazards, which involve rapid mass movement of granular solids, water and air as a single-phase system. The dynamics of a granular flow involve at least three distinct scales: the micro-scale, meso-scale, and the macro-scale. This study aims to understand the ability of continuum models to capture the micro-mechanics of dry granular collapse. Material Point Method (MPM), a hybrid Lagrangian and Eulerian approach, with Mohr-Coulomb failure criterion is used to describe the continuum behaviour of granular column collapse, while the micromechanics is captured using Discrete Element Method (DEM) with tangential contact force model. The run-out profile predicted by the continuum simulations matches with DEM simulations for columns with small aspect ratios ('h/r' < 2), however MPM predicts larger run-out distances for columns with higher aspect ratios ('h/r' > 2). Energy evolution studies in DEM simulations reveal higher collisional dissipation in the initial free-fall regime for tall columns. The lack of a collisional energy dissipation mechanism in MPM simulations results in larger run-out distances. Micro-structural effects, such as shear band formations, were observed both in DEM and MPM simulations. A sliding flow regime is observed above the distinct passive zone at the core of the column. Velocity profiles obtained from both the scales are compared to understand the reason for a slow flow run-out mobilization in MPM simulations. © 2013 AIP Publishing LLC.