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

Analysis of plate bending by means of high order finite hyperelements

The possibilities and limitations of high order hyperelements in plate bending analysis are discussed. Explicit shape functions for some types of triangular elements are given. These elements are applied to simple cases to assess their computational efficiency.

Application of the dynamic models to the detection of imperfections in plate Structures

This paper is part of a set of publications related with the development of mathematical models aimed to simulate the dynamic input and output of experimental nondestructive tests in order to detect structural imperfections. The structures to be considered are composed by steel plates of thin thickness. The imperfections in these cases are cracks and they can penetrate either a significant part of the plate thickness or be micro cracks or superficial imperfections. The first class of cracks is related with structural safety and the second one is more connected to the structural protection to the environment, particularly if protective paintings can be deteriorated. Two mathematical groups of models have been developed. The first group tries to locate the position and extension of the imperfection of the first class of imperfections, i.e. cracks and it is the object of the present paper. Bending Kirchoff thin plate models belong to this first group and they are used to this respect. The another group of models is dealt with membrane structures under the superficial Rayleigh waves excitation. With this group of models the micro cracks detection is intended. In the application of the first group of models to the detection of cracks, it has been observed that the differences between the natural frequencies of the non cracked and the cracked structures are very small. However, geometry and crack position can be identified quite accurately if this comparison is carried out between first derivatives (mode rotations) of the natural modes are used instead. Finally, in relation with the analysis of the superficial crack existence the use of Rayleigh waves is very promising. The geometry and the penetration of the micro crack can be detected very accurately. The mathematical and numerical treatment of the generation of these Rayleigh waves present and a numerical application has been shown.

A C1 finite element for Kirchhoff plate bending

After a short introduction the possibilities and limitations of polynomial simple elements with C1 continuity are discussed with reference to plate bending analysis. A family of this kind of elements is presented.. These elements are applied to simple cases in order to assess their computational efficiency. Finally some conclusions are shown, and future research is also proposed.

Structural Analysis of Thin Tile Vaults and Domes: The Inner Oval Dome of the Basilica de los Desamparados in Valencia

The inner oval dome of the Basílica de la Virgen los Desamparados, built in 1701, is one of the most slender masonry vaults ever built. It is a tile dome with a total thickness of 80 mm and a main span of 18.50 m. It was built without centering with great ingenuity and economy of means, thirty three years after the termination of the building in 1667. The dome is in contact with the external dome only in the inferior part with the projecting ribs of the intrados, the lunettes of the windows, and, in the upper part, through 126 inclined iron bars. This unique construction was revealed in the 1990's in the studies previous to the restoration of the Basílica, and has given rise to different theories about the mode of construction and the structural behaviour and safety of the dome. The present contribution aims to provide a plausible hypothesis about the mode of construction and to explain the safety of the inner dome which has stood, without need of repairs or reinforcement, for 300 hundred years.

A computer formulation for the analysis on continuous non prismatic folded plate structures of arbitrary cross-section

A computer solution to analyze nonprismatic folded plate structures is shown. Arbitrary cross-sections (simple and multiple), continuity over intermediate supports and general loading and longitudinal boundary conditions are dealt with. The folded plates are assumed to be straight and long (beam like structures) and some simplifications are introduced in order to reduce the computational effort. The formulation here presented may be very suitable to be used in the bridge deck analysis.

Vibro-acoustic analysis of spacecraft structures with thin air layers

