Elasticity, yielding and episodicity in simple models of mantle convection

We explore the implications of refinements in the mechanical description of planetary constituents on the convection modes predicted by finite-element simulations. The refinements consist in the inclusion of incremental elasticity, plasticity (yielding) and multiple simultaneous creep mechanisms in addition to the usual visco-plastic models employed in the context of unified plate-mantle models. The main emphasis of this paper rests on the constitutive and computational formulation of the model. We apply a consistent incremental formulation of the non-linear governing equations avoiding the computationally expensive iterations that are otherwise necessary to handle the onset of plastic yield. In connection with episodic convection simulations, we point out the strong dependency of the results on the choice of the initial temperature distribution. Our results also indicate that the inclusion of elasticity in the constitutive relationships lowers the mechanical energy associated with subduction events.

Deformation and failure of welded steels used in offshore constructions

Hydrocarbons are the most common form of energy used to date. The activities involving exploration and exploitation of large oil and gas fields are constantly in operation and have extended to such hostile environments as the North Sea. This enforces much greater demands on the materials which are used, and the need for enhancing the endurance of the existing ones which must continue parallel to the explorations. Due to their ease in fabrication, relatively high mechanical properties and low costs, steels are the most widely favoured material for the construction of offshore platforms. The most critical part of an offshore structure prone to failure are the welded nodal joints, particulary those which are used within the vicinity of the splash zones. This is an area of high complex stress concentrations, varying mechanical and metallurgical properties in addition to severe North Sea environmental conditions. The main are of this work has been concerned with the durability studies of this type of steel, based on the concept of the worst case analysis, consisting of combinations of welds of varying qualities, various degrees of stress concentrations and the environmental conditions of stress corrosion and hydrogen embrittlement. The experiments have been designed to reveal significance of defects as sites of crack initiation in the welded steels and the extent to which stress corrosion and hydrogen embrittlement will limit their durability. This has been done for various heat treatments and in some experiments deformation has been forced through the welded zone of the specimens to reveal the mechanical properties of the welds themselves to provide data for finite element simulations. A comparison of the results of these simulations with the actual deformation and fracture behaviour has been done to reveal the extent to which both mechanical and metallurgical factors control behaviour of the steels in the hostile environments of high stress, corrosion, and hydrogen embrittlement at their surface.

Pattern formation in volumetrically heated fluids

Finite element simulations have been performed along side normal mode analysis on the linear stability that examined the development of volumetrically heated flow patterns in a horizontal layer controlled by the Prandtl number, Pr, and the Grashof number, Gr. The fluid was bounded by an isothermal plane above an adiabatic plane. In the simulations performed here, a number of convective polygonal planforms occurred, as Gr increased above the critical Grashof number, Grc at Pr = 7, while roll structures were observed for Pr < 1 at 2Grc.

Pattern competition in volumetricaly heated fluids

Finite element simulations have been performed along side Galerkin-type calculations that examined the development of volumetrically heated flow patterns in a horizontal layer controlled by the Prandtl number, Pr, and the Grashof number, Gr. The fluid was bounded by an isothermal plane above an adiabatic plane. In the simulations performed here, a number of convective polygonal planforms occurred, as Gr increased above the critical Grashof number, Grc at Pr = 7, while roll structures were observed for Pr < 1 at 2Grc.

Using overlay test to evaluate fracture properties of field-aged asphalt concrete

