Load Transfer from a Smooth Elastic Pin to a Large Sheet

The plane problem of load transfer from an elastic interference or clearance fit pin to a large elastic sheet with a perfectly smooth interface is solved. As the load on the pin is monotonically increased, the pin-hole interface is in partial contact above certain critical load in interference fit and throughout the loading range in clearance fit.Such situations result in mixed boundary-value problems with moving boundaries and the arc of contact varies nonlinearly with applied load. These problems are analyzed by an inverse technique in which the arcs of contact/separation are prescribed and the causative loads are evaluated. A direct method of analysis is adopted using biharmonic polar trigonometric stress functions and a simple collocation method for satisfying the boundary conditions. A unified analytical formulation is achieved for interference and clearance fits. The solutions for the linear problem of push fits are inherent in the unified analysis. Numerical results highlighting the effects of pin and sheet elasticity parameters are presented.

Wet and dry model granulates under mechanical load : a confocal microscopy study

Granular matter, also known as bulk solids, consists of discrete particles with sizes between micrometers and meters. They are present in many industrial applications as well as daily life, like in food processing, pharmaceutics or in the oil and mining industry. When handling granular matter the bulk solids are stored, mixed, conveyed or filtered. These techniques are based on observations in macroscopic experiments, i.e. rheological examinations of the bulk properties. Despite the amply investigations of bulk mechanics, the relation between single particle motion and macroscopic behavior is still not well understood. For exploring the microscopic properties on a single particle level, 3D imaging techniques are required.rnThe objective of this work was the investigation of single particle motions in a bulk system in 3D under an external mechanical load, i.e. compression and shear. During the mechanical load the structural and dynamical properties of these systems were examined with confocal microscopy. Therefor new granular model systems in the wet and dry state were designed and prepared. As the particles are solid bodies, their motion is described by six degrees of freedom. To explore their entire motion with all degrees of freedom, a technique to visualize the rotation of spherical micrometer sized particles in 3D was developed. rnOne of the foci during this dissertation was a model system for dry cohesive granular matter. In such systems the particle motion during a compression of the granular matter was investigated. In general the rotation of single particles was the more sensitive parameter compared to the translation. In regions with large structural changes the rotation had an earlier onset than the translation. In granular systems under shear, shear dilatation and shear zone formation were observed. Globally the granular sediments showed a shear behavior, which was known already from classical shear experiments, for example with Jenike cells. Locally the shear zone formation was enhanced, when near the applied load a pre-diluted region existed. In regions with constant volume fraction a mixing between the different particle layers occurred. In particular an exchange of particles between the current flowing region and the non-flowing region was observed. rnThe second focus was on model systems for wet granular matter, where an additional binding liquid is added to the particle suspension. To examine the 3D structure of the binding liquid on the micrometer scale independently from the particles, a second illumination and detection beam path was implemented. In shear and compression experiments of wet clusters and bulk systems completely different dynamics compared to dry cohesive models systems occured. In a Pickering emulsion-like system large structural changes predominantly occurred in the local environment of binding liquid droplets. These large local structural changes were due to an energy interplay between the energy stored in the binding droplet during its deformation and the binding energy of particles at the droplet interface. rnConfocal microscopy in combination with nanoindentation gave new insights into the single particle motions and dynamics of granular systems under a mechanical load. These novel experimental results can help to improve the understanding of the relationship between bulk properties of granular matter, such as volume fraction or yield stress and the dynamics on a single particle level.rnrn

Multifunctional cement composites strain and damage sensors applied on reinforced concrete (RC) structural elements

