Numerical modelling of axially loaded piled foundations with physical tests validation

L'oggetto della tesi è la modellazione numerica di pali di fondazione per turbine eoliche in ambiente offshore. Il metodo di modellazione comprende l'interpretazione di prove CPT per l'estrapolazione dei parametri geotecnici del terreno, la creazione di un modello dell'interazione tra struttura e terreno tramite il software agli elementi finiti Abaqus, la sua verifica, l'ottimizzazione dell'estensione del modello del terreno e della mesh e la simulazione numerica del comportamento di pali caricati assialmente in compressione. Grazie al confronto con i risultati di test su modelli fisici, eseguiti nel laboratorio dell'istituto Fraunhofer IWES di Hannover, di tre pali aventi la stessa geometria ma installati in punti con diverse condizioni di compattazione della sabbia, è stata possibile la messa a punto di una strategia di simulazione più accurata possibile e la sua validazione. A conclusione del lavoro è stato eseguito un esempio di design di una fondazione jacket per una turbina eolica offshore.

Numerical prediction of the lateral displacement capacity of flexural reinforced concrete shear walls

Previous earthquakes showed that shear wall damage could lead to catastrophic failures of the reinforced concrete building. The lateral load capacity of shear walls needs to be estimated to minimize associated losses during catastrophic events; hence it is necessary to develop and validate reliable and stable numerical methods able to converge to reasonable estimations with minimum computational effort. The beam-column 1-D line element with fiber-type cross-section model is a practical option that yields results in agreement with experimental data. However, shortcomings of using this model to predict the local damage response may come from the fact that the model requires fine calibration of material properties to overcome regularization and size effects. To reduce the mesh-dependency of the numerical model, a regularization method based on the concept of post-yield energy is applied in this work to both the concrete and the steel material constitutive laws to predict the nonlinear cyclic response and failure mechanism of concrete shear walls. Different categories of wall specimens known to produce a different response under in plane cyclic loading for their varied geometric and detailing characteristics are considered in this study, namely: 1) scaled wall specimens designed according to the European seismic design code and 2) unique full-scale wall specimens detailed according to the U.S. design code to develop a ductile behavior under cyclic loading. To test the boundaries of application of the proposed method, two full-scale walls with a mixed shear-flexure response and different values of applied axial load are also considered. The results of this study show that the use of regularized constitutive models considerably enhances the response predictions capabilities of the model with regards to global force-drift response and failure mode. The simulations presented in this thesis demonstrate the proposed model to be a valuable tool for researchers and engineers.

Modelagem matemática e simulação numérica do resfriamento rápido de morango com ar forçado

This work approaches the forced air cooling of strawberry by numerical simulation. The mathematical model that was used describes the process of heat transfer, based on the Fourier's law, in spherical coordinates and simplified to describe the one-dimensional process. For the resolution of the equation expressed for the mathematical model, an algorithm was developed based on the explicit scheme of the numerical method of the finite differences and implemented in the scientific computation program MATLAB 6.1. The validation of the mathematical model was made by the comparison between theoretical and experimental data, where strawberries had been cooled with forced air. The results showed to be possible the determination of the convective heat transfer coefficient by fitting the numerical and experimental data. The methodology of the numerical simulations was showed like a promising tool in the support of the decision to use or to develop equipment in the area of cooling process with forced air of spherical fruits.

Direct Computation of Operational Matrices for Polynomial Bases

Several numerical methods for boundary value problems use integral and differential operational matrices, expressed in polynomial bases in a Hilbert space of functions. This work presents a sequence of matrix operations allowing a direct computation of operational matrices for polynomial bases, orthogonal or not, starting with any previously known reference matrix. Furthermore, it shows how to obtain the reference matrix for a chosen polynomial base. The results presented here can be applied not only for integration and differentiation, but also for any linear operation.

Looking at the Gregory-Laflamme instability through quasinormal modes

We study evolution of gravitational perturbations of black strings. It is well known that for all wave numbers less than some threshold value, the black string is unstable against the scalar type of gravitational perturbations, which is named the Gregory-Laflamme instability. Using numerical methods, we find the quasinormal modes and time-domain profiles of the black string perturbations in the stable sector and also show the appearance of the Gregory-Laflamme instability in the time domain. The dependence of the black string quasinormal spectrum and late-time tails on such parameters as the wave vector and the number of extra dimensions is discussed. There is numerical evidence that at the threshold point of instability, the static solution of the wave equation is dominant. For wave numbers slightly larger than the threshold value, in the region of stability, we see tiny oscillations with very small damping rate. While, for wave numbers slightly smaller than the threshold value, in the region of the Gregory-Laflamme instability, we observe tiny oscillations with very small growth rate. We also find the level crossing of imaginary part of quasinormal modes between the fundamental mode and the first overtone mode, which accounts for the peculiar time domain profiles.

Dynamic behavior analysis of drill-threading process when machining AISI Al-Si-Cu(4) alloy

