Weak form equation–based finite-element modeling of viscoelastic asphalt mixtures

The objective of this study is to demonstrate using weak form partial differential equation (PDE) method for a finite-element (FE) modeling of a new constitutive relation without the need of user subroutine programming. The viscoelastic asphalt mixtures were modeled by the weak form PDE-based FE method as the examples in the paper. A solid-like generalized Maxwell model was used to represent the deforming mechanism of a viscoelastic material, the constitutive relations of which were derived and implemented in the weak form PDE module of Comsol Multiphysics, a commercial FE program. The weak form PDE modeling of viscoelasticity was verified by comparing Comsol and Abaqus simulations, which employed the same loading configurations and material property inputs in virtual laboratory test simulations. Both produced identical results in terms of axial and radial strain responses. The weak form PDE modeling of viscoelasticity was further validated by comparing the weak form PDE predictions with real laboratory test results of six types of asphalt mixtures with two air void contents and three aging periods. The viscoelastic material properties such as the coefficients of a Prony series model for the relaxation modulus were obtained by converting from the master curves of dynamic modulus and phase angle. Strain responses of compressive creep tests at three temperatures and cyclic load tests were predicted using the weak form PDE modeling and found to be comparable with the measurements of the real laboratory tests. It was demonstrated that the weak form PDE-based FE modeling can serve as an efficient method to implement new constitutive models and can free engineers from user subroutine programming.

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.

Multiobjective design optimization of IGBT power modules considering power cycling and thermal cycling

Insulated-gate bipolar transistor (IGBT) power modules find widespread use in numerous power conversion applications where their reliability is of significant concern. Standard IGBT modules are fabricated for general-purpose applications while little has been designed for bespoke applications. However, conventional design of IGBTs can be improved by the multiobjective optimization technique. This paper proposes a novel design method to consider die-attachment solder failures induced by short power cycling and baseplate solder fatigue induced by the thermal cycling which are among major failure mechanisms of IGBTs. Thermal resistance is calculated analytically and the plastic work design is obtained with a high-fidelity finite-element model, which has been validated experimentally. The objective of minimizing the plastic work and constrain functions is formulated by the surrogate model. The nondominated sorting genetic algorithm-II is used to search for the Pareto-optimal solutions and the best design. The result of this combination generates an effective approach to optimize the physical structure of power electronic modules, taking account of historical environmental and operational conditions in the field.

State detection of bond wires in IGBT modules using eddy current pulsed thermography

Insulated gate bipolar transistor (IGBT) modules are important safety critical components in electrical power systems. Bond wire lift-off, a plastic deformation between wire bond and adjacent layers of a device caused by repeated power/thermal cycles, is the most common failure mechanism in IGBT modules. For the early detection and characterization of such failures, it is important to constantly detect or monitor the health state of IGBT modules, and the state of bond wires in particular. This paper introduces eddy current pulsed thermography (ECPT), a nondestructive evaluation technique, for the state detection and characterization of bond wire lift-off in IGBT modules. After the introduction of the experimental ECPT system, numerical simulation work is reported. The presented simulations are based on the 3-D electromagnetic-thermal coupling finite-element method and analyze transient temperature distribution within the bond wires. This paper illustrates the thermal patterns of bond wires using inductive heating with different wire statuses (lifted-off or well bonded) under two excitation conditions: nonuniform and uniform magnetic field excitations. Experimental results show that uniform excitation of healthy bonding wires, using a Helmholtz coil, provides the same eddy currents on each, while different eddy currents are seen on faulty wires. Both experimental and numerical results show that ECPT can be used for the detection and characterization of bond wires in power semiconductors through the analysis of the transient heating patterns of the wires. The main impact of this paper is that it is the first time electromagnetic induction thermography, so-called ECPT, has been employed on power/electronic devices. Because of its capability of contactless inspection of multiple wires in a single pass, and as such it opens a wide field of investigation in power/electronic devices for failure detection, performance characterization, and health monitoring.

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.

Lattice Boltzmann method for flow and heat transfer in microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Design optimization of modern machine drive systems for maximum fault tolerant and optimal operation

