Iber: herramienta de simulación numérica del flujo en ríos

The recent requirements of Spanish regulations and directives, on their turn based on European directives, have led to the development of a new two dimensional open channel flow modelling tool. The tool, named Iber, combines a hydrodynamic module, a turbulence module and a sediment transport module, and is based in the finite volume method to solve the involved equations. The simulation code has been integrated in a pre-process and post-process interface based on GiD software, developed by CIMNE. The result is a flow and sediment modelling system for rivers and estuaries that uses advanced numerical schemes, robust and stable, which are especially suitable for discontinuous flows taking place in torrential and hydrologically irregular rivers.     

Evaluación numérica del efecto del espesor de la placa de contacto en la acción de palanca en conexión de acero tipo «T»

This paper presents a numerical study of «T-Stub» steel connections using FEM analyses. In such connections prying action phenomenon may take place. Prying action effect on the bolts has been studied before but the location of the prying action forces has always been simplified. The effect of the thickness of «T-Stub» flanges on contact areas between flanges and support base is investigated in this paper. A 3D finite element model is used and interface elements are employed for the investigation. Nonlinear FE analyses are undertaken on connection with different flange thickness and bolt preload with two bolts. Nonlinear gap elements are used as interfaces. Discussion and conclusions on contact areas and stresses, prying action distributions and resultant of loads on bolts in the «T-Stub» connections are presented.     

Modelación numérica de deslizamientos de ladera en embalses mediante el Método de Partículas y Elementos Finitos (PFEM)

The paper presents the results of the application of the Particle Finite Element Method (PFEM) to the analysis of landslides in reservoirs. This is a complex phenomenon, because of the interaction between the landslide, the still water in the reservoir and the dam. PFEM combines a Lagrangian approach with the solution of the governing equations of the problem using the FEM. A mesh connecting the initial set of particles (nodes) is re-generated in every time step. Some validation cases are presented, in which PFEM results are compared with experimental data. More complex calculations have been made over the actual geometry of reservoirs taken from the cartographic information of the sites. In these cases the wave generation, its propagation and dam overtopping are reproduced. Finally, Lituya bay rock slide in which 90 × 10⁶ tons of rocks fell on the bay, generating a huge wave that caused a maximum run-up of 524 m on the opposite shore, has been simulated in 3D. The results show that PFEM is a useful tool in risk assessment related with landslides in reservoirs as it gives a good approximation to the potential affections, thus allowing the appropriate design of protection measures.     

Simulación numérica del flujo turbulento de aire con gotas dispersas de agua a través de separadores de torres de refrigeración

This work presents a numerical study on the turbulent flow of air with dispersed water droplets in separators of mechanical cooling towers. The averaged Navier-Stokes equations are discretised through a finite volume method, using the Fluent and Phoenics codes, and alternatively employing the turbulence models k ? ?, k ? ω and the Reynolds stress model, with low-Re version and wall enhanced treatment refinements. The results obtained are compared with numerical and experimental results taken from the literature. The degree of accuracy obtained with each of the considered models of turbulence is stated. The influence of considering whether or not the simulation of the turbulent dispersion of droplets is analyzed, as well as the effects of other relevant parameters on the collection efficiency and the coefficient of pressure drop. Focusing on four specific eliminators (‘Belgian wave’, ‘H1-V’, ‘L-shaped’ and ‘Zig-zag’), the following ranges of parameters are outlined: 1 ≤ U_e ≤ 5 m/s for the entrance velocity, 2 ≤ D_p ≤ 50 μm for the droplet diameter, 650 ≤ Re ≤ 8.500 for Reynolds number, and 0.05 ≤ P_i ≤ 5 for the inertial parameter. Results reached alternately with Fluent and Phoenics codes are compared. The best results correspond to the simulations performed with Fluent, using the SST k ? ω turbulence model, with values of the dimensionless scaled distance to wall y⁺ in the range 0.2 to 0.5. Finally, correlations are presented to predict the conditions for maximum collection efficiency (100 %), depending on the geometric parameter of removal efficiency of each of the separators, which is introduced in this work.     

Modelación numérica de la descarga térmica de la Central Nucleoeléctrica Laguna Verde

