System identification of linear vibrating structures

Methods of dynamic modelling and analysis of structures, for example the finite element method, are well developed. However, it is generally agreed that accurate modelling of complex structures is difficult and for critical applications it is necessary to validate or update the theoretical models using data measured from actual structures. The techniques of identifying the parameters of linear dynamic models using Vibration test data have attracted considerable interest recently. However, no method has received a general acceptance due to a number of difficulties. These difficulties are mainly due to (i) Incomplete number of Vibration modes that can be excited and measured, (ii) Incomplete number of coordinates that can be measured, (iii) Inaccuracy in the experimental data (iv) Inaccuracy in the model structure. This thesis reports on a new approach to update the parameters of a finite element model as well as a lumped parameter model with a diagonal mass matrix. The structure and its theoretical model are equally perturbed by adding mass or stiffness and the incomplete number of eigen-data is measured. The parameters are then identified by an iterative updating of the initial estimates, by sensitivity analysis, using eigenvalues or both eigenvalues and eigenvectors of the structure before and after perturbation. It is shown that with a suitable choice of the perturbing coordinates exact parameters can be identified if the data and the model structure are exact. The theoretical basis of the technique is presented. To cope with measurement errors and possible inaccuracies in the model structure, a well known Bayesian approach is used to minimize the least squares difference between the updated and the initial parameters. The eigen-data of the structure with added mass or stiffness is also determined using the frequency response data of the unmodified structure by a structural modification technique. Thus, mass or stiffness do not have to be added physically. The mass-stiffness addition technique is demonstrated by simulation examples and Laboratory experiments on beams and an H-frame.

Deformation of sand in plane strain and oxisymmetric compression

Small scale laboratory experiments, in which the specimen is considered to represent an element of soil in the soil mass, are essential to the evolution of fundamental theories of mechanical behaviour. In this thesis, plane strain and axisymmetric compression tests, performed on a fine sand, are reported and the results are compared with various theoretical predictions. A new apparatus is described in which cuboidal samples can be tested in either axisymmetric compression or plane strain. The plane strain condition is simulated either by rigid side platens, in the conventional manner, or by flexible side platens which also measure the intermediate principal stress. Close control of the initial porosity of the specimens is achieved by a vibratory method of sample preparation. The strength of sand is higher in plane strain than in axisymmetric compression, and the strains required to mobilize peak strength are much smaller. The difference between plane strain and axisymmetric compression behaviour is attributed to the restrictions on particle movement enforced by the plane strain condition; this results in an increase in the frictional component of shear strength. The stress conditions at failure in plane strain, including the intermediate principal stress, are accurately predicted by a theory based on the stress- dilatancy interpretation of Mohr's circles. Detailed observations of rupture modes are presented and measured rupture plane inclinations are predicted by the stress-dilatancy theory. Although good correlation with the stress-dilatancy theory is obtained during virgin loading, in both axisymmetric compression and plane strain, the stress-dilatancy rule is only obeyed during reloading if the specimen has been unloaded to approximate ambient stress conditions. The shape of the stress-strain curves during pre-peak deformation, in both plane strain and axisymmetric compression, is accurately described bv a combined parabolic-hyperbolic specification.

The impact on agriculture of the drawdown of shallow watertables

