Application of the finite element method to problems in fracture mechanics

Numerical techniques have been finding increasing use in all aspects of fracture mechanics, and often provide the only means for analyzing fracture problems. The work presented here, is concerned with the application of the finite element method to cracked structures. The present work was directed towards the establishment of a comprehensive two-dimensional finite element, linear elastic, fracture analysis package. Significant progress has been made to this end, and features which can now be studied include multi-crack tip mixed-mode problems, involving partial crack closure. The crack tip core element was refined and special local crack tip elements were employed to reduce the element density in the neighbourhood of the core region. The work builds upon experience gained by previous research workers and, as part of the general development, the program was modified to incorporate the eight-node isoparametric quadrilateral element. Also. a more flexible solving routine was developed, and provided a very compact method of solving large sets of simultaneous equations, stored in a segmented form. To complement the finite element analysis programs, an automatic mesh generation program has been developed, which enables complex problems. involving fine element detail, to be investigated with a minimum of input data. The scheme has proven to be versati Ie and reasonably easy to implement. Numerous examples are given to demonstrate the accuracy and flexibility of the finite element technique.

Mechanics of longitudinal rolling of tube through three-grooved rolls

This investigation examined the process of the longitudinal rolling of tubes through a set of three driven grooved rolls. Tubes were rolled with or without internal support i.e. under mandrel rolling or sinking conditions. Knowledge was required of the way in which the roll separating force and rolling torque vary for different conditions of rolling. The objective of this work being to obtain a better understanding and optimization of the mechanics of the process. The design and instrumentation of a complete experimental three-roll mill for the rolling of lead tube as an analogue material for hot steel, with the measurement of the individual roll force and torque is described. A novel type of roll load cell was incorporated and its design and testing discussed. Employing three roll sizes of 170 mm, 255 mm and 340 mm shroud diameter, precise tube specimens of various tube diameter to thickness ratios were rolled under sinking and mandrel rolling conditions. To obtain an indication of the tube-roll contact areas some of the specimens were partially rolled. For comparative purposes the remaining tubes were completely rolled as a single pass. The roll forces, torques and tube parameters e.g. reduction of area, D/t ratio, were collated and compared for each of the three roll diameters considered. The influence of friction, particularly in the mandrel rolling process, was commented upon. Theoretical studies utilising the equilibrium and energy methods were applied to both the sinking and mandrel rolling processes. In general, the energy approach gave better comparison with experiment, especially for mandrel rolling. The influence of the tube deformation zones on the two processes was observed and on the subsequent modification of the tube-roll arc contact length. A rudimentary attempt was made in the theoretical sinking analysis to allow for the deformation zone prior to roll contact; some success was noted. A general survey of the available tube rolling literature, for both the sinking and mandrel processes has been carried out.

The mechanics of ultrasonic tube bending

The technology of precision bending of tubes has recently increased in importance and is widely demanded for many industrial applications. However, whilst attention has been concentrated on automation and increasing the production rate of the bending machines, it seems that with one exception very little work has been done in order to understand and therefore fundamentally improve the bending process. A new development for the process of draw-bending of tubes, in which the supporting mandrel is axially vibrated at an ultrasonic frequency, has been perfected. A research programme was undertaken to study the mechanics of tube• bending under both vibratory and non-vibratory conditions. For this purpose, a conventional tube-bending machine was modified and equipped with an oscillatory system. Thin-walled mild steel tubes of different diameter to thickness ratios were bent to mean bend radii having various values from 1.5 to 2.0 times the tube diameter. It was found that the application of ultrasonic vibration reduces the process forces and that the force reduction increases with increasing the vibration amplitude. A reduction in the bending torque of up to 30 per cent was recorded and a reduction in the maximum tube-wall thinning of about 15 per cent was observed. The friction vector reversal mechanism as well as a reduction in friction account for the changes of the forces and the strains. Monitoring the mandrel friction during bending showed, in some cases, that the axial vibration reverses the mandrel .mean force from tension to compression and, thus, the mandrel is assisting the tube motion instead of resisting it. A theory has been proposed to describe the mechanics of deformation during draw-bending of tubes, which embodies the conditions of both "with" and "without" mandrel axial vibration. A theoretical analysis, based on the equilibrium of forces approach, has been developed in which the basic process parameters were taken into consideration. The stresses, the strains and the bending torque were calculated utilising this new solution, and a specially written computer programme was used to perform the computations. It was shown that the theory is in good agreement with the measured values of the strains under vibratory and non-vibratory conditions. Also, the predicted bending 'torque showed a similar trend to that recorded experimentally.

A study of the micro-mechanics of granular media

The development of more realistic constitutive models for granular media, such as sand, requires ingredients which take into account the internal micro-mechanical response to deformation. Unfortunately, at present, very little is known about these mechanisms and therefore it is instructive to find out more about the internal nature of granular samples by conducting suitable tests. In contrast to physical testing the method of investigation used in this study employs the Distinct Element Method. This is a computer based, iterative, time-dependent technique that allows the deformation of granular assemblies to be numerically simulated. By making assumptions regarding contact stiffnesses each individual contact force can be measured and by resolution particle centroid forces can be calculated. Then by dividing particle forces by their respective mass, particle centroid velocities and displacements are obtained by numerical integration. The Distinct Element Method is incorporated into a computer program 'Ball'. This program is effectively a numerical apparatus which forms a logical housing for this method and allows data input and output, and also provides testing control. By using this numerical apparatus tests have been carried out on disc assemblies and many new interesting observations regarding the micromechanical behaviour are revealed. In order to relate the observed microscopic mechanisms of deformation to the flow of the granular system two separate approaches have been used. Firstly a constitutive model has been developed which describes the yield function, flow rule and translation rule for regular assemblies of spheres and discs when subjected to coaxial deformation. Secondly statistical analyses have been carried out using data which was extracted from the simulation tests. These analyses define and quantify granular structure and then show how the force and velocity distributions use the structure to produce the corresponding stress and strain-rate tensors.

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.