Analysis of transient eddy currents in MRI using a cylindrical FDTD method

Most magnetic resonance imaging (MRI) spatial encoding techniques employ low-frequency pulsed magnetic field gradients that undesirably induce multiexponentially decaying eddy currents in nearby conducting structures of the MRI system. The eddy currents degrade the switching performance of the gradient system, distort the MRI image, and introduce thermal loads in the cryostat vessel and superconducting MRI components. Heating of superconducting magnets due to induced eddy currents is particularly problematic as it offsets the superconducting operating point, which can cause a system quench. A numerical characterization of transient eddy current effects is vital for their compensation/control and further advancement of the MRI technology as a whole. However, transient eddy current calculations are particularly computationally intensive. In large-scale problems, such as gradient switching in MRI, conventional finite-element method (FEM)-based routines impose very large computational loads during generation/solving of the system equations. Therefore, other computational alternatives need to be explored. This paper outlines a three-dimensional finite-difference time-domain (FDTD) method in cylindrical coordinates for the modeling of low-frequency transient eddy currents in MRI, as an extension to the recently proposed time-harmonic scheme. The weakly coupled Maxwell's equations are adapted to the low-frequency regime by downscaling the speed of light constant, which permits the use of larger FDTD time steps while maintaining the validity of the Courant-Friedrich-Levy stability condition. The principal hypothesis of this work is that the modified FDTD routine can be employed to analyze pulsed-gradient-induced, transient eddy currents in superconducting MRI system models. The hypothesis is supported through a verification of the numerical scheme on a canonical problem and by analyzing undesired temporal eddy current effects such as the B-0-shift caused by actively shielded symmetric/asymmetric transverse x-gradient head and unshielded z-gradient whole-body coils operating in proximity to a superconducting MRI magnet.

Damage analysis of bridge structures using vibrational techniques