Field material testing provides firsthand information on pavement conditions which are most helpful in evaluating performance and identifying preventive maintenance or overlay strategies. High variability of field asphalt concrete due to construction raises the demand for accuracy of the test. Accordingly, the objective of this study is to propose a reliable and repeatable methodology to evaluate the fracture properties of field-aged asphalt concrete using the overlay test (OT). The OT is selected because of its efficiency and feasibility for asphalt field cores with diverse dimensions. The fracture properties refer to the Paris’ law parameters based on the pseudo J-integral (A and n) because of the sound physical significance of the pseudo J-integral with respect to characterizing the cracking process. In order to determine A and n, a two-step OT protocol is designed to characterize the undamaged and damaged behaviors of asphalt field cores. To ensure the accuracy of determined undamaged and fracture properties, a new analysis method is then developed for data processing, which combines the finite element simulations and mechanical analysis of viscoelastic force equilibrium and evolution of pseudo displacement work in the OT specimen. Finally, theoretical equations are derived to calculate A and n directly from the OT test data. The accuracy of the determined fracture properties is verified. The proposed methodology is applied to a total of 27 asphalt field cores obtained from a field project in Texas, including the control Hot Mix Asphalt (HMA) and two types of warm mix asphalt (WMA). The results demonstrate a high linear correlation between n and −log A for all the tested field cores. Investigations of the effect of field aging on the fracture properties confirm that n is a good indicator to quantify the cracking resistance of asphalt concrete. It is also indicated that summer climatic condition clearly accelerates the rate of aging. The impact of the WMA technologies on fracture properties of asphalt concrete is visualized by comparing the n-values. It shows that the Evotherm WMA technology slightly improves the cracking resistance, while the foaming WMA technology provides the comparable fracture properties with the HMA. After 15 months aging in the field, the cracking resistance does not exhibit significant difference between HMA and WMAs, which is confirmed by the observations of field distresses.

Improved methodology to evaluate fracture properties of warm-mix asphalt using overlay test

This study developed a reliable and repeatable methodology to evaluate the fracture properties of asphalt mixtures with an overlay test (OT). In the proposed methodology, first, a two-step OT protocol was used to characterize the undamaged and damaged behaviors of asphalt mixtures. Second, a new methodology combining the mechanical analysis of viscoelastic force equilibrium in the OT specimen and finite element simulations was used to determine the undamaged properties and crack growth function of asphalt mixtures. Third, a modified Paris's law replacing the stress intensity factor by the pseudo J-integral was employed to characterize the fracture behavior of asphalt mixtures. Theoretical equations were derived to calculate the parameters A and n (defined as the fracture properties) in the modified Paris's law. The study used a detailed example to calculate A and n from the OT data. The proposed methodology was successfully applied to evaluate the impact of warm-mix asphalt (WMA) technologies on fracture properties. The results of the tested specimens showed that Evotherm WMA technology slightly improved the cracking resistance of asphalt mixtures, while foaming WMA technology provided comparable fracture properties. In addition, the study found that A decreased with the increase in n in general. A linear relationship between 2log(A) and n was established.

Multimode Analysis of Nanoscale Biomolecular Interactions

Biomolecular interactions, including protein-protein, protein-DNA, and protein-ligand interactions, are of special importance in all biological systems. These interactions may occer during the loading of biomolecules to interfaces, the translocation of biomolecules through transmembrane protein pores, and the movement of biomolecules in a crowded intracellular environment. The molecular interaction of a protein with its binding partners is crucial in fundamental biological processes such as electron transfer, intracellular signal transmission and regulation, neuroprotective mechanisms, and regulation of DNA topology. In this dissertation, a customized surface plasmon resonance (SPR) has been optimized and new theoretical and label free experimental methods with related analytical calculations have been developed for the analysis of biomolecular interactions. Human neuroglobin (hNgb) and cytochrome c from equine heart (Cyt c) proteins have been used to optimize the customized SPR instrument. The obtained Kd value (~13 µM), from SPR results, for Cyt c-hNgb molecular interactions is in general agreement with a previously published result. The SPR results also confirmed no significant impact of the internal disulfide bridge between Cys 46 and Cys 55 on hNgb binding to Cyt c. Using SPR, E. coli topoisomerase I enzyme turnover during plasmid DNA relaxation was found to be enhanced in the presence of Mg2+. In addition, a new theoretical approach of analyzing biphasic SPR data has been introduced based on analytical solutions of the biphasic rate equations. In order to develop a new label free method to quantitatively study protein-protein interactions, quartz nanopipettes were chemically modified. The derived Kd (~20 µM) value for the Cyt c-hNgb complex formations matched very well with SPR measurements (Kd ~16 µM). The finite element numerical simulation results were similar to the nanopipette experimental results. These results demonstrate that nanopipettes can potentially be used as a new class of a label-free analytical method to quantitatively characterize protein-protein interactions in attoliter sensing volumes, based on a charge sensing mechanism. Moreover, the molecule-based selective nature of hydrophobic and nanometer sized carbon nanotube (CNT) pores was observed. This result might be helpful to understand the selective nature of cellular transport through transmembrane protein pores.

