Refine
Document Type
- Conference Proceeding (24)
- Article (17)
Institute
- Professur Informatik im Bauwesen (41) (remove)
Keywords
- Finite-Elemente-Methode (41) (remove)
Im vorliegenden Beitrag wird ein in das FE-Programmsystem ANSYS implementiertes elastoplastisches Berechnungsmodell zur nichtlinearen, räumlichen Untersuchung von Mauerwerkstrukturen vorgestellt. Die Modellierung des heterogenen Baustoffs Mauerwerk erfolgt mit Hilfe eines verschmierten Ersatzkontinuums. Das anisotrope Materialverhalten wird sowohl hinsichtlich der Spannungs-Dehnungsbeziehung als auch bei der Beschreibung der Festigkeit berücksichtigt. Durch die Verwendung einer zusammengesetzten Fließbedingung ist es möglich, das Versagen der einzelnen Mauerwerkkomponenten Stein und Mörtelfugen und des Verbundes zu berücksichtigen. Dadurch ist die Anwendbarkeit des Modells für mehrere Mauerwerksarten gegeben. Die hierfür verwendeten Materialparameter sind aus einfachen Kleinkörperversuchen bestimmbar oder innerhalb gewisser Grenzen aus empirischen Formeln berechenbar. Die notwendige Beschränkung der Anzahl der Materialparameter sichert die praktische Anwendbarkeit des entwickelten Berechnungsmodells. Die numerische Umsetzung des hier verwendeten impliziten Berechnungsverfahrens lässt sich in eine lokale und eine globale Iterationsebene gliedern. Die lokale Iteration am Integrationspunkt dient der Spannungsrückführung. Dabei sind die Besonderheiten der Verarbeitung mehrflächiger Fließfiguren zu beachten. Die globale Iteration auf Systemebene sichert die Umlagerung des Residuums. Mit der Nachrechnung von Versuchsergebnissen soll das entwickelte Modell verifiziert und seine physikalische Leistungsfähigkeit eingeschätzt werden.
Numerische Approximation makroskopischer Verkehrsmodelle mit der Methode der Finiten Elemente
(2000)
Makroskopische Verkehrsmodelle sind ein wesentliches Hilfsmittel bei der Beurteilung und Steuerung von Verkehrsflüssen auf Hauptverkehrsadern. Für die notwendige Beeinflussung des Verkehrsablaufs werden Online-Messungen und prognostische numerische Simulationen benötigt. Für die Simulationen bieten sich makroskopische Verkehrsmodelle an, die den Verkehr als kontinuierliche Fahrzeugströmeabbilden. Aufgrund der Analogie zu den Modellen der Strömungsmechanik lassen sich die numerischen Verfahren aus diesem Bereich auch zur Lösung makroskopischer Verkehrsmodelle verwenden. Es wird eine Finite-Elemente-Approximation für die numerische Umsetzung makroskopischer Verkehrsmodelle vorgestellt. Exemplarisch wird sie am Verkehrsmodell von Kerner und Konhäuser erläutert. Dieses und andere makroskopische Verkehrsmodelle wurden bisher mit der Methode der Finiten Differenzen gelöst. Die vorgestellte Approximation entspricht einem Petrov-Galerkin-Verfahren, bei dem der Fehler eines Standard-Galerkin-Verfahrens mit Hilfe eines Upwinding-Koeffizienten minimiert wird. Die Wahl des Upwinding-Koeffizienten ist übertragbar und basiert ausschließlich auf dem Charakter der zugrundeliegenden Gleichungen. Die Ergebnisse zeigen typische Phänomene eines Verkehrsablaufs wie die Entstehung von Stop-and-Go-Wellen oder Staus. Die Finite-Elemente-Methode erweist sich für unter-schiedlichste Verkehrsmodelle als ausgesprochen stabil.
The dynamic behaviour of shells, which are widely used in construction and mechanical engineering as critical components of machinery and 3-D structures, under static and dynamic loadings is described by system of deep nonlinear differential equations. Solution of these equations can be received with assistance of technique basing on a modern numerical algorithms and computer modeling.. The system of nonlinear differential equations of vibration of the shells is proposed taking into account the inertia forces in the tangential and normal directions. Its solution is based on combination of parameter prolongation method, finite-difference method and the Newton-Kantorovich iterative algorithm that allows plotting the loading trajectories and determination of bifurcation points on them. Package of Applied Programs >SEVSOR< is a computation means to be used in research of deformation, stability and vibration in thin axically-symmetric shells of complicated shape Input data include information on shell geometry, physical and mechanical properties, bearing conditions, types of loadings and load application. Frame output of motion forms in real time or either in decelerated or accelerated time scales for creating cartoons or video films is used for analysis of the compound dynamic processes in shell-type structures.
The worldwide growth of communication networks and associated technologies provide the basic infrastructure for new ways of executing the engineering process. Collaboration amongst team members seperated in time and location is of particular importance. Two broad themes can be recognized in research pertaining to distributed collaboration. One theme focusses on the technical and technological aspects of distributed work, while the other emphasises human aspects thereof. The case of finite element structural analysis in a distributed collaboratory is examined in this paper. An approach is taken which has its roots in human aspects of the structural analysis task. Based on experience of how structural engineers currently approach and execute this task while utilising standard software designed for use on local workstations only, criteria are stated for a software architechture that could support collaborative structural analysis. Aspects of a pilot application and the results of qualitative performance measurements are discussed.
Dynamic testing for damage assessment as non-destructive method has attracted growing in-terest for systematic inspections and maintenance of civil engineering structures. In this con-text the paper presents the Stochastic Finite Element (SFE) Modeling of the static and dy-namic results of own four point bending experiments with R/C beams. The beams are dam-aged by an increasing load. Between the load levels the dynamic properties are determined. Calculated stiffness loss factors for the displacements and the natural frequencies show differ-ent histories. A FE Model for the beams is developed with a discrete crack formulation. Cor-related random fields are used for structural parameters stiffness and tension strength. The idea is to simulate different crack evolutions. The beams have the same design parameters, but because of the stochastic material properties their undamaged state isn't yet the same. As the structure is loaded a stochastic first crack occurs on the weakest place of the structure. The further crack evolution is also stochastic. These is a great advantage compared with de-terministic formulations. To reduce the computational effort of the Monte Carlo simulation of this nonlinear problem the Latin-Hypercube sampling technique is applied. From the results functions of mean value and standard deviation of displacements and frequencies are calcu-lated. Compared with the experimental results some qualitative phenomena are good de-scribed by the model. Differences occurs especially in the dynamic behavior of the higher load levels. Aim of the investigations is to assess the possibilities of dynamic testing under consideration of effects from stochastic material properties
In this paper we consider modelling of composite material with inclusions where the elastic material properties of both matrix and inclusions are uncertain and vary within prescribed bounds. Such mechanical systems, involving interval uncertainties and modelled by finite element method, can be described by parameter dependent systems of linear interval equations and process variables depending on the system solution. A newly developed hybrid interval approach for solving parametric interval linear systems is applied to the considered model and the results are compared to other interval methods. The hybrid approach provides very sharp bounds for the process variables - element strains and stresses. The sources for overestimation when dealing with interval computations are demonstrated. Based on the element strains and stresses, we introduce a definition for the values of nodal strains and stresses by using a set-theoretic approach.
The influence of vortex-induces vibrations on vertical tie rods has been proved as a determinant load factor in the lifetime-oriented dimensioning of arched steel bridges. Particularly, the welded connection plates between the suspenders and the arches often exhibit cracks induced primarily rods. In this context, the synchronization of the vortex-shedding to the rod motion in a critical wind velocity range, the so-called lock-in effect, is of essential interest.
A large-scale computer modeling and simulation method is presented for environmental flows in urban area. Several GIS and CAD data were used for the preparation of shape model and an automatic mesh generation method based on Delaunay method was developed. Parallel finite element method based on domain decomposition method was employed for the numerical simulation of natural phenomena. The present method was applied to the simulation of flood flow and wind flow in urban area. The present method is shown to be a useful planning and design tool for the natural disasters and the change of environments.
There was suggested a phenomenological modified quadratic condition of the beginning of plasticity for plastic and quasifragile orthotropic materials. Limiting surface in the shape of a paraboloid with an axis bend over hydrostatic axis corresponds to the condition. The equations of theory of current with the isotropic and anisotropic hardenings, associated with the suggested yield condition, modified into the version of determining equations of strain theory of plasticity are received. These defining equations formed the basis of highlyprecise non-classic continual (along thickness) theory of non-linear deformation of thick sandwich plates and sloping shells. In the approximations along the cross coordinate the specificity of flexural and non-flexural deformations is taken into account. The necessity of introducing the approximations of higher order, as well as accounting for the cross compression while decreasing of the relatively cross normal and shear layer rigidness is shown. The specifications, obtained in comparison with the known physically nonlinear specified model of the bending of plates with orthotropic layers are distinguished. An effective procedure of linearization of the solving equations and getting the solutions in frames of the discrete-continual scheme of the finite-element method is suggested. The approximations of higher order let to model the appearance of the cracs of layers being split by the introducing of slightly hard thin layers into the finite element, not violating the idea of continuality of theory. Calculation of a threelayer plate with rigid face diaphragms on the contour is considered
Die Methode der Finiten Elemente ist ein numerisches Verfahren zur Interpolation vorgegebener Werte und zur numerischen Approximation von Lösungen stationärer oder instationärer partieller Differentialgleichungen bzw. Systemen partieller Differentialgleichungen. Grundlage dieser Verfahren ist die Formulierung geeigneter Finiter Elemente und Finiter Element Zerlegungen. Finite Elemente besitzen in der Regel eine geometrische Basis bestehend aus Strecken im eindimensionalen, Drei- oder Vierecken im zweidimensionalen und Tetra- oder Hexaedern im dreidimensionalen euklidischen Raum, eine Menge von Freiheitsgraden und eine Basis von Funktionen. Die geometrische Basis eines Finiten Elements wird verallgemeinert als geometrische Zelle formuliert. Diese geschlossene geometrische Formulierung führt zu einer geometrieunabhängigen Definition der Basisfunktionen eines Finiten Elements in den Zellkoordinaten der geometrischen Zelle. Finite Elemente auf der Basis geometrischer Zellen werden als Bestandteile Finiter Element Zerlegungen in Finiten Element Interpolationen und Finiten Element Approximationen verwendet. Die Finiten Element Approximationen werden am Beispiel der 2-dimensionalen Diffusionsgleichung über das Standard-Galerkin-Verfahren ermittelt.
The method of the finite elements is an adaptable numerical procedure for interpolation as well as for the numerical approximation of solutions of partial differential equations. The basis of these procedure is the formulation of suitable finite elements and element decompositions of the solution space. Classical finite elements are based on triangles or quadrangles in the two-dimensional space and tetrahedron or hexahedron in the threedimensional space. The use of arbitrary-dimensional convex and non-convex polyhedrons as the geometrical basis of finite elements increases the flexibility of generating finite element decompositions substantially and is sometimes the only way to get a clear decomposition...
Detailuntersuchungen an Tragwerken führen bei FE-Berechnungen immer wieder auf das Problem einer geeigneten Netzgestaltung. Während in weiten Bereichen ein grobes Netz ausreicht, muß an kritischen Stellen ein sehr feines Netz gewählt werden, um gerade dort hinreichend genaue Ergebnisse zu erhalten. Bei der Realisierung lokaler Netzverdichtungen stellt die Gestaltung des Übergangs vom groben zum feinen Netz das Hauptproblem dar. Im Beitrag wird hierzu eine Familie von FE-Übergangselementen vorgestellt, mit denen sich eine voll-kompatible Kopplung von wenigen großen Elementen mit vielen kleinen Elementen bereits über nur eine Stufe erzielen läßt. Diese neu entwickelten sogenannten pNh-Elemente ermöglichen an einer oder mehreren Seiten den Anschluß von N kleineren Elementen (Elementseiten für h-Verfeinerung). Das wird durch N stückweise definierte Ansatzfunktionen an den entsprechenden Seiten erreicht, wobei die Teilung nicht äquidistant sein braucht. Darüber hinaus ist es möglich, Elemente unterschiedlichen Polynomgrades p an den Standardseiten und den Verfeinerungsseiten anzuschließen. Der praktische Einsatz der Übergangselemente setzt geeignete automatische oder halbautomatische Netzgeneratoren voraus, die diese Elemente einbeziehen. Im Rahmen einer substrukturorientierten Modellierung läßt sich dies besonders günstig realisieren. Im Beitrag wird gezeigt, wie durch Zerlegung des Gesamtmodells in Bereiche mit grobem Netz, mit Übergangsnetz und mit feinem Netz, eine effektive Generierung der Netzverdichtungen zu erreichen ist. An einem praktischen Beispiel aus dem Bauingenieurwesen werden die Vorteile des vorgestellten Übergangselementkonzeptes umfassend demonstriert.
In der vorliegenden Arbeit werden dickwandige Schalentragwerke unter statischen Belastungen betrachtet. Die Schale besteht aus verschiedenen Zonen und in jeder Zone wird die Schalenmittelflaeche mittels eines eigenen Geometriegleichungssystems definiert. Das Verzerrungsfeld hat allen 6 unabhaengigen Komponenten unter der Annahme, dass die Querschnittsfasern, die normal zu der Mittelflaeche der unbelasteten Schale sind, geradelinig bleiben. Ein dreidimensionales isoparametrisches finites Element wird vorgeschlagen. Die Berechnung wird mit der Hilfe der Makroelemententechnik durchgefuehrt. In der Arbeit werden die wesentliche Parameter der Schalengeometrie, sowie auch entsprechendes Anteil von Klassen des konstruktiven Modells, definiert. Ein konstruktives Informationsmodell und ein FEM-Informationsmodell, werden entwickelt. Die Informationsverbindungen zwischen den beiden Modellen werden definiert. Diese objektorientierten Modelle werden in Programmiersprache Microsoft Visual C++ v.4.0 unter Windows 95 implementiert. Als numerisches Beispiel wird ein Bogenmauertragwerk betrachtet.
In this paper, systematic analyses for the shoring systems installed to support the applied loads during construction are performed on the basis of the numerical approach. On the basis of a rigorous time-dependent analysis, structural behaviors of reinforced concrete (RC) frame structures according to the changes in design variables such as the types of shoring systems, shore stiffness and shore spacing are analyzed and discussed. The time-dependent deformations of concrete such as creep and shrinkage and construction sequences of frame structures are also taken into account to minimize the structural instability and to reach to an improved design of shoring system because these effects may increase the axial forces delivered to the shores. In advance, the influence of the column shortening effect, generally mentioned in a tall building structure, is analyzed. From many parametric studies, it has been finally concluded that the most effective shoring system in RC frame structures is 2S1R (two shores and one reshore) regardless of the changes in design variables.
The paper deals with the simulation of the non-linear and time dependent behaviour of complex structures in engineering. Such simulations have to provide high accuracy in the prediction of deformations and stability, by taking into account the long term influences of the non-linear behaviour of the material as well as the large deformation and contact conditions. The limiting factors of the computer simulation are the computer run time and the memory requirement during solving large scale problems. To overcome these problems we use a dynamic-explicit time integration procedure for the solution of the semi-discrete equations of motion, which is very suited for parallel processing. In the paper at first we give a brief review of the theoretical background of the mechanical modelling and the dynamic-explicit technique for the solution of the semi-discrete equations of motion. Then the concept of parallel processing will be discussed . A test example concludes the paper.
Bei komplexen Gründungskonstruktionen sind Planungsfehler durch eine konsistente Modellierung vermeidbar. Manuelle Berechnungsmethoden ermöglichen im allgemeinen ein dreidimensionales Vorgehen nicht. Numerische Berechnungsmethoden, wie z.B. die Finite-Element-Methode, sind ein optimales Werkzeug zur ganzheitlichen Simulation des Problems. Die für die Finite-Element-Analyse notwendige Diskretisierung komplexer Bau- grundstrukturen ist manuell nicht zu bewältigen. Der vorliegende Beitrag zeigt wie ein Finite-Element-Modell automatisch aus einem geotechnischen Modell unter Berücksichtigung der spezifischen Anforderungen der Baugrund-Tragwerk-Struktur und des Bauablaufes erzeugt werden kann. Hierbei wird die Berücksichtigung der geometrischen und der mechanischen Besonderheiten bei der Netzgenerierung dargestellt.
Discrete-continual Finite Element Method of Analysis for Three-dimensional Curvilinear Structures
(2003)
This paper is devoted to discrete-continual finite element method (DCFEM) of analysis for three-dimensional curvilinear structures. Operational and variational formulations of the problem in the ring coordinate system are presented. The discrete-continual design model for structures with constant physical and geometrical parameters in longitudinal direction is offered on the basis of so-called curvilinear discrete-continual finite elements. Element coordinate system, approximation of nodal unknowns, construction of element nodal load vector are under consideration. Element system of differential equations is formulated with use of special generalized block-structured stiffness matrix of discrete-continual finite element. Local differential relations are formulated. Resultant multipoint boundary problem for system of ordinary differential equations is given. Method of analytical solution of multipoint boundary problems in structural analysis is offered as well. Its major peculiarities include universality, computer-oriented algorithm involving theory of distributions, computational stability, optimal conditionality of resultant systems, partial Jordan decomposition of matrix of coefficients, eliminating necessity of calculation of root vectors. Brief information concerning developed software is provided.
Die Eisenbahnbrücken des Lehrter Bahnhofs in Berlin - Ein ganzheitliches FE-Berechnungskonzept
(2000)
Der Komplexität moderner Brückenbauwerke scheinen die verwendeten Berechungsmodelle oft nicht angemessen. Tragwerksberechnungen basieren in vielen Fällen noch auf der Vorgehensweise, das Brückenbauwerk in Einzelbauteile zu zerlegen und mit unterschiedlichen Teilmodellen zu behandeln. Das erscheint, auch vor dem Hintergrund ständig wachsender Rechnerleistung, nicht mehr zeitgemäß. Dies gilt zum Beispiel auch für die gängige Praxis, flächenhafte Brückenüberbauten mit Balkenmodellen zu berechnen. Der vorliegende Beitrag stellt ein ganzheitliches Berech-nungskonzept vor, welches auf der Basis eines einzigen FE-Modells die Berechnung des Gesamtbauwerks erlaubt. Damit wird für alle Bauteile neben der Zustandsgrößenberechnung auch die Bemessung von Stahl- und Spannbetonbauteilen bis hin zu Nachweisen wie zur Beschränkung der Rissbreite geführt. Die Anwendung dieses Berechnungskonzeptes wird am Beispiel der Eisenbahnüberführung des neuen Lehrter Bahn-hofs in Berlin gezeigt. Das verwendete FE-Modell umfasst Baugrund, Fundamente, Stahl- bzw. Gußstahlunterkonstruktion sowie den Stahl- bzw. Spannbetonüberbau. Besonderheiten sind unter anderem die Modellierung des plattenbalkenartigen Überbaus durch exzentrische, vorspannbare Schalenelemente und das getrennte Vorhalten von tragwerks- und lastbezogenen Eingabefiles. Damit gelingt die sequentielle Erfassung unterschiedlicher Bettungsmoduli zur Simulation statischer und dynamischer Beanspruchungen, die Berücksichtigung des Anspannens und der Interaktion zwischen vorgespannten Stahlverbänden zur Aufnahme von Horizontallasten sowie die Berücksichtigung unterschiedlicher statischer Systeme bei der Herstellung des Spannbetonüberbaus.
Development and Analysis of Sparse Matrix Concepts for Finite Element Approximation on general Cells
(2004)
In engineering and computing, the finite element approximation is one of the most well-known computational solution techniques. It is a great tool to find solutions for mechanic, fluid mechanic and ecological problems. Whoever works with the finite element method will need to solve a large system of linear equations. There are different ways to find a solution. One way is to use a matrix decomposition technique such as LU or QR. The other possibility is to use an iterative solution algorithm like Conjugate Gradients, Gauß-Seidel, Multigrid Methods, etc. This paper will focus on iterative solvers and the needed storage techniques...
For modeling of singular fields of stresses and deformations in elasters with a crack is offered to use of three-dimesional a special finite element. Weak compessible of elasters is taken into account on the basis of threefold approximation of fields of displacements, deformations and function of volume change. At intensive cyclic loading of the elastomer constructions with a crack it is necessary to take into account warming and large deformations at the crack top. The stress-deformed state elasters with a crack is determined from the decision of a nonlinear problem by a modified method Newton-Kantorovich. Account stress intensity factors for a rectangular plate with a various arrangement of a through crack is executed. Process of development of a surface crack and dissipative warming in prismatic a element of shift is investigated.
The paper investigates accuracy of deflection predictions made by the finite element package ATENA and design code methods ACI and EC2. Deflections have been calculated for a large number of experimental reinforced concrete beams reported by three investigators. Statistical parameters have been established for each of the technique at different load levels, separately for the beams with small and moderate reinforcement ratio.
Creation of hierarchical sequence of the plastic and viscoplastic models according to different levels of structure approximations is considered. Developed strategy of multimodel analysis, which consists of creation of the inelastic models library, determination of selection criteria system and caring out of multivariant sequential clarifying computations, is described. Application of the multimodel approach in numerical computations has demonstrated possibility of reliable prediction of stress-strain response under wide variety of combined nonproportional loading.
This paper presents the combination of two different parallelization environments, OpenMP and MPI, in one numerical simulation tool. The computation of the system matrices and vectors is parallelized with OpenMP and the solution of the system of equations is done with the MPIbased solver MUMPS. The efficiency of both algorithms is shown on several linear and nonlinear examples using the Finite Element Method and a meshless discretization technique.
The steel structure design codes require to check up the member strength when evaluating plastic deformations. The model of perfectly plastic material is accepted. The strength criteria for simple cross-sections (I section, etc.) of steel members are given in design codes. The analytical strength criteria for steel cross-sections and numerical approaches based on stepwise procedure are investigated in many articles. Another way for checking the carrying capacity of cross-sections is the use of methods that are applied for defining strain-deformed state of elastic perfectly plastic systems. In this paper non-iterative methods are suggested for checking strength of cross-sections. Carrying capacity of cross section is verified according to extremum principle of plastic fail under monotonically loading and the strain-deformed state of cross-section is defined according to extremum energy principals of elastic potential of residual stresses and complementary work of residual displacements. The mathematical expressions of these principals for discrete cross-section are formulated as problems of convex mathematical programming. The cross-section of steel member using finite element method is divided into free form plane elements. The constant distribution of stresses along the finite element is accepted. The relationships of finite elements for static formulation of the problem are formed so, that kinematics formulation relationships could be obtained in a formal way using the theory of duality. Numerical examples of determination of cross-section strength, composition of interactive curves and composition of moment-curvature curves for different axial force levels are presented.
Framed-tube system with multiple internal tubes is analysed using an orthotropic box beam analogy approach in which each tube is individually modelled by a box beam that accounts for the flexural and shear deformations, as well as the shear-lag effects. A simple numerical modeling technique is proposed for estimating the shear-lag phenomenon in tube structures with multiple internal tubes. The proposed method idealizes the framed-tube structures with multiple internal tubes as equivalent multiple tubes, each composed of four equivalent orthotropic plate panels. The numerical analysis is based on the minimum potential energy principle in conjunction with the variational approach. The shear-lag phenomenon of such structures is studied taking into account the additional bending moments in the tubes. A detailed work is carried out through the numerical analysis of the additional bending moment. The moment factor is further introduced to identify the shear lag phenomenon along with the additional moment.
To fulfil safety requirements the changes in the static and/or dynamic behaviour of the structure must be analysed with great care. These changes are often caused by local reduction of the stiffness of the structure caused by the irregularities in the structure, as for example cracks. In simple structures such analysis can be performed directly, by solving equations of motion, but for more complex structures a different approach, usually numerical, must be applied. The problem of crack implementation into the structure behaviour has been studied by many authors who have usually modelled the crack as a massless rotational spring of suitable stiffness placed at the beam at the location where the crack occurs. Recently, the numerical procedure for the computation of the stiffness matrix for a beam element with a single transverse crack has been replaced with the element stiffness matrix written in fully symbolic form. A detailed comparison of the results obtained by using 200 2D finite elements with those obtained with a single cracked beam element has confirmed the usefulness of such element.
The displacements and stresses in arch dams and their abutments are frequently determined with 20-node brick elements. The elements are distorted near the contact plane between the wall and the abutment. A cantilever beam testbed has been developed to investigate the consequences of this distortion. It is shown that the deterioration of the accuracy in the computed stresses is significant. A compatible 18-node wedge element with linear stress variation is developed as an alternative to the brick element. The shape of this element type is readily adapted to the shape of the contact plane. It is shown that the accuracy of the computed stresses in the vicinity of the contact plane is improved significantly by the use of wedge elements.
For the analysis of arbitrary, by Finite Elements discretized shell structures, an efficient numerical simulation strategy with quadratic convergence including geometrically and physically nonlinear effects will be presented. In the beginning, a Finite-Rotation shell theory allowing constant shear deformations across the shell thickness is given in an isoparametric formulation. The assumed-strain concept enables the derivation of a locking-free finite element. The Layered Approach will be applied to ensure a sufficiently precise prediction of the propagation of plastic zones even throughout the shell thickness. The Riks-Wempner-Wessels global iteration scheme will be enhanced by a Line-Search procedure to ensure the tracing of nonlinear deformation paths with rather great load steps even in the post-peak range. The elastic-plastic material model includes isotropic hardening. A new Operator-Split return algorithm ensures considerably exact solution of the initial-value problem even for greater load steps. The combination with consistently linearized constitutive equations ensures quadratic convergence in a close neighbourhood to the exact solution. Finally, several examples will demonstrate accuracy and numerical efficiency of the developed algorithm.
A geometrical inclusion-matrix model for the finite element analysis of concrete at multiple scales
(2003)
This paper introduces a method to generate adequate inclusion-matrix geometries of concrete in two and three dimensions, which are independent of any specific numerical discretization. The article starts with an analysis on shapes of natural aggregates and discusses corresponding mathematical realizations. As a first prototype a two-dimensional generation of a mesoscale model is introduced. Particle size distribution functions are analysed and prepared for simulating an adequate three-dimensional representation of the aggregates within a concrete structure. A sample geometry of a three-dimensional test cube is generated and the finite element analysis of its heterogeneous geometry by a uniform mesh is presented. Concluding, aspects of a multiscale analysis are discussed and possible enhancements are proposed.
Iso-parametric finite elements with linear shape functions show in general a too stiff element behavior, called locking. By the investigation of structural parts under bending loading the so-called shear locking appears, because these elements can not reproduce pure bending modes. Many studies dealt with the locking problem and a number of methods to avoid the undesirable effects have been developed. Two well known methods are the >Assumed Natural Strain< (ANS) method and the >Enhanced Assumed Strain< (EAS) method. In this study the EAS method is applied to a four-node plane element with four EAS-parameters. The paper will describe the well-known linear formulation, its extension to nonlinear materials and the modeling of material uncertainties with random fields. For nonlinear material behavior the EAS parameters can not be determined directly. Here the problem is solved by using an internal iteration at the element level, which is much more efficient and stable than the determination via a global iteration. To verify the deterministic element behavior the results of common test examples are presented for linear and nonlinear materials. The modeling of material uncertainties is done by point-discretized random fields. To show the applicability of the element for stochastic finite element calculations Latin Hypercube Sampling was applied to investigate the stochastic hardening behavior of a cantilever beam with nonlinear material. The enhanced linear element can be applied as an alternative to higher-order finite elements where more nodes are necessary. The presented element formulation can be used in a similar manner to improve stochastic linear solid elements.
The primary objective of initial shape analysis of a cable stayed bridge is to calculate initial installation cable tension forces and to evaluate fabrication camber of main span and pylon providing the final longitudinal profile of the bridge at the end of construction. In addition, the initial cable forces depending on the alternation of the bridge’s shape can be obtained from the analysis, and will be used to provide construction safety during construction. In this research, we conducted numerical experiments for initial shape of Ko-ha bridge, which will be constructed in the near future, using three different typical methods such as continuous beam method, linear truss method, and IIMF (Introducing Initial Member Force) method