A Finite Element Numerical Algorithm for Modelling and Data Fitting in Complex Systems

Numerical modelling methodologies are important by their application to engineering and scientific problems, because there are processes where analytical mathematical expressions cannot be obtained to model them. When the only available information is a set of experimental values for the variables that determine the state of the system, the modelling problem is equivalent to determining the hyper-surface that best fits the data. This paper presents a methodology based on the Galerkin formulation of the finite elements method to obtain representations of relationships that are defined a priori, between a set of variables: y = z(x1, x2,...., xd). These representations are generated from the values of the variables in the experimental data. The approximation, piecewise, is an element of a Sobolev space and has derivatives defined in a general sense into this space. The using of this approach results in the need of inverting a linear system with a structure that allows a fast solver algorithm. The algorithm can be used in a variety of fields, being a multidisciplinary tool. The validity of the methodology is studied considering two real applications: a problem in hydrodynamics and a problem of engineering related to fluids, heat and transport in an energy generation plant. Also a test of the predictive capacity of the methodology is performed using a cross-validation method.

Barandillas de protección personal con solicitación dinámica

Las barandillas de seguridad con solicitación dinámica, tipos B y C según UNE-EN 13374, son el conjunto de componentes destinados a la protección colectiva de las personas contra las caídas a un nivel inferior y a retener materiales. Frecuentemente incorporan una red de seguridad. En este artículo se analizan aspectos fundamentales como su diseño, sus requisitos mecánicos y la solución a la problemática que suponen los soportes o puntos duros capaces de producir lesiones al accidentado. Para este fin, se ha estudiado el fenómeno de la caída de un lastre sobre tales redes con distintos modelos numéricos de elementos finitos en régimen dinámico y mediante consideraciones analíticas cinemáticas, obteniendo conclusiones útiles como la forma geométrica necesaria, la escuadría de los perfiles que constituyen el bastidor de la red, y los factores de impacto que sufren los accidentados. Ello permite sugerir mejoras al texto normativo.

Redes de seguridad para cierre vertical de fachadas

Entre otros medios de protección colectiva, en los últimos años se vienen utilizando en España redes verticales de cierre de fachadas, a paño completo entre forjados en fase de estructura, y también para huecos menores como ventanas o puertas, en fases posteriores. Este tipo de protección no cuenta hasta el momento con ningún tipo de regulación oficial, si bien en el grupo de trabajo AEN/CTN81/SC2/GT7 se desarrolla actualmente el borrador de la norma Pr UNE 81651 a tal fin. Desde el Dtº. de Ing. de la Construcción de la Universidad de Alicante, varios miembros del grupo de trabajo han desarrollado, contando con la experiencia en otros sistemas de protección, algunos modelos numéricos con elementos finitos para valorar y cuantificar variables mecánicas asociadas al episodio de una eventual caída de una persona sobre este tipo de redes, caracterizado por una menor energía cinética que para otros tipos de ellas. En el artículo se exponen los resultados obtenidos para ambos tipos de paños, grandes y pequeños, donde se han valorado la separación de anclajes de la red, las aberturas generadas en su perímetro durante la retención, los esfuerzos en anclajes y elementos textiles, el factor de impacto que sufriría el accidentado y la repercusión de la disposición geométrica de la red, al cuadro o al rombo.

Generation of representation models for complex systems using Lagrangian functions

In this article, a new methodology is presented to obtain representation models for a priori relation z = u(x1, x2, . . . ,xn) (1), with a known an experimental dataset zi; x1i ; x2i ; x3i ; . . . ; xni i=1;2;...;p· In this methodology, a potential energy is initially defined over each possible model for the relationship (1), what allows the application of the Lagrangian mechanics to the derived system. The solution of the Euler–Lagrange in this system allows obtaining the optimal solution according to the minimal action principle. The defined Lagrangian, corresponds to a continuous medium, where a n-dimensional finite elements model has been applied, so it is possible to get a solution for the problem solving a compatible and determined linear symmetric equation system. The computational implementation of the methodology has resulted in an improvement in the process of get representation models obtained and published previously by the authors.

