Estudios de confort de vehículos ferroviarios mediante técnicas combinadas de simulación FEM y MBS

La elasticidad de la caja de viajeros de los vehículos ferroviarios tiene una gran influencia sobre el confort. Por esta razón, cuando se desea simular el comportamiento dinámico del vehículo para estudios de confort, resulta conveniente construir un modelo elástico de caja, a fin de obtener resultados más precisos. La construcción de este tipo de modelos pasa por el desarrollo de dos etapas fundamentales, que comprenden la generación de un modelo de caja mediante la técnica de los elementos finitos (FEM) para su caracterización dinámica como cuerpo elástico y la definición de un modelo de sistema multicuerpo (MBS) que englobe los restantes componentes del vehículo. En este artículo se presentan los resultados obtenidos en un estudio comparativo llevado a cabo con modelos de caja rígida y elástica, en los que se ha valorado el nivel de confort obtenido con ambas configuraciones. Para ello, se ha simulado el comportamiento del vehículo a dos velocidades distintas, de 70 y 110km/h, y con dos niveles de irregularidades. Se han analizado las aceleraciones de la caja, que se han procesado de acuerdo a las especificaciones de la norma EN12299, a fin de obtener el índice de comodidad. Este parámetro se ha utilizado para comparar el nivel de confort obtenido con ambos modelos, habiéndose encontrado una gran diferencia en los índices calculados con caja rígida y con caja elástica, lo que confirma la gran influencia de la elasticidad de la caja en los estudios de confort llevados a cabo mediante técnicas de simulación dinámica.

Development of FEM/BEM and SEA models from experimental results for structural elements with attached equipment

This work focuses on the analysis of a structural element of MetOP-A satellite. Given the special interest in the influence of equipment installed on structural elements, the paper studies one of the lateral faces on which the Advanced SCATterometer (ASCAT) is installed. The work is oriented towards the modal characterization of the specimen, describing the experimental set-up and the application of results to the development of a Finite Element Method (FEM) model to study the vibro-acoustic response. For the high frequency range, characterized by a high modal density, a Statistical Energy Analysis (SEA) model is considered, and the FEM model is used when modal density is low. The methodology for developing the SEA model and a compound FEM and Boundary Element Method (BEM) model to provide continuity in the medium frequency range is presented, as well as the necessary updating, characterization and coupling between models required to achieve numerical models that match experimental results.

A FEM-BEM coupling procedure through the Steklov-Poincarè operator

Many advantages can be got in combining finite and boundary elements.It is the case, for example, of unbounded field problems where boundary elements can provide the appropriate conditions to represent the infinite domain while finite elements are suitable for more complex properties in the near domain. However, in spite of it, other disadvantages can appear. It would be, for instance, the loss of symmetry in the finite elements stiffness matrix, when the combination is made. On the other hand, in our days, with the strong irruption of the parallel proccessing the techniques of decomposition of domains are getting the interest of numerous scientists. With their application it is possible to separate the resolution of a problem into several subproblems. That would be beneficial in the combinations BEM-FEM as the loss of symmetry would be avoided and every technique would be applicated separately. Evidently for the correct application of these techniques it is necessary to establish the suitable transmission conditions in the interface between BEM domain and FEM domain. In this paper, one parallel method is presented which is based in the interface operator of Steklov Poincarè.

A comparative analysis between FEM and Boundary Methods in bending of Plates

Two different methods of analysis of plate bending, FEM and BM are discussed in this paper. The plate behaviour is assumed to be represented by using the linear thin plate theory where the Poisson-Kirchoff assumption holds. The BM based in a weighted mean square error technique produced good results for the problem of plate bending. The computational effort demanded in the BM is smaller than the one needed in a FEM analysis for the same level of accuracy. The general application of the FEM cannot be matched by the BM. Particularly, different types of geometry (plates of arbitrary geometry) need a similar but not identical treatment in the BM. However, this loss of generality is counterbalanced by the computational efficiency gained in the BM in the solution achievement

Calculation Methodology and Fabrication Procedures for Particle Accelerator Strip-Line Kickers: Application to the CTF3 Combiner Ring Extraction Kicker and TL2 Tail Clippers

