Direct numerical simulation of the dynamics of sliding rough surfaces

The noise generated by the friction of two rough surfaces under weak contact pressure is usually called roughness noise. The underlying vibration which produces the noise stems from numerous instantaneous shocks (in the microsecond range) between surface micro-asperities. The numerical simulation of this problem using classical mechanics requires a fine discretization in both space and time. This is why the finite element method takes much CPU time. In this study, we propose an alternative numerical approach which is based on a truncated modal decomposition of the vibration, a central difference integration scheme and two algorithms for contact: The penalty algorithm and the Lagrange multiplier algorithm. Not only does it reproduce the empirical laws of vibration level versus roughness and sliding speed found experimentally but it also provides the statistical properties of local events which are not accessible by experiment. The CPU time reduction is typically a factor of 10.Comment: 16 pages, 16 figures, accepted versio

Estimation of the normal contact stiffness for frictional interface in sticking and sliding conditions

Modeling of frictional contact systems with high accuracy needs the knowledge of several contact parameters, which are mainly related to the local phenomena at the contact interfaces and affect the complex dynamics of mechanical systems in a prominent way. This work presents a newer approach for identifying reliable values of the normal contact stiffness between surfaces in contact, in both sliding and sticking conditions. The combination of experimental tests, on a dedicated set-up, with finite element modeling, allowed for an indirect determination of the normal contact stiffness. The stiffness was found to increase with increasing contact pressure and decreasing roughness, while the evolution of surface topography and third-body rheology affected the contact stiffness when sliding

Packing Characteristics of Different Shaped Proppants for use with Hydrofracing - A Numerical Investigation using 3D FEMDEM

Imperial Users onl

Friction Variability in Planar Pushing Data: Anisotropic Friction and Data-collection Bias

Friction plays a key role in manipulating objects. Most of what we do with our hands, and most of what robots do with their grippers, is based on the ability to control frictional forces. This paper aims to better understand the variability and predictability of planar friction. In particular, we focus on the analysis of a recent dataset on planar pushing by Yu et al. [1] devised to create a data-driven footprint of planar friction. We show in this paper how we can explain a significant fraction of the observed unconventional phenomena, e.g., stochasticity and multi-modality, by combining the effects of material non-homogeneity, anisotropy of friction and biases due to data collection dynamics, hinting that the variability is explainable but inevitable in practice. We introduce an anisotropic friction model and conduct simulation experiments comparing with more standard isotropic friction models. The anisotropic friction between object and supporting surface results in convergence of initial condition during the automated data collection. Numerical results confirm that the anisotropic friction model explains the bias in the dataset and the apparent stochasticity in the outcome of a push. The fact that the data collection process itself can originate biases in the collected datasets, resulting in deterioration of trained models, calls attention to the data collection dynamics.Comment: 8 pages, 13 figure

Fully implicit frictional dynamics with soft constraints

Dynamics simulation with frictional contacts is important for a wide range of applications, from cloth simulation to object manipulation. Recent methods using smoothed friction forces have enabled robust and differentiable simulation of elastodynamics with friction. However, the resulting frictional behaviors can be qualitatively inaccurate and may not converge to analytic solutions. Here we propose an alternative, fully implicit, formulation for simulating elastodynamics subject to frictional contacts with realistic friction behavior. Furthermore, we demonstrate how higher-order time integration can be used in our method, as well as in incremental potential methods. We develop an inexact Newton method with forward-mode automatic differentiation that simplifies the implementation and improves performance. Finally, we show how our method can be extended to respond to volume changes using a unified penalty function derived from first principles and capable of emulating compressible as well as nearly incompressible media

A novel approach to modelling and simulating the contact behaviour between a human hand model and a deformable object

A deeper understanding of biomechanical behaviour of human hands becomes fundamental for any human hand-operated Q2 activities. The integration of biomechanical knowledge of human hands into product design process starts to play an increasingly important role in developing an ergonomic product-to-user interface for products and systems requiring high level of comfortable and responsive interactions. Generation of such precise and dynamic models can provide scientific evaluation tools to support product and system development through simulation. This type of support is urgently required in many applications such as hand skill training for surgical operations, ergonomic study of a product or system developed and so forth. The aim of this work is to study the contact behaviour between the operators’ hand and a hand-held tool or other similar contacts, by developing a novel and precise nonlinear 3D finite element model of the hand and by investigating the contact behaviour through simulation. The contact behaviour is externalised by solving the problem using the bi-potential method. The human body’s biomechanical characteristics, such as hand deformity and structural behaviour, have been fully modelled by implementing anisotropic hyperelastic laws. A case study is given to illustrate the effectiveness of the approac

A finite element framework for modeling internal frictional contact in three-dimensional fractured media using unstructured tetrahedral meshes

AbstractThis paper introduces a three-dimensional finite element (FE) formulation to accurately model the linear elastic deformation of fractured media under compressive loading. The presented method applies the classic Augmented Lagrangian(AL)-Uzawa method, to evaluate the growth of multiple interacting and intersecting discrete fractures. The volume and surfaces are discretized by unstructured quadratic triangle-tetrahedral meshes; quarter-point triangles and tetrahedra are placed around fracture tips. Frictional contact between crack faces for high contact precisions is modeled using isoparametric integration point-to-integration point contact discretization, and a gap-based augmentation procedure. Contact forces are updated by interpolating tractions over elements that are adjacent to fracture tips, and have boundaries that are excluded from the contact region. Stress intensity factors are computed numerically using the methods of displacement correlation and disk-shaped domain integral. A novel square-root singular variation of the penalty parameter near the crack front is proposed to accurately model the contact tractions near the crack front. Tractions and compressive stress intensity factors are validated against analytical solutions. Numerical examples of cubes containing one, two, twenty four and seventy interacting and intersecting fractures are presented