Esta Tesis presenta un estudio sobre el comportamiento vibroacústico de estructuras espaciales que incluyen capas de aire delgadas, así como sobre su modelización numérica. Las capas de aire pueden constituir un elemento fundamental en estos sistemas, como paneles solares plegados, que se consideran el caso de estudio en este trabajo. Para evaluar la influencia de las capas de aire en la respuesta dinámica del sistema se presenta el uso de modelos unidimensionales. La modelización de estos sistemas se estudia para los rangos de baja y alta frecuencia. En el rango de baja frecuencia se propone un conjunto de estrategias de simulación basadas en técnicas numéricas que se utilizan habitualmente en la industria aeroespacial para facilitar la aplicación de los resultados de la Tesis en los modelos numéricos actuales. Los resultados muestran el importante papel de las capas de aire en la respuesta del sistema. El uso de modelos basados en elementos finitos o de contorno para estos elementos proporciona resultados equivalentes aunque la aplicabilidad de estos últimos puede estar condicionada por la geometría del problema. Se estudia asimismo el uso del Análisis Estadístico de la Energía (SEA) para estos elementos. Una de las estrategias de simulación propuestas, que incluye una formulación energética para el aire que rodea a la estructura, se propone como estimador preliminar de la respuesta del sistema y sus frecuencias propias. Para el rango de alta frecuencia, se estudia la influencia de la definición del propio modelo SEA. Se presenta el uso de técnicas de reducción para determinar una matriz de pérdidas SEA reducida para definiciones incompletas del sistema (si algún elemento que interactúa con el resto no se incluye en el modelo). Esta nueva matriz tiene en cuenta la contribución de las subestructuras que no se consideran parte del modelo y que suelen ignorarse en el procedimiento habitual para reducir el tamaño del mismo. Esta matriz permite también analizar sistemas que incluyen algún componente con problemas de accesibilidad para medir su respuesta. Respecto a la determinación de los factores de pérdidas del sistema, se presenta una metodología que permite abordar casos en los que el método usual, el Método de Inyección de Potencia (PIM), no puede usarse. Se presenta un conjunto de métodos basados en la técnicas de optimización y de actualización de modelos para casos en los que no se puede medir la respuesta de todos los elementos del sistema y también para casos en los que no todos los elementos pueden ser excitados, abarcando un conjunto de casos más amplio que el abordable con el PIM. Para ambos rangos de frecuencia se presentan diferentes casos de análisis: modelos numéricos para validar los métodos propuestos y un panel solar plegado como caso experimental que pone de manifiesto la aplicación práctica de los métodos presentados en la Tesis. ABSTRACT This Thesis presents an study on the vibro-acoustic behaviour of spacecraft structures with thin air layers and their numerical modelling. The air layers can play a key role in these systems as solar wings in folded configuration that constitute the study case for this Thesis. A method based on one-dimensional models is presented to assess the influence of the air layers in the dynamic response of the system. The modelling of such systems is studied for low and high frequency ranges. In the low frequency range a set of modelling strategies are proposed based on numerical techniques used in the industry to facilitate the application of the results in the current numerical models. Results show the active role of the air layers in the system response and their great level of influence. The modelling of these elements by means of Finite Elements (FE) and Boundary Elements (BE) provide equivalent results although the applicability of BE models can be conditioned by the geometry of the problem. The use of Statistical Energy Analysis (SEA) for these systems is also presented. Good results on the system response are found for models involving SEA beyond the usual applicability limit. A simulation strategy, involving energetic formulation for the surrounding fluid is proposed as fast preliminary approach for the system response and the coupled eigenfrequencies. For the high frequency range, the influence of the definition of the SEA model is presented. Reduction techniques are used to determine a Reduced SEA Loss Matrix if the system definition is not complete and some elements, which interact with the rest, are not included. This new matrix takes into account the contribution of the subsystems not considered that are neglected in the usual approach for decreasing the size of the model. It also allows the analysis of systems with accessibility restrictions on some element in order to measure its response. Regarding the determination of the loss factors of a system, a methodology is presented for cases in which the usual Power Injection Method (PIM) can not be applied. A set of methods are presented for cases in which not all the subsystem responses can be measured or not all the subsystems can be excited, as solar wings in folded configuration. These methods, based on error minimising and model updating techniques can be used to calculate the system loss factors in a set of cases wider than the PIM’s. For both frequency ranges, different test problems are analysed: Numerical models are studied to validate the methods proposed; an experimental case consisting in an actual solar wing is studied on both frequency ranges to highlight the industrial application of the new methods presented in the Thesis.

A plastic analysis of nonprismatic folded plate structures

A computer method for the plastic analysis of folded plate structures is presented. The method considers the specific characteristics of the folded plate structural model using a simplified one-dimensional theory. and it can be applied to the analysis of any type of folded pIates, either prismatic or nonprismatic, with arbitrary cross-section. A simple example is analyzed in order to show the possibilities of the procedure and some results of interest are presented

The folded plate model in structural analysis. Elastic, plastic and dynamic formulations for nonprismatic structures

