Robust Procedures for Obtaining Assembly Contact State Extremal Configurations

Two important components in the selection of an admittance that facilitates force-guided assembly are the identification of: 1) the set of feasible contact states, and 2) the set of configurations that span each contact state, i.e., the extremal configurations. We present a procedure to automatically generate both sets from CAD models of the assembly parts. In the procedure, all possible combinations of principle contacts are considered when generating hypothesized contact states. The feasibility of each is then evaluated in a genetic algorithm based optimization procedure. The maximum and minimum value of each of the 6 configuration variables spanning each contact state are obtained by again using genetic algorithms. Together, the genetic algorithm approach, the hierarchical data structure containing the states, the relationships among the states, and the extremals within each state are used to provide a reliable means of identifying all feasible contact states and their associated extremal configurations

The Double Hierarchy Method: a parallel 3D contact method for the interaction of spherical particles with rigid FE boundaries using the DEM

The final publication is available at Springer via http://dx.doi.org/10.1007/s40571-016-0109-4In this work, we present a new methodology for the treatment of the contact interaction between rigid boundaries and spherical discrete elements (DE). Rigid body parts are present in most of large-scale simulations. The surfaces of the rigid parts are commonly meshed with a finite element-like (FE) discretization. The contact detection and calculation between those DE and the discretized boundaries is not straightforward and has been addressed by different approaches. The algorithm presented in this paper considers the contact of the DEs with the geometric primitives of a FE mesh, i.e. facet, edge or vertex. To do so, the original hierarchical method presented by Horner et al. (J Eng Mech 127(10):1027–1032, 2001) is extended with a new insight leading to a robust, fast and accurate 3D contact algorithm which is fully parallelizable. The implementation of the method has been developed in order to deal ideally with triangles and quadrilaterals. If the boundaries are discretized with another type of geometries, the method can be easily extended to higher order planar convex polyhedra. A detailed description of the procedure followed to treat a wide range of cases is presented. The description of the developed algorithm and its validation is verified with several practical examples. The parallelization capabilities and the obtained performance are presented with the study of an industrial application example.Peer ReviewedPostprint (author's final draft

Articular human joint modelling

Copyright @ Cambridge University Press 2009.The work reported in this paper encapsulates the theories and algorithms developed to drive the core analysis modules of the software which has been developed to model a musculoskeletal structure of anatomic joints. Due to local bone surface and contact geometry based joint kinematics, newly developed algorithms make the proposed modeller different from currently available modellers. There are many modellers that are capable of modelling gross human body motion. Nevertheless, none of the available modellers offer complete elements of joint modelling. It appears that joint modelling is an extension of their core analysis capability, which, in every case, appears to be musculoskeletal motion dynamics. It is felt that an analysis framework that is focused on human joints would have significant benefit and potential to be used in many orthopaedic applications. The local mobility of joints has a significant influence in human motion analysis, in understanding of joint loading, tissue behaviour and contact forces. However, in order to develop a bone surface based joint modeller, there are a number of major problems, from tissue idealizations to surface geometry discretization and non-linear motion analysis. This paper presents the following: (a) The physical deformation of biological tissues as linear or non-linear viscoelastic deformation, based on spring-dashpot elements. (b) The linear dynamic multibody modelling, where the linear formulation is established for small motions and is particularly useful for calculating the equilibrium position of the joint. This model can also be used for finding small motion behaviour or loading under static conditions. It also has the potential of quantifying the joint laxity. (c) The non-linear dynamic multibody modelling, where a non-matrix and algorithmic formulation is presented. The approach allows handling complex material and geometrical nonlinearity easily. (d) Shortest path algorithms for calculating soft tissue line of action geometries. The developed algorithms are based on calculating minimum ‘surface mass’ and ‘surface covariance’. An improved version of the ‘surface covariance’ algorithm is described as ‘residual covariance’. The resulting path is used to establish the direction of forces and moments acting on joints. This information is needed for linear or non-linear treatment of the joint motion. (e) The final contribution of the paper is the treatment of the collision. In the virtual world, the difficulty in analysing bodies in motion arises due to body interpenetrations. The collision algorithm proposed in the paper involves finding the shortest projected ray from one body to the other. The projection of the body is determined by the resultant forces acting on it due to soft tissue connections under tension. This enables the calculation of collision condition of non-convex objects accurately. After the initial collision detection, the analysis involves attaching special springs (stiffness only normal to the surfaces) at the ‘potentially colliding points’ and motion of bodies is recalculated. The collision algorithm incorporates the rotation as well as translation. The algorithm continues until the joint equilibrium is achieved. Finally, the results obtained based on the software are compared with experimental results obtained using cadaveric joints

A New Framework Based on a Discrete Element Method to Model the Fracture Behavior for Brittle Polycrystalline Materials

This work aims to develop and implement a linear elastic grain-level micromechanical model based on the discrete element method using bonded contacts and an improved fracture criteria to capture both intergranular and transgranular microcrack initiation and evolution in polycrystalline ceramics materials. Gaining a better understanding of the underlying mechanics and micromechanics of the fracture process of brittle polycrystalline materials will aid in high performance material design. Continuum mechanics approaches cannot accurately simulate the crack propagation during fracture due to the discontinuous nature of the problem. In this work we distinguish between predominately intergranular failure (along the grain boundaries) versus predominately transgranular failure (across the grains) based on grain orientation and microstructural parameters to describe the contact interfaces and present the first approach at fracturing discrete elements. Specifically, the influence of grain boundary strength and stiffness on the fracture behavior of an idealized ceramic material is studied under three different loading conditions: uniaxial compression, brazilian, and four-point bending. Digital representations of the sample microstructures for the test cases are composed of hexagonal, prismatic, honeycomb-packed grains represented by rigid, discrete elements. The principle of virtual work is used to develop a microscale fracture criteria for brittle polycrystalline materials for tensile, shear, torsional and rolling modes of intergranular motion. The interactions between discrete elements within each grain are governed by traction displacement relationships

Grasping and Assembling with Modular Robots

A wide variety of problems, from manufacturing to disaster response and space exploration, can benefit from robotic systems that can firmly grasp objects or assemble various structures, particularly in difficult, dangerous environments. In this thesis, we study the two problems, robotic grasping and assembly, with a modular robotic approach that can facilitate the problems with versatility and robustness. First, this thesis develops a theoretical framework for grasping objects with customized effectors that have curved contact surfaces, with applications to modular robots. We present a collection of grasps and cages that can effectively restrain the mobility of a wide range of objects including polyhedra. Each of the grasps or cages is formed by at most three effectors. A stable grasp is obtained by simple motion planning and control. Based on the theory, we create a robotic system comprised of a modular manipulator equipped with customized end-effectors and a software suite for planning and control of the manipulator. Second, this thesis presents efficient assembly planning algorithms for constructing planar target structures collectively with a collection of homogeneous mobile modular robots. The algorithms are provably correct and address arbitrary target structures that may include internal holes. The resultant assembly plan supports parallel assembly and guarantees easy accessibility in the sense that a robot does not have to pass through a narrow gap while approaching its target position. Finally, we extend the algorithms to address various symmetric patterns formed by a collection of congruent rectangles on the plane. The basic ideas in this thesis have broad applications to manufacturing (restraint), humanitarian missions (forming airfields on the high seas), and service robotics (grasping and manipulation)

The Stability of Heavy Objects with Multiple Contacts

In both robot grasping and robot locomotion, we wish to hold objects stably in the presence of gravity. We present a derivation of second-order stability conditions for a supported heavy object, employing the tool of Stratified Morse theory. We then apply these general results to the case of objects in the plane