Modern electric machine drives, particularly three phase permanent magnet machine drive systems represent an indispensable part of high power density products. Such products include; hybrid electric vehicles, large propulsion systems, and automation products. Reliability and cost of these products are directly related to the reliability and cost of these systems. The compatibility of the electric machine and its drive system for optimal cost and operation has been a large challenge in industrial applications. The main objective of this dissertation is to find a design and control scheme for the best compromise between the reliability and optimality of the electric machine-drive system. The effort presented here is motivated by the need to find new techniques to connect the design and control of electric machines and drive systems. ^ A highly accurate and computationally efficient modeling process was developed to monitor the magnetic, thermal, and electrical aspects of the electric machine in its operational environments. The modeling process was also utilized in the design process in form finite element based optimization process. It was also used in hardware in the loop finite element based optimization process. The modeling process was later employed in the design of a very accurate and highly efficient physics-based customized observers that are required for the fault diagnosis as well the sensorless rotor position estimation. Two test setups with different ratings and topologies were numerically and experimentally tested to verify the effectiveness of the proposed techniques. ^ The modeling process was also employed in the real-time demagnetization control of the machine. Various real-time scenarios were successfully verified. It was shown that this process gives the potential to optimally redefine the assumptions in sizing the permanent magnets of the machine and DC bus voltage of the drive for the worst operating conditions. ^ The mathematical development and stability criteria of the physics-based modeling of the machine, design optimization, and the physics-based fault diagnosis and the physics-based sensorless technique are described in detail. ^ To investigate the performance of the developed design test-bed, software and hardware setups were constructed first. Several topologies of the permanent magnet machine were optimized inside the optimization test-bed. To investigate the performance of the developed sensorless control, a test-bed including a 0.25 (kW) surface mounted permanent magnet synchronous machine example was created. The verification of the proposed technique in a range from medium to very low speed, effectively show the intelligent design capability of the proposed system. Additionally, to investigate the performance of the developed fault diagnosis system, a test-bed including a 0.8 (kW) surface mounted permanent magnet synchronous machine example with trapezoidal back electromotive force was created. The results verify the use of the proposed technique under dynamic eccentricity, DC bus voltage variations, and harmonic loading condition make the system an ideal case for propulsion systems.^

Physics-based modeling of power system components for the evaluation of low-frequency radiated electromagnetic fields

The low-frequency electromagnetic compatibility (EMC) is an increasingly important aspect in the design of practical systems to ensure the functional safety and reliability of complex products. The opportunities for using numerical techniques to predict and analyze system's EMC are therefore of considerable interest in many industries. As the first phase of study, a proper model, including all the details of the component, was required. Therefore, the advances in EMC modeling were studied with classifying analytical and numerical models. The selected model was finite element (FE) modeling, coupled with the distributed network method, to generate the model of the converter's components and obtain the frequency behavioral model of the converter. The method has the ability to reveal the behavior of parasitic elements and higher resonances, which have critical impacts in studying EMI problems. For the EMC and signature studies of the machine drives, the equivalent source modeling was studied. Considering the details of the multi-machine environment, including actual models, some innovation in equivalent source modeling was performed to decrease the simulation time dramatically. Several models were designed in this study and the voltage current cube model and wire model have the best result. The GA-based PSO method is used as the optimization process. Superposition and suppression of the fields in coupling the components were also studied and verified. The simulation time of the equivalent model is 80-100 times lower than the detailed model. All tests were verified experimentally. As the application of EMC and signature study, the fault diagnosis and condition monitoring of an induction motor drive was developed using radiated fields. In addition to experimental tests, the 3DFE analysis was coupled with circuit-based software to implement the incipient fault cases. The identification was implemented using ANN for seventy various faulty cases. The simulation results were verified experimentally. Finally, the identification of the types of power components were implemented. The results show that it is possible to identify the type of components, as well as the faulty components, by comparing the amplitudes of their stray field harmonics. The identification using the stray fields is nondestructive and can be used for the setups that cannot go offline and be dismantled

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.

Lattice Boltzmann Method for Flow and Heat Transfer in Microgeometries

Recent technological developments have made it possible to design various microdevices where fluid flow and heat transfer are involved. For the proper design of such systems, the governing physics needs to be investigated. Due to the difficulty to study complex geometries in micro scales using experimental techniques, computational tools are developed to analyze and simulate flow and heat transfer in microgeometries. However, conventional numerical methods using the Navier-Stokes equations fail to predict some aspects of microflows such as nonlinear pressure distribution, increase mass flow rate, slip flow and temperature jump at the solid boundaries. This necessitates the development of new computational methods which depend on the kinetic theory that are both accurate and computationally efficient. In this study, lattice Boltzmann method (LBM) was used to investigate the flow and heat transfer in micro sized geometries. The LBM depends on the Boltzmann equation which is valid in the whole rarefaction regime that can be observed in micro flows. Results were obtained for isothermal channel flows at Knudsen numbers higher than 0.01 at different pressure ratios. LBM solutions for micro-Couette and micro-Poiseuille flow were found to be in good agreement with the analytical solutions valid in the slip flow regime (0.01 < Kn < 0.1) and direct simulation Monte Carlo solutions that are valid in the transition regime (0.1 < Kn < 10) for pressure distribution and velocity field. The isothermal LBM was further extended to simulate flows including heat transfer. The method was first validated for continuum channel flows with and without constrictions by comparing the thermal LBM results against accurate solutions obtained from analytical equations and finite element method. Finally, the capability of thermal LBM was improved by adding the effect of rarefaction and the method was used to analyze the behavior of gas flow in microchannels. The major finding of this research is that, the newly developed particle-based method described here can be used as an alternative numerical tool in order to study non-continuum effects observed in micro-electro-mechanical-systems (MEMS).