Logic-Based Outer Approximation for the Design of Discrete-Continuous Dynamic Systems with Implicit Discontinuities

We address the optimization of discrete-continuous dynamic optimization problems using a disjunctive multistage modeling framework, with implicit discontinuities, which increases the problem complexity since the number of continuous phases and discrete events is not known a-priori. After setting a fixed alternative sequence of modes, we convert the infinite-dimensional continuous mixed-logic dynamic (MLDO) problem into a finite dimensional discretized GDP problem by orthogonal collocation on finite elements. We use the Logic-based Outer Approximation algorithm to fully exploit the structure of the GDP representation of the problem. This modelling framework is illustrated with an optimization problem with implicit discontinuities (diver problem).

Análisis numérico del comportamiento vibratorio del conjunto móvil de un altavoz dinámico

En este trabajo se describe el proceso llevado a cabo para analizar numéricamente mediante el método de los elementos finitos (MEF) el comportamiento vibratorio del conjunto móvil de un altavoz dinámico de bobina con doble suspensión inferior y en ausencia de la superior. El estudio se centra en el rango de baja frecuencia. El calibrado del modelo se realiza en base a medidas experimentales de la frecuencia de resonancia y desplazamiento del diafragma. Se hace énfasis en la importancia de los factores de participación asociados a la fuerza de excitación y en los cambios que se producen en estos al introducir fuerzas no equilibradas. Así mismo, el análisis proporciona datos para decidir sobre la distancia entre suspensiones, su número óptimo de pliegues y la ubicación de las trencillas, siendo estos parámetros de gran interés en el diseño de este tipo de altavoces.

Sistemas de protección individual anticaídas sometidos a impacto. Simulaciones numéricas

Los sistemas anticaída se establecen en UNE-EN-363 para evitar o retener caídas de personas sin que las fuerzas de retención causen daños relevantes. Constituyen el último recurso para prevención de caídas tras considerar otras medidas. Este artículo analiza las variables que afectan al comportamiento del sistema: altura de caída, longitud de cuerda de retención, rigidez y amortiguamiento del material de retención, y peso del operario. A tal fin, el fenómeno ha sido simulado con modelos dinámicos de elementos finitos mediante una cuerda elástica y un lastre rígido. Ello permite obtener valores verosímiles de las máximas fuerzas sobre el accidentado y el sistema durante el impacto. Los resultados demuestran el papel fundamental de la proporción entre altura de caída y longitud de cuerda (factor de caída), más determinante que la propia altura de caída. Las elevadas fuerzas máximas obtenidas indican la necesidad de incluir en el modelo absorbedores de energía elastoplásticos, en futuras investigaciones, con el fin de disminuir el impacto a valores razonables.

A systematic study on the influence of carbon on the behavior of hard-to-ionize elements in inductively coupled plasma–mass spectrometry

A systematic study on the influence of carbon on the signal of a large number of hard-to-ionize elements (i.e. B, Be, P, S, Zn, As, Se, Pd, Cd, Sb, I, Te, Os, Ir, Pt, Au, and Hg) in inductively coupled plasma–mass spectrometry has been carried out. To this end, carbon matrix effects have been evaluated considering different plasma parameters (i.e. nebulizer gas flow rate, r.f. power and sample uptake rate), sample introduction systems, concentration and type of carbon matrix (i.e. glycerol, citric acid, potassium citrate and ammonium carbonate) and type of mass spectrometer (i.e. quadrupole filter vs. double-focusing sector field mass spectrometer). Experimental results show that P, As, Se, Sb, Te, I, Au and Hg sensitivities are always higher for carbon-containing solutions than those obtained without carbon. The other hard-to-ionize elements (Be, B, S, Zn, Pd, Cd, Os, Ir and Pt) show no matrix effect, signal enhancement or signal suppression depending on the experimental conditions selected. The matrix effects caused by the presence of carbon are explained by changes in the plasma characteristics and the corresponding changes in ion distribution in the plasma (as reflected in the signal behavior plot, i.e. the signal intensity as a function of the nebulizer gas flow rate). However, the matrix effects for P, As, Se, Sb, Te, I, Au and Hg are also related to an increase in analyte ion population caused as a result of charge transfer reactions involving carbon-containing charged species in the plasma. The predominant specie is C+, but other species such as CO+, CO2+, C2+ and ArC+ could also play a role. Theoretical data suggest that B, Be, S, Pd, Cd, Os, Ir and Pt could also be involved in carbon based charge transfer reactions, but no experimental evidence substantiating this view has been found.