Conventional threading operations involve two distinct machining processes: drilling and threading. Therefore, it is time consuming for the tools must be changed and the workpiece has to be moved to another machine. This paper presents an analysis of the combined process (drilling followed by threading) using a single tool for both operations: the tap-milling tool. Before presenting the methodology used to evaluate this hybrid tool, the ODS (operating deflection shapes) basics is shortly described. ODS and finite element modeling (FEM) were used during this research to optimize the process aiming to achieve higher stable machining conditions and increasing the tool life. Both methods allowed the determination of the natural frequencies and displacements of the machining center and optimize the workpiece fixture system. The results showed that there is an excellent correlation between the dynamic stability of the machining center-tool holder and the tool life, avoiding a tool premature catastrophic failure. Nevertheless, evidence showed that the tool is very sensitive to work conditions. Undoubtedly, the use of ODS and FEM eliminate empiric decisions concerning the optimization of machining conditions and increase drastically the tool life. After the ODS and FEM studies, it was possible to optimize the process and work material fixture system and machine more than 30,000 threaded holes without reaching the tool life limit and catastrophic fail.

A coupled boundary element/differential equation method formulation for plate-beam interaction analysis

In this work, a new boundary element formulation for the analysis of plate-beam interaction is presented. This formulation uses a three nodal value boundary elements and each beam element is replaced by its actions on the plate, i.e., a distributed load and end of element forces. From the solution of the differential equation of a beam with linearly distributed load the plate-beam interaction tractions can be written as a function of the nodal values of the beam. With this transformation a final system of equation in the nodal values of displacements of plate boundary and beam nodes is obtained and from it, all unknowns of the plate-beam system are obtained. Many examples are analyzed and the results show an excellent agreement with those from the analytical solution and other numerical methods. (C) 2009 Elsevier Ltd. All rights reserved.

A discussion on volume change in the plastic phase