Implementation of X-ray diffraction in the measurement of residual thermal strains

Microelectronic systems are multi-material, multi-layer structures, fabricated and exposed to environmental stresses over a wide range of temperatures. Thermal and residual stresses created by thermal mismatches in films and interconnections are a major cause of failure in microelectronic devices. Due to new device materials, increasing die size and the introduction of new materials for enhanced thermal management, differences in thermal expansions of various packaging materials have become exceedingly important and can no longer be neglected. X-ray diffraction is an analytical method using a monochromatic characteristic X-ray beam to characterize the crystal structure of various materials, by measuring the distances between planes in atomic crystalline lattice structures. As a material is strained, this interplanar spacing is correspondingly altered, and this microscopic strain is used to determine the macroscopic strain. This thesis investigates and describes the theory and implementation of X-ray diffraction in the measurement of residual thermal strains. The design of a computer controlled stress attachment stage fully compatible with an Anton Paar heat stage will be detailed. The stress determined by the diffraction method will be compared with bimetallic strip theory and finite element models.

Composition of sediments from the Arabian Sea

Thirty seven deep-sea sediment cores from the Arabian Sea were studied geochemically (49 major and trace elements) for four time slices during the Holocene and the last glacial, and in one high sedimentation rate core (century scale resolution) to detect tracers of past variations in the intensity of the atmospheric monsoon circulation and its hydrographic expression in the ocean surface. This geochemical multi-tracer approach, coupled with additional information on the grain size composition of the clastic fraction, the bulk carbonate and biogenic opal contents makes it possible to characterize the sedimentological regime in detail. Sediments characterized by a specific elemental composition (enrichment) originated from the following sources: river suspensions from the Tapti and Narbada, draining the Indian Deccan traps (Ti, Sr); Indus sediments and dust from Rajasthan and Pakistan (Rb, Cs); dust from Iran and the Persian Gulf (Al, Cr); dust from central Arabia (Mg); dust from East Africa and the Red Sea (Zr/Hf, Ti/Al). Corg, Cd, Zn, Ba, Pb, U, and the HREE are associated with the intensity of upwelling in the western Arabian Sea, but only those patterns that are consistently reproduced by all of these elements can be directly linked with the intensity of the southwest monsoon. Relying on information from a single element can be misleading, as each element is affected by various other processes than upwelling intensity and nutrient content of surface water alone. The application of the geochemical multi-tracer approach indicates that the intensity of the southwest monsoon was low during the LGM, declined to a minimum from 15,000-13,000 14C year BP, intensified slightly at the end of this interval, was almost stable during the Bölling, Alleröd and the Younger Dryas, but then intensified in two abrupt successions at the end of the Younger Dryas (9900 14C year BP) and especially in a second event during the early Holocene (8800 14C year BP). Dust discharge by northwesterly winds from Arabia exhibited a similar evolution, but followed an opposite course: high during the LGM with two primary sources-the central Arabian desert and the dry Persian Gulf region. Dust discharge from both regions reached a pronounced maximum at 15,000-13,000 14C year. At the end of this interval, however, the dust plumes from the Persian Gulf area ceased dramatically, whereas dust discharge from central Arabia decreased only slightly. Dust discharge from East Africa and the Red Sea increased synchronously with the two major events of southwest monsoon intensification as recorded in the nutrient content of surface waters. In addition to the tracers of past dust flux and surface water nutrient content, the geochemical multi-tracer approach provides information on the history of deep sea ventilation (Mo, S), which was much lower during the last glacial maximum than during the Holocene. The multi-tracer approach-i.e. a few sedimentological parameters plus a set of geochemical tracers widely available from various multi-element analysis techniques-is a highly applicable technique for studying the complex sedimentation patterns of an ocean basin, and, specifically in the case of the Arabian Sea, can even reveal the seasonal structure of climate change.

Análise do fluxo d’água na barragem Engenheiro Armando Ribeiro Gonçalves - Açu/RN