In this research, strain-sensing and damage-sensing functional properties of cement composites have been studied on a conventional reinforced concrete (RC) beam. Carbon nanofiber (CNFCC) and fiber (CFCC) cement composites were used as sensors on a 4 m long RC beam. Different casting conditions (in situ or attached), service location (under tension or compression) and electrical contacts (embedded or superficial) were compared. Both CNFCC and CFCC were suitable as strain sensors in reversible (elastic) sensing condition testing. CNFCC showed higher sensitivities (gage factor up to 191.8), while CFCC only reached gage factors values of 178.9 (tension) or 49.5 (compression). Furthermore, damage-sensing tests were run, increasing the applied load progressively up to the RC beam failure. In these conditions, CNFCC sensors were also strain sensitive, but no damage sensing mechanism was detected for the strain levels achieved during the tests. Hence, these cement composites could act as strain sensors, even for severe damaged structures near to their collapse.

Sensitivity of proximal femoral stiffness and areal bone mineral density to changes in bone geometry and density

Areal bone mineral density (aBMD) is the most common surrogate measurement for assessing the bone strength of the proximal femur associated with osteoporosis. Additional factors, however, contribute to the overall strength of the proximal femur, primarily the anatomical geometry. Finite element analysis (FEA) is an effective and widely used computerbased simulation technique for modeling mechanical loading of various engineering structures, providing predictions of displacement and induced stress distribution due to the applied load. FEA is therefore inherently dependent upon both density and anatomical geometry. FEA may be performed on both three-dimensional and two-dimensional models of the proximal femur derived from radiographic images, from which the mechanical stiffness may be redicted. It is examined whether the outcome measures of two-dimensional FEA, two-dimensional, finite element analysis of X-ray images (FEXI), and three-dimensional FEA computed stiffness of the proximal femur were more sensitive than aBMD to changes in trabecular bone density and femur geometry. It is assumed that if an outcome measure follows known trends with changes in density and geometric parameters, then an increased sensitivity will be indicative of an improved prediction of bone strength. All three outcome measures increased non-linearly with trabecular bone density, increased linearly with cortical shell thickness and neck width, decreased linearly with neck length, and were relatively insensitive to neck-shaft angle. For femoral head radius, aBMD was relatively insensitive, with two-dimensional FEXI and threedimensional FEA demonstrating a non-linear increase and decrease in sensitivity, respectively. For neck anteversion, aBMD decreased non-linearly, whereas both two-dimensional FEXI and three dimensional FEA demonstrated a parabolic-type relationship, with maximum stiffness achieved at an angle of approximately 15o. Multi-parameter analysis showed that all three outcome measures demonstrated their highest sensitivity to a change in cortical thickness. When changes in all input parameters were considered simultaneously, three and twodimensional FEA had statistically equal sensitivities (0.41±0.20 and 0.42±0.16 respectively, p = ns) that were significantly higher than the sensitivity of aBMD (0.24±0.07; p = 0.014 and 0.002 for three-dimensional and two-dimensional FEA respectively). This simulation study suggests that since mechanical integrity and FEA are inherently dependent upon anatomical geometry, FEXI stiffness, being derived from conventional two-dimensional radiographic images, may provide an improvement in the prediction of bone strength of the proximal femur than currently provided by aBMD.

Design and evaluation of an MRI compatible axial compression device for 3D assessment of spinal deformity and flexibility

Magnetic Resonance Imaging (MRI) offers a valuable research tool for the assessment of 3D spinal deformity in AIS, however the horizontal patient position imposed by conventional scanners removes the axial compressive loading on the spine which is an important determinant of deformity shape and magnitude in standing scoliosis patients. The objective of this study was to design, construct and test an MRI compatible compression device for research into the effect of axial loading on spinal deformity using supine MRI scans. The compression device was designed and constructed, consisting of a vest worn by the patient, which was attached via straps to a pneumatically actuated footplate. An applied load of 0.5 x bodyweight was remotely controlled by a unit in the scanner operator’s console. The entire device was constructed using non-metallic components for MRI compatibility. The device was evaluated by performing unloaded and loaded supine MRI scans on a series of 10 AIS patients. The study concluded that an MRI compatible compression device had been successfully designed and constructed, providing a research tool for studies into the effect of axial loading on 3D spinal deformity in scoliosis. The 3D axially loaded MR imaging capability developed in this study will allow future research investigations of the effect of axial loading on spinal rotation, and for imaging the response of scoliotic spinal tissues to axial loading.