A particle accelerator is any device that, using electromagnetic fields, is able to communicate energy to charged particles (typically electrons or ionized atoms), accelerating and/or energizing them up to the required level for its purpose. The applications of particle accelerators are countless, beginning in a common TV CRT, passing through medical X-ray devices, and ending in large ion colliders utilized to find the smallest details of the matter. Among the other engineering applications, the ion implantation devices to obtain better semiconductors and materials of amazing properties are included. Materials supporting irradiation for future nuclear fusion plants are also benefited from particle accelerators. There are many devices in a particle accelerator required for its correct operation. The most important are the particle sources, the guiding, focalizing and correcting magnets, the radiofrequency accelerating cavities, the fast deflection devices, the beam diagnostic mechanisms and the particle detectors. Most of the fast particle deflection devices have been built historically by using copper coils and ferrite cores which could effectuate a relatively fast magnetic deflection, but needed large voltages and currents to counteract the high coil inductance in a response in the microseconds range. Various beam stability considerations and the new range of energies and sizes of present time accelerators and their rings require new devices featuring an improved wakefield behaviour and faster response (in the nanoseconds range). This can only be achieved by an electromagnetic deflection device based on a transmission line. The electromagnetic deflection device (strip-line kicker) produces a transverse displacement on the particle beam travelling close to the speed of light, in order to extract the particles to another experiment or to inject them into a different accelerator. The deflection is carried out by the means of two short, opposite phase pulses. The diversion of the particles is exerted by the integrated Lorentz force of the electromagnetic field travelling along the kicker. This Thesis deals with a detailed calculation, manufacturing and test methodology for strip-line kicker devices. The methodology is then applied to two real cases which are fully designed, built, tested and finally installed in the CTF3 accelerator facility at CERN (Geneva). Analytical and numerical calculations, both in 2D and 3D, are detailed starting from the basic specifications in order to obtain a conceptual design. Time domain and frequency domain calculations are developed in the process using different FDM and FEM codes. The following concepts among others are analyzed: scattering parameters, resonating high order modes, the wakefields, etc. Several contributions are presented in the calculation process dealing specifically with strip-line kicker devices fed by electromagnetic pulses. Materials and components typically used for the fabrication of these devices are analyzed in the manufacturing section. Mechanical supports and connexions of electrodes are also detailed, presenting some interesting contributions on these concepts. The electromagnetic and vacuum tests are then analyzed. These tests are required to ensure that the manufactured devices fulfil the specifications. Finally, and only from the analytical point of view, the strip-line kickers are studied together with a pulsed power supply based on solid state power switches (MOSFETs). The solid state technology applied to pulsed power supplies is introduced and several circuit topologies are modelled and simulated to obtain fast and good flat-top pulses.

Análisis del estado de tensiones en uniones carpinteras de empalme de llave por el método de los elementos finitos

El objeto del presente trabajo es la determinación de la distribución de tensiones en la unión carpintera de empalme de llave mediante el método de los elementos finitos y su comparación con los valores obtenidos mediante la teoría de Resistencia de Materiales. Se analizan las zonas donde se produce concentración de tensiones y se estudia la influencia del refino de la malla sobre los resultados con el fin de determinar el tamaño que proporciona los valores de tensión más ajustados a la teoría. En las zonas donde la concentración de tensiones es menor, distintos tamaños de la malla muestran valores de tensión similares. Los resultados muestran una simetría central de la distribución de isotensiones donde el centro de simetría se corresponde con el centro geométrico de la unión. En las zonas donde se produce concentración de tensiones, los valores de las mismas aumentan considerablemente con el refino de la malla. La comparación de los niveles de tensiones normales obtenidos por el MEF y la teoría clásica muestran diferencias reducidas, excepto en los puntos de concentración de tensiones. The purpose of this study is to determine the stress distribution in the joint carpenter of halved and tabled joint with the finite element method and its comparison with the values obtained using the theory of Strength of Materials. The stress concentration areas where analyzed and the influence of mesh refinement was studied on the results in order to determine the mesh size that provides the stress values more consistent with the theory. In areas where stress concentration is lower, different mesh sizes show similar stress values. The results show a central symmetry of the isobar lines distribution where the centre of symmetry corresponds to the geometric centre of the joint. In areas where stress concentration occurs, the same values increase considerably with the refinement of the mesh. Comparison of normal stress levels obtained by the FEM and the classical theory shows small differences, except at points of stress concentration.

Ethics and materials: Some Spanish Cases

The problem is general: modern architects and engineers are trying to understand historic structures using the wrong theoretical frame, the classic (elastic) thery of structures developed in the 19th Century for iron and stell, and in the 20th century for reinforced concrete, disguised with "modern" computer packages, mainly FEM, but also others. Masonry is an essentially different material, and the structural equations must be adapted accordingly. It is not a matter of "taste" or "opinion", and the consequences are before us. Since, say 1920s, historic monuments have suffered the aggression of generations of archietcts and engineers, trying to transform masonry in reinfored concrete or steel. The damage to the monuments and the expense has been, and is, enormous. However, as we have an adequate theory (modern limit analysis of masonry structures, Heyman 1966) which encompasses the "old theory" used successfully by the 18th and 19th Century practical engineers (from Perronet to Sejourné), it is a matter of "Ethics" not to use the wrong approach. It is also "contra natura" to modify the material masonry with indiscriminate injections, stitchings, etc. It is insane to consider, suddenly, that buildings which are Centuries or milennia old, are suddenly in danger of collapse. Maintenance is necessary but not the actual destruction of the constructive essence of the monument. A cocktail of "ignorance, fear and greed" is acting under the best of intentions.

