Aplicación de un algoritmo evolutivo para la optimización del modelado computacional de la pirólisis de materiales

The pyrolysis models undergoes a fast development due to the enhancement of both the computational power and the new test used to characterize the behaviour of materials under thermal stresses, which were widely used in the field of the chemical engineering to obtain the reaction rates. Thus, these models allow us to characterize either the transient heat of material (thermal inertia) or the complete chemical scheme of mass loss processes (kinetic triplet). The pyrolysis model needs a several number of parameters what does the optimization of a suitable set of parameters a difficult task. Two kinds of materials have been investigated; the first one was a real material which mass loss process was characterized as “one-step” reaction and the second one as a “two-steps” process. Further, it has been analyzed the influence of some algorithm features (initial population number, parameter range, crossover influence) in the optimization time and also in the accuracy of results.     

Propuesta y validación de un algoritmo Simulated annealing modificado para la solución de problemas de optimización

Over the last decades, heuristic optimization methods based on imitating natural, biological, social or cultural processes on a computational level have aroused the interest of the scientific community due to their ability to effectively explore multimodal and multidimensional solution spaces. Despite all the papers published in the international literature, most heuristic algorithms still have low precision and accuracy. In this context, a modified Simulated annealing algorithm (MSAA) is proposed and validated for solving optimization problems. Performance evaluation was performed on test functions (benchmark functions) with and without restrictions reported in the international literature and practical design problems in civil engineering. In all cases analyzed MSAA obtained equal or better results than those reported by other authors, illustrating the applicability of the proposed algorithm.     

Desarrollo de un algoritmo de optimización global en colonias de hormigas con selección de región factible para espacios continuos

This study introduces a new algorithm for the ant colony optimization (ACO) method, which has been proposed to solve global optimization problems with continuous decision variables. This algorithm, namely ACO-FRS, involves a strategy for the selection of feasible regions during optimization search and it performs the exploration of the search space using a similar approach to that used by the ants during the search of food. Four variants of this algorithm have been tested in several benchmark problems and the results of this study have been compared with those reported in literature for other ACO-type methods for continuous spaces. Overall, the results show that the incorporation of the selection of feasible regions allows the performing of a global search to explore those regions with low level of pheromone, thus increasing the feasibility of ACO for finding the global optimal solution.     

Optimización mixta de estructuras de transporte de energía: aplicación del algoritmo de recocido simulado

A general methodology to optimize the weight of power transmission structures is presented in this article. This methodology is based on the simulated annealing algorithm defined by Kirkpatrick in the early ‘80s. This algorithm consists of a stochastic approach that allows to explore and analyze solutions that do not improve the objective function in order to develop a better exploration of the design region and to obtain the global optimum. The proposed algorithm allows to consider the discrete behavior of the sectional variables for each element and the continuous behavior of the general geometry variables. Thus, an optimization methodology that can deal with a mixed optimization problem and includes both continuum and discrete design variables is developed. In addition, it does not require to study all the possible design combinations defined by discrete design variables. The algorithm proposed usually requires to develop a large number of simulations (structural analysis in this case) in practical applications. Thus, the authors have developed first order Taylor expansions and the first order sensitivity analysis involved in order to reduce the CPU time required. Exterior penalty functions have been also included to deal with the design constraints. Thus, the general methodology proposed allows to optimize real power transmission structures in acceptable CPU time.     

Modelo de un proceso de cristalización continua de un sorbete por medio de la metodología de momentos

Freezing is an important step in the manufacturing process of ice-cream and sorbet, since the operating conditions have a strong influence on the micro-structure, and consequently on the sensorial attributes of the final product. This steep of freezing is carried out by a scraped surface heat exchanger (SSHE) where the product quality is conditioned by process conditions as the evaporation temperature of a refrigerant fluid, the mix flow rate, the dasher speed and the cylinder pressure due to the air introduction. In order to study the relevance of a control system based on the influence of process variables on product quality, this paper presents a model for a continuous crystallization of a sorbet using the method of moments, which is validated by experimental data.The model created by this methodology has been able to represent the influence of the process conditions during the crystallization of the sorbet on the final product characteristics such as crystal size and the draw temperature in the outlet of the SSHE in absence of air. The model based in moments is studied as a reduced model of the population balance equation and includes the phenomena of heterogeneous nucleation and growth. This model developed represents minimal computational requirements and is highly adapted for optimization and/or process control tasks.     

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.     