Development and evaluation of hydrocarbon sorbent materials with the aid of chemometrics

Hydrocarbon spills on roads are a major safety concern for the driving public and can have severe cost impacts both on pavement maintenance and to the economy through disruption to services. The time taken to clean-up spills and re-open roads in a safe driving condition is an issue of increasing concern given traffic levels on major urban arterials. Thus, the primary aim of the research was to develop a sorbent material that facilitates rapid clean-up of road spills. The methodology involved extensive research into a range of materials (organic, inorganic and synthetic sorbents), comprehensive testing in the laboratory, scale-up and field, and product design (i.e. concept to prototype). The study also applied chemometrics to provide consistent, comparative methods of sorbent evaluation and performance. In addition, sorbent materials at every stage were compared against a commercial benchmark. For the first time, the impact of diesel on asphalt pavement has been quantified and assessed in a systematic way. Contrary to conventional thinking and anecdotal observations, the study determined that the action of diesel on asphalt was quite rapid (i.e. hours rather than weeks or months). This significant finding demonstrates the need to minimise the impact of hydrocarbon spills and the potential application of the sorbent option. To better understand the adsorption phenomenon, surface characterisation techniques were applied to selected sorbent materials (i.e. sand, organo-clay and cotton fibre). Brunauer Emmett Teller (BET) and thermal analysis indicated that the main adsorption mechanism for the sorbents occurred on the external surface of the material in the diffusion region (sand and organo-clay) and/or capillaries (cotton fibre). Using environmental scanning electron microscopy (ESEM), it was observed that adsorption by the interfibre capillaries contributed to the high uptake of hydrocarbons by the cotton fibre. Understanding the adsorption mechanism for these sorbents provided some guidance and scientific basis for the selection of materials. The study determined that non-woven cotton mats were ideal sorbent materials for clean-up of hydrocarbon spills. The prototype sorbent was found to perform significantly better than the commercial benchmark, displaying the following key properties: • superior hydrocarbon pick-up from the road pavement; • high hydrocarbon retention capacity under an applied load; • adequate field skid resistance post treatment; • functional and easy to use in the field (e.g. routine handling, transportation, application and recovery); • relatively inexpensive to produce due to the use of raw cotton fibre and simple production process; • environmentally friendly (e.g. renewable materials, non-toxic to environment and operators, and biodegradable); and • rapid response time (e.g. two minutes total clean-up time compared with thirty minutes for reference sorbents). The major outcomes of the research project include: a) development of a specifically designed sorbent material suitable for cleaning up hydrocarbon spills on roads; b) submission of patent application (serial number AU2005905850) for the prototype product; and c) preparation of Commercialisation Strategy to advance the sorbent product to the next phase (i.e. R&D to product commercialisation).

The feasibility of vibration analysis as a technique to detect osseointegration of transfemoral implants

