Rotor retaining sleeve design for a 1.12-MW high-speed PM machine

Permanent-magnet (PM) synchronous machines (PMSMs) can provide excellent performance in terms of torque density, energy efficiency, and controllability. However, PMs on the rotor are prone to centrifugal force, which may break their physical integrity, particularly at high-speed operation. Typically, PMs are bound with carbon fiber or retained by alloy sleeves on the rotor surface. This paper is concerned with the design of a rotor retaining sleeve for a 1.12-MW 18-kr/min PM machine; its electromagnetic performance is investigated by the 2-D finite-element method (FEM). Theoretical and numerical analyses of the rotor stress are carried out. For the carbon fiber protective measure, the stresses of three PM configurations and three pole filler materials are compared in terms of operating temperature, rotor speed, retaining sleeve thickness, and interference fit. Then, a new hybrid protective measure is proposed and analyzed by the 2-D FEM for operational speeds up to 22 kr/min (1.2 times the rated speed). The rotor losses and machine temperatures with the carbon fiber retaining sleeve and the hybrid retaining sleeve are compared, and the sleeve design is refined. Two rotors using both designs are prototyped and experimentally tested to validate the effectiveness of the developed techniques for PM machines. The developed retaining sleeve makes it possible to operate megawatt PM machines at high speeds of 22 kr/min. This opens doors for many high-power high-speed applications such as turbo-generator, aerospace, and submarine motor drives.

Design and flux-weakening control of an interior permanent magnet synchronous motor for electric vehicles

Permanent magnet synchronous motors (PMSMs) provide a competitive technology for EV traction drives owing to their high power density and high efficiency. In this paper, three types of interior PMSMs with different PM arrangements are modeled by the finite element method (FEM). For a given amount of permanent magnet materials, the V shape interior PMSM is found better than the U-shape and the conventional rotor topologies for EV traction drives. Then the V shape interior PMSM is further analyzed with the effects of stator slot opening and the permanent magnet pole chamfering on cogging torque and output torque performance. A vector-controlled flux-weakening method is developed and simulated in matlab to expand the motor speed range for EV drive system. The results show good dynamic and steady-state performance with a capability of expanding speed up to 4 times of the rated. A prototype of the V shape interior PMSM is also manufactured and tested to validate the numerical models built by the finite element method.

Lateral load distribution for curved steel I -girder bridges

Lateral load distribution factor is a key factor for designing and analyzing curved steel I-girder bridges. In this dissertation, the effects of various parameters on moment and shear distribution for curved steel I-girder bridges were studied using the Finite Element Method (FEM). The parameters considered in the study were: radius of curvature, girder spacing, overhang, span length, number of girders, ratio of girder stiffness to overall bridge stiffness, slab thickness, girder longitudinal stiffness, cross frame spacing, and girder torsional inertia. The variations of these parameters were based on the statistical analysis of the real bridge database, which was created by extracting data from existing or newly designed curved steel I-girder bridge plans collected all over the nation. A hypothetical bridge superstructure model that was made of all the mean values of the data was created and used for the parameter study. ^ The study showed that cross frame spacing and girder torsional inertia had negligible effects. Other parameters had been identified as key parameters. Regression analysis was conducted based on the FEM analysis results and simplified formulas for predicting positive moment, negative moment, and shear distribution factors were developed. Thirty-three real bridges were analyzed using FEM to verify the formulas. The ratio of the distribution factor obtained from the formula to the one obtained from the FEM analysis, which was referred to as the g-ratio, was examined. The results showed that the standard deviation of the g-ratios was within 0.04 to 0.06 and the mean value of the g-ratios was greater than unity by one standard deviation. This indicates that the formulas are conservative in most cases but not overly conservative. The final formulas are similar in format to the current American Association of State Highway and Transportation Officials (AASHTO) Load Resistance and Factor Design (LRFD) specifications. ^ The developed formulas were compared with other simplified methods. The outcomes showed that the proposed formulas had the most accurate results among all methods. ^ The formulas developed in this study will assist bridge engineers and researchers in predicting the actual live load distribution in horizontally curved steel I-girder bridges. ^

Durability Assessment of Polymer Trileaflet Heart Valves

