Search

The optimal design of a vertical cantilever beam is presented in this paper. The beam is assumed immersed in an elastic Winkler soil and subjected to several loads: a point force at the tip section, its self weight and a uniform distributed load along its length. lbe optimal design problem is to find the beam of a given length and minimum volume, such that the resultant compressive stresses are admisible. This prohlem is analyzed according to linear elasticity theory and within different alternative structural models: column, Navier-Bernoulli beam-column, Timoshenko beamcolumn (i.e. with shear strain) under conservative loads, typically, constant direction loads. Results obtained in each case are compared, in order to evaluate the sensitivity of model on the numerical results. The beam optimal design is described by the section distribution layout (area, second moment, shear area etc.) along the beam span and the corresponding beam total volume. Other situations, some of them very interesting from a theoretical point of view, with follower loads (Beck and Leipholz problems) are also discussed, leaving for future work numerical details and results.

This paper is a preliminary version of Chapter 3 of a State-of-the-Art Report by the IASS Working Group 5: Concrete Shell Roofs. The intention of this chapter is to set forth for those who intend to design concrete shell roofs information and advice about the selection, verification and utilization of commercial computer tools for analysis and design tasks.The computer analysis and design steps for a concrete shell roof are described. Advice follows on the aspects to be considered in the application of commercial finite element (FE)computer programs to concrete shell analysis, starting with recommendations on how novices can gain confidence and competence in the use of software. To establish vocabulary and provide background references, brief surveys are presented of, first,element types and formulations for shells and, second, challenges presented by advanced analyses of shells. The final section of the chapter indicates what capabilities to seek in selecting commercial FE software for the analysis and design of concrete shell roofs. Brief concluding remarks summarize advice regarding judicious use of computer analysis in design practice.

The theory of guided wave propagation in N-layered media is presented. The derivation of the dispersion equation is obtained from the application of appropriate boundary and continuity conditions to the solution of the general wave propagation. The resulting dispersion equation is given in the form of a determinant of a 4Nx4N coefficient matrix. A numerical procedure is proposed to represent and solve the implicit equation resulting. The validity and efficiency of the proposed numerical model is discussed. Dispersion curves characterizing the N-layered material are obtained and compared to published results. In order to illustrate the use of the model to different practical applications, such as coating problems of plasma spray on a turbine blade, aircraft multiple layers, ice detection, etc. a sensitivity analysis of the dispersion curves respect to depth imperfections is given.

The principal goal in the design of any general engineering structure is to achieve stability. A structure must remain stable and therefore, fit for purpose, throughout its useful life. This requirement can be achieved with the help of structural analysis techniques, which nowadays can be performed using computers. In this paper, the authors consider the three distinct phases of structural analysis, namely Idealisation, Analysis and Interpretation. The idealisation phase describes the process of translating the proposed engineering works into a model that represents the physical response of the structure to applied loading. This model is a simplification of the real structure and only considers the characteristics of the components that have the greatest influence on its stability. In the following phase, the analysis phase, the aforementioned model is used to study and understand how the structure might behave in reality. The framing of this model may draw upon knowledge obtained in many other disciplines, including Applied Mathematics and Materials Science. The significance of the results are determined during the interpretation phase and the engineer either decides that stability can be achieved or works through the process again, modifying the initial idealisation. As computers continue to grow in capability, it is the purpose of this paper to consider the role of computational methods and how they may empower present and future structural engineers.

The optimal design of a column subjected to compression forces is studied. Different support conditions are considered for the column. Also, a point load at the top section, self weight and unifonnly distributed vertical loading may be included in the analysis. The beam is assumed immersed in an elastic Winkler-type soil. Conservative loading has been assumed, i.e. with constant direction and therefore follower loads situation has been disregarded. The optimal column design is described by both cross-section distribution (area; second moment, shear area etc.) along the column span and the corresponding total volume. The optimal design problem consists of detennining the cross-section distribution layout for a beam of given length, so that the column volume becomes minimum, and resultant compressive stresses are admisible. This problem is analyzed assuming a linear elastic material behavior. The column is analyzed by comparing the results obtained from the Navier-Bemoulli beam-column model, and from Timoshenko beam- column model. Special shape functions are introduced in each case.

