The application of contact mechanics to the numerical simulation of particulate material

Particulate solids are complex redundant systems which consist of discrete particles. The interactions between the particles are complex and have been the subject of many theoretical and experimental investigations. Invetigations of particulate material have been restricted by the lack of quantitative information on the mechanisms occurring within an assembly. Laboratory experimentation is limited as information on the internal behaviour can only be inferred from measurements on the assembly boundary, or the use of intrusive measuring devices. In addition comparisons between test data are uncertain due to the difficulty in reproducing exact replicas of physical systems. Nevertheless, theoretical and technological advances require more detailed material information. However, numerical simulation affords access to information on every particle and hence the micro-mechanical behaviour within an assembly, and can replicate desired systems. To use a computer program to numerically simulate material behaviour accurately it is necessary to incorporte realistic interaction laws. This research programme used the finite difference simulation program `BALL', developed by Cundall (1971), which employed linear spring force-displacement laws. It was thus necessary to incorporate more realistic interaction laws. Therefore, this research programme was primarily concerned with the implementation of the normal force-displacement law of Hertz (1882) and the tangential force-displacement laws of Mindlin and Deresiewicz (1953). Within this thesis the contact mechanics theories employed in the program are developed and the adaptations which were necessary to incorporate these laws are detailed. Verification of the new contact force-displacement laws was achieved by simulating a quasi-static oblique contact and single particle oblique impact. Applications of the program to the simulation of large assemblies of particles is given, and the problems in undertaking quasi-static shear tests along with the results from two successful shear tests are described.

Computer simulation of moist agglomerate collisions

This thesis considers the computer simulation of moist agglomerate collisions using the discrete element method (DEM). The study is confined to pendular state moist agglomerates, at which liquid is presented as either absorbed immobile films or pendular liquid bridges and the interparticle force is modelled as the adhesive contact force and interstitial liquid bridge force. Algorithms used to model the contact force due to surface adhesion, tangential friction and particle deformation have been derived by other researchers and are briefly described in the thesis. A theoretical study of the pendular liquid bridge force between spherical particles has been made and the algorithms for the modelling of the pendular liquid bridge force between spherical particles have been developed and incorporated into the Aston version of the DEM program TRUBAL. It has been found that, for static liquid bridges, the more explicit criterion for specifying the stable solution and critical separation is provided by the total free energy. The critical separation is given by the cube root of liquid bridge volume to a good approximation and the 'gorge method' of evaluation based on the toroidal approximation leads to errors in the calculated force of less than 10%. Three dimensional computer simulations of an agglomerate impacting orthogonally with a wall are reported. The results demonstrate the effectiveness of adding viscous binder to prevent attrition, a common practice in process engineering. Results of simulated agglomerate-agglomerate collisions show that, for colinear agglomerate impacts, there is an optimum velocity which results in a near spherical shape of the coalesced agglomerate and, hence, minimises attrition due to subsequent collisions. The relationship between the optimum impact velocity and the liquid viscosity and surface tension is illustrated. The effect of varying the angle of impact on the coalescence/attrition behaviour is also reported. (DX 187, 340).

Consensus on recording of gas permeable contact lens fit

Purpose: To develop a new schematic scheme for efficiently recording the key parameters of gas permeable contact lens (GP) fits based on current consensus. Methods: Over 100 established GP fitters and educators met to discuss the parameters proposed in educational material for evaluating GP fit and concluded on the key parameters that should be recorded. The accuracy and variability of evaluating the fluorescein pattern of GP fit was determined by having 35 experienced contact lens practitioners from across the world, grading 5 images of a range of fits and the topographer simulation of the same fits, in random, order using the proposed scheme. The accuracy of the grading was compared to objective image analysis of the fluorescein intensity of the same images. Results: The key information to record to adequately describe the fit of an GP was agreed as: the manufacturer, brand and lens parameters; settling time; comfort on a 5 point scale; centration; movement on blink on a ±2 scale; and the Primary Fluorescein Pattern in the central, mid-peripheral and edge regions of the lens averaged along the horizontal and vertical lens axes, on a ±2 scale. On average 50-60% of practitioners selected the median grade when subjectively rating fluorescein intensity and this was correlated to objective quantification (r= 0.602, p< 0.001). Objective grading suggesting horizontal median fluorescein intensity was generally symmetrical, as was the vertical meridian, but this was not the case for subjective grading. Simulated fluorescein patterns were subjectively and objectively graded as being less intense than real photographs (p< 0.01). Conclusion: GP fit recording can be standardised and simplified to enhance GP practice. © 2013 British Contact Lens Association.