Cálculo de la función distancia para el método Level Set usando la formulación estabilizada USFEM/Rothe

In this work we use the Unusual Stabilized Finite Element Method (USFEM) associated to Rothe's method for solving the redistancing problem in the Level Set Method. Rothe's method is used first for advancing the solution in (pseudo)time and USFEM for solving the resulting steady advective–reaction problem in each time step. Several 2D problems are solved and results compared with SUPG scheme supplemented with a nonlinear discontinuity–capturing operator.     

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.     

Sobre la utilización de códigos de elementos finitos basados en mallados cartesianos en optimización estructural

The structural shape optimization is an iterative process built up by a higher level, which proposes the geometries to analyze, and a lower level which is in charge of analyzing, numerically, their structural response, usually by means of the Finite Element Method (FEM). These techniques normally report notorious advantages in an industrial environment, but their high computational cost is the main drawback. The efficiency of the global process requires the efficiency of both levels. This work focuses on the improvement of the efficiency of the lower level by using a methodology that uses a 2D linear elasticity code based on geometry-independent Cartesian grids, combined with FEM solution and recovery techniques, adapted to this framework. This mesh type simplifies the mesh generation and, in combination with a hierarchical data structure, reuses a great calculus amount. The recovery technique plays a double role: a) it is used in the Zienkiewicz-Zhu type error estimators allowing to quantify the FEM solution quality to guide the h-adaptive refinement process which minimizes the computational cost for a given accuracy; and b) it provides a solution, more accurate than the FEM one, that can be used. The numerical results, which include a comparative with a commercial code, show the effect of the proposed methodology improving the efficiency in the optimization process and in the solution quality.     

Comportamiento dinámico de viaductos cortos considerando la interacción vehículo-vía-estructura-suelo

The dynamic vehicle-track-bridge-soil interaction is studied on high speed lines. The analysis is carried out using a general and fully three dimensional multi-body-finite element-boundary element model, formulated in the time domain to predict vibrations due to the train passage over the bridge. The vehicle is modelled as a multi-body system, the track and the bridge are modelled using finite elements and the soil is considered as a homogeneous half-space by the boundary element method. Usually, moving force model and moving mass model are employed to study the dynamic response of bridges. In this work, the multi-body system allows one to consider the quasi-static and dynamic excitation mechanisms. Soil-structure interaction is taken into account on the dynamic structure behaviour on simply-supported short span bridges. The influence of soil-structure interaction is analysed in both resonant and non-resonant regimes.     

El concepto de hiperviga como generalización de los modelos unidimensionales en la estática de piezas alargadas

This paper is concerned with the general formulation of linear problems in rod elastostatics, and with the identification of their common formal and structural features, valid for every kinematical hypothesis. The generalized variables (section forces and generalized deformations) defining the 1-dimensional model are introduced in a consistent and natural way, through a convenient factorization of the density of complementary potential energy, for every kinematical constraint which can be expressed as a linear combination of the generalized displacements. The identification of this complementary energy function with the Hamiltonian functional of analytical mechanics allows a systematic procedure to construct the equations which rule this class of problems. In this frame, the main target is to establish the required conditions to write the rod equilibrium equations in a purely statical form (with no interplay of kinematical variables). We primarily conclude that this is possible when the kinematics of the cross-sections is constrained to a rigid body movement. As a consequence, the concept of hyperbeam is proposed in order to define those models with deformation modes beyond rigid body movements of the cross section: in them, section forces and generalized displacements are coupled in the equilibrium equations. This is related to the idea of local static indeterminacy (hyperstaticity) and justifies the new name.     

Sobre la necesidad de controlar el error de discretización de elementos finitos en optimización de forma estructural con algoritmos evolutivos