The promoters of the large groundwater developments implemented in the 1970's paid little attention to the effects of pumping on soil moisture. A field study, conducted in 1979 in the Tern Area of the Shropshire Groundwater Scheme, revealed that significant quantities of the available moisture could be removed from the root zone of vegetation when drawdown of shallow watertables occurred. Arguments to this effect, supported by the field study evidence, were successfully presented at the Shropshire Groundwater Scheme public inquiry. The aim of this study has been to expand the work which was undertaken in connection with the Shropshire Groundwater Scheme, and to develop a method whereby the effects of groundwater pumping on vegetation can be assessed, and hence the impacts minimised. Two concepts, the critical height and the soil sensitivity depth, formulated during the initial work are at the core of the Environmental Impact Assessment method whose development is described. A programme of laboratory experiments on soil columns is described, as is the derivation of relationships for determining critical heights and field capacity moisture profiles. These relationships are subsequently employed in evaluating the effects of groundwater drawdown. In employing the environmental assessment technique, digitised maps of relevant features of the Tern Area are combined to produce composite maps delineating the extent of the areas which are potentially sensitive to groundwater drawdown. A series of crop yield/moisture loss functions are then employed to estimate the impact of simulated pumping events on the agricultural community of the Tern Area. Finally, guidelines, based on experience gained through evaluation of the Tern Area case study, are presented for use in the design of soil moisture monitoring systems and in the siting of boreholes. In addition recommendations are made for development of the EIA technique, and further research needs are identified.

Die load and stresses in press forging

Several axi-symmetric EN3B steel components differing in shape and size were forged on a 100 ton joint knuckle press. A load cell fitted under the lower die inserts recorded the total deformation forces. Job parameters were measured off the billets and the forged parts. Slug temperatures were varied and two lubricants - aqueous colloidal graphite and oil - were used. An industrial study was also conducted to check the results of the laboratory experiments. Loads were measured (with calibrated extensometers attached to the press frames) when adequately heated mild steel slugs were being forged in finishing dies. Geometric parameters relating to the jobs and the dies were obtained from works drawings. All the variables considered in the laboratory study could not, however, be investigated without disrupting production. In spite of this obvious limitation, the study confirmed that parting area is the most significant geometric factor influencing the forging load. Multiple regression analyses of the laboratory and industrial results showed that die loads increase significantly with the weights and parting areas of press forged components, and with the width to thickness ratios of the flashes formed, but diminish with increasing slug temperatures and higher billet diameter to height ratios. The analyses also showed that more complicated parts require greater loads to forge them. Die stresses, due to applied axial loads, were investigated by the photoelastic method. The three dimensional frozen stress technique was employed. Model dies were machined from cast araldite cylinders, and the slug material was simulated with plasticene. Test samples were cut from the centres of the dies after the stress freezing. Examination of the samples, and subsequent calculations, showed that the highest stresses were developed in die outer corners. This observation partly explains why corner cracking occurs frequently in industrial forging dies. Investigation of die contact during the forging operation revealed the development of very high stresses.

A dynamic model for the activated sludge treatment of coke oven effluents

This is an Inter-Disciplinary Higher Degree (IHD) thesis about Water Pollution Control in the Iron and Steel Industry. After examining the compositions, and various treatment methods, for the major effluent streams from a typical Integrated Iron and Steel works, it was decided to concentrate investigative work on the activated-sludge treatment of coke-oven effluents. A mathematical model of this process was developed in an attempt to provide a tool for plant management that would enable improved performance, and enhanced control of Works Units. The model differs from conventional models in that allowance is made for the presence of two genera of microorganisms, each of which utilises a particular type of substrate as its energy source. Allowance is also made for the inhibitive effect of phenol on thiocyanate biodegradation, and for the self-toxicity of the bacteria when present in a high substrate concentration environment. The enumeration of the kinetic characteristics of the two groups of micro-organisms was shown to be of major importance. Laboratory experiments were instigated in an attempt to determine accurate values of these coefficients. The use of the Suspended Solids concentration was found to be too insensitive a measure of viable active mass. Other measures were investigated, and Adenosine Triphosphate concentration was chosen as the most effective measure of bacterial populations. Using this measure, a model was developed for phenol biodegradation from experimental results which implicated the possibility of storage of substate prior to metabolism. A model for thiocyanate biodegradation was also developed, although the experimental results indicate that much work is still required in this area.

Development of an integrated reverse osmosis-greenhouse system driven by solar photovoltaic generators