One of the main causes of above knee or transfemoral amputation (TFA) in the developed world is trauma to the limb. The number of people undergoing TFA due to limb trauma, particularly due to war injuries, has been increasing. Typically the trauma amputee population, including war-related amputees, are otherwise healthy, active and desire to return to employment and their usual lifestyle. Consequently there is a growing need to restore long-term mobility and limb function to this population. Traditionally transfemoral amputees are provided with an artificial or prosthetic leg that consists of a fabricated socket, knee joint mechanism and a prosthetic foot. Amputees have reported several problems related to the socket of their prosthetic limb. These include pain in the residual limb, poor socket fit, discomfort and poor mobility. Removing the socket from the prosthetic limb could eliminate or reduce these problems. A solution to this is the direct attachment of the prosthesis to the residual bone (femur) inside the residual limb. This technique has been used on a small population of transfemoral amputees since 1990. A threaded titanium implant is screwed in to the shaft of the femur and a second component connects between the implant and the prosthesis. A period of time is required to allow the implant to become fully attached to the bone, called osseointegration (OI), and be able to withstand applied load; then the prosthesis can be attached. The advantages of transfemoral osseointegration (TFOI) over conventional prosthetic sockets include better hip mobility, sitting comfort and prosthetic retention and fewer skin problems on the residual limb. However, due to the length of time required for OI to progress and to complete the rehabilitation exercises, it can take up to twelve months after implant insertion for an amputee to be able to load bear and to walk unaided. The long rehabilitation time is a significant disadvantage of TFOI and may be impeding the wider adoption of the technique. There is a need for a non-invasive method of assessing the degree of osseointegration between the bone and the implant. If such a method was capable of determining the progression of TFOI and assessing when the implant was able to withstand physiological load it could reduce the overall rehabilitation time. Vibration analysis has been suggested as a potential technique: it is a non destructive method of assessing the dynamic properties of a structure. Changes in the physical properties of a structure can be identified from changes in its dynamic properties. Consequently vibration analysis, both experimental and computational, has been used to assess bone fracture healing, prosthetic hip loosening and dental implant OI with varying degrees of success. More recently experimental vibration analysis has been used in TFOI. However further work is needed to assess the potential of the technique and fully characterise the femur-implant system. The overall aim of this study was to develop physical and computational models of the TFOI femur-implant system and use these models to investigate the feasibility of vibration analysis to detect the process of OI. Femur-implant physical models were developed and manufactured using synthetic materials to represent four key stages of OI development (identified from a physiological model), simulated using different interface conditions between the implant and femur. Experimental vibration analysis (modal analysis) was then conducted using the physical models. The femur-implant models, representing stage one to stage four of OI development, were excited and the modal parameters obtained over the range 0-5kHz. The results indicated the technique had limited capability in distinguishing between different interface conditions. The fundamental bending mode did not alter with interfacial changes. However higher modes were able to track chronological changes in interface condition by the change in natural frequency, although no one modal parameter could uniquely distinguish between each interface condition. The importance of the model boundary condition (how the model is constrained) was the key finding; variations in the boundary condition altered the modal parameters obtained. Therefore the boundary conditions need to be held constant between tests in order for the detected modal parameter changes to be attributed to interface condition changes. A three dimensional Finite Element (FE) model of the femur-implant model was then developed and used to explore the sensitivity of the modal parameters to more subtle interfacial and boundary condition changes. The FE model was created using the synthetic femur geometry and an approximation of the implant geometry. The natural frequencies of the FE model were found to match the experimental frequencies within 20% and the FE and experimental mode shapes were similar. Therefore the FE model was shown to successfully capture the dynamic response of the physical system. As was found with the experimental modal analysis, the fundamental bending mode of the FE model did not alter due to changes in interface elastic modulus. Axial and torsional modes were identified by the FE model that were not detected experimentally; the torsional mode exhibited the largest frequency change due to interfacial changes (103% between the lower and upper limits of the interface modulus range). Therefore the FE model provided additional information on the dynamic response of the system and was complementary to the experimental model. The small changes in natural frequency over a large range of interface region elastic moduli indicated the method may only be able to distinguish between early and late OI progression. The boundary conditions applied to the FE model influenced the modal parameters to a far greater extent than the interface condition variations. Therefore the FE model, as well as the experimental modal analysis, indicated that the boundary conditions need to be held constant between tests in order for the detected changes in modal parameters to be attributed to interface condition changes alone. The results of this study suggest that in a clinical setting it is unlikely that the in vivo boundary conditions of the amputated femur could be adequately controlled or replicated over time and consequently it is unlikely that any longitudinal change in frequency detected by the modal analysis technique could be attributed exclusively to changes at the femur-implant interface. Therefore further development of the modal analysis technique would require significant consideration of the clinical boundary conditions and investigation of modes other than the bending modes.

Mechanical considerations in fracture fixation