A specific numerical procedure for the analysis of arbitrary nonprismatic folded plate structures is presented. An elastic model is studied and compared with a harmonic solution for a prismatic structure. An extension to the plastic analysis is developed, and the influence of the structural geometry and loading pattern is analyzed. Nonprismatic practical cases, with arbitrary geometry and loading are shown, as well in the elastic range as in the plastic one. Finally, a dynamic formulation is outlined

Analysis by finite element calculations of light scattering in laser-textured AZO films for PV thin-film solar cells

In the thin-film photovoltaic industry, to achieve a high light scattering in one or more of the cell interfaces is one of the strategies that allow an enhancement of light absorption inside the cell and, therefore, a better device behavior and efficiency. Although chemical etching is the standard method to texture surfaces for that scattering improvement, laser light has shown as a new way for texturizing different materials, maintaining a good control of the final topography with a unique, clean, and quite precise process. In this work AZO films with different texture parameters are fabricated. The typical parameters used to characterize them, as the root mean square roughness or the haze factor, are discussed and, for deeper understanding of the scattering mechanisms, the light behavior in the films is simulated using a finite element method code. This method gives information about the light intensity in each point of the system, allowing the precise characterization of the scattering behavior near the film surface, and it can be used as well to calculate a simulated haze factor that can be compared with experimental measurements. A discussion of the validation of the numerical code, based in a comprehensive comparison with experimental data is included.

Performance Analysis of 10,000 Residential PV Systems in France and Belgium

The main objective of this paper is to review the state of the art of residential PV systems in France and Belgium. This is done analyzing the operational data of 10650 PV systems (9657 located in France and 993 in Belgium). Three main questions are posed. How much energy do they produce? What level of performance is associated to their production? Which are the key parameters that most influence their quality? During the year 2010, the PV systems in France have produced a mean annual energy of 1163 kWh/kWp in France and 852 kWh/kWp in Belgium. As a whole, the orientation of PV generators causes energy productions to be some 7% inferior to optimally oriented PV systems. The mean Performance Ratio is 76% in France and 78% in Belgium, and the mean Performance Index is 85% in both countries. On average, the real power of the PV modules falls 4.9% below its corresponding nominal power announced on the manufacturer?s datasheet. A brief analysis by PV modules technology has lead to relevant observations about two technologies in particular. On the one hand, the PV systems equipped with Heterojunction with Intrinsic. Thin layer (HIT) modules show performances higher than average. On the other hand, the systems equipped with Copper Indium (di)Selenide (CIS) modules show a real power that is 16 % lower than their nominal value.

The PSE-3D instability analysis methodology for flows depending strongly on two and weakly on the third spatial dimension

The present contribution discusses the development of a PSE-3D instability analysis algorithm, in which a matrix forming and storing approach is followed. Alternatively to the typically used in stability calculations spectral methods, new stable high-order finitedifference-based numerical schemes for spatial discretization 1 are employed. Attention is paid to the issue of efficiency, which is critical for the success of the overall algorithm. To this end, use is made of a parallelizable sparse matrix linear algebra package which takes advantage of the sparsity offered by the finite-difference scheme and, as expected, is shown to perform substantially more efficiently than when spectral collocation methods are used. The building blocks of the algorithm have been implemented and extensively validated, focusing on classic PSE analysis of instability on the flow-plate boundary layer, temporal and spatial BiGlobal EVP solutions (the latter necessary for the initialization of the PSE-3D), as well as standard PSE in a cylindrical coordinates using the nonparallel Batchelor vortex basic flow model, such that comparisons between PSE and PSE-3D be possible; excellent agreement is shown in all aforementioned comparisons. Finally, the linear PSE-3D instability analysis is applied to a fully three-dimensional flow composed of a counter-rotating pair of nonparallel Batchelor vortices.

Wedges and plate-bandes: mechanical theories after De la Hire