The durability of a polymer trileaflet valve is dependent on leaflet stress concentrations, so valve designs that reduce stress can, hypothetically, increase durability. Design aspects that are believed to contribute to reduced leaflet stress include stent flexibility, parabolic coaptation curvature, and leaflet anisotropy. With this in mind, the purpose of this investigation was to elucidate what specific combinations of these parameters promote optimal acute and long-term valve function. A combination of four stent designs, seven leaflet reinforcement materials, and three coaptation geometries were evaluated through a combination of experimentation and modeling. Static tensile and Poisson’s ratio tests and dynamic tensile fatigue testing were used to evaluate the individual leaflet components; and hydrodynamic testing and accelerated valve fatigue was used to assess complete valve prototypes. The two most successful designs included a 0.40 mm thick knit-reinforced valve with a fatigue life of 10.35 years, and a 0.20 mm thick knit-reinforced valve with a 28.9 mmHg decrease in pressure drop over the former. A finite element model was incorporated to verify the impact of the above-mentioned parameters on leaflet stress concentrations. Leaflet anisotropy had a large impact on stress concentrations, and matching the circumferential modulus to that of the natural valve showed the greatest benefit. Varying the radial modulus had minimal impact. Varying coaptation geometry had no impact, but stent flexibility did have a marked effect on the stress at the top of the commissure, where a completely rigid stent resulted in a higher peak stress than a flexible stent (E = 385 MPa). In conclusion, stent flexibility and leaflet anisotropy do effect stress concentrations in the SIBS trileaflet valve, but coaptation geometry does not. Regions of high stress concentrations were linked to failure locations in vitro, so a fatigue prediction model was developed from the S/N curves generated during dynamic tensile testing of the 0.20 mm knit-reinforced leaflets. Failure was predicted at approximately 400 million cycles (10 years) at the top of the commissure. In vitro fatigue of this valve showed failure initiation after approximately 167 million cycles (4.18 years), but it was related to a design defect that is subsequently being changed.

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.

Modelagem numérica de uma estrutura de contenção de estacas espaçadas atirantadas em areia

Retaining walls design involves factors such as plastification, loading and unloading, pre-stressing, excessive displacements and earth and water thrust. Furthermore, the interaction between the retained soil and the structure is rather complex and hard to predict. Despite the advances in numerical simulation techniques and monitoring of forces and displacements with field instrumentation, design projects are still based on classical methods, whose simplifying assumptions may overestimate structural elements of the retaining wall. This dissertation involves a three-dimensional numerical study on the behavior of a retaining wall using the finite element method (FEM). The retaining wall structure is a contiguous bored pile wall with tie-back anchors. The numerical results were compared to data obtained from field instrumentation. The influence of the position of one or two layers of anchors and the effects of the construction of a slab bounded at the top of the retaining wall was evaluated. Furthermore, this study aimed at investigating the phenomenon of arching in the soil behind the wall. Arching was evaluated by analyzing the effects of pile spacing on horizontal stresses and displacements. Parametric analysis with one layers of anchors showed that the smallest horizontal displacements of the structure were achieved for between 0.3 and 0.5 times the excavation depth. Parametric analyses with two anchor layers showed that the smallest horizontal displacements were achieve for anchors positioned in depths of 0.4H and 0.7H. The construction of a slab at the top of the retaining wall decreased the horizontal displacements by 0.14% times the excavation depth as compared to analyses without the slab. With regard to the arching , analyzes showed an optimal range of spacing between the faces of the piles between 0.4 and 0.6 times the diameter of the pile

Contribuição ao estudo de treliças não lineares considerando efeitos térmicos

This work presents the numerical analysis of nonlinear trusses summited to thermomechanical actions with Finite Element Method (FEM). The proposed formulation is so-called positional FEM and it is based on the minimum potential energy theorem written according to nodal positions, instead of displacements. The study herein presented considers the effects of geometric and material nonlinearities. Related to dynamic problems, a comparison between different time integration algorithms is performed. The formulation is extended to impact problems between trusses and rigid wall, where the nodal positions are constrained considering nullpenetration condition. In addition, it is presented a thermodynamically consistent formulation, based on the first and second law of thermodynamics and the Helmholtz free-energy for analyzing dynamic problems of truss structures with thermoelastic and thermoplastic behavior. The numerical results of the proposed formulation are compared with examples found in the literature.

Novel double-sided linear generator for the wave energy conversion