In this article, a model for the determination of displacements, deformations and tensions of a submarine pipeline during the construction is presented. The process is carried out from an initial floating situation to the final laying position on the seabed. The existence of currents and small waves are also considered. Firstly, this technique, usually applied to polyethylene pipelines, is described in this paper as well as some real world examples, as well as the variables that can be modified to control the behavior of the structure. A detailed description of the actions in this process is considered, specially the ones related to marine environment, as Archimedes force, current and sea waves. The behavior of the pipeline is modeled with a non linear elasto dynamic model where geometric non linearities are taken into account. A 3-D beam model, without cross section deformation effects, is developed. Special care is taken in the numerical analysis, developed within an updated lagrangian formulation framework, with the sea bed contact, the follower forces due to the external water pressures and the dynamic actions. Finally, some subroutines are implemented into ANSYS to simulate the two dimensional case, where the whole construction process is achieved. With this software, a sensibility analysis of the bending moments, axial forces and stresses obtained with different values of the control variables in order to optimize the construction steps. These control variables are, the axial load in the pipe, the inundated inner length and the distance of the control barge from the coast.

A Mindlin plate with periodically distributed ribs patterns is analyzed by using homogenization techniques based on asymptotic expansion methods. The stiffness matrix of the homogenized plate is found to be dependent on the geometrical characteristics of the periodical cell, i.e. its skewness, plan shape, thickness variation etc. and on the plate material elastic constants. The computation of this plate stiffness matrix is carried out by averaging over the cell domain some solutions of different periodical boundary value problems. These boundary value problems are defined in variational form by linear first order differential operators on the cell domain and the boundary conditions of the variational equation correspond to a periodic structural problem. The elements of the stiffness matrix of homogenized plate are obtained by linear combinations of the averaged solution functions of the above mentioned boundary value problems. Finally, an illustrative example of application of this homogenization technique to hollowed plates and plate structures with ribs patterns regularly arranged over its area is shown. The possibility of using in the profesional practice the present procedure to the actual analysis of floors of typical buildings is also emphasized.

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.

In a Finite Element (FE) analysis of elastic solids several items are usually considered, namely, type and shape of the elements, number of nodes per element, node positions, FE mesh, total number of degrees of freedom (dot) among others. In this paper a method to improve a given FE mesh used for a particular analysis is described. For the improvement criterion different objective functions have been chosen (Total potential energy and Average quadratic error) and the number of nodes and dof's of the new mesh remain constant and equal to the initial FE mesh. In order to find the mesh producing the minimum of the selected objective function the steepest descent gradient technique has been applied as optimization algorithm. However this efficient technique has the drawback that demands a large computation power. Extensive application of this methodology to different 2-D elasticity problems leads to the conclusion that isometric isostatic meshes (ii-meshes) produce better results than the standard reasonably initial regular meshes used in practice. This conclusion seems to be independent on the objective function used for comparison. These ii-meshes are obtained by placing FE nodes along the isostatic lines, i.e. curves tangent at each point to the principal direction lines of the elastic problem to be solved and they should be regularly spaced in order to build regular elements. That means ii-meshes are usually obtained by iteration, i.e. with the initial FE mesh the elastic analysis is carried out. By using the obtained results of this analysis the net of isostatic lines can be drawn and in a first trial an ii-mesh can be built. This first ii-mesh can be improved, if it necessary, by analyzing again the problem and generate after the FE analysis the new and improved ii-mesh. Typically, after two first tentative ii-meshes it is sufficient to produce good FE results from the elastic analysis. Several example of this procedure are presented.

In this work we introduce a new mathematical model for elastic beams with a cross-section of constant width and longitudinally variable depth, obtained from the classical three-dimensional linear elasticity problem by using an asymptotic expansion method. We characterize the first- and second-order displacements and the first-order stress field, giving results related to existence, uniqueness and convergence for the limit model solution. Finally, we present the computations for a classical example.

El libro Apología de un matemático de G. H. Hardy (1877-1947), uno de los mejores matemáticos de este siglo, ha vuelto a ser publicado recientemente en español, dentro del año de la Matemática que estamos celebrando (…)

Examples of global solutions of the shell equations are presented, such as the ones based on the well known Levy series expansion. Also discussed are some natural extensions of the Levy method as well as the inherent limitations of these methods concerning the shell model assumptions, boundary conditions and geometric regularity. Finally, some open additional design questions are noted mainly related to the simultaneous use in analysis of these global techniques and the local methods (like the finite elements) to finding the optimal shell shape, and to determining the reinforcement layout.