This work contributes to the study of nuclear plant thermal discharges in coastal areas by using a numerical model which solves the Navier-Stokes-Reynolds equations for shallow waters and the energy equation for computing temperature variations. The numerical model takes into account the heat flux given in the upper layer, where the free surface and the atmosphere interact. In this study, the thermal plume dispersion from the nuclear power plant Laguna Verde, Veracruz, Mexico, is analyzed. Bathymetry, oceanographic, meteorological, hydrologic and plant operating data are used to run numerical simulations. The results are compared against observed data showing good agreement. The Nash-Suffle's criterion is also applied to verify the quality of the numerical solution obtaining suitable results.     

Caracterización de la delaminación en materiales compuestos mediante la teoría de mezclas serie/paralelo

This paper presents a new procedure to deal with the delamination problem found in laminated composites, based in a continuum mechanics formulation. The procedure proposed obtains the composite constitutive performance with the Serial/Parallel mixing theory, developed by F. Rastellini. This theory characterizes composite materials by coupling the constitutive behaviour of the composite components, imposing an iso–strain relation among the components in the fibre (or parallel) direction and an iso-stress relation in the remaining directions (serial directions). The proposed procedure uses a damage formulation to characterize the constitutive behaviour of matrix component in order to obtain the stress-strain performance of this material.With these two formulations, the delamination phenomenon is characterized naturally by the numerical simulation, being unnecessary the definition of special elements or computationally expensive techniques like the definition of contact elements or mesh separation. Matrix failure, as a result of the stress state found in it, leads to a reduction of the stiffness and strength capacity of the composite in its serial directions, among them, the shear component. This stiffness reduction provides a composite performance equivalent to what is found in a delaminated material.To prove the ability of the formulation proposed to solve delamination problems, the End Notch Failure test is numerically simulated and the results obtained are compared with experimental ones. The agreement found in the results with both simulations, numerical and experimental, validate the proposed methodology to solve the delamination problem.     

Validación numérica del comportamiento cíclico de conexiones interiores en estructuras con forjados reticulares

Several approaches are investigated to model interior reinforced concrete waffle-flat-plate-column connections. A model is proposed that provides very good results with reasonable low computational cost. The proposed model is validated with the experimental results obtained on a 3/5 scale specimen, subjected to quasi-static in cyclic loads up to collapse. To this end, the non-linear advanced theory of reinforced concrete is applied on a three-dimensional finite element model and non-linear analysis are conducted. Both fiber and layer elements are used for the one-dimensional and bi-dimensional components respectively. The main results of the simulation were: (i) the capacity curve obtained through out a push-over analysis with displacement control, (ii) the hysteretic curves of the slab, and (iii) the crack patterns. A very good agreement is found between numerical and experimental results.     

Mejora de la solución fuertemente acoplada de problemas FSI mediante una aproximación de la matriz tangente de presión

The interaction of fluids with surrounding structures constitutes a classical challenge for the different numerical techniques. The aim of current work is twofold: first we provide a simple theoretical explanation of the problems to be faced in incompressible FSI. Then we introduce and justify a new procedure for the solution of complex fluid-structure interaction problems. Such a new strategy is based on the introduction of an «interface Laplacian» at the coupling boundary. The idea is to consider the dependence between fluid pressure and structural velocity as a non linear problem for which a Quasi-Newton scheme is sought. The new interface term is then proved to be an approximation of the tangent matrix for such non-linear problem. In the derivation of this result we make use exclusively of discrete linear algebra. Finally, we prove the efficiency of the new approach showing its ability to tackle standard benchmark problems.     

Modelización numérica del comportamiento estructural cíclico de barras esbeltas de acero con pandeo restringido

This work presents a numerical model of the cyclic structural behavior of dissipative buckling-restrained braces, commonly used as an alternative to classical concentric braces for seismic protection of building frames and other structures. Such devices are usually composed of a slender steel core embedded in a stockiest casing that is intended to prevent its buckling when it is under compression. The casing is made either of mortar or steel, and a sliding interface is interposed between the core and the casing to prevent excessive shear stress transfer. The behavior of the steel core is described by a damage and plasticity model; the behavior of the mortar casing is described by an isotropic damage model and the sliding behavior of the interface is described by a contact penalty model. These 3 models are implemented in the Abaqus software package following an explicit formulation. In a previous article (published in an earthquake engineering journal) the model was briefly described, its ability to reproduce the cyclical behavior of buckling-restrained braces was preliminarily pointed out and their results were satisfactorily compared with those of experimental tests. The aim of this paper is to describe the model thoroughly and to present new judgments about its usefulness.     