Plate-bandes are straight masonry arches (they are called, also, flat arches or lintel arches). Ideally they have the surfaces of extrados and intrados plane and horizontal. The stones or bricks have radial joints converging usually in one centre. The voussoirs have the form of wedges and in French they are called "claveaux". A plate-bande is, in fact, a lintel made of several stones and the proportions of lintels and plate-bandes are similar. Proportions of plate-bandes, that is the relationship between the thickness t and the span s (t/s)varies, typically between 1/4–1/3 in thick plate-bandes, and is less than 1/20 in the most slender ones. A ratio of circa 1/8 was usual in the 18th Century and follows a simple geometrical rule: the centre form with the intrados an equilateral triangle and the plate-bande should contain an arc of circle. The joints are usually plane, but in some cases present a «rebated» or «stepped» form. Plate-bandes exert an inclined thrust as any masonry arch. This thrust is usually very high and it requires either massive buttresses, or to be built in the middle of thick walls. Master builders and architects have tried since antiquity to calculate the abutment necessary for any arch. A modern architect or engineer will measure the arch thrust in units of force, kN or tons. Traditionally, the thrust has been measured as the size of the buttresses to resist it safely. Old structural rules, then, addressed the design problem establishing a relationship between the span and the depth of the buttress. These were empirical rules, particular for every type of arch or structure in every epoch. Thus, the typical gothic buttress is 1/4 of the vault span, but a Renaissance or baroque barrel vault will need more than 1/3 of the span. A plate-bande would require more than one half of the span; this is precisely the rule cited by the French engineer Gautier, who tried unsuccessfully to justify it by static reasons. They were used, typically, to form the lintels of windows or doors (1-2 m, typically); in Antiquity they were used, also, though rarely, at the gates of city walls or in niches (ca. 2 m, reaching 5.2 m). Plate-bandes may show particular problems: it is not unusual that some sliding of the voussoirs can be observed, particularly in thick plate-bandes. The stepped joints on Fig. 1, left, were used to avoid this problem. There are other «hidden» methods, like iron cramps or the use of stone wedges, etc. In seismic zones these devices were usual. Another problem relates to the deformation; a slight yielding of the abutments, or even the compression of the mortar joints, may lead to some cracking and the descent of the central keystone. Even a tiny descent will convert the original straight line of the intrados in a broken line with a visible «kink» or angle in the middle. Of course, both problems should be avoided. Finally, the wedge form of the voussoirs lead to acute angles in the stones and this can produce partial fractures; this occurs usually at the inferior border of the springers at the abutments. It follows, that to build a successful plate-bande is not an easy matter. Also, the structural study of plate-bandes is far from simple, and mechanics and geometry are related in a particular way. In the present paper we will concentrate on the structural aspects and their constructive consequences, with a historical approach. We will outline the development of structural analysis of plate-bandes from ca. 1700 until today. This brief history has a more than purely academic interest. Different approaches and theories pointed to particular problem, and though the solution given may have been incorrect, the question posed was often pertinent. The paper ends with the application of modern Limit Analysis of Masonry Structures, developed mainly by professor Heyman in the last fifty years. The work aims, also, to give some clues for the actual architect and engineer involved in the analysis or restoration of masonry buildings.

Electromagnetic analysis of periodic structures under a waveguide viewpoint