Sub-wavelength resonators as seismic Rayleigh waves shield

We explore the thesis that tall structures can be protected by means of seismic metamaterials. Seismic metamaterials can be built as some elements are created over soil layer with different shapes, dimensions, patterns and from different materials. Resonances in these elements are acting as locally resonant metamaterials for Rayleigh surface waves in the geophysics context. Analytically we proved that if we put infinite chain of SDOF resonator over the soil layer as an elastic, homogeneous and isotropic material, vertical component of Rayleigh wave, longitudinal resonance of oscillators will couple with each other, they would create a Rayleigh bandgap frequency, and wave will experience attenuation before it reaches the structure. As it is impossible to use infinite chain of resonators over soil layer, we considered finite number of resonators throughout our simulations. Analytical work is interpreted using finite element simulations that demonstrates the observed attenuation is due to bandgaps when oscillators are arranged at sub-wavelength scale with respect to the incident Rayleigh wave. For wavelength less than 5 meters, the resulting bandgaps are remarkably large and strongly attenuating when impedance of oscillators matches impedance of soil. Since longitudinal resonance of SDOF resonator are proportional to its length inversely, a formed array of resonators that attenuates Rayleigh waves at frequency ≤10 Hz could be designed starting from vertical pillars coupled to the ground. Optimum number of vertical pillars and their interval spacing called effective area of resonators are investigated. For 10 pillars with effective area of 1 meter and resonance frequency of 4.9 Hz, bandgap frequency causes attenuation and a sinusoidal impulsive force illustrate wave steering down phenomena. Simulation results proved analytical findings of this work.

Application of Numerical Methods to Study Arrangement and Fracture of Lithium-Ion Microstructure

The focus of this work is to develop and employ numerical methods that provide characterization of granular microstructures, dynamic fragmentation of brittle materials, and dynamic fracture of three-dimensional bodies.

We first propose the fabric tensor formalism to describe the structure and evolution of lithium-ion electrode microstructure during the calendaring process. Fabric tensors are directional measures of particulate assemblies based on inter-particle connectivity, relating to the structural and transport properties of the electrode. Applying this technique to X-ray computed tomography of cathode microstructure, we show that fabric tensors capture the evolution of the inter-particle contact distribution and are therefore good measures for the internal state of and electronic transport within the electrode.

We then shift focus to the development and analysis of fracture models within finite element simulations. A difficult problem to characterize in the realm of fracture modeling is that of fragmentation, wherein brittle materials subjected to a uniform tensile loading break apart into a large number of smaller pieces. We explore the effect of numerical precision in the results of dynamic fragmentation simulations using the cohesive element approach on a one-dimensional domain. By introducing random and non-random field variations, we discern that round-off error plays a significant role in establishing a mesh-convergent solution for uniform fragmentation problems. Further, by using differing magnitudes of randomized material properties and mesh discretizations, we find that employing randomness can improve convergence behavior and provide a computational savings.

The Thick Level-Set model is implemented to describe brittle media undergoing dynamic fragmentation as an alternative to the cohesive element approach. This non-local damage model features a level-set function that defines the extent and severity of degradation and uses a length scale to limit the damage gradient. In terms of energy dissipated by fracture and mean fragment size, we find that the proposed model reproduces the rate-dependent observations of analytical approaches, cohesive element simulations, and experimental studies.

Lastly, the Thick Level-Set model is implemented in three dimensions to describe the dynamic failure of brittle media, such as the active material particles in the battery cathode during manufacturing. The proposed model matches expected behavior from physical experiments, analytical approaches, and numerical models, and mesh convergence is established. We find that the use of an asymmetrical damage model to represent tensile damage is important to producing the expected results for brittle fracture problems.