Bone’s capacity to repair following trauma is both unique and astounding. However, fractures sometimes fail to heal. Hence, the goal of fracture treatment is the restoration of bone’s structure, composition and function. Fracture fixation devices should provide a favourable mechanical and biological environment for healing to occur. The use of internal fixation is increasing as these devices may be applied with less invasive techniques. Recent studies suggest however that, internal fixation devices may be overly stiff and suppresses callus formation. The degree of mechanical stability influences the healing outcome. This is determined by the stiffness of the fixation device and the degree of limb loading. This project aims to characterise the fixation stability of an internal plate fixation device and the influence of modifications to its configuration on implant stability. As there are no standardised methods for the determination of fixation stiffness, the first part of this project aims to compares different methodologies and determines the most appropriate method to characterise the stiffness of internal plate fixators. The stiffness of a fixation device also influences the physiological loads experienced by the healing bone. Since bone adapts to this applied load by undergoing changes through a remodelling process, undesirable changes could occur during the period of treatment with an implant. The second part of this project aims to develop a methodology to quantify remodelling changes. This quantification is expected to aid our understanding of the changes in pattern due to implant related remodelling and on the factors driving the remodelling process. Knowledge gained in this project is useful to understand how the configuration of internal fixation devices can promote timely healing and prevent undesirable bone loss.

Integrated health prediction of bridge systems using dynamic object oriented Bayesian networks (DOOBNS)

The serviceability and safety of bridges are crucial to people’s daily lives and to the national economy. Every effort should be taken to make sure that bridges function safely and properly as any damage or fault during the service life can lead to transport paralysis, catastrophic loss of property or even casualties. Nonetheless, aggressive environmental conditions, ever-increasing and changing traffic loads and aging can all contribute to bridge deterioration. With often constrained budget, it is of significance to identify bridges and bridge elements that should be given higher priority for maintenance, rehabilitation or replacement, and to select optimal strategy. Bridge health prediction is an essential underpinning science to bridge maintenance optimization, since the effectiveness of optimal maintenance decision is largely dependent on the forecasting accuracy of bridge health performance. The current approaches for bridge health prediction can be categorised into two groups: condition ratings based and structural reliability based. A comprehensive literature review has revealed the following limitations of the current modelling approaches: (1) it is not evident in literature to date that any integrated approaches exist for modelling both serviceability and safety aspects so that both performance criteria can be evaluated coherently; (2) complex system modelling approaches have not been successfully applied to bridge deterioration modelling though a bridge is a complex system composed of many inter-related bridge elements; (3) multiple bridge deterioration factors, such as deterioration dependencies among different bridge elements, observed information, maintenance actions and environmental effects have not been considered jointly; (4) the existing approaches are lacking in Bayesian updating ability to incorporate a variety of event information; (5) the assumption of series and/or parallel relationship for bridge level reliability is always held in all structural reliability estimation of bridge systems. To address the deficiencies listed above, this research proposes three novel models based on the Dynamic Object Oriented Bayesian Networks (DOOBNs) approach. Model I aims to address bridge deterioration in serviceability using condition ratings as the health index. The bridge deterioration is represented in a hierarchical relationship, in accordance with the physical structure, so that the contribution of each bridge element to bridge deterioration can be tracked. A discrete-time Markov process is employed to model deterioration of bridge elements over time. In Model II, bridge deterioration in terms of safety is addressed. The structural reliability of bridge systems is estimated from bridge elements to the entire bridge. By means of conditional probability tables (CPTs), not only series-parallel relationship but also complex probabilistic relationship in bridge systems can be effectively modelled. The structural reliability of each bridge element is evaluated from its limit state functions, considering the probability distributions of resistance and applied load. Both Models I and II are designed in three steps: modelling consideration, DOOBN development and parameters estimation. Model III integrates Models I and II to address bridge health performance in both serviceability and safety aspects jointly. The modelling of bridge ratings is modified so that every basic modelling unit denotes one physical bridge element. According to the specific materials used, the integration of condition ratings and structural reliability is implemented through critical failure modes. Three case studies have been conducted to validate the proposed models, respectively. Carefully selected data and knowledge from bridge experts, the National Bridge Inventory (NBI) and existing literature were utilised for model validation. In addition, event information was generated using simulation to demonstrate the Bayesian updating ability of the proposed models. The prediction results of condition ratings and structural reliability were presented and interpreted for basic bridge elements and the whole bridge system. The results obtained from Model II were compared with the ones obtained from traditional structural reliability methods. Overall, the prediction results demonstrate the feasibility of the proposed modelling approach for bridge health prediction and underpin the assertion that the three models can be used separately or integrated and are more effective than the current bridge deterioration modelling approaches. The primary contribution of this work is to enhance the knowledge in the field of bridge health prediction, where more comprehensive health performance in both serviceability and safety aspects are addressed jointly. The proposed models, characterised by probabilistic representation of bridge deterioration in hierarchical ways, demonstrated the effectiveness and pledge of DOOBNs approach to bridge health management. Additionally, the proposed models have significant potential for bridge maintenance optimization. Working together with advanced monitoring and inspection techniques, and a comprehensive bridge inventory, the proposed models can be used by bridge practitioners to achieve increased serviceability and safety as well as maintenance cost effectiveness.