The option value problem with two costs is written as a variational inequality. The advantage of this formulation is that it takes place in a fixed domain. Thus no front tracking is needed for numerical approximation of the free boundary. An iterative algorithm is proposed which can be used to solve the nonlinear system obtained by finite differences or finite elements procedures. Especial care has to be taken in the design of differences finites schemes o finite elements due to the degeneracy of the differential operator. These schemes can be absortion or convection dominated nearly to the axis. This is a preliminary note to the study of this kind of problems.

Se plantea el problema de conseguir un método que, partiendo de una malla inicial en elementos finitos, genere una malla próxima a la óptima, conservando el mismo número de grados de libertad y, por consiguiente, sin penalizar los tiempos de análisis por ordenador. La bondad de una malla se mide por un funcional determinado (Energía Potencial Total, Error Cuadrático Medio, etc.) La Técnica de gradiente descendente, aplicada al funcional de Energía Potencial Total (1) permite, a partir de una inicial dada, obtener una malla mejorada. No obstante, este método requiere un esfuerzo de computación muy grande. Como consecuencia de la aplicación del método del gradiente descendente a numerosos casos, se ha observado que la geometría que adopta la malla mejorada, se aproxima a la de una malla tal que sus nudos se sitúen sobre las líneas isostáticas (envolventes de las tensiones principales) y generen elementos regulares (de lados iguales) o cuasirregulares. La conclusión más importante es, precisamente, que a partir de una malla inicial, razonablemente regular, se puede realizar un único análisis y definir las isostáticas correspondientes a este modelo. Ajustando una malla, con el mismo número de nudos que la inicial, a las líneas isostáticas con nudos situados de forma regular, se obtiene otra que está próxima a la óptima. A la malla así obtenida, se la denomina isostática isométrica. Esta conclusión se ha probado que es válida para los diferentes funcionales que se utilizan para evaluar la bondad de la malla en elementos finitos.

The design of shell and spatial structures represents an important challenge even with the use of the modern computer technology.If we concentrate in the concrete shell structures many problems must be faced,such as the conceptual and structural disposition, optimal shape design, analysis, construction methods, details etc. and all these problems are interconnected among them. As an example the shape optimization requires the use of several disciplines like structural analysis, sensitivity analysis, optimization strategies and geometrical design concepts. Similar comments can be applied to other space structures such as steel trusses with single or double shape and tension structures. In relation to the analysis the Finite Element Method appears to be the most extended and versatile technique used in the practice. In the application of this method several issues arise. First the derivation of the pertinent shell theory or alternatively the degenerated 3-D solid approach should be chosen. According to the previous election the suitable FE model has to be adopted i.e. the displacement,stress or mixed formulated element. The good behavior of the shell structures under dead loads that are carried out towards the supports by mainly compressive stresses is impaired by the high imperfection sensitivity usually exhibited by these structures. This last effect is important particularly if large deformation and material nonlinearities of the shell may interact unfavorably, as can be the case for thin reinforced shells. In this respect the study of the stability of the shell represents a compulsory step in the analysis. Therefore there are currently very active fields of research such as the different descriptions of consistent nonlinear shell models given by Simo, Fox and Rifai, Mantzenmiller and Buchter and Ramm among others, the consistent formulation of efficient tangent stiffness as the one presented by Ortiz and Schweizerhof and Wringgers, with application to concrete shells exhibiting creep behavior given by Scordelis and coworkers; and finally the development of numerical techniques needed to trace the nonlinear response of the structure. The objective of this paper is concentrated in the last research aspect i.e. in the presentation of a state-of-the-art on the existing solution techniques for nonlinear analysis of structures. In this presentation the following excellent reviews on this subject will be mainly used.

The aim of this study is to explain the changes in the real estate prices as well as in the real estate stock market prices, using some macro-economic explanatory variables, such as the gross domestic product (GDP), the real interest rate and the unemployment rate. Several regressions have been carried out in order to express some types of incremental and absolute deflated real estate lock market indexes in terms of the macro-economic variables. The analyses are applied to the Swedish economy. The period under study is 1984-1994. Time series on monthly data are used. i.e. the number of data-points is 132. If time leads/lags are introduced in the e regressions, significant improvements in the already high correlations are achieved. The signs of the coefficients for IR, UE and GDP are all what one would expect to see from an economic point of view: those for GDP are all positive, those for both IR and UE are negative. All the regressions have high R2 values. Both markets anticipate change in the unemployment rate by 6 to 9 months, which seems reasonable because such change can be forecast quite reliably. But, on the contrary, there is no reason why they should anticipate by 3-6 months changes in the interest rate that can hardly be reliably forecast so far in advance.