The impact of this work is that designers of lithium-ion battery components can employ the numerical methods presented herein to analyze the evolving electrode microstructure during manufacturing, operational, and extraordinary loadings. This allows for enhanced designs and manufacturing methods that advance the state of battery technology. Further, these numerical tools have applicability in a broad range of fields, from geotechnical analysis to ice-sheet modeling to armor design to hydraulic fracturing.

Experimental study of forces and energies during shearing of steel sheet with angled tools

Shearing is a fast and inexpensive method to cut sheet metal that has been used since the beginning of the industrialism. Consequently, published experimental studies of shearing can be found from over a century back in time. Recent studies, however, are due to the availability of low cost digital computation power, mostly based on finite element simulations that guarantees quick results. Still, for validation of models and simulations, accurate experimental data is a requisite. When applicable, 2D models are in general desirable over 3D models because of advantages like low computation time and easy model formulation. Shearing of sheet metal with parallel tools is successfully modelled in 2D with a plane strain approximation, but with angled tools the approximation is less obvious. Therefore, plane strain approximations for shearing with angled tools were evaluated by shear experiments of high accuracy. Tool angle, tool clearance, and clamping of the sheet were varied in the experiments. The results showed that the measured forces in shearing with angled tools can be approximately calculated using force measurements from shearing with parallel tools. Shearing energy was introduced as a quantifiable measure of suitable tool clearance range. The effects of the shearing parameters on forces were in agreement with previous studies. Based on the agreement between calculations and experiments, analysis based on a plane strain assumption is considered applicable for angled tools with a small (up to 2 degrees) rake angle.

Quality balancing for parallel adaptive FEM

We present a dynamic distributed load balancing algorithm for parallel, adaptive finite element simulations using preconditioned conjugate gradient solvers based on domain-decomposition. The load balancer is designed to maintain good partition aspect ratios. It can calculate a balancing flow using different versions of diffusion and a variant of breadth first search. Elements to be migrated are chosen according to a cost function aiming at the optimization of subdomain shapes. We show how to use information from the second step to guide the first. Experimental results using Bramble's preconditioner and comparisons to existing state-ot-the-art load balancers show the benefits of the construction.

Load balancing for parallel adaptive FEM

We present a dynamic distributed load balancing algorithm for parallel, adaptive finite element simulations using preconditioned conjugate gradient solvers based on domain-decomposition. The load balancer is designed to maintain good partition aspect ratios. It calculates a balancing flow using different versions of diffusion and a variant of breadth first search. Elements to be migrated are chosen according to a cost function aiming at the optimization of subdomain shapes. We show how to use information from the second step to guide the first. Experimental results using Bramble's preconditioner and comparisons to existing state-of-the-art balancers show the benefits of the construction.

Análise da fluência em vigas de concreto armado simplesmente apoiadas usando o método dos elementos finitos

Reinforced concrete creep is a phenomenon of great importance. Despite being appointed as the main cause of several pathologies, its effects are yet considered in a simplified way by the structural designers. In addition to studying the phenomenon in reinforced concrete structures and its current account used in the structural analysis, this paper compares creep strains at simply supported reinforced concrete beams in analytical and in experimental forms with the finite element method (FEM) simulation results. The strains and deflections obtained through the analytical form were calculated with the Brazilian code NBR 6118 (2014) recommendations and the simplified method from CEB-FIP 90 and the experimental results were extracted from tests available in the literature. Finite element simulations are performed using ANSYS Workbench software, using its 3D SOLID 186 elements and the structure symmetry. Analyzes of convergence using 2D PLANE 183 elements are held as well. At the end, it is concluded that FEM analyses are quantitative and qualitative efficient for the estimation of this non-linearity and that the method utilized to obtain the creep coefficients values is sufficiently accurate.

Estudo comparativo de pás para aerogeradores de grande porte fabricadas em materiais compósitos reforçadas com fibra de carbono ou fibra de vidro