This technical note discusses the possibility of using a more simplified scheme to estimate the plastic multiplier when some material shows volume changes, e.g. soil, balsa wood foam and other similar materials. Two procedures regarding volume changes during the plastic phase are discussed here. The first one is the classic procedure applied to non-associative plasticity, for which a Drucker-Prager-like surface is adopted to represent the plastic potential. For the second procedure, the plastic potential is not explicitly known, however, its orthogonal direction is chosen respecting a plastic volume change parameter similar to Poisson`s ratio. Copyright (C) 2007 John Wiley & Sons, Ltd.

A simple way to introduce fibers into FEM models

This communication proposes a simple way to introduce fibers into finite element modelling. This is a promising formulation to deal with fiber-reinforced composites by the finite element method (FEM), as it allows the consideration of short or long fibers placed arbitrarily inside a continuum domain (matrix). The most important feature of the formulation is that no additional degree of freedom is introduced into the pre-existent finite element numerical system to consider any distribution of fiber inclusions. In other words, the size of the system of equations used to solve a non-reinforced medium is the same as the one used to solve the reinforced counterpart. Another important characteristic is the reduced work required by the user to introduce fibers, avoiding `rebar` elements, node-by-node geometrical definitions or even complex mesh generation. An additional characteristic of the technique is the possibility of representing unbounded stresses at the end of fibers using a finite number of degrees of freedom. Further studies are required for non-linear applications in which localization may occur. Along the text the linear formulation is presented and the bounded connection between fibers and continuum is considered. Four examples are presented, including non-linear analysis, to validate and show the capabilities of the formulation. Copyright (c) 2007 John Wiley & Sons, Ltd.

A hybrid-Trefftz formulation for plane elasticity with selective enrichment of the approximations

This work is related to the so-called non-conventional finite element formulations. Essentially, a methodology for the enrichment of the initial approximation which is typical of the meshless methods and based on the clouds concept is introduced in the hybrid-Trefftz formulation for plane elasticity. The formulation presented allows for the approximation and direct enrichment of two independent fields: stresses in the domains and displacements on the boundaries of the elements. Defined by a set of elements and interior boundaries sharing a common node, the cloud notion is employed to select the enrichment support for the approximation fields. The numerical analysis performed reveals an excellent performance of the resulting formulation, characterized by the good approximation ability and a reduced computational effort. Copyright (C) 2009 John Wiley & Sons, Ltd.

Solution of a heterogeneous modeling of a horizontal-flow anaerobic immobilized biomass (HAIB) reactor by the sequencing method

A modeling study was completed to develop a methodology that combines the sequencing and finite difference methods for the simulation of a heterogeneous model of a tubular reactor applied in the treatment of wastewater. The system included a liquid phase (convection diffusion transport) and a solid phase (diffusion reaction) that was obtained by completing a mass balance in the reactor and in the particle, respectively. The model was solved using a pilot-scale horizontal-flow anaerobic immobilized biomass (HAIB) reactor to treat domestic sewage, with the concentration results compared with the experimental data. A comparison of the behavior of the liquid phase concentration profile and the experimental results indicated that both the numerical methods offer a good description of the behavior of the concentration along the reactor. The advantage of the sequencing method over the finite difference method is that it is easier to apply and requires less computational time to model the dynamic simulation of outlet response of HAIB.

Verification of a mixed high-order accurate DNS code for laminar turbulent transition by the method of manufactured solutions

This paper presents results on a verification test of a Direct Numerical Simulation code of mixed high-order of accuracy using the method of manufactured solutions (MMS). This test is based on the formulation of an analytical solution for the Navier-Stokes equations modified by the addition of a source term. The present numerical code was aimed at simulating the temporal evolution of instability waves in a plane Poiseuille flow. The governing equations were solved in a vorticity-velocity formulation for a two-dimensional incompressible flow. The code employed two different numerical schemes. One used mixed high-order compact and non-compact finite-differences from fourth-order to sixth-order of accuracy. The other scheme used spectral methods instead of finite-difference methods for the streamwise direction, which was periodic. In the present test, particular attention was paid to the boundary conditions of the physical problem of interest. Indeed, the verification procedure using MMS can be more demanding than the often used comparison with Linear Stability Theory. That is particularly because in the latter test no attention is paid to the nonlinear terms. For the present verification test, it was possible to manufacture an analytical solution that reproduced some aspects of an instability wave in a nonlinear stage. Although the results of the verification by MMS for this mixed-order numerical scheme had to be interpreted with care, the test was very useful as it gave confidence that the code was free of programming errors. Copyright (C) 2009 John Wiley & Sons, Ltd.

Shell curvature as an initial deformation: A geometrically exact finite element approach

An alternative approach for the analysis of arbitrarily curved shells is developed in this paper based on the idea of initial deformations. By `alternative` we mean that neither differential geometry nor the concept of degeneration is invoked here to describe the shell surface. We begin with a flat reference configuration for the shell mid-surface, after which the initial (curved) geometry is mapped as a stress-free deformation from the plane position. The actual motion of the shell takes place only after this initial mapping. In contrast to classical works in the literature, this strategy enables the use of only orthogonal frames within the theory and therefore objects such as Christoffel symbols, the second fundamental form or three-dimensional degenerated solids do not enter the formulation. Furthermore, the issue of physical components of tensors does not appear. Another important aspect (but not exclusive of our scheme) is the possibility to describe exactly the initial geometry. The model is kinematically exact, encompasses finite strains in a totally consistent manner and is here discretized under the light of the finite element method (although implementation via mesh-free techniques is also possible). Assessment is made by means of several numerical simulations. Copyright (C) 2009 John Wiley & Sons, Ltd.

Multi-actuated functionally graded piezoelectric micro-tools design: A multiphysics topology optimization approach

Micro-tools offer significant promise in a wide range of applications Such as cell Manipulation, microsurgery, and micro/nanotechnology processes. Such special micro-tools consist of multi-flexible structures actuated by two or more piezoceramic devices that must generate output displacements and forces lit different specified points of the domain and at different directions. The micro-tool Structure acts as a mechanical transformer by amplifying and changing the direction of the piezoceramics Output displacements. The design of these micro-tools involves minimization of the coupling among movements generated by various piezoceramics. To obtain enhanced micro-tool performance, the concept of multifunctional and functionally graded materials is extended by, tailoring elastic and piezoelectric properties Of the piezoceramics while simultaneously optimizing the multi-flexible structural configuration using multiphysics topology optimization. The design process considers the influence of piezoceramic property gradation and also its polarization sign. The method is implemented considering continuum material distribution with special interpolation of fictitious densities in the design domain. As examples, designs of a single piezoactuator, an XY nano-positioner actuated by two graded piezoceramics, and a micro-gripper actuated by three graded piezoceramics are considered. The results show that material gradation plays an important role to improve actuator performance, which may also lead to optimal displacements and coupling ratios with reduced amount of piezoelectric material. The present examples are limited to two-dimensional models because many of the applications for Such micro-tools are planar devices. Copyright (c) 2008 John Wiley & Sons, Ltd.

Mechanical strength evaluation for Nd-YAG laser and electric resistance spot weld (ERSW) joint under multiaxial loading

This work presents a comparison between laser weld (LBW) and electric resistance spot weld (ERSW) processes used for assemblies of components in a body-in-white (BIW) at a world class automotive industry. It is carried out by evaluating the mechanical strength modeled both by experimental and numerical methods. An ""Arcan"" multiaxial test was designed and manufactured in order to enable 0 degrees, 45 degrees and 90 degrees directional loadings. The welded specimens were uncoated low carbon steel sheets (S-y = 170 MPa) used currently at the automotive industry, with two different thicknesses: 0.80 and 1.20 mm. A numerical analysis was carried out using the finite element method (FEM) through LS-DYNA code. (c) 2007 Elsevier B.V. All rights reserved.