The development of a system that integrates reverse osmosis (RO) with a horticultural greenhouse has been advanced through laboratory experiments. In this concept, intended for the inland desalination of brackish groundwater in dry areas, the RO concentrate will be reduced in volume by passing it through the evaporative cooling pads of the greenhouse. The system will be powered by solar photovoltaics (PV). Using a solar array simulator, we have verified that the RO can operate with varying power input and recovery rates to meet the water demands for irrigation and cooling of a greenhouse in north-west India. Cooling requires ventilation by a fan which has also been built, tested and optimised with a PV module outdoors. Results from the experiments with these two subsystems (RO and fan) are compared to theoretical predictions to reach conclusions about energy usage, sizing and cost. For example, the optimal sizing for the RO system is 0.12–1.3 m2 of PV module per m2 of membrane, depending on feed salinity. For the fan, the PV module area equals that of the fan aperture. The fan consumes <30 J of electrical energy per m3 of air moved which is 3 times less than that of standard fans. The specific energy consumption of the RO, at 1–2.3 kWh ?m-3, is comparable to that reported by others. Now that the subsystems have been verifi ed, the next step will be to integrate and test the whole system in the field.

Mean zonal flows generated by librations of a rotating spherical cavity

Longitudinal librations represent oscillations about the axis of a rotating axisymmetric fluid-filled cavity. An analytical theory is developed for the case of a spherical cavity in the limit when the libration frequency is small in comparison with the rotation rate, but large in comparison with the inverse of the spin-up time. It is shown that longitudinal librations create a steady zonal flow through the nonlinear advection in the Ekman layers. The theory can be applied to laboratory experiments as well as to solid planets and satellites with a liquid core. Copyright © Cambridge University Press 2010.

Performance analysis of 20 Gb/s RZ-DPSK non-slope matched transoceanic submarine links

Direct computation of the bit-error rate (BER) and laboratory experiments are used to assess the performance of a non-slope matched transoceanic submarine transmission link operating at 20Gb/s channel rate and employing return-to-zero differential-phase shift keying (RZ-DPSK) signal modulation. Using this system as an example, we compare the accuracies of the existing theoretical approaches to the BER estimation for the RZ-DPSK format.

Robustness of 40 Gb/s ASK modulation formats in the practical system infrastructure

In this work, we analyzed by means of numerical and laboratory experiments the resilience of 40 Gb/s amplitude shift keying modulation formats to transmission impairments in standard single-mode fiber lines as well as to optical filtering introduced by the optical add/drop multiplexer cascade. Our study is a pre-requisite to assess the implementation of cost-effective 40 Gb/s modulation technology in next generation high bit-rate robust optical transport networks.

Robustness of 40 Gb/s ASK modulation formats in the practical system infrastructure

In this work, we analyzed by means of numerical and laboratory experiments the resilience of 40 Gb/s amplitude shift keying modulation formats to transmission impairments in standard single-mode fiber lines as well as to optical filtering introduced by the optical add/drop multiplexer cascade. Our study is a pre-requisite to assess the implementation of cost-effective 40 Gb/s modulation technology in next generation high bit-rate robust optical transport networks. © 2006 Optical Society of America.

Performance analysis of 20Gb/s RZ-DPSK non-slope matched transoceanic submarine links

Direct computation of the bit-error rate (BER) and laboratory experiments are used to assess the performance of a non-slope matched transoceanic submarine transmission link operating at 20Gb/s channel rate and employing return-to-zero differential-phase shift keying (RZ-DPSK) signal modulation. Using this system as an example, we compare the accuracies of the existing theoretical approaches to the BER estimation for the RZ-DPSK format. © 2007 Optical Society of America.

Failure mechanism of the all-polyethylene glenoid implant