The research and development of wind turbine blades are essential to keep pace with worldwide growth in the renewable energy sector. Although currently blades are typically produced using glass fiber reinforced composite materials, the tendency for larger size blades, particularly for offshore applications, has increased the interest on carbon fiber reinforced composites because of the potential for increased stiffness and weight reduction. In this study a model of blade designed for large generators (5 MW) was studied on a small scale. A numerical simulation was performed to determine the aerodynamic loading using a Computational Fluid Dynamics (CFD) software. Two blades were then designed and manufactured using epoxy matrix composites: one reinforced with glass fibers and the other with carbon fibers. For the structural calculations, maximum stress failure criterion was adopted. The blades were manufactured by Vacuum Assisted Resin Transfer Molding (VARTM), typical for this type of component. A weight comparison of the two blades was performed and the weight of the carbon fiber blade was approximately 45% of the weight of the fiberglass reinforced blade. Static bending tests were carried out on the blades for various percentages of the design load and deflections measurements were compared with the values obtained from finite element simulations. A good agreement was observed between the measured and calculated deflections. In summary, the results of this study confirm that the low density combined with high mechanical properties of carbon fibers are particularly attractive for the production of large size wind turbine blades

Design and application of experimental methods for steel sheet shearing

Shearing is the process where sheet metal is mechanically cut between two tools. Various shearing technologies are commonly used in the sheet metal industry, for example, in cut to length lines, slitting lines, end cropping etc. Shearing has speed and cost advantages over competing cutting methods like laser and plasma cutting, but involves large forces on the equipment and large strains in the sheet material. The constant development of sheet metals toward higher strength and formability leads to increased forces on the shearing equipment and tools. Shearing of new sheet materials imply new suitable shearing parameters. Investigations of the shearing parameters through live tests in the production are expensive and separate experiments are time consuming and requires specialized equipment. Studies involving a large number of parameters and coupled effects are therefore preferably performed by finite element based simulations. Accurate experimental data is still a prerequisite to validate such simulations. There is, however, a shortage of accurate experimental data to validate such simulations. In industrial shearing processes, measured forces are always larger than the actual forces acting on the sheet, due to friction losses. Shearing also generates a force that attempts to separate the two tools with changed shearing conditions through increased clearance between the tools as result. Tool clearance is also the most common shearing parameter to adjust, depending on material grade and sheet thickness, to moderate the required force and to control the final sheared edge geometry. In this work, an experimental procedure that provides a stable tool clearance together with accurate measurements of tool forces and tool displacements, was designed, built and evaluated. Important shearing parameters and demands on the experimental set-up were identified in a sensitivity analysis performed with finite element simulations under the assumption of plane strain. With respect to large tool clearance stability and accurate force measurements, a symmetric experiment with two simultaneous shears and internal balancing of forces attempting to separate the tools was constructed. Steel sheets of different strength levels were sheared using the above mentioned experimental set-up, with various tool clearances, sheet clamping and rake angles. Results showed that tool penetration before fracture decreased with increased material strength. When one side of the sheet was left unclamped and free to move, the required shearing force decreased but instead the force attempting to separate the two tools increased. Further, the maximum shearing force decreased and the rollover increased with increased tool clearance. Digital image correlation was applied to measure strains on the sheet surface. The obtained strain fields, together with a material model, were used to compute the stress state in the sheet. A comparison, up to crack initiation, of these experimental results with corresponding results from finite element simulations in three dimensions and at a plane strain approximation showed that effective strains on the surface are representative also for the bulk material. A simple model was successfully applied to calculate the tool forces in shearing with angled tools from forces measured with parallel tools. These results suggest that, with respect to tool forces, a plane strain approximation is valid also at angled tools, at least for small rake angles. In general terms, this study provide a stable symmetric experimental set-up with internal balancing of lateral forces, for accurate measurements of tool forces, tool displacements, and sheet deformations, to study the effects of important shearing parameters. The results give further insight to the strain and stress conditions at crack initiation during shearing, and can also be used to validate models of the shearing process.