This work analyzes the influence of the discretization error contained in the Finite Element (FE) analyses of each design configuration proposed by the structural shape optimization algorithms over the behavior of the algorithm. The paper clearly shows that if FE analyses are not accurate enough, the final solution provided by the optimization algorithm will neither be optimal nor satisfy the constraints. The need for the use of adaptive FE analysis techniques in shape optimum design will be shown. The paper proposes the combination of two strategies to reduce the computational cost related to the use of mesh adaptivity in evolutionary optimization algorithms: (a) the use of the algorithm described by Bugeda et al. [1] which reduces the computational cost associated to the adaptive FE analysis of each geometrical configuration and, (b) the successive increase of the required accuracy of the FE analyses in order to obtain a considerable reduction of the computational cost in the early stages of the optimization process.     

Relación entre el parámetro flecha-luz y la carga crítica de pandeo lateral para arcos parabólicos biempotrados bajo diferentes hipótesis de cálculo

Curved graphs that relate the f/L ratio and the λ = qL³/EI_y parameter for clamped arches with the different calculation hypothesis are going to be proposed in this article.It has been necessary to determine the real critical lateral strength of a series of clamped arches considering different calculation hypothesis respect to geometrical linearity and a steel material characteristic.The elastic-plastic material hypothesis with non geometrical linearity considers the existence of residual stresses and initial lateral deflection.     

Integración temporal explícita con grandes pasos de tiempo de la ecuación de transmisión del calor

We present a novel fully explicit time integration method that remains stable for large time steps, requires neither matrix inversions nor solving a system of equations and therefore allows for nearly effort-less parallelization. In this first paper the proposed approach is applied to solve conduction heat transfer problems, showing that it is stable for any time step as is the case with implicit methods but with a much lower computation time.     

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.     

Determinación de la presión sonora radiada por pistones circulares y anillos no planos usando un método numérico simplificado

Although the study of the sound pressure radiation from membranes and plates is not new, current and future applications have produced a large body of recent research in the field. Several works have been published on the radiation from general plane surfaces and some particular geometries such as rectangular, circular, elliptic and annular. However, the case of sound radiation from non-planar axisymmetric rings that could be applied to the design of coaxial loudspeakers has not received much attention. In this article, a simplified numerical approach for determining the sound pressure radiated from symmetric non-planar pistons and rings is presented. The method can also include those cases having a radially-symmetric velocity distribution.     

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.     

Optimización robusta en aplicaciones aeronáuticas con la combinación de cálculo estocástico y algoritmos evolutivos

Uncertainties are a daily issue to deal with in aerospace engineering and applications. Robust optimization methods commonly use a random generation of the inputs and take advantage of multi-point criteria to look for robust solutions accounting with uncertainty definition. From the computational point of view, the application to coupled problems, like fluid-dynamics (CFD) or fluid-structure interaction (FSI), can be extremely expensive. This work presents a coupling between stochastic analysis techniques and evolutionary optimization algorithms for the definition of a stochastic robust optimization procedure. At first, a stochastic procedure is proposed to be applied into optimization problems. The proposed method has been applied to both CFD and FSI problems for the reduction of drag and flutter, respectively.     

Evaluación probabilista de la capacidad,fragilidad y daño sísmico de edificios de hormigón armado

Currently, many structures existing in seismic areas are highly vulnerable because they have been built without the use of seismic design codes or by using outdated codes. Often, methods for assessing the vulnerability of the structures do not take into account that their seismic behavior is dynamic and highly nonlinear and, moreover, that the structural characteristics and action have large uncertainties. This article aims to assess the vulnerability of structures taking into account that the mechanical properties of materials and the seismic action are random variables, by using advanced techniques based on the Monte Carlo method and on the nonlinear stochastic dynamics. The results obtained with these techniques are compared with those corresponding to a standard vulnerability assessment, based on deterministic models, in order to highlight the differences between both approaches. The main conclusion of this work is the need to address the vulnerability assessment problem from a probabilistic perspective which, combined with advanced nonlinear static and dynamic structural analysis techniques, provides a powerful tool giving information impossible to be captured by means of deterministic models. Finally, detailed results obtained for a building with waffle slabs, which is a structural typology widely used in Spain, are included and discussed.