Un modelo numérico para la simulación de disolución de precipitados en aleaciones de aluminio con endurecimiento utilizando redes neuronales

The motivation of this work is the modeling of the hardening precipitate and hardness evolutions of fully hardened heat treatable aluminium alloys during friction stir welding (FSW) and/or heat treatment processes. The model used is based on the kinetics of dissolution of precipitates model for hardened aluminium alloys given by Myhr and Grong (1991). This model contains a single independent variable, the time, and a single state variable, the volume fraction of hardening precipitates. A key point of this model is the identification of the effective activation energy for precipitates dissolution and the master curve defining the model, which was given by a look-up table. The goal of this work is to find an estimation of the effective activation energy and to model the dissolution rate of hardening precipitate in aluminium alloys using neural networks, avoiding the use of look-up tables. For this purpose a new and more convenient parametrization of the master curve is defined, a neural networks class is proposed, an objective functional is defined and a variational problem including independent parameters is solved. The novel methodology has been applied to different aluminium alloys, including the AA 6005A T6, AA 7449 T79 and AA 2198 T8. Experimental tests have been carried out for those aluminium alloys in order to get the HV1 hardness after isothermal heat treatments at different temperatures and for different treatment time durations. The effective activation energy for hardening precipitates dissolution and the master curve of the model have been obtained, using different network architectures, for the aluminium alloys considered in this work.     

Método numérico para la evaluación holística del riesgo sísmico utilizando la teoría de conjuntos difusos

The fuzzy set theory offers a bridge between the symbolic and numerical processing, allowing managing qualitative concepts useful in the decision-making process related to the seismic risk management and, in general, to the disaster risk management. Its use in the seismic risk evaluation is necessary in the cases where the data required to apply a conventional method of assessing risk are not available or are insufficient. One possible solution, considered in this article, is to replace the missing information by expert opinions and to process the resulting qualitative variables and linguistic qualifications instead of numerical values. This process is based on the fuzzy set theory. In order to achieve an effective management, the risk must be defined as the potential physical, economic, social and environmental consequences which occur due to hazards in a given period of time. From this holistic perspective and using the fuzzy set theory, the proposed numerical method calculates a level of the physical risk and level of the aggravating conditions related to social fragility and to the lack of resilience, to determine a total risk level. In the article are included two examples of application of the proposed method and the obtained results are compared with those corresponding to a conventional method of holistic evaluation.     

Resolución numérica de las ecuaciones del movimiento de edificios de varias plantas con no linealidades severas

This paper presents a new algorithm for solving the equations of motion of multi-storey buildings that incorporate frictional energy dissipators as seismic protection. The behavior of the dissipators is represented by Coulomb dry friction models; they introduce severe nonlinearities in the dynamic behavior of the structure every time that the contact conditions (stick or slip) change in the dissipators. These nonlinearities complicate the resolution of the equations of motion as it usually is described by lumped masses models whose degrees of freedom are the displacements of the floors and, as the stick or slip conditions change, the degrees of freedom must be modified: for blocking conditions they are only the displacements of the storeys while under sliding conditions the displacements of the dissipators have to be also considered. In previous articles the accuracy of the proposed algorithm has been verified by comparison with experimental results; as well, the computational efficiency of the algorithm has been confirmed by comparing the required resources (in terms of computation time and of memory allocation) with those of other algorithms. The objectives of this paper are to describe in detail the numerical solution of the equations of motion and present representative examples confirming the ability of the algorithm to reproduce the dynamic behavior of buildings with friction dissipators and reporting preliminarily about the usefulness of such devices to reduce the oscillations of the structure to be protected.     

Discusión sobre el artículo «Evaluación numérica del efecto del espesor de la placa de contacto en la acción de palanca en conexión de acero tipo T» de L.M. Bezerra,C.S. de Freitas,W.T. Matias y J.E. Carmona

A discussion is proposed on the paper “Evaluación numérica del efecto del espesor de la placa de contacto en la acción de palanca en conexión de acero tipo T” by L.M. Bezerra, C.S. de Freitas, W.T. Matias and J.E. Carmona. The topics to be discussed are the definition of the T-stub model, the bibliographic references used in the work and the finite element analysis assumptions adopted in the modeling of the components.     

Modelado numérico del proceso de soldadura FSW incorporando una técnica de estimación de parámetros