Fixation failure of glenoid components is the main cause of unsuccessful total shoulder arthroplasties. The characteristics of these failures are still not well understood, hence, attempts at improving the implant fixation are somewhat blind and the failure rate remains high. This lack of understanding is largely due to the fundamental problem that direct observations of failure are impossible as the fixation is inherently embedded within the bone. Twenty custom made implants, reflecting various common fixation designs, and a specimen set-up was prepared to enable direct observation of failure when the specimens were exposed to cyclic superior loads during laboratory experiments. Finite element analyses of the laboratory tests were also carried out to explain the observed failure scenarios. All implants, irrespective of the particular fixation design, failed at the implant-cement interface and failure initiated at the inferior part of the component fixation. Finite element analyses indicated that this failure scenario was caused by a weak and brittle implant-cement interface and tensile stresses in the inferior region possibly worsened by a stress raiser effect at the inferior rim. The results of this study indicate that glenoid failure can be delayed or prevented by improving the implant/cement interface strength. Also any design features that reduce the geometrical stress raiser and the inferior tensile stresses in general should delay implant loosening.

Oribatid mites (Acari : Oribatida) in microcosms-A review

Oribatid mites are one of the most abundant and species rich group in soil mesofauna. In spite of this, we have limited knowledge on the role they play in ecological processes. Since they are small like other mesofauna members and live a hidden life, their observation is difficult. Their life style, interactions with other organisms and role in soil decomposition processes can be investigated with laboratory experiments more exactly and effectively. While the literature of microcosm studies is very extensive, relatively few review papers have been written about the methodology of microcosm studies. The present review summarizes details of techniques that have been applied laboratory microcosms involving oribatid mites. It is shown what one should keep in mind in planning, composing and setting up a microcosm. Our comparative evaluation reveals how the laboratory experiments were maintained and manipulated and what kind of information was extracted. This methodological review can be useful in preparing microcosm experiments applied to other animal groups as well.

The spectral distribution of biologically active solar radiation at Miami, Florida, USA

The spectral distribution of solar radiation was studied under different sky conditions during a 15- month period in Miami, Florida (USA), and over a latitudinal gradient at solar maximum. Spectroradiometric scans were characterized for total irradiance (300- 3000 nm) and the relative energetic and photon contributions of the following wavelength regions: UV-B (300-320nm); UV-A (320-400nm); B (400-500rim); PAR (400-700 nm); R (600-700 nm); and FR (728- 732 rim). Notable results include: (i) significantly higher UV-A energy fluxes than currently in use for laboratory experiments involving the biological effects of this bandwidth (values ranged from 33.6 to 55.4 W/m 2 in Miami over the year); (ii) marked diurnal shifts in B:R and R:FR, with elevated R:FR values in early morning: (iii) a strong correlation between R: FR and atmospheric water content; and (iv) unusually high PAR values under direct sunlight with cloudy skies (2484 ~tmot/2 per s).

Eulerian-Lagrangian two phase debris flow model

The main objective of this work is to develop a quasi three-dimensional numerical model to simulate stony debris flows, considering a continuum fluid phase, composed by water and fine sediments, and a non-continuum phase including large particles, such as pebbles and boulders. Large particles are treated in a Lagrangian frame of reference using the Discrete Element Method, the fluid phase is based on the Eulerian approach, using the Finite Element Method to solve the depth-averaged Navier-Stokes equations in two horizontal dimensions. The particle’s equations of motion are in three dimensions. The model simulates particle-particle collisions and wall-particle collisions, taking into account that particles are immersed in a fluid. Bingham and Cross rheological models are used for the continuum phase. Both formulations provide very stable results, even in the range of very low shear rates. Bingham formulation is better able to simulate the stopping stage of the fluid when applied shear stresses are low. Results of numerical simulations have been compared with data from laboratory experiments on a flume-fan prototype. Results show that the model is capable of simulating the motion of big particles moving in the fluid flow, handling dense particulate flows and avoiding overlap among particles. An application to simulate debris flow events that occurred in Northern Venezuela in 1999 shows that the model could replicate the main boulder accumulation areas that were surveyed by the USGS. Uniqueness of this research is the integration of mud flow and stony debris movement in a single modeling tool that can be used for planning and management of debris flow prone areas.