Se expone brevemente el problema matemático de homogeneización, en sus facetas: homogeneización en el dominio y en el contorno. Con respecto a esta última, se introducen los conceptos de talla crítica y problema límite. Los resultados obtenidos se aplican a un caso particular representado por un cilindro con condiciones de borde periódicas. Se comprueba la validez del desarrollo teórico mediante una experimentación numérica utilizando elementos finitos 3-D observándose la velocidad de convergencia que se obtiene en relación con el caso límite. Finalmente se comenta la posibilidad de utilizar la técnica de homogeneización para obtener unos criterios hacia un diseño eficiente y óptimo de los resultados de las barras corrugadas en el hormigón armado.

La solución al problema de encontrar la malla óptima en Elementos Finitos (EF), con un determinado número de grados de libertad, presenta un indudable interés en la aplicación del método. En la actualidad, el problema se plantea en términos de un proceso que permite obtener una mejor malla de elementos finitos a partir de una inicial. La nueva malla se diseña matemáticamente (remallado) de forma que el error del método sea lo más uniforme posible en todo el dominio de cálculo. Sin embargo, esta técnica de indudable interés y aplicación, al aumentar el número de grados de libertad (gdl) de la aproximación, no permite deducir de un modo directo el problema de la malla óptima condicionada a un número fijo de gdl. Con la solución de este problema se podrán deducir algunos criterios y recomendaciones para el diseño de una malla de elementos finitos, que exigirá, en general, en un proceso de remallado, modificaciones menores. Para problemas unidimensionales (barras y pilares simples), se pueden encontrar soluciones analíticas. Para problemas 2-D más complicados (tensión y deformación plana), se han utilizado métodos numéricos para obtener la malla óptima. Existen varios criterios de optimización, aquí se utiliza el del mínimo de la energía potencial total (EPT). Algunos ejemplos ilustrativos del método de optimización se presentan, indicándose algunas conclusiones.

Since the advent of the computer into the engineering field, the application of the numerical methods to the solution of engineering problems has grown very rapidly. Among the different computer methods of structural analysis the Finite Element (FEM) has been predominantly used. Shells and space structures are very attractive and have been constructed to solve a large variety of functional problems (roofs, industrial building, aqueducts, reservoirs, footings etc). In this type of structures aesthetics, structural efficiency and concept play a very important role. This class of structures can be divided into three main groups, namely continuous (concrete) shells, space frames and tension (fabric, pneumatic, cable etc )structures. In the following only the current applications of the FEM to the analysis of continuous shell structures will be discussed. However, some of the comments on this class of shells can be also applied to some extend to the others, but obviously specific computational problems will be restricted to the continuous shells. Different aspects, such as, the type of elements,input-output computational techniques etc, of the analysis of shells by the FEM will be described below. Clearly, the improvements and developments occurring in general for the FEM since its first appearance in the fifties have had a significative impact on the particular class of structures under discussion.

Se ha desarrollado un modelo implícito no lineal 2-D en EF para la resolución de las ecuaciones de aguas poco profundas. La discretización espacial se ha realizado por medio de elementos lagrangianos isoparamétricos. Se ha aplicado la integración numérica de Simpson para obtener las matrices elementales, y para la integración temporal se han utilizado diferentes esquemas en diferencias finitas, comprobándose el modelo con diferentes ejemplos.

El Cálculo matricial de estructuras corresponde a un planteamiento moderno y original del Análisis de Estructuras, en su acepción más genuina del término, de vuelta a los orígenes. Su desarrollo, y por lo tanto aprendizaje más eficaz, se realiza con la ayuda de un computador. Sin embargo, relegar al estudiante al simple papel de usuario de programas generales de cálculo de estructuras por computador, que puede utilizar como una misteriosa caja negra, sin comprensión de sus fundamentos, parece una penosa y poco formativa perspectiva educacional. Por otra parte, el desarrollo por el alumno de pequeños programas de cálculo, permite una adecuada y natural asimilación de algunos conceptos del cálculo matricial de estructuras, si bien exige un previo conocimiento de procedimientos numéricos e informativos y, particularmente, un esfuerzo, a veces desproporcionado, para un estudiante que sólo desea comprender pero no desarrollar nuevos programas de cálculo. En esta publicación se intenta un camino intermedio, en donde se da especial énfasis a los conceptos fundamentales de matriz de rigidez y cargas equivalentes de elementos y subestructuras. Conceptos que son útiles y prácticos, tanto para el futuro usuario como para el realizador de programas. Este hecho se puede comprobar en los modernos programas generales donde existen facilidades de subestructuración a diferentes niveles. Además, se incluyen aquí simples problemas de análisis matricial de estructuras que se plantean de un modo manual y que no precisan para su resolución, necesariamente, un computador. Por último se muestran ejemplos de uso de un programa general de computador, implementado en el Centro de Cálculo de la Universidad de Santander.