Numerical models of heat transfer and fluid flow used in the simulation of the friction-stir welding (FSW) process have contributed to the understanding of the process. However, there are some input model parameters that cannot be easily determined from fundamental principles or the welding conditions. As a result, the model predictions are not always in agreement with experimental results. In this work, the Levenberg-Marquardt (LM) method is used in order to perform a non-linear estimation of the unknown parameters present in the heat transfer and fluid flow models, by adjusting the temperatures results obtained with the models to temperature experimental measurements. These models are implemented in a general-purpose software that uses a numerical formulation developed from the finite element method (FEM). The unknown parameters are: the friction coefficient and the amount of adhesion of material to the surface of the tool, the heat transfer coefficient on the bottom surface and the amount of viscous dissipation converted into heat. The obtained results show an improvement in the numerical model predictions from the incorporation of parameter estimation techniques.     

Aplicación de métodos computacionales en la evaluación de la respuesta aeroelástica de puentes soportados por cables

The possibilities of computational methods for assessing the response of cable supported bridges under wind action are considered in this work. The main objective is to study the possibilities of substituting wind tunnel campaigns by computer based analyses, particularly at the early design stage. The preliminary proposed design for a continuous cable-stayed bridge with two main spans of 650 m and a single box girder deck has been considered as a case study. The force coefficients of the deck cross-section have been computed and the unsteady response associated to vortex-shedding has been simulated using CFD commercial software. Furthermore, an in-house piece of software has been employed to obtain the response for flutter and buffeting phenomena adopting the hybrid approach; with that purpose the experimental flutter functions of a similar box girder deck were adopted. The computational results have been validated by comparison with similar experimental results published by other researchers. It has been verified that the set of adopted methods offers reliable results with moderate costs; therefore, the proposed approach is very suitable at the early design stage of long span bridges or at conceptual design works.     

Una formulación de mínimo peso con restricciones en tensión en optimización topológica de estructuras

Optimum design of structures has been traditionally focused on the analysis of shape and dimensions optimization problems. However, more recently a new discipline has emerged: the topology optimization of the structures. This discipline states innovative models that allow to obtain optimal solutions without a previous definition of the type of structure being considered. These formulations obtain the optimal topology and the optimal shape and size of the resulting elements. The most usual formulations of the topology optimization problem try to obtain the structure of maximum stiffness. These approaches maximize the stiffness for a given amount of material to be used. These formulations have been widely analyzed and applied in engineering but they present considerable drawbacks from a numerical and from a practical point of view. In this paper the author propose a different formulation, as an alternative to maximum stiffness approaches, that minimizes the weight and includes stress constraints. The advantages of this kind of formulations are crucial since the cost of the structure is minimized, which is the most frequent objective in engineering, and they guarantee the structural feasibility since stresses are constrained. In addition, this approach allows to avoid some of the drawbacks and numerical instabilities related to maximum stiffness approaches. Finally, some practical examples have been solved in order to verify the validity of the results obtained and the advantages of the proposed formulation.     

Descomposición en pistones circulares para la obtención de la radiación acústica de un pistón rectangular

Obtaining pressure radiated by flat surfaces is not a new problem. This problem has studied by the complexity of the topic and its application to design flat speakers. These kind of speakers are the speakers that we use in televisions, ceilings, cinema screens, panels, etc. in this cases usually we have rectangular speakers.The single source model is the simplest model to convert the surface vibration to radiated pressure in a point. This is an easy model but it is very slow, especially when we want calculate at high frequencies. For rectangular surfaces there are models that use relatively complex auxiliary functions. In this case the calculation is accelerated but its implementation is more complicated and is necessary to particularize each situation.This paper presents the decomposition of a rectangular surface in several circular surfaces, by means of area associations, seeking a rapid method based on circular pistons whose behavior is known with an acceptable error in the allocation.     

Conmutación temporal entre técnicas numéricas implícita y explícita para simulación de estructuras no lineales inestables

Simulation problems involving non-linear materials imply in numerous cases divergence of the implicit method which use return mapping algorithms for modelling of the nonlinear response. A switching implicit-explicit numerical technique in the context of Finite Element Methods is presented in this paper. Implicit/explicit mesh partitions are not considered whatsoever. Formulation for application to nonlinear hyperelastic materials and nonlinear elastic-plastic materials is provided. Furthermore, the response of the solid subjected to large deformations is presented and is embedded in the proposed technique. Numerical tests for nonlinear problems (geometric and/or material) showed the accurateness of the technique.