Safety assessment of underground vehicles passing over highly resilient straight tracks in the presence of a broken rail

In this paper, vehicle-track interaction for a new slab track design, conceived to reduce noise and vibration levels has been analyzed, assessing the derailment risk for trains running on curved track when encountering a broken rail. Two different types of rail fastening systems with different elasticities have been analysed and compared. Numerical methods were used in order to simulate the dynamic behaviour of the train-track interaction. Multibody system (MBS) modelling techniques were combined with techniques based on the finite element method (FEM). MBS modelling was used for modelling the vehicle and FEM for simulating the elastic track. The simulation model was validated by comparing simulated results to experimental data obtained in field testing. During the simulations various safety indices, characteristic of derailment risk, were analysed. The simulations realised at the maximum running velocity of 110 km/h showed a similar behaviour for several track types. When reducing the running speed, the safety indices worsened for both cases. Although the worst behaviour was observed for the track with a greater elasticity, in none of the simulations did a derailment occur when running over the broken rail.

MicroRaman Spectroscopy of Si Nanowires: Influence of Size

Si Nanowires (NWs) were studied by Raman microspectroscopy. The Raman spectrum of the NWs reveals important thermal effects, which broaden and shift the one phonon Raman bands. The low thermal conductivity of the NWs and the low thermal dissipation are responsible for the temperature enhancement in the NW under the excitation with the laser beam. We have modeled, using finite element methods, the interaction between the laser beam and the NWs. The Raman spectrum of Si NWs is interpreted in terms of the temperature induced by the laser beam excitation, in correlation with finite element methods (fem) for studying the interaction between the laser beam and the NWs.

Coupled model of initiation and propagation of corrosion in reinforced concrete

Corrosion of reinforcing steel in concrete due to chloride ingress is one of the main causes of the deterioration of reinforced concrete structures. Structures most affected by such a corrosion are marine zone buildings and structures exposed to de-icing salts like highways and bridges. Such process is accompanied by an increase in volume of the corrosión products on the rebarsconcrete interface. Depending on the level of oxidation, iron can expand as much as six times its original volume. This increase in volume exerts tensile stresses in the surrounding concrete which result in cracking and spalling of the concrete cover if the concrete tensile strength is exceeded. The mechanism by which steel embedded in concrete corrodes in presence of chloride is the local breakdown of the passive layer formed in the highly alkaline condition of the concrete. It is assumed that corrosion initiates when a critical chloride content reaches the rebar surface. The mathematical formulation idealized the corrosion sequence as a two-stage process: an initiation stage, during which chloride ions penetrate to the reinforcing steel surface and depassivate it, and a propagation stage, in which active corrosion takes place until cracking of the concrete cover has occurred. The aim of this research is to develop computer tools to evaluate the duration of the service life of reinforced concrete structures, considering both the initiation and propagation periods. Such tools must offer a friendly interface to facilitate its use by the researchers even though their background is not in numerical simulation. For the evaluation of the initiation period different tools have been developed: Program TavProbabilidade: provides means to carry out a probability analysis of a chloride ingress model. Such a tool is necessary due to the lack of data and general uncertainties associated with the phenomenon of the chloride diffusion. It differs from the deterministic approach because it computes not just a chloride profile at a certain age, but a range of chloride profiles for each probability or occurrence. Program TavProbabilidade_Fiabilidade: carries out reliability analyses of the initiation period. It takes into account the critical value of the chloride concentration on the steel that causes breakdown of the passive layer and the beginning of the propagation stage. It differs from the deterministic analysis in that it does not predict if the corrosion is going to begin or not, but to quantifies the probability of corrosion initiation. Program TavDif_1D: was created to do a one dimension deterministic analysis of the chloride diffusion process by the finite element method (FEM) which numerically solves Fick’second Law. Despite of the different FEM solver already developed in one dimension, the decision to create a new code (TavDif_1D) was taken because of the need to have a solver with friendly interface for pre- and post-process according to the need of IETCC. An innovative tool was also developed with a systematic method devised to compare the ability of the different 1D models to predict the actual evolution of chloride ingress based on experimental measurements, and also to quantify the degree of agreement of the models with each others. For the evaluation of the entire service life of the structure: a computer program has been developed using finite elements method to do the coupling of both service life periods: initiation and propagation. The program for 2D (TavDif_2D) allows the complementary use of two external programs in a unique friendly interface: • GMSH - an finite element mesh generator and post-processing viewer • OOFEM – a finite element solver. This program (TavDif_2D) is responsible to decide in each time step when and where to start applying the boundary conditions of fracture mechanics module in function of the amount of chloride concentration and corrosion parameters (Icorr, etc). This program is also responsible to verify the presence and the degree of fracture in each element to send the Information of diffusion coefficient variation with the crack width. • GMSH - an finite element mesh generator and post-processing viewer • OOFEM – a finite element solver. The advantages of the FEM with the interface provided by the tool are: • the flexibility to input the data such as material property and boundary conditions as time dependent function. • the flexibility to predict the chloride concentration profile for different geometries. • the possibility to couple chloride diffusion (initiation stage) with chemical and mechanical behavior (propagation stage). The OOFEM code had to be modified to accept temperature, humidity and the time dependent values for the material properties, which is necessary to adequately describe the environmental variations. A 3-D simulation has been performed to simulate the behavior of the beam on both, action of the external load and the internal load caused by the corrosion products, using elements of imbedded fracture in order to plot the curve of the deflection of the central region of the beam versus the external load to compare with the experimental data.