Se hace una pequeña introducción y después un estudio sobre las posibilidades y limitaciones en análisis de placas delgadas de elementos simples polinómicos de clase C.1. Se expone una familia jerárquica de dichos elementos, que se aplica a varios casos particulares. En base a estos se deducen algunas conclusiones, especialmente en lo que se refiere a eficacia computacional. Al final se proponen trabajos a realizar a partir de los datos existentes.

Los métodos semianalíticos poseen gran idoneidad para abordar problemas estructurales reales, obteniendo los niveles de precisión necesarios con esfuerzos de ordenador menores que los métodos exclusivamente numéricos.

The computational advantages of the use of different approaches -numerical and analytical ones- to the analysis of different parts of the same shell structure are discussed. Examples of large size problems that can be reduced to those more suitable to be handled by a personal related axisyrometric finite elements, local unaxisymmetric shells, geometric quasi-regular shells, infinite elements and homogenization techniques are described

En esta comunicación se pretende dar una breve introducción a las inecuaciones variacionales, y exponer las líneas de trabajo que sobre este tema se desarrollan en el Departamento de Ecuaciones Funcionales de la Universidad de Santiago.Asimismo se resumen algunas contribu ciones publicadas o en curso de realización.

En la actualidad el estudio del comportamiento de los tableros de puentes bajo ciertas acciones de carácter dinámico puede constituir un importante factor a considerar en el curso del proyecto. Ello se debe al progresivo incremento de la esbeltez de los tableros de puentes Y ello es debido a un progresivo incremento de la esbeltez de debido a causas múltiples, entre las que pueden citarse las de índole económica y las resultantes de un mejor conocimiento de la fenomenología estructural. Las acciones que son susceptibles de provocar una respuesta dinámica en un puente,son de origen vario: eólicas, sísmicas, impacto, circulación, etc. El carácter errático en su actuación,constituye una característica común a la mayoría de las excitaciones dinámicas y dificulta extraordinariamente una definición determinista adecuada de las mismas. Hasta la fecha se ha utilizado el métódo normal de análisis dinámico que corresponde al denominado cálculo quasi-estático.Es decir, la respuesta dinámica de una estructura se evalúa como el resultado obtenido de un cálculo estático (sin consideración de las fuerzas de inercia) multiplicado por un coeficiente de mayoración dinámica. De esta forma, se evita el cálculo dinámico de la estructura ,en general más complejo, pero exige el conocimiento de los coeficientes de mayoración dinámica. Por ello, se observa la necesidad de un cálculo dinámico completo en aquellos casos en los que la imprecisión del análisis quasi-estático es excesiva y se desea, por lo tanto, bien determinar el coeficiente de mayoración o bien el cálculo dinámico directo de la estructura.

En el estudio de determinados elementos estructurales de clmentacl6n -vigas y losas flotantes, pilotes, etc. - se han desarrollado modelos simplificados para simular el omportamiento del terreno, asomilándolo a un medio continuo y elástico. Esta idealización, claramente alejada de la realidad, ha sido y aún sigue siendo utilizada en numerosos casos en los que el grado de fiabilidad requerido no hace recomendable el recurrir a técnicas de análisis más sofisticadas y económicamente más costosas, como podrfa ser, por ejemplo, la discretización del conjunto "elemento estructural-terreno de cimentación" mediante elementos finitos.

Se presenta una aplicación muy concreta del Método de Elementos Finitos al trazado de carreteras, que puede ser de interés en algunos casos específicos como en el tratamiento del problema de optimización.