Much research is currently centred on the detection of damage in structures using vibrational data. The work presented here examined several areas of interest in support of a practical technique for identifying and locating damage within bridge structures using apparent changes in their vibrational response to known excitation. The proposed goals of such a technique included the need for the measurement system to be operated on site by a minimum number of staff and that the procedure should be as non-invasive to the bridge traffic-flow as possible. Initially the research investigated changes in the vibrational bending characteristics of two series of large-scale model bridge-beams in the laboratory and these included ordinary-reinforced and post-tensioned, prestressed designs. Each beam was progressively damaged at predetermined positions and its vibrational response to impact excitation was analysed. For the load-regime utilised the results suggested that the infuced damage manifested itself as a function of the span of a beam rather than a localised area. A power-law relating apparent damage with the applied loading and prestress levels was then proposed, together with a qualitative vibrational measure of structural damage. In parallel with the laboratory experiments a series of tests were undertaken at the sites of a number of highway bridges. The bridges selected had differing types of construction and geometric design including composite-concrete, concrete slab-and-beam, concrete-slab with supporting steel-troughing constructions together with regular-rectangular, skewed and heavily-skewed geometries. Initial investigations were made of the feasibility and reliability of various methods of structure excitation including traffic and impulse methods. It was found that localised impact using a sledge-hammer was ideal for the purposes of this work and that a cartridge `bolt-gun' could be used in some specific cases.

Microstructural analysis of surface and interface zones in concrete

Several of OPC paste and concrete specimens, with different mix proportions, were cast against CPF and impermeable formwork (IF) and the profiles of pore structure, microhardness and scratch hardness of the cover zone were established. The chloride ingress and the depth of carbonation of the surface zone of concrete cast against CPF and IF were investigated. The main mechanisms controlling the ECR processes and the factors affecting such treatment were critically reviewed. Subsequently, as a means of restoring passivation of steel embedded in carbonated concrete, such HCP specimens were subjected to ECR. The influence of ECR on the chemistry of the pore solution and the microstructure of the surface and the steel/cement past interface zones were also studied. The main findings of this investigation were as follows: (a) The thickness of the microstructure gradient of cover concrete is significantly decreased with increasing period of water curing but is relatively unaffected by curing temperature, w/e ratio and the use of cement replacement materials. (b) The scratch hardness technique was shown to be potentially useful for characterising the microstructure and microhardness gradients of the surface zone. (c) A relationship between the microstructure gradient and mass transport properties of the surface zone was established. (d) The use of CPF resulted in a significant reduction in porosity of both the cement paste matrix and the aggregate/cement paste transition zone, and a marked improvement in the resistance of the surface zone to carbonation and the ingress of chloride ions. (e) The ECR treatment resulted in a marked densification of the pore structure and in changes to the pore solution chemistry and the cement phases of near-surface and steel/cement paste transition zones. This effect was more pronounced with current density, period of treatment and particularly with the use of sodium phosphate as an electrolyte.

Effect of poor curing conditions and remedies on the durability of steel in concrete

An investigation was undertaken to study the effect of poor curing simulating hot climatic conditions and remedies on the durability of steel in concrete. Three different curing environments were used i.e. (1) Saturated Ca(OH)2 solution at 20°C, (2) Saturated Ca(OH)2 solution at 50°C and (3) Air at 50°C at 30% relative humidity. The third curing condition corresponding to the temperature and relative humidity typical of Middle Eastern Countries. The nature of the hardened cement paste matrix, cured under the above conditions was studied by means of Mercury Intrusion Porosimetry for measuring pore size distribution. The results were represented as total pore volume and initial pore entry diameter. The Scanning Electron Microscope was used to look at morphological changes during hydration, which were compared to the Mercury Intrusion Porosimetry results. X-ray defraction and Differential Thermal Analysis techniques were also employed for looking at any phase transformations. Polymer impregnation was used to reduce the porosity of the hardened cement pastes, especially in the case of the poorly cured samples. Carbonation rates of unimpregnated and impregnated cements were determined. Chloride diffusion studies were also undertaken to establish the effect of polymer impregnation and blending of the cements. Finally the corrosion behaviour of embedded steel bars was determined by the technique of Linear Polarisation. The steel was embedded in both untreated and polymer impregnated hardened cement pastes placed in either a solution containing NaCl or an environmental cabinet which provided carbonation at 40°C and 50% relative humidity.

Design of sway frames

The aim of the work presented in this thesis is to produce a direct method to design structures subject to deflection constraints at the working loads. The work carried out can be divided into four main parts. In the first part, a direct design procedure for plane steel frames subjected to sway limitations is proposed. The stiffness equations are modified so that the sway in each storey is equal to some specified values. The modified equations are then solved by iteration to calculate the cross-sectional properties of the columns as well as the other joint displacements. The beam sections are selected initially and then altered in an effort to reduce the total material cost of the frame. A linear extrapolation technique is used to reduce this cost. In this design, stability functions are used so that the effect of axial loads in the members are taken into consideration. The final reduced cost design is checked for strength requirements and the members are altered accordingly. In the second part, the design method is applied to the design of reinforced concrete frames in which the sway in the columns play an active part in the design criteria. The second moment of area of each column is obtained by solving the modified stiffness equations and then used to calculate the mlnlmum column depth required. Again the frame has to be checked for all the ultimate limit state load cases. In the third part, the method is generalised to design pin-jointed space frames for deflection limitatlions. In these the member areas are calculated so that the deflection at a specified joint is equal to its specified value. In the final part, the Lagrange multiplier technique is employed to obtain an optimum design for plane rigidly jointed steel frames. The iteration technique is used here to solve the modified stiffness equations as well as derivative equations obtained in accordance to the requirements of the optimisation method.

Kineto-elastodynamic analysis of high-speed four-bar mechanism

This thesis addresses the kineto-elastodynamic analysis of a four-bar mechanism running at high-speed where all links are assumed to be flexible. First, the mechanism, at static configurations, is considered as structure. Two methods are used to model the system, namely the finite element method (FEM) and the dynamic stiffness method. The natural frequencies and mode shapes at different positions from both methods are calculated and compared. The FEM is used to model the mechanism running at high-speed. The governing equations of motion are derived using Hamilton's principle. The equations obtained are a set of stiff ordinary differential equations with periodic coefficients. A model is developed whereby the FEM and the dynamic stiffness method are used conjointly to provide high-precision results with only one element per link. The principal concern of the mechanism designer is the behaviour of the mechanism at steady-state. Few algorithms have been developed to deliver the steady-state solution without resorting to costly time marching simulation. In this study two algorithms are developed to overcome the limitations of the existing algorithms. The superiority of the new algorithms is demonstrated. The notion of critical speeds is clarified and a distinction is drawn between "critical speeds", where stresses are at a local maximum, and "unstable bands" where the mechanism deflections will grow boundlessly. Floquet theory is used to assess the stability of the system. A simple method to locate the critical speeds is derived. It is shown that the critical speeds of the mechanism coincide with the local maxima of the eigenvalues of the transition matrix with respect to the rotational speed of the mechanism.

Design and analysis of a 2-DoF split-stator induction motor

A two degrees of freedom (2-DOF) actuator capable of producing linear translation, rotary motion, or helical motion would be a desirable asset to the fields of machine tools, robotics, and various apparatuses. In this paper, a novel 2-DOF split-stator induction motor was proposed and electromagnetic structure pa- rameters of the motor were designed and optimized. The feature of the direct-drive 2-DOF induction motor lies in its solid mover ar- rangement. In order to study the complex distribution of the eddy current field on the ferromagnetic cylinder mover and the motor’s operating characteristics, the mathematical model of the proposed motor was established, and characteristics of the motor were ana- lyzed by adopting the permeation depth method (PDM) and finite element method (FEM). The analytical and numerical results from motor simulation clearly show a correlation between the PDM and FEM models. This may be considered as a fair justification for the proposed machine and design tools.

Cross-country and inter-temporal analysis of Supply Chain Integration (SCI): China, Hungary and the Netherlands - changes in SCI from 1996 to 2008

Using plant level data from a global survey with multiple time frames, one begun in the late 1990s, this paper introduces measures of supply chain integration and discusses the dynamic relationship between the level of integration and a set of internal and external performance measurements. Specifically, data from Hungary, The Netherlands and The People’s Republic of China are used in the analyses. The time frames considered range from the late 1990s till 2009, encompassing major changes and transitions. Our results seem to indicate that SCI has an underlying structure of four sets of indicators, namely: (1) delivery frequency from the supplier or to the customer; (2) sharing internal processes with suppliers; (3) sharing internal processes with buyers and (4) joint facility location with partners. The differences between groups in terms of several performance measures proved to be small, being mostly statistically insignificant - but looking at the ANOVA table we can conclude that in this sample of companies those having joint location with their partners seem to outperform others.

A frame analysis of the leadership styles of dental and medical school deans

The purpose of this study was to determine how dental and medical school deans perceived their own leadership styles and how many management/leadership perspectives they used when approaching a problem. A descriptive and correlational study was conducted which examined the leadership orientations or "frames" utilized by these deans. Four frames--the structural, the human resource, the political, and the symbolic--emerged from leadership studies which described the behavior of leaders in organizational actions and problem solving.^ Employing the repeated measures ANCOVA, it was found that there were no significant interactions between type of dean and perceptions or perspectives. However, the main effects of both leadership style perceptions and the use of perspectives were significant. This indicated that over the total sample of deans, both perceptions and perspectives were used differently; the deans perceived themselves and used perspectives the highest for the human resource frame and the lowest for the political frame. Also, dental deans' mean use overall of the four perspective frames was higher than the medical deans. In addition to the ANCOVA, t-tests performed on perceptions revealed that the deans perceived the frames differently in three human resource items. Pearson's correlations indicated that for the combined sample of deans, significant negative correlations existed when the human resource or symbolic frame was compared to the structural frame. All of the deans used multiple frames, allowing them to re-frame or combine frames according to the situation. ^

Advanced thermal analysis of microelectronics using spreading resistance models

Thermal analysis of electronic devices is one of the most important steps for designing of modern devices. Precise thermal analysis is essential for designing an effective thermal management system of modern electronic devices such as batteries, LEDs, microelectronics, ICs, circuit boards, semiconductors and heat spreaders. For having a precise thermal analysis, the temperature profile and thermal spreading resistance of the device should be calculated by considering the geometry, property and boundary conditions. Thermal spreading resistance occurs when heat enters through a portion of a surface and flows by conduction. It is the primary source of thermal resistance when heat flows from a tiny heat source to a thin and wide heat spreader. In this thesis, analytical models for modeling the temperature behavior and thermal resistance in some common geometries of microelectronic devices such as heat channels and heat tubes are investigated. Different boundary conditions for the system are considered. Along the source plane, a combination of discretely specified heat flux, specified temperatures and adiabatic condition are studied. Along the walls of the system, adiabatic or convective cooling boundary conditions are assumed. Along the sink plane, convective cooling with constant or variable heat transfer coefficient are considered. Also, the effect of orthotropic properties is discussed. This thesis contains nine chapters. Chapter one is the introduction and shows the concepts of thermal spreading resistance besides the originality and importance of the work. Chapter two reviews the literatures on the thermal spreading resistance in the past fifty years with a focus on the recent advances. In chapters three and four, thermal resistance of a twodimensional flux channel with non-uniform convection coefficient in the heat sink plane is studied. The non-uniform convection is modeled by using two functions than can simulate a wide variety of different heat sink configurations. In chapter five, a non-symmetrical flux channel with different heat transfer coefficient along the right and left edges and sink plane is analytically modeled. Due to the edge cooling and non-symmetry, the eigenvalues of the system are defined using the heat transfer coefficient on both edges and for satisfying the orthogonality condition, a normalized function is calculated. In chapter six, thermal behavior of two-dimensional rectangular flux channel with arbitrary boundary conditions on the source plane is presented. The boundary condition along the source plane can be a combination of the first kind boundary condition (Dirichlet or prescribed temperature) and the second kind boundary condition (Neumann or prescribed heat flux). The proposed solution can be used for modeling the flux channels with numerous different source plane boundary conditions without any limitations in the number and position of heat sources. In chapter seven, temperature profile of a circular flux tube with discretely specified boundary conditions along the source plane is presented. Also, the effect of orthotropic properties are discussed. In chapter 8, a three-dimensional rectangular flux channel with a non-uniform heat convection along the heat sink plane is analytically modeled. In chapter nine, a summary of the achievements is presented and some systems are proposed for the future studies. It is worth mentioning that all the models and case studies in the thesis are compared with the Finite Element Method (FEM).

Seismic Response Analysis of a Full-Scale Base-Isolated Structure via Measurements and Modeling

The full-scale base-isolated structure studied in this dissertation is the only base-isolated building in South Island of New Zealand. It sustained hundreds of earthquake ground motions from September 2010 and well into 2012. Several large earthquake responses were recorded in December 2011 by NEES@UCLA and by GeoNet recording station nearby Christchurch Women's Hospital. The primary focus of this dissertation is to advance the state-of-the art of the methods to evaluate performance of seismic-isolated structures and the effects of soil-structure interaction by developing new data processing methodologies to overcome current limitations and by implementing advanced numerical modeling in OpenSees for direct analysis of soil-structure interaction.

This dissertation presents a novel method for recovering force-displacement relations within the isolators of building structures with unknown nonlinearities from sparse seismic-response measurements of floor accelerations. The method requires only direct matrix calculations (factorizations and multiplications); no iterative trial-and-error methods are required. The method requires a mass matrix, or at least an estimate of the floor masses. A stiffness matrix may be used, but is not necessary. Essentially, the method operates on a matrix of incomplete measurements of floor accelerations. In the special case of complete floor measurements of systems with linear dynamics, real modes, and equal floor masses, the principal components of this matrix are the modal responses. In the more general case of partial measurements and nonlinear dynamics, the method extracts a number of linearly-dependent components from Hankel matrices of measured horizontal response accelerations, assembles these components row-wise and extracts principal components from the singular value decomposition of this large matrix of linearly-dependent components. These principal components are then interpolated between floors in a way that minimizes the curvature energy of the interpolation. This interpolation step can make use of a reduced-order stiffness matrix, a backward difference matrix or a central difference matrix. The measured and interpolated floor acceleration components at all floors are then assembled and multiplied by a mass matrix. The recovered in-service force-displacement relations are then incorporated into the OpenSees soil structure interaction model.

Numerical simulations of soil-structure interaction involving non-uniform soil behavior are conducted following the development of the complete soil-structure interaction model of Christchurch Women's Hospital in OpenSees. In these 2D OpenSees models, the superstructure is modeled as two-dimensional frames in short span and long span respectively. The lead rubber bearings are modeled as elastomeric bearing (Bouc Wen) elements. The soil underlying the concrete raft foundation is modeled with linear elastic plane strain quadrilateral element. The non-uniformity of the soil profile is incorporated by extraction and interpolation of shear wave velocity profile from the Canterbury Geotechnical Database. The validity of the complete two-dimensional soil-structure interaction OpenSees model for the hospital is checked by comparing the results of peak floor responses and force-displacement relations within the isolation system achieved from OpenSees simulations to the recorded measurements. General explanations and implications, supported by displacement drifts, floor acceleration and displacement responses, force-displacement relations are described to address the effects of soil-structure interaction.

RELIBILITY ANALYSIS OF SIMPLE SLOPES AND SOIL-STRUCTURES WITH LINEAR LIMIT STATES

The first objective of this research was to develop closed-form and numerical probabilistic methods of analysis that can be applied to otherwise conventional methods of unreinforced and geosynthetic reinforced slopes and walls. These probabilistic methods explicitly include random variability of soil and reinforcement, spatial variability of the soil, and cross-correlation between soil input parameters on probability of failure. The quantitative impact of simultaneously considering the influence of random and/or spatial variability in soil properties in combination with cross-correlation in soil properties is investigated for the first time in the research literature. Depending on the magnitude of these statistical descriptors, margins of safety based on conventional notions of safety may be very different from margins of safety expressed in terms of probability of failure (or reliability index). The thesis work also shows that intuitive notions of margin of safety using conventional factor of safety and probability of failure can be brought into alignment when cross-correlation between soil properties is considered in a rigorous manner. The second objective of this thesis work was to develop a general closed-form solution to compute the true probability of failure (or reliability index) of a simple linear limit state function with one load term and one resistance term expressed first in general probabilistic terms and then migrated to a LRFD format for the purpose of LRFD calibration. The formulation considers contributions to probability of failure due to model type, uncertainty in bias values, bias dependencies, uncertainty in estimates of nominal values for correlated and uncorrelated load and resistance terms, and average margin of safety expressed as the operational factor of safety (OFS). Bias is defined as the ratio of measured to predicted value. Parametric analyses were carried out to show that ignoring possible correlations between random variables can lead to conservative (safe) values of resistance factor in some cases and in other cases to non-conservative (unsafe) values. Example LRFD calibrations were carried out using different load and resistance models for the pullout internal stability limit state of steel strip and geosynthetic reinforced soil walls together with matching bias data reported in the literature.

Simulation and analysis of an emergency lowering system for crane applications

An emergency lowering system for use in safety critical crane applications is discussed. The system is used to safely lower the payload of a crane in case of an electric blackout. The system is based on a backup power source, which is used to operate the crane while the regular supply is not available. The system enables both horizontal and vertical movements of the crane. Two different configurations for building the system are described, one with an uninterruptible power source (UPS) or a diesel generator connected in parallel to the crane’s power supply and one with a customized energy storage connected to the intermediate DC-link in the crane. In order to be able to size the backup power source, the power required during emergency lowering needs to be understood. A simulation model is used to study and optimize the power used during emergency lowering. The simulation model and optimizations are verified in a test hoist. Simulation results are presented with non-optimized and optimized controls for two example applications: a paper roll crane and a steel mill ladle crane. The optimizations are found to significantly reduce the required power for the crane movements during emergency lowering.

Finite element analysis of hot rolling in the blooming mill

During this thesis work a coupled thermo-mechanical finite element model (FEM) was builtto simulate hot rolling in the blooming mill at Sandvik Materials Technology (SMT) inSandviken. The blooming mill is the first in a long line of processes that continuously or ingotcast ingots are subjected to before becoming finished products. The aim of this thesis work was twofold. The first was to create a parameterized finiteelement (FE) model of the blooming mill. The commercial FE software package MSCMarc/Mentat was used to create this model and the programing language Python was used toparameterize it. Second, two different pass schedules (A and B) were studied and comparedusing the model. The two pass series were evaluated with focus on their ability to healcentreline porosity, i.e. to close voids in the centre of the ingot. This evaluation was made by studying the hydrostatic stress (σm), the von Mises stress (σeq)and the plastic strain (εp) in the centre of the ingot. From these parameters the stress triaxiality(Tx) and the hydrostatic integration parameter (Gm) were calculated for each pass in bothseries using two different transportation times (30 and 150 s) from the furnace. The relationbetween Gm and an analytical parameter (Δ) was also studied. This parameter is the ratiobetween the mean height of the ingot and the contact length between the rolls and the ingot,which is useful as a rule of thumb to determine the homogeneity or penetration of strain for aspecific pass. The pass series designed with fewer passes (B), many with greater reduction, was shown toachieve better void closure theoretically. It was also shown that a temperature gradient, whichis the result of a longer holding time between the furnace and the blooming mill leads toimproved void closure.