The direct drive point absorber is a robust and efficient system for wave energy harvesting, where the linear generator represents the most complex part of the system. Therefore, its design and optimization are crucial tasks. The tubular shape of a linear generator’s magnetic circuit offers better permanent magnet flux encapsulation and reduction in radial forces on the translator due to its symmetry. A double stator topology can improve the power density of the linear tubular machine. Common designs employ a set of aligned stators on each side of a translator with radially magnetized permanent magnets. Such designs require doubling the amount of permanent magnet material and lead to an increase in the cogging force. The design presented in this thesis utilizes a translator with buried axially magnetized magnets and axially shifted positioning of the two stators such that no additional magnetic material, compared to single side machine, is required. In addition to the conservation of magnetic material, a significant improvement in the cogging force occurs in the two phase topology, while the double sided three phase system produces more power at the cost of a small increase in the cogging force. The analytical and the FEM models of the generator are described and their results compared to the experimental results. In general, the experimental results compare favourably with theoretical predictions. However, the experimentally observed permanent magnet flux leakage in the double sided machine is larger than predicted theoretically, which can be justified by the limitations in the prototype fabrication and resulting deviations from the theoretical analysis.

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).

Semi-engineered earthquake-resistant structures: one-storey buildings built up with gabion-box walls

This thesis studies the static and seismic behavior of simple structures made with gabion box walls. The analysis was performed considering a one-story building with standard dimensions in plan (6m x 5m) and a lightweight timber roof. The main focus of the present investigation is to find the principals aspects of the seismic behavior of a one story building made with gabion box walls, in order to prevent a failure due to seismic actions and in this way help to reduce the seismic risk of developing countries where this natural disaster have a significant intensity. Regarding the gabion box wall, it has been performed some calculations and analysis in order to understand the static and dynamic behavior. From the static point of view, it has been performed a verification of the normal stress computing the normal stress that arrives at the base of the gabion wall and the corresponding capacity of the ground. Moreover, regarding the seismic analysis, it has been studied the in-plane and out-of-plane behavior. The most critical aspect was discovered to be the out-of-plane behavior, for which have been developed models considering the “rigid- no tension model” for masonry, finding a kinematically admissible multiplier that will create a collapse mechanism for the structure. Furthermore, it has been performed a FEM and DEM models to find the maximum displacement at the center of the wall, maximum tension stresses needed for calculating the steel connectors for joining consecutive gabions and the dimensions (length of the wall and distance between orthogonal walls or buttresses) of a geometrical configuration for the standard modulus of the structure, in order to ensure an adequate safety margin for earthquakes with a PGA around 0.4-0.5g. Using the results obtained before, it has been created some rules of thumb, that have to be satisfy in order to ensure a good behavior of these structure.

Biomechanical modeling of parotid glands morphing in head & neck radiation therapy treatments

Durante i trattamenti radioterapici dei pazienti oncologici testa-collo, le ghiandole parotidee (PGs) possono essere indebitamente irradiate a seguito di modificazioni volumetriche-spaziali inter/intra-frazione causate da fattori quali il dimagrimento, l’esposizione a radiazioni ionizzanti ed il morphing anatomico degli organi coinvolti nelle aree d’irraggiamento. Il presente lavoro svolto presso la struttura di Fisica Medica e di Radioterapia Oncologica dell’A.O.U di Modena, quale parte del progetto di ricerca del Ministero della Salute (MoH2010, GR-2010-2318757) “ Dose warping methods for IGRT and Adaptive RT: dose accumulation based on organ motion and anatomical variations of the patients during radiation therapy treatments ”, sviluppa un modello biomeccanico in grado di rappresentare il processo di deformazione delle PGs, considerandone la geometria, le proprietà elastiche e l'evoluzione durante il ciclo terapeutico. Il modello di deformazione d’organo è stato realizzato attraverso l’utilizzo di un software agli elementi finiti (FEM). Molteplici superfici mesh, rappresentanti la geometria e l’evoluzione delle parotidi durante le sedute di trattamento, sono state create a partire dai contorni dell’organo definiti dal medico radioterapista sull’immagine tomografica di pianificazione e generati automaticamente sulle immagini di setup e re-positioning giornaliere mediante algoritmi di registrazione rigida/deformabile. I constraints anatomici e il campo di forze del modello sono stati definiti sulla base di ipotesi semplificative considerando l’alterazione strutturale (perdita di cellule acinari) e le barriere anatomiche dovute a strutture circostanti. L’analisi delle mesh ha consentito di studiare la dinamica della deformazione e di individuare le regioni maggiormente soggette a cambiamento. Le previsioni di morphing prodotte dal modello proposto potrebbero essere integrate in un treatment planning system per metodiche di Adaptive Radiation Therapy.