Existen situaciones para las que los resultados de interés corresponden a una zona pequeña de la estructura en comparación con su totalidad. La utilización de elementos infinitos constituye la posibilidad que evita un esfuerzo excesivo de cálculo. Esta idea se aplica en este artículo al análisis estático de la flexión de placas, desarrollando un elemento finito rectangular con un lado situado a distancia infinita. La función de decrecimiento de los movimientos en el interior del elemento es exponencial con un parámetro cuyo valor se determina mediante minimización de la energía potencial total. Se presentan resultados obtenidos con el elemento, contrastando su bondad y eficacia en comparación con un cálculo convencional. Se comentan finalmente distintas extensiones posibles del elemento.

Hasta fecha relativamente reciente, el cálculo estructural de placas con geometría y condiciones de borde arbitrarias ha constituido un problema complejo, frecuentemente inabordable sin el recurso de importantes y drásticas simplificaciones. Los métodos entonces existentes, -analíticos, semianalíticos (desarrollos en series) o numéricos (diferencias finitas)- eran incapaces o exigían un elevado grado de laboriosidad en su resolución. Con la aparición del método de los elementos finitos (l), (MEF), la situación se modificó de un modo considerable, al existir la posibilidad de un tratamiento unificado dentro de la teoría general del análisis matricial de estructuras -del cálculo de placas- con un mínimo de aproximaciones. No obstante se comprobó que la versión en movimientos del MEF encontraba mayores dificultades en su aplicación a problemas de la clase C1 –como son los de la flexión de placas delgadas- en comparación con los problemas C0- correspondientes a casos de elasticidad (tensión y deformación plana, por ejemplo).

Existen numerosas situaciones en la Técnica para las que es preciso resolver problemas de condiciones iniciales y contorno con una elevada exigencia de continuidad en las soluciones. Algunos ejemplos no exhaustivos se citan a continuación: flexión de vigas y placas en el análisis de las estructuras, problemas de láminas, topografía, trazado de vías de comunicación, definición geométrica de estructuras o reconocimiento caligráfico. Para ello se utilia el método de los elementos finitos que constituye en la actualidad una técnica matemática de discretización bien establecida . Su eficiencia se manifiesta en particular en la resolución de problemas de contorno planteados en su formulación débil y en los problemas de determinación de extrema les de funcionales.

Se desarrolla un elemento finito especial que permite la conexión entre un elemento viga (con flexión) y un elemento cuadrilátero de extensión (tensión o deformación plana). Asimismo, este elemento especial sirve para modelizar la discontinuidad existente en el contacto entre diferentes materiales. La formulación en movimientos del elemento, permite su inserción directa en un programa general de elementos finitos y, de esta forma, calcular estructuras de hormigón en contacto con el suelo (muros de contención, cimentación, estructuras enterradas como túneles y tuberías, etc.), tanto en el rango elástico como en el elastoplástico, utilizando una potencia limitada en medios computacionales.

Es una aplicación del método de los elementos finitos (M.E.F) al cálculo de losas delgadas isotrópicas. Es pues un desarrollo de la función solución en suma de funciones a trozos. Dentro del M.E.F se utilizan como funciones interpolantes polinomios (integración numérica sencilla). La continuidad conseguida es C elevado a 1 (para el caso planteado representa convergencia monotónica ). Son elementos simples (de fácil extensión a láminas) y que forman una familia jerárquica (distintos grados de aproximación sin cambiar la malla). El primer elemento de la familia es el clough- felippa. Al final se dan resultados comparativos de algunas placas con otro tipo de elementos y la solución exacta.

Curso monográfico realizado del 26 al 30 de mayo de 1980, en el que se resuelven una serie de problemas de elementos finitos. Los temas tratados son: 1. Problemas generales. Ancho de banda. Diseño de mallas. 2.-Términos numéricos de ecuaciones diferenciales. 3.-Funciones de forma. 4 y 5 - Matrices de rigidez y cargas equivalentes. 6.- Aplicaciones. 7.- Ejercicios de cálculo numérico. 8.-Ejercicios de programación

El método de los elementos finitos ha encontrado numerosas aplicaciones en la resolución de problemas de tipo estructural. No obstante sus posibilidades son mucho más amplias, como se pone de manifiesto en la aplicación que en el siguiente artículo se hace de dicho método al trazado de carreteras, de interés en casos específicos, como el tratamiento de problemas de optimización.