The solution of partial differential equation of seepage problems is difficult to find analytically, especially for situations that involve great complexity. To overcome this problem, software based on finite differences and finite elements are usually used. This work presents the use of a finite element software, the GEO5, to solve the seepage problem at a dam of very complex section, the dam Eng. Armando Ribeiro Gonçalves, which at the end of its construction suffered rupture of the upstream slope at the central dam and then went through a process of reconstruction and auscultation. The analyses were performed for the operating condition of the reservoir, with an established flow. A numerical model was developed based on the level readings of the reservoir water and their piezometric readings as a proposal for the evaluation and future behavior prediction of the dam on established flow conditions. The use of constitutive models with the aid of computer systems is reflected in a way to predict future risk situations so they can be prevented

Comportamento mecânico e biológico de implante dentário com interface cônica interna

The general aim of this study was to evaluate the conical interface of pilar/implant. The specific aims were to evaluate the influence of hexagonal internal index in the microleakage and mechanical strength of Morse taper implants; the effect of axial loading on the deformation in cervical region of Morse taper implants of different diameters through strain gauge; the effect of axial loading in cervical deformation and sliding of abutment into the implant by tridimensional measurements; the integrity of conical interface before and after dynamic loading by microscopy and microleakage; and the stress distribution in tridimensional finite element models of Morse taper implants assembled with 2 pieces abutment. According to the obtained results, could be concluded that the diameter had influence in the cervical deformation of Morse taper implants; the presence of internal hexagonal index in the end of internal cone of implant didn´t influenced the bacterial microleakage under static loading neither reduced the mechanical strength of implants; one million cycles of vertical and off-center load had no negative influence in Morse taper implant integrity.

Avaliação comparativa da remodelação óssea periimplantar em implantes com junção cone-morse e hexágono externo: estudo clínico randomizado controlado, duplo cego, boca dividida

Background: Several theories, such as the biological width formation, the inflammatory reactions due to the implant-abutment microgap contamination, and the periimplant stress/strain concentration causing bone microdamage accumulation, have been suggested to explain early periimplant bone loss. However, it is yet not well understood to which extent the implant-abutment connection type may influence the remodeling process around dental implants. Aim: to evaluate clinical, bacteriological, and biomechanical parameters related to periimplant bone loss at the crestal region, comparing external hexagon (EH) and Morse-taper (MT) connections. Materials and methods: Twelve patients with totally edentulous mandibles received four custom made Ø 3.8 x 13 mm implants in the interforaminal region of the mandible, with the same design, but different prosthetic connections (two of them EH or MT, randomly placed based on a split-mouth design), and a immediate implant- supported prosthesis. Clinical parameters (periimplant probing pocket depth, modified gingival index and mucosal thickness) were evaluated at 6 sites around the implants, at a 12 month follow-up. The distance from the top of the implant to the first bone-to-implant contact – IT-FBIC was evaluated on standardized digital peri-apical radiographs acquired at 1, 3, 6 and 12 months follow-up. Samples of the subgingival microbiota were collected 1, 3 and 6 months after implant loading. DNA were extracted and used for the quantification of Tanerella forsythia, Porphyromonas gingivalis, Aggragatibacter actinomycetemcomitans, Prevotella intermedia and Fusobacterium nucleatum. Comparison among multiple periods of observation were performed using repeated-measures Analysis of Variance (ANOVA), followed by a Tukey post-hoc test, while two-period based comparisons were made using paired t- test. Further, 36 computer-tomographic based finite element (FE) models were accomplished, simulating each patient in 3 loading conditions. The results for the peak EQV strain in periimplant bone were interpreted by means of a general linear model (ANOVA). Results: The variation in periimplant bone loss assessed by means of radiographs was significantly different between the connection types (P<0.001). Mean IT-FBIC was 1.17±0.44 mm for EH, and 0.17±0.54 mm for MT, considering all evaluated time periods. All clinical parameters presented not significant differences. No significant microbiological differences could be observed between both connection types. Most of the collected samples had very few pathogens, meaning that these regions were healthy from a microbiological point of view. In FE analysis, a significantly higher peak of EQV strain (P=0.005) was found for EH (mean 3438.65 µ∑) compared to MT (mean 840.98 µ∑) connection. Conclusions: Varying implant-abutment connection type will result in diverse periimplant bone remodeling, regardless of clinical and microbiological conditions. This fact is more likely attributed to the singular loading transmission through different implant-abutment connections to the periimplant bone. The present findings suggest that Morse-taper connection is more efficient to prevent periimplant bone loss, compared to an external hexagon connection.