Piezoelectric energy harvesting using blood-flow-like excitation for implantable devices

The goal of this research is to produce a system for powering medical implants to increase the lifetime of the implanted devices and reduce the battery size. The system consists of a number of elements – the piezoelectric material for generating power, the device design, the circuit for rectification and energy storage. The piezoelectric material is analysed and a process for producing a repeatable high quality piezoelectric material is described. A full width half maximum (FWHM) of the rocking curve X-Ray diffraction (XRD) scan of between ~1.5° to ~1.7° for test wafers was achieved. This is state of the art for AlN on silicon and means devices with good piezoelectric constants can be fabricated. Finite element modelling FEM) was used to design the structures for energy harvesting. The models developed in this work were established to have an accuracy better than 5% in terms of the difference between measured and modelled results. Devices made from this material were analysed for power harvesting ability as well as the effect that they have on the flow of liquid which is an important consideration for implantable devices. The FEM results are compared to experimental results from laser Doppler vibrometry (LDV), magnetic shaker and perfusion machine tests. The rectifying circuitry for the energy harvester was also investigated. The final solution uses multiple devices to provide the power to augment the battery and so this was a key feature to be considered. Many circuits were examined and a solution based on a fully autonomous circuit was advanced. This circuit was analysed for use with multiple low power inputs similar to the results from previous investigations into the energy harvesting devices. Polymer materials were also studied for use as a substitute for the piezoelectric material as well as the substrate because silicon is more brittle.

Modeling of a Heat-Induced Buckling of Plates Using the Mesh-free Method

In the process of engineering design of structural shapes, the flat plate analysis results can be generalized to predict behaviors of complete structural shapes. In this case, the purpose of this project is to analyze a thin flat plate under conductive heat transfer and to simulate the temperature distribution, thermal stresses, total displacements, and buckling deformations. The current approach in these cases has been using the Finite Element Method (FEM), whose basis is the construction of a conforming mesh. In contrast, this project uses the mesh-free Scan Solve Method. This method eliminates the meshing limitation using a non-conforming mesh. I implemented this modeling process developing numerical algorithms and software tools to model thermally induced buckling. In addition, convergence analysis was achieved, and the results were compared with FEM. In conclusion, the results demonstrate that the method gives similar solutions to FEM in quality, but it is computationally less time consuming.

He isotope, U/Th isotope, Calcium carbonate, biogenic opal, rare earth element concentrations, and grain size distribution measured on core-top sediments during SONNE cruise SO202 (INOPEX)

Eolian dust is a significant source of iron and other nutrients that are essential for the health of marine ecosystems and potentially a controlling factor of the high nutrient-low chlorophyll status of the Subarctic North Pacific. We map the spatial distribution of dust input using three different geochemical tracers of eolian dust, 4He, 232Th and rare earth elements, in combination with grain size distribution data, from a set of core-top sediments covering the entire Subarctic North Pacific. Using the suite of geochemical proxies to fingerprint different lithogenic components, we deconvolve eolian dust input from other lithogenic inputs such as volcanic ash, ice-rafted debris, riverine and hemipelagic input. While the open ocean sites far away from the volcanic arcs are dominantly composed of pure eolian dust, lithogenic components other than eolian dust play a more crucial role along the arcs. In sites dominated by dust, eolian dust input appears to be characterized by a nearly uniform grain size mode at ~4 µm. Applying the 230Th-normalization technique, our proxies yield a consistent pattern of uniform dust fluxes of 1-2 g/m**2/yr across the Subarctic North Pacific. Elevated eolian dust fluxes of 2-4 g/m**2/yr characterize the westernmost region off Japan and the southern Kurile Islands south of 45° N and west of 165° E along the main pathway of the westerly winds. The core-top based dust flux reconstruction is consistent with recent estimates based on dissolved thorium isotope concentrations in seawater from the Subarctic North Pacific. The dust flux pattern compares well with state-of-the-art dust model predictions in the western and central Subarctic North Pacific, but we find that dust fluxes are higher than modeled fluxes by 0.5-1 g/m**2/yr in the northwest, northeast and eastern Subarctic North Pacific. Our results provide an important benchmark for biogeochemical models and a robust approach for downcore studies testing dust-induced iron fertilization of past changes in biological productivity in the Subarctic North Pacific.