Estas notas que se publican a continuación corresponden a un curso de postgrado impartido durante el primer semestre del año 1983. El interés mostrado por los asistentes a dicho curso nos ha animado a escribir un resumen de las clases. Este libro supone un conocimiento te6rico de las ideas básicas del método de los elementos finitos. No obstante en una primera lección se resumen y ordenan aquellos aspectos mas importantes, que serán utilizados en lecciones sucesivas. En estas se desarrolla un programa de computador muy sencillo -sin complicaciones informáticas que obscurezcan la simplicidad del método- y se analiza de un modo detallado -en forma de organigramas y listados comentados- las distintas rutinas en lenguaje FORTRAN de este programa. Asimismo, y respetando el carácter elemental de la exposición se abren algunas posibilidades de ampliación y nuevos desarrollos del método. Algunos ejercicios y ejemplos al final de cada capítulo se espera permitan clarificar los puntos mas conflictivos del método. Finalmente se reúne en un apéndice, los distintos programas que se han mostrado en las sucesivas lecciones y que con objeto de que puedan ser procesados en microcomputadores se han traducido al lenguaje BASIC. Creemos y la experiencia del curso así nos la ha confirmado, que el método de elementos finitos se debe enseñar y aprender mediante la praxis y presentar los sucesivos desarrollos del método de un modo motivado como solución a problemas numéricos e informáticos que aparecen en su desarrollo. Si las lecciones que aquí se presentan permiten transmitir mejor estas ideas, los autores se sentirán más que recompensados por el trabajo que ha supuesto dar a luz a esta publicación.

El problema de la optimización en el diseño de presas - bóveda, planteado en toda su generalidad, corresponde a la definición de este diseño, es decir, de los siguientes parámetros dela presa bóveda: 1) Forma. 2) Altura. 3) Localización. 4) Propiedades. 5) Procedimiento de construcción. 6) Tratamiento de la cimentación y juntas con el terreno. 7) Localización y tipo de obras secundarias (aliviadero, galerías, etc.). 8) Uso de armado o pretensado de modo que se minimice una cierta función mérito,que exprese el coste de la estructura en función de estos parámetros.

The Department of Structural Analysis of the University of Santander has been for a longtime involved in the solution of the country´s practical engineering problems. Some of these have required the use of non-conventional methods of analysis, in order to achieve adequate engineering answers. As an example of the increasing application of non-linear computer codes in the nowadays engineering practice, some cases will be briefly presented. In each case, only the main features of the problem involved and the solution used to solve it will be shown

El objetivo de este artículo es fundamentalmente didáctico. Corresponde a un resumen de una monografía escrita sobre el tema, como report interno del Departamento de Análisis de las Estructuras de la Escuela Técnica Superior de Ingenieros de Caminos, Canales y Puertos de Santander, e intenta presentar como varía la exactitud -o velocidad de la convergencia- en los resultados obtenidos mediante el método de los elenentos finitos(MEF).Con objeto de reducir al máximo el esfuerzo computacional se ha considerado un problema de continuidad Cº, y únicamente un elemento aislado, concretamente, el caso de una barra aislada de sección variable bajo acción axil. Este ejemplo sencillo no permite deducir conclusiones definitivas y extensibles a otros casos pero sí probablenente, obtener algunas indicaciones sobre la importancia de algunos aspectos de MEF en relación con la exactitud de los resultados.

Existen diferentes métodos para el estudio de la calidad del agua en los estuarios, ríos, lagos, en zonas afectadas por los vertidos de sustancias contaminantes; que van desde los métodos experimentales tradicionales, hasta los más recientes modelos matemáticos. El modelo desarrollado permite resolver la ecuación de dispersión 2—D mediante técnicas en elementos finitos y por lo tanto obtener la evolución espacio-temporal de la concentración de constituyente. Se ha observado su exactitud en diferentes situaciones prácticas, y en particular se ha aplicado a la Bahía de Santander, analizándose el efecto de diferentes tipos de vertidos. El modelo es susceptible de utilizarse conjuntamente con otro modelo hidrodinámico capaz de simular la evolución del campo de velocidades.

En el presente artículo se plantea el problema de las condiciones de borde que deben ser consideradas, para los límites laterales de un modelo de elementos finitos, en el caso de que se pretenda analizar un problema.geotécnico en condiciones estáticas. Partiendo de las soluciones existentes para el caso de problemas dinámicos de interacción suelo-estructura, los autores demuestran que puede deducirse un nuevo borde aplicable a las condiciones estáticas. Dicho borde se materializa en una matriz de rigidez que interrelaciona entre si todos los nudos del lado que se considera. La bondad de la solución propuesta, así como la correspondiente a otras condiciones de borde, es estudiada en un problema concreto que dispone de una solución, estableciéndose unas conclusiones particulares respecto al comportamiento de las diferentes condiciones de borde analizadas.

