Refine
Has Fulltext
- yes (73) (remove)
Document Type
- Article (27)
- Conference Proceeding (24)
- Doctoral Thesis (17)
- Master's Thesis (4)
- Diploma Thesis (1)
Institute
- Professur Informatik im Bauwesen (41)
- Institut für Strukturmechanik (ISM) (17)
- Professur Baumechanik (3)
- Professur Stahlbau (3)
- Bauhaus-Institut für zukunftsweisende Infrastruktursysteme (b.is) (2)
- Graduiertenkolleg 1462 (2)
- Professur Stahl- und Hybridbau (2)
- Professur Angewandte Mathematik (1)
- Professur Baustatik und Bauteilfestigkeit (1)
- Professur Modellierung und Simulation - Mechanik (1)
Keywords
- Finite-Elemente-Methode (73) (remove)
Piezoelectric materials are used in several applications as sensors and actuators where they experience high stress and electric field concentrations as a result of which they may fail due to fracture. Though there are many analytical and experimental works on piezoelectric fracture mechanics. There are very few studies about damage detection, which is an interesting way to prevent the failure of these ceramics.
An iterative method to treat the inverse problem of detecting cracks and voids in piezoelectric structures is proposed. Extended finite element method (XFEM) is employed for solving the inverse problem as it allows the use of a single regular mesh for large number of iterations with different flaw geometries.
Firstly, minimization of cost function is performed by Multilevel Coordinate Search (MCS) method. The XFEM-MCS methodology is applied to two dimensional electromechanical problems where flaws considered are straight cracks and elliptical voids. Then a numerical method based on combination of classical shape derivative and level set method for front propagation used in structural optimization is utilized to minimize the cost function. The results obtained show that the XFEM-level set methodology is effectively able to determine the number of voids in a piezoelectric structure and its corresponding locations.
The XFEM-level set methodology is improved to solve the inverse problem of detecting inclusion interfaces in a piezoelectric structure. The material interfaces are implicitly represented by level sets which are identified by applying regularisation using total variation penalty terms. The formulation is presented for three dimensional structures and inclusions made of different materials are detected by using multiple level sets. The results obtained prove that the iterative procedure proposed can determine the location and approximate shape of material subdomains in the presence of higher noise levels.
Piezoelectric nanostructures exhibit size dependent properties because of surface elasticity and surface piezoelectricity. Initially a study to understand the influence of surface elasticity on optimization of nano elastic beams is performed. The boundary of the nano structure is implicitly represented by a level set function, which is considered as the design variable in the optimization process. Two objective functions, minimizing the total potential energy of a nanostructure subjected to a material volume constraint and minimizing the least square error compared to a target
displacement, are chosen for the numerical examples. The numerical examples demonstrate the importance of size and aspect ratio in determining how surface effects impact the optimized topology of nanobeams.
Finally a conventional cantilever energy harvester with a piezoelectric nano layer is analysed. The presence of surface piezoelectricity in nano beams and nano plates leads to increase in electromechanical coupling coefficient. Topology optimization of these piezoelectric structures in an energy harvesting device to further increase energy conversion using appropriately modified XFEM-level set algorithm is performed .
Phase Field Modeling for Fracture with Applications to Homogeneous and Heterogeneous Materials
(2017)
The thesis presents an implementation including different applications of a variational-based approach for gradient type standard dissipative solids. Phase field model for brittle fracture is an application of the variational-based framework for gradient type solids. This model allows the prediction of different crack topologies and states. Of significant concern is the application of theoretical and numerical formulation of the phase field modeling into the commercial finite element software Abaqus in 2D and 3D. The fully coupled incremental variational formulation of phase field method is implemented by using the UEL and UMAT subroutines of Abaqus. The phase field method
considerably reduces the implementation complexity of fracture problems as it removes the need for numerical tracking of discontinuities in the displacement field that are characteristic of discrete crack methods. This is accomplished by replacing the sharp discontinuities with a scalar damage phase field representing the diffuse crack topology wherein the amount of diffusion is controlled by a regularization parameter. The nonlinear coupled system consisting of the linear momentum equation and a diffusion type equation governing the phase field evolution is solved simultaneously via a Newton-
Raphson approach. Post-processing of simulation results to be used as visualization
module is performed via an additional UMAT subroutine implemented in the standard Abaqus viewer.
In the same context, we propose a simple yet effective algorithm to initiate and propagate cracks in 2D geometries which is independent of both particular constitutive laws and specific element technology and dimension. It consists of a localization limiter in the form of the screened Poisson equation with, optionally, local mesh refinement. A staggered scheme for standard equilibrium and screened Cauchy equations is used. The remeshing part of the algorithm consists of a sequence of mesh subdivision and element erosion steps. Element subdivision is based on edge split operations using a
given constitutive quantity (either damage or void fraction). Mesh smoothing makes use of edge contraction as function of a given constitutive quantity such as the principal stress or void fraction. To assess the robustness and accuracy of this algorithm, we use both quasi-brittle benchmarks and ductile tests.
Furthermore, we introduce a computational approach regarding mechanical loading in microscale on an inelastically deforming composite material. The nanocomposites material of fully exfoliated clay/epoxy is shaped to predict macroscopic elastic and fracture related material parameters based on their fine–scale features. Two different configurations of polymer nanocomposites material (PNCs) have been studied. These configurations are fully bonded PNCs and PNCs with an interphase zone formation between the matrix and the clay reinforcement. The representative volume element of PNCs specimens with different clay weight contents, different aspect ratios, and different
interphase zone thicknesses are generated by adopting Python scripting. Different constitutive models are employed for the matrix, the clay platelets, and the interphase zones. The brittle fracture behavior of the epoxy matrix and the interphase zones material are modeled using the phase field approach, whereas the stiff silicate clay platelets of the composite are designated as a linear elastic material. The comprehensive study investigates the elastic and fracture behavior of PNCs composites, in addition to predict Young’s modulus, tensile strength, fracture toughness, surface energy dissipation, and cracks surface area in the composite for different material parameters, geometry, and interphase zones properties and thicknesses.
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.
We conducted extensive molecular dynamics simulations to investigate the thermal conductivity of polycrystalline hexagonal boron-nitride (h-BN) films. To this aim, we constructed large atomistic models of polycrystalline h-BN sheets with random and uniform grain configuration. By performing equilibrium molecular dynamics (EMD) simulations, we investigated the influence of the average grain size on the thermal conductivity of polycrystalline h-BN films at various temperatures. Using the EMD results, we constructed finite element models of polycrystalline h-BN sheets to probe the thermal conductivity of samples with larger grain sizes. Our multiscale investigations not only provide a general viewpoint regarding the heat conduction in h-BN films but also propose that polycrystalline h-BN sheets present high thermal conductivity comparable to monocrystalline sheets.
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.
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.
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...
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.
Hydro- und morphodynamischen Prozesse in Binnengewässern und im Küstennahbereich erzeugen hochkomplexe Phänomene. Zur Beurteilung der Entwicklung von Küstenzohnen, von Flussbetten sowie von Eingriffen des Menschen in Form von Schutzbauwerken sind geeignete numerische Modellwerkzeuge notwendig. Es wird ein holistischer Modellansatz zur Approximation gekoppelter Seegangs-, Strömungs- und Morphodynamischer Prozesse auf der Basis stabilisierter Finiter Elemente vorgestellt. Der Großteil der Modellgleichungen der Hydro- und Morphodynamik sind Transportgleichungen. Dem Transportcharakter dieser Gleichungen entsprechend wird ein stabilisiertes Finites Element Verfahren auf Dreiecken vorgestellt. Die vorgestellte Approximation entspricht einem streamline upwinding Petrov-Galerkin-Verfahrens für vektorwertige mehrdimensionale Probleme, bei dem der Fehler eines Standard-Galerkin-Verfahrens mit Hilfe eines Upwinding-Koeffizienten minimiert wird. Die Wahl des Upwinding-Koeffizienten ist übertragbar auf andere Problemklassen und basiert ausschließlich auf dem Charakter der zugrundeliegene Das Modell wurde für Seegangs- und Strömungs-Untersuchungen im Jade-Weser-Ästuar an der deutschen Nordseeküste eingesetzt.
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.
The Finite Element Method (FEM) is widely used in engineering for solving Partial Differential Equations (PDEs) over complex geometries. To this end, it is required to provide the FEM software with a geometric model that is typically constructed in a Computer-Aided Design (CAD) software. However, FEM and CAD use different approaches for the mathematical description of the geometry. Thus, it is required to generate a mesh, which is suitable for FEM, based on the CAD model. Nonetheless, this procedure is not a trivial task and it can be time consuming. This issue becomes more significant for solving shape and topology optimization problems, which consist in evolving the geometry iteratively. Therefore, the computational cost associated to the mesh generation process is increased exponentially for this type of applications.
The main goal of this work is to investigate the integration of CAD and CAE in shape and topology optimization. To this end, numerical tools that close the gap between design and analysis are presented. The specific objectives of this work are listed below:
• Automatize the sensitivity analysis in an isogeometric framework for applications in shape optimization. Applications for linear elasticity are considered.
• A methodology is developed for providing a direct link between the CAD model and the analysis mesh. In consequence, the sensitivity analysis can be performed in terms of the design variables located in the design model.
• The last objective is to develop an isogeometric method for shape and topological optimization. This method should take advantage of using Non-Uniform Rational B-Splines (NURBS) with higher continuity as basis functions.
Isogeometric Analysis (IGA) is a framework designed to integrate the design and analysis in engineering problems. The fundamental idea of IGA is to use the same basis functions for modeling the geometry, usually NURBS, for the approximation of the solution fields. The advantage of integrating design and analysis is two-fold. First, the analysis stage is more accurate since the system of PDEs is not solved using an approximated geometry, but the exact CAD model. Moreover, providing a direct link between the design and analysis discretizations makes possible the implementation of efficient sensitivity analysis methods. Second, the computational time is significantly reduced because the mesh generation process can be avoided.
Sensitivity analysis is essential for solving optimization problems when gradient-based optimization algorithms are employed. Automatic differentiation can compute exact gradients, automatically by tracking the algebraic operations performed on the design variables. For the automation of the sensitivity analysis, an isogeometric framework is used. Here, the analysis mesh is obtained after carrying out successive refinements, while retaining the coarse geometry for the domain design. An automatic differentiation (AD) toolbox is used to perform the sensitivity analysis. The AD toolbox takes the code for computing the objective and constraint functions as input. Then, using a source code transformation approach, it outputs a code for computing the objective and constraint functions, and their sensitivities as well. The sensitivities obtained from the sensitivity propagation method are compared with analytical sensitivities, which are computed using a full isogeometric approach.
The computational efficiency of AD is comparable to that of analytical sensitivities. However, the memory requirements are larger for AD. Therefore, AD is preferable if the memory requirements are satisfied. Automatic sensitivity analysis demonstrates its practicality since it simplifies the work of engineers and designers.
Complex geometries with sharp edges and/or holes cannot easily be described with NURBS. One solution is the use of unstructured meshes. Simplex-elements (triangles and tetrahedra for two and three dimensions respectively) are particularly useful since they can automatically parameterize a wide variety of domains. In this regard, unstructured Bézier elements, commonly used in CAD, can be employed for the exact modelling of CAD boundary representations. In two dimensions, the domain enclosed by NURBS curves is parameterized with Bézier triangles. To describe exactly the boundary of a two-dimensional CAD model, the continuity of a NURBS boundary representation is reduced to C^0. Then, the control points are used to generate a triangulation such that the boundary of the domain is identical to the initial CAD boundary representation. Thus, a direct link between the design and analysis discretizations is provided and the sensitivities can be propagated to the design domain.
In three dimensions, the initial CAD boundary representation is given as a collection of NURBS surfaces that enclose a volume. Using a mesh generator (Gmsh), a tetrahedral mesh is obtained. The original surface is reconstructed by modifying the location of the control points of the tetrahedral mesh using Bézier tetrahedral elements and a point inversion algorithm. This method offers the possibility of computing the sensitivity analysis using the analysis mesh. Then, the sensitivities can be propagated into the design discretization. To reuse the mesh originally generated, a moving Bézier tetrahedral mesh approach was implemented.
A gradient-based optimization algorithm is employed together with a sensitivity propagation procedure for the shape optimization cases. The proposed shape optimization approaches are used to solve some standard benchmark problems in structural mechanics. The results obtained show that the proposed approach can compute accurate gradients and evolve the geometry towards optimal solutions. In three dimensions, the moving mesh approach results in faster convergence in terms of computational time and avoids remeshing at each optimization step.
For considering topological changes in a CAD-based framework, an isogeometric phase-field based shape and topology optimization is developed. In this case, the diffuse interface of a phase-field variable over a design domain implicitly describes the boundaries of the geometry. The design variables are the local values of the phase-field variable. The descent direction to minimize the objective function is found by using the sensitivities of the objective function with respect to the design variables. The evolution of the phase-field is determined by solving the time dependent Allen-Cahn equation.
Especially for topology optimization problems that require C^1 continuity, such as for flexoelectric structures, the isogeometric phase field method is of great advantage. NURBS can achieve the desired continuity more efficiently than the traditional employed functions. The robustness of the method is demonstrated when applied to different geometries, boundary conditions, and material configurations. The applications illustrate that compared to piezoelectricity, the electrical performance of flexoelectric microbeams is larger under bending. In contrast, the electrical power for a structure under compression becomes larger with piezoelectricity.
Parallele Netzgenerierung
(1997)
Bei der Berechnung von statischen oder dynamischen Problemen mit Hilfe der Methode der Finiten Elemente ist eine Diskretisierung des zu berechnenden Gebietes notwendig. Bei einer sinnvollen Modellierung des Gebietes ist die Elementgröße meist nicht konstant, sondern ist an kritischen Stellen kleiner. Die Vorgaben hierfür können einerseits aus Erfahrungen des Anwenders, andererseits aus einer Fehlerabschätzung einer vorangegangenen FE-Berechnung resultieren [5]. Soll die FE-Berechnung auf einem Parallelrechner geschehen, ist eine Partitionierung des Gebietes, d.h. eine Zuordnung der Elemente zu den Prozessoren, notwendig. Bei dem hier beschriebenen Ansatz werden nun im Gegensatz zu den üblichen Verfahren erst die Eingangsdaten für den Netzgenerator umgewandelt und dann das Elementnetz direkt auf dem Parallelrecher gleichzeitig auf allen Prozessoren erzeugt. Eine Aufteilung der Elemente auf die Prozessoren entsteht als Nebenprodukt der Netzaufteilung. Die entstehenden Teilgebietsgrenzen werden geometrisch minimiert. Die Lastbalance der Netzaufteilung sowie der FE-Rechnung wird durch ein annähernd gleiche Anzahl der Elemente je Partition gewährleistet. Als Eingabedaten wird eine Beschreibung des Gebietes durch Polygonzüge, sowie einer Netzdichtefunktion, z.B. durch Punkte mit Angaben über die angestrebte Elementgröße, benötigt.
Auf der Basis der Literaturrecherche wird in dieser Arbeit eine 5-lagige MAG-geschweißte Stumpfnaht an austenitisch-ferritischen Stahl X2CrNiMoN22-5-3 (Duplex-Stahl 1.4462) mit dem FE-Programm „SYSWELD®“ simuliert. Die Berech-nungen der Temperaturfelder werden unter der Berücksichtigung sowohl von tempe-raturunabhängigen als auch temperaturabhängigen thermophysikalischen Material-eigenschaften am drei-dimensionalen und zwei-dimensionalen Modell durchgeführt. Die berechneten Temperatur-Zeit-Verläufe und Gefügeumwandlungen beim MAG-Schweißen der Stumpfnaht werden hinsichtlich der Einflüsse und Veränderun-gen analysiert und die ermittelten Abkühlzeiten t12/8 werden für jede Schweißlage bewertet. Anschließend werden die Berechnungen des Eigenspannungszustandes für einzelne Schweißlagen untersucht.
In der täglichen Ingenieurpraxis werden in zunehmenden Maße numerische Analysen im Rahmen der Finite-Elemente-Methode auch zur Untersuchung stabilitätsgefährdeter Strukturen eingesetzt. Für die aktuelle Praxis, insbesondere im konstruktiven Stahlbau, ist jedoch festzustellen, dass zwischen der fortgeschrittenen Theorie und dem Niveau der praktischen Anwendung numerischer Stabilitätsanalysen eine große Kluft besteht. Aus praktischer Sicht erscheint es unumgänglich, die weiter wachsende Diskrepanz zwischen den umfangreichen theoretischen Möglichkeiten und der gegenwärtigen Praxis abzubauen. Damit steht der praktisch tätige Ingenieur vor der Aufgabe, sein Wissen auf dem Gebiet numerischer Stabilitätsanalysen zu vertiefen und bereits vorhandene FE-Programme um Berechnungsalgorithmen für umfassende numerische Stabilitätsanalysen zu erweitern. Dafür werden in der Arbeit die Grundlagen einer FEM- orientierten modernen Stabilitätstheorie einheitlich und aus Sicht einer praktischen Anwendung aufbereitet. Die Darstellung von realisierten programmtechnischen Umsetzungen für erweiterte Analysenmethoden wie Nachbeulanalysen, Pfadwechsel und Approximationen imperfekter Pfade ermöglicht eine Erweiterung des Methodenvorrates. Die innerhalb der Arbeit untersuchten Beispiele zeigen, dass durch die Anwendung der behandelten Verfahren das Tragverhalten einer stabilitätsgefährdeten Struktur wesentlich besser eingeschätzt werden kann als bei Beschränkung auf die herkömmlichen Analysemethoden.
In this paper we present a theoretical background for a coupled analytical–numerical approach to model a crack propagation process in two-dimensional bounded domains. The goal of the coupled analytical–numerical approach is to obtain the correct solution behaviour near the crack tip by help of the analytical solution constructed by using tools of complex function theory and couple it continuously with the finite element solution in the region far from the singularity. In this way, crack propagation could be modelled without using remeshing. Possible directions of crack growth can be calculated through the minimization of the total energy composed of the potential energy and the dissipated energy based on the energy release rate. Within this setting, an analytical solution of a mixed boundary value problem based on complex analysis and conformal mapping techniques is presented in a circular region containing an arbitrary crack path. More precisely, the linear elastic problem is transformed into a Riemann–Hilbert problem in the unit disk for holomorphic functions. Utilising advantages of the analytical solution in the region near the crack tip, the total energy could be evaluated within short computation times for various crack kink angles and lengths leading to a potentially efficient way of computing the minimization procedure. To this end, the paper presents a general strategy of the new coupled approach for crack propagation modelling. Additionally, we also discuss obstacles in the way of practical realisation of this strategy.
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.
We propose an enhanced iterative scheme for the precise reconstruction of piezoelectric material parameters from electric impedance and mechanical displacement measurements. It is based on finite-element simulations of the full three-dimensional piezoelectric equations, combined with an inexact Newton or nonlinear Landweber iterative inversion scheme. We apply our method to two piezoelectric materials and test its performance. For the first material, the manufacturer provides a full data set; for the second one, no material data set is available. For both cases, our inverse scheme, using electric impedance measurements as input data, performs well.
Wirklichkeitsnahe Erfassung und Beschreibung des Trag- und Verformungsverhaltens von Strukturen baulicher Anlagen hat in den letzten Jahrzehnten ständig an Bedeutung gewonnen. Konstruktionen im Hoch- und Industriebau werden zunehmend multifunktional genutzt - die >Grenzen< zwischen Bauwerk und Tragwerk, zwischen Hüll- und Tragkonstruktion lösen sich auf. Werden raumabschließende Elemente (Wände, Decken, Dächer) gleichzeitig als Tragelemente und wärme- und schalldämmende Konstruktionen ausgeführt, so entstehen beispielsweise Sandwichplatten, deren Schichten sehr stark differierende Materialeigenschaften aufweisen. Beim Aufbau des FEM-Modells für vielschichtige Schalen können die Formänderungshypothesen für jede Schicht einzeln als auch für die Schale insgesamt gegeben werden. Im ersten Fall ist der Knotenfreiheitsgrad von der Schichtenzahl abhängig, im zweiten Fall nicht. Im weiteren wird eine Formänderungshypothese für das Schichtenpaket angenommen. Ausgegangen wird von den Gleichungen der 3D-Elastizitätstheorie. Die Berücksichtigung der Querkraftschubverformungen ergibt die Möglichkeit einer adäquaten Beschreibung der Verformungen sowohl dünner Schalen als auch von Schalen mittlerer Dicke; die Berechnung der Krümmungen und der LAMEschen Parameter der Bezugsfläche zu umgehen, was für komplizierte Schalenformen eine selbständige Aufgabe ist; eines natürlichen Übergangs von homogenen zu geschichteten Schalen. Das vielschichtige isoparametrische Schalen-FE wird vorgestellt, seine Implementierung in das in Entwicklung befindliche Programmsystem SLANG wird vorbereitet.
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.