Direct measurement as a tool for assessing the accuracy of inverse dynamics

The understanding of the loads generated within the prosthetic leg can aid engineers in the design of components and clinicians in the process of rehabilitation. Traditional methods to assess these loads have relied on inverse dynamics. This indirect method estimates the applied load using video recordings and force-plates located at a distance from the region of interest, such as the base of the residuum. The well-known limitations of this method are related to the accuracy of this recursive model and the experimental conditions required (Frossard et al., 2003). Recent developments in sensors (Frossard et al., 2003) and prosthetic fixation (Brånemark et al., 2000) permit the direct measurement of the loads applied on the residuum of transfemoral amputees. In principle, direct measurement should be an appropriate tool for assessing the accuracy of inverse dynamics. The purpose of this paper is to determine the validity of this assumption. The comparative variable used in this study is the velocity of the relative body center of mass (VCOM(t)). The relativity is used to align the static (w.r.t. position) force plate measurement with the dynamic load cell measurement.

Dynamic initiation toughness measu~ment using double Grooved discs

Among the multitude of test specimen geometries used for dynamic fiacture toughness evaluation, the most widely uscd specimen is lhc Chavpy specimen due its simple geomclry and availability of testing machines. The standard Chatpy specimen dimensions may llOl always give plane st~ain condilions and hence, it may be necessary Io coilduct lcs/s using specimens of dillEvcnt thicknesses to establish the plane strain K~a. An axisymmct/ic specimen, on the otlaev hand would always give flow constraints l~n a nominal specimen thickness i~rcspcctive of the test matctial. The notched disk specimen pVOl)oscd by Bcrn:ud ctal. [1] for static and dynamic initiation toughness measurement although p~ovicles plain-strain conditions, the crack plopagatcs at an angle to the direction of applied load. This makes inteq~retation of the test results difficult us it ~Ccluivcs ~actial slices to be cut fiom the fractured specimen to ascertain the angle o1 crack growth and a linite element model l~)r tl);t{ pa~ticulat ctack o~icntalion.

Fasteners in Composite Plates - Effects of Interfacial Friction

The effects of tangential friction at pin—hole interfaces are appropriately modelled for the analysis of fasteners in large composite (orthotropic) plate loaded along its edges. The pin—hole contact could be of interference, clearance or neat fit. When the plate load is monotonically increased, interference fits give rise to receding contact, whereas clearance fits result in advancing contact. In either case, the changing contact situations lead to non-linear moving boundary value problems. The neat fit comes out as a special case in which the contact and separation regions are invariant with the applied load level and so the problem remains linear. The description of boundary conditions in the presence of tangential friction, will depend on whether the problem is one of advancing or receding contact, advancing contact presenting a special problem. A model is developed for the limiting case of a rigid pin and an ideally rough interface (infinitely large friction coefficient). The non-linearity resulting from the continuously varying proportions of contact and separation at the interface, is handled by an “Inverse Formulation” which was successfully applied earlier by the authors for smooth (zero friction) interfacial conditions. The additional difficulty introduced by advancing contact is handled by adopting a “Marching Solution”. The modelling and the procedure are illustrated in respect of symmetric plate load cases. Numerical results are presented bringing out the effects of interfacial friction and plate orthotropy on load-contact relations and plate stresses.

Analysis of a finite composite plate with smooth rigid pin

A continuum method of analysis is presented in this paper for the problem of a smooth rigid pin in a finite composite plate subjected to uniaxial loading. The pin could be of interference, push or clearance fit. The plate is idealized to an orthotropic sheet. As the load on the plate is progressively increased, the contact along the pin-hole interface is partial above certain load levels in all three types of fit. In misfit pins (interference or clearance), such situations result in mixed boundary value problems with moving boundaries and in all of them the arc of contact and the stress and displacement fields vary nonlinearly with the applied load. In infinite domains similar problems were analysed earlier by ‘inverse formulation’ and, now, the same approach is selected for finite plates. Finite outer domains introduce analytical complexities in the satisfaction of boundary conditions. These problems are circumvented by adopting a method in which the successive integrals of boundary error functions are equated to zero. Numerical results are presented which bring out the effects of the rectangular geometry and the orthotropic property of the plate. The present solutions are the first step towards the development of special finite elements for fastener joints.

Achilles tendon loading patterns during barefoot walking and slow running on a treadmill: An ultrasonic propagation study

Measurement of tendon loading patterns during gait is important for understanding the pathogenesis of tendon "overuse" injury. Given that the speed of propagation of ultrasound in tendon is proportional to the applied load, this study used a noninvasive ultrasonic transmission technique to measure axial ultrasonic velocity in the right Achilles tendon of 27 healthy adults (11 females and 16 males; age, 26 ± 9 years; height, 1.73 ± 0.07 m; weight, 70.6 ± 21.2 kg), walking at self-selected speed (1.1 ± 0.1 m/s), and running at fixed slow speed (2 m/s) on a treadmill. Synchronous measures of ankle kinematics, spatiotemporal gait parameters, and vertical ground reaction forces were simultaneously measured. Slow running was associated with significantly higher cadence, shorter step length, but greater range of ankle movement, higher magnitude and rate of vertical ground reaction force, and higher ultrasonic velocity in the tendon than walking (P < 0.05). Ultrasonic velocity in the Achilles tendon was highly reproducible during walking and slow running (mean within-subject coefficient of variation < 2%). Ultrasonic maxima (P1, P2) and minima (M1, M2) were significantly higher and occurred earlier in the gait cycle (P1, M1, and M2) during running than walking (P < 0.05). Slow running was associated with higher and earlier peaks in loading of the Achilles tendon than walking.

Development of special fastener elements

A special finite element (FASNEL) is developed for the analysis of a neat or misfit fastener in a two-dimensional metallic/composite (orthotropic) plate subjected to biaxial loading. The misfit fasteners could be of interference or clearance type. These fasteners, which are common in engineering structures, cause stress concentrations and are potential sources of failure. Such cases of stress concentration present considerable numerical problems for analysis with conventional finite elements. In FASNEL the shape functions for displacements are derived from series stress function solutions satisfying the governing difffferential equation of the plate and some of the boundary conditions on the hole boundary. The region of the plate outside FASNEL is filled with CST or quadrilateral elements. When a plate with a fastener is gradually loaded the fastener-plate interface exhibits a state of partial contact/separation above a certain load level. In misfit fastener, the extent of contact/separation changes with applied load, leading to a nonlinear moving boundary problem and this is handled by FASNEL using an inverse formulation. The analysis is developed at present for a filled hole in a finite elastic plate providing two axes of symmetry. Numerical studies are conducted on a smooth rigid fastener in a finite elastic plate subjected to uniaxial loading to demonstrate the capability of FASNEL.