A simulation method to estimate closing forces in car-sealing rubber elements

The door-closing process can reinforce the impression of a solid, rock-proof, car body or of a rather cheap, flimsy vehicle. As there are no real prototypes during rubber profile bidding-out stages, engineers need to carry out non-linear numerical simulations that involve complex phenomena as well as static and dynamic loads for several profile candidates. This paper presents a structured virtual design tool based on FEM, including constitutive laws and incompressibility constraints allowing to predict more realistically the final closing forces and even to estimate sealing overpressure as an additional guarantee of noise insulation. Comparisons with results of physical tests are performed.

Estructuras de hormigón armado bajo carga dinámica severa. Parte III: ejemplo de simulación numérica

Publicamos la tercera parte y última de este artículo sobre estructuras de hormigón armado bajo carga dinámica severa en el que, los autores, nos han ido exponiendo diferentes ejemplos. En esta ocasión, se pone patente que con el desarrollo de los ordenadores durante las últimas décadas se ha dado la posibilidad de utilizar el método de los elementos finitos (FEM: Finite Elements Method) para estudiar los efectos de una carga dinámica severa, como en el caso de las explosiones. A continuación se describe el uso del método para este tipo de situaciones de carga.

3D hp-adaptive finite element simulations of a magic-T electromagnetic waveguide structure

This paper employs a 3D hp self-adaptive grid-refinement finite element strategy for the solution of a particular electromagnetic waveguide structure known as Magic-T. This structure is utilized as a power divider/combiner in communication systems as well as in other applications. It often incorporates dielectrics, metallic screws, round corners, and so on, which may facilitate its construction or improve its design, but significantly difficult its modeling when employing semi-analytical techniques. The hp-adaptive finite element method enables accurate modeling of a Magic-T structure even in the presence of these undesired materials/geometries. Numerical results demonstrate the suitability of the hp-adaptive method for modeling a Magic-T rectangular waveguide structure, delivering errors below 0.5% with a limited number of unknowns. Solutions of waveguide problems delivered by the self-adaptive hp-FEM are comparable to those obtained with semi-analytical techniques such as the Mode Matching method, for problems where the latest methods can be applied. At the same time, the hp-adaptive FEM enables accurate modeling of more complex waveguide structures.

A three-dimensional self-adaptive hp finite element method for the characterization of waveguide discontinuities

We propose the use of a highly-accurate three-dimensional (3D) fully automatic hp-adaptive finite element method (FEM) for the characterization of rectangular waveguide discontinuities. These discontinuities are either the unavoidable result of mechanical/electrical transitions or deliberately introduced in order to perform certain electrical functions in modern communication systems. The proposed numerical method combines the geometrical flexibility of finite elements with an accuracy that is often superior to that provided by semi-analytical methods. It supports anisotropic refinements on irregular meshes with hanging nodes, and isoparametric elements. It makes use of hexahedral elements compatible with high-order H(curl)H(curl) discretizations. The 3D hp-adaptive FEM is applied for the first time to solve a wide range of 3D waveguide discontinuity problems of microwave communication systems in which exponential convergence of the error is observed.

Analysis of the stress state of a halved and tabled traditional timber scarf joint with the Finite Element Method

The purpose of this study is to determine the stress distribution in the carpentry joint of halved and tabled scarf joint with the finite element method (FEM) and its comparison with the values obtained using the theory of Strength of Materials. The stress concentration areas where analyzed and the influence of mesh refinement was studied on the results in order to determine the mesh size that provides the stress values more consistent with the theory. In areas where stress concentration is lower, different mesh sizes show similar stress values. In areas where stress concentration occurs, the same values increase considerably with the refinement of the mesh. The results show a central symmetry of the isobar lines distribution where the centre of symmetry corresponds to the geometric centre of the joint. Comparison of normal stress levels obtained by the FEM and the classical theory shows small differences, except at points of stress concentration.