A finite element model of perforated panel absorbers including viscothermal effects

Most of the analytical models devoted to determine the acoustic properties of a rigid perforated panel consider the acoustic impedance of a single hole and then use the porosity to determine the impedance for the whole panel. However, in the case of not homogeneous hole distribution or more complex configurations this approach is no longer valid. This work explores some of these limitations and proposes a finite element methodology that implements the linearized Navier Stokes equations in the frequency domain to analyse the acoustic performance under normal incidence of perforated panel absorbers. Some preliminary results for a homogenous perforated panel show that the sound absorption coefficient derived from the Maa analytical model does not match those from the simulations. These differences are mainly attributed to the finite geometry effect and to the spatial distribution of the perforations for the numerical case. In order to confirm these statements, the acoustic field in the vicinities of the perforations is analysed for a more complex configuration of perforated panel. Additionally, experimental studies are carried out in an impedance tube for the same configuration and then compared to previous methods. The proposed methodology is shown to be in better agreement with the laboratorial measurements than the analytical approach.

New methodology for describing the equilibrium beach profile applied to the Valencia's beaches

Mathematical models used for the understanding of coastal seabed morphology play a key role in beach nourishment projects. These projects have become the fundamental strategy for coastal maintenance during the last few years. Accordingly, the accuracy of these models is vital to optimize the costs of coastal regeneration projects. Planning of such interventions requires methodologies that do not generate uncertainties in their interpretation. A study and comparison of mathematical simulation models of the coastline is carried out in this paper, as well as elements that are part of the model that are a source of uncertainty. The equilibrium profile (EP) and the offshore limit corresponding to the depth of closure (DoC) have been analyzed taking into account different timescale ranges. The results have thus been compared using data sets from three different periods which are identified as present, past and future. Accuracy in data collection for the beach profiles and the definition of the median grain size calculation using collected samples are the two main factors that have been taken into account in this paper. These data can generate high uncertainties and can produce a lack of accuracy in nourishment projects. Together they can generate excessive costs due to possible excess or shortage of sand used for the nourishment. The main goal of this paper is the development of a new methodology to increase the accuracy of the existing equilibrium beach profile models, providing an improvement to the inputs used in such models and in the fitting of the formulae used to obtain seabed shape. This new methodology has been applied and tested on Valencia's beaches.

Quantum theory of spin waves in finite chiral spin chains

We calculate the effect of spin waves on the properties of finite-size spin chains with a chiral spin ground state observed on biatomic Fe chains deposited on iridium(001). The system is described with a Heisenberg model supplemented with a Dzyaloshinskii-Moriya coupling and a uniaxial single ion anisotropy that presents a chiral spin ground state. Spin waves are studied using the Holstein-Primakoff boson representation of spin operators. Both the renormalized ground state and the elementary excitations are found by means of Bogoliubov transformation, as a function of the two variables that can be controlled experimentally, the applied magnetic field and the chain length. Three main results are found. First, because of the noncollinear nature of the classical ground state, there is a significant zero-point reduction of the ground-state magnetization of the spin spiral. Second, there is a critical external field from which the ground state changes from chiral spin ground state to collinear ferromagnetic order. The character of the two lowest-energy spin waves changes from edge modes to confined bulk modes over this critical field. Third, in the spin-spiral state, the spin-wave spectrum exhibits oscillatory behavior as function of the chain length with the same period of the spin helix.