The fundamental objective of this Ph. D. dissertation is to demonstrate that, under particular circumstances which cover most of the structures with practical interest, periodic structures can be understood and analyzed by means of closed waveguide theories and techniques. To that aim, in the first place a transversely periodic cylindrical structure is considered and the wave equation, under a combination of perfectly conducting and periodic boundary conditions, is studied. This theoretical study runs parallel to the classic analysis of perfectly conducting closed waveguides. Under the light shed by the aforementioned study it is clear that, under certain very common periodicity conditions, transversely periodic cylindrical structures share a lot of properties with closed waveguides. Particularly, they can be characterized by a complete set of TEM, TE and TM modes. As a result, this Ph. D. dissertation introduces the transversely periodic waveguide concept. Once the analogies between the modes of a transversely periodic waveguide and the ones of a closed waveguide have been established, a generalization of a well-known closed waveguide characterization method, the generalized Transverse Resonance Technique, is developed for the obtention of transversely periodic modes. At this point, all the necessary elements for the consideration of discontinuities between two different transversely periodic waveguides are at our disposal. The analysis of this type of discontinuities will be carried out by means of another well known closed waveguide method, the Mode Matching technique. This Ph. D. dissertation contains a sufficient number of examples, including the analysis of a wire-medium slab, a cross-shaped patches periodic surface and a parallel plate waveguide with a textured surface, that demonstrate that the Transverse Resonance Technique - Mode Matching hybrid is highly precise, efficient and versatile. Thus, the initial statement: ”periodic structures can be understood and analyzed by means of closed waveguide theories and techniques”, will be corroborated. Finally, this Ph. D. dissertation contains an adaptation of the aforementioned generalized Transverse Resonance Technique by means of which the analysis of laterally open periodic waveguides, such as the well known Substrate Integrated Waveguides, can be carried out without any approximation. The analysis of this type of structures has suscitated a lot of interest in the recent past and the previous analysis techniques proposed always resorted to some kind of fictitious wall to close the structure. vii Resumen El principal objetivo de esta tesis doctoral es demostrar que, bajo ciertas circunstancias que se cumplen para la gran mayoría de estructuras con interés práctico, las estructuras periódicas se pueden analizar y entender con conceptos y técnicas propias de las guías de onda cerradas. Para ello, en un primer lugar se considera una estructura cilíndrical transversalmente periódica y se estudia la ecuación de onda bajo una combinación de condiciones de contorno periódicas y de conductor perfecto. Este estudio teórico y de caracter general, sigue el análisis clásico de las guías de onda cerradas por conductor eléctrico perfecto. A la luz de los resultados queda claro que, bajo ciertas condiciones de periodicidad (muy comunes en la práctica) las estructuras cilíndricas transversalmente periódicas guardan multitud de analogías con las guías de onda cerradas. En particular, pueden ser descritas mediante un conjunto completo de modos TEM, TE y TM. Por ello, ésta tesis introduce el concepto de guía de onda transversalmente periódica. Una vez establecidas las similitudes entre las soluciones de la ecuación de onda, bajo una combinación de condiciones de contorno periódicas y de conductor perfecto, y los modos de guías de onda cerradas, se lleva a cabo, con éxito, la adaptación de un conocido método de caracterización de guías de onda cerradas, la técnica de la Resonancia Transversal Generalizada, para la obtención de los modos de guías transversalmente periódicas. En este punto, se tienen todos los elementos necesarios para considerar discontinuidades entre guías de onda transversalmente periódicas. El analisis de este tipo de discontinuidades se llevará a cabo mediante otro conocido método de análisis de estructuras cerradas, el Ajuste Modal. Esta tesis muestra multitud de ejemplos, como por ejemplo el análisis de un wire-medium slab, una superficie de parches con forma de cruz o una guía de placas paralelas donde una de dichas placas tiene cierta textura, en los que se demuestra que el método híbrido formado por la Resonancia Transversal Generalizada y el Ajuste Modal, es tremendamente preciso, eficiente y versátil y confirmará la validez de el enunciado inicial: ”las estructuras periódicas se pueden analizar y entender con conceptos y técnicas propias de las guías de onda cerradas” Para terminar, esta tésis doctoral incluye también una modificación de la técnica de la Resonancia Transversal Generalizada mediante la cual es posible abordar el análisis de estructuras periódica abiertas en los laterales, como por ejemplo las famosas guías de onda integradas en sustrato, sin ninguna aproximación. El análisis de este tipo de estructuras ha despertado mucho interés en los últimos años y las técnicas de análisis propuestas hasta ix el momento acostumbran a recurrir a algún tipo de pared ficticia para simular el carácter abierto de la estructura.

Selective ablation with UV lasers of a-Si:H thin film solar cells in direct scribing configuration

Monolithical series connection of silicon thin-film solar cells modules performed by laser scribing plays a very important role in the entire production of these devices. In the current laser process interconnection the two last steps are developed for a configuration of modules where the glass is essential as transparent substrate. In addition, the change of wavelength in the employed laser sources is sometimes enforced due to the nature of the different materials of the multilayer structure which make up the device. The aim of this work is to characterize the laser patterning involved in the monolithic interconnection process in a different configurations of processing than the usually performed with visible laser sources. To carry out this study, we use nanosecond and picosecond laser sources working at 355nm of wavelength in order to achieve the selective ablation of the material from the film side. To assess this selective removal of material has been used EDX (energy dispersive using X-ray) analysis