Se presenta el cálculo de esfuerzos de una losa ortótropa rectangular apoyada en dos bordes opuestos y en los otros dos con condiciones de contorno muy generales. La formulación se realiza en forma compacta al utilizar la notación matricial, lo que permite una programación directa en un computador electrónico. La solución obtenida se compara con la que se deduce mediante el método de Guyon-Massonet-Rowe, que constituye una simplificación muy conocida, y se comentan las diferencias existentes. Se muestra a continuación un estudio alternativo mediante la teoría de láminas plegadas prismáticas, que representa un modelo matemático más adecuado para este tipode estructuras de tableros de puentes. Se comparan finalmente resultados numéricos entre ambas teorías.

Este trabajo estudia la aportación que los métodos de agregación de juicios de expertos pueden realizar en el cálculo de la peligrosidad sísmica de emplazamientos. Se han realizado cálculos en dos emplazamientos de la Península Ibérica: Mugardos (La Coruña) y Cofrentes (Valencia) que están sometidos a regímenes tectónicos distintos y que, además, alojan instalaciones industriales de gran responsabilidad. Las zonas de estudio, de 320 Km de radio, son independientes. Se ha aplicado un planteamiento probabilista a la estimación de la tasa anual de superación de valores de la aceleración horizontal de pico y se ha utilizado el Método de Montecarlo para incorporar a los resultados la incertidumbre presente en los datos relativos a la definición de cada fuente sismogenética y de su sismicidad. Los cálculos se han operado mediante un programa de ordenador, desarrollado para este trabajo, que utiliza la metodología propuesta por el Senior Seismic Hazard Analysis Commitee (1997) para la NRC. La primera conclusión de los resultados ha sido que la Atenuación es la fuente principal de incertidumbre en las estimaciones de peligrosidad en ambos casos. Dada la dificultad de completar los datos históricos disponibles de esta variable se ha estudiado el comportamiento de cuatro métodos matemáticos de agregación de juicios de expertos a la hora de estimar una ley de atenuación en un emplazamiento. Los datos de partida se han obtenido del Catálogo de Isosistas del IGN. Los sismos utilizados como variables raíz se han elegido con el criterio de cubrir uniformemente la serie histórica disponible y los valores de magnitud observados. Se ha asignado un panel de expertos particular a cada uno de los dos emplazamientos y se han aplicado a sus juicios los métodos de Cooke, equipesos, Apostolakis_Mosleh y Morris. Sus propuestas se han comparado con los datos reales para juzgar su eficacia y su facilidad de operación. A partir de los resultados se ha concluido que el método de Cooke ha mostrado el comportamiento más eficiente y robusto para ambos emplazamientos. Este método, además, ha permitido identificar, razonadamente, a aquellos expertos que no deberían haberse introducido en un panel. The present work analyses the possible contribution of the mathematical methods of aggregation in the assessment of Seismic Hazzard. Two sites, in the Iberian Peninsula, have been considered: Mugardos ( La Coruña) and Cofrentes (Valencia).Both of them are subjected to different tectonic regimes an both accommodate high value industrial plants. Their areas of concern, with radius of 320 Km, are not overlapping. A probabilistic approach has been applied in the assessment the annual probability of exceedence of the horizontal peak acceleration. The Montecarlo Method has allowed to transfer the uncertainty in the models and parameters to the final results. A computer program has been developed for this purpose. The methodology proposed by the Senior Seismic Analysis Committee (1997) for the NRC has been considered. Attenuation in Ground motion has been proved to be the main source of uncertainty in seismic hazard for both sites. Taking into account the difficulties to complete existing historical data in this subject the performance of four mathematical methods of aggregation has been studied. Original data have been obtained from the catalogs of the Spanish National Institute of Geography. The seismic events considered were chosen to cover evenly the historical records and the observed values of magnitude. A panel of experts have been applied to each site and four aggregation methods have been developed : equal weights, Cooke, Apostolakis-Mosleh and Morris The four proposals have been compaired with the actual data to judge their performance and ease of application. The results have shown that the Method of Cooke have proved the most efficient and robust for both sites. This method, besides, allow the reasoned identification of those experts who should be rejected from the panel