This paper is concerned with the investigation of the nonlinear eigenvalue problem describing the natural oscillations of a plate with a load that elastically attached to it. We study properties of eigenvalues and eigenfunctions of this eigenvalue problem and prove the existence theorem for eigensolutions. Theoretical results are illustrated by numerical experiments.

In the paper we construct exact three-point discretizations of linear and nonlinear two-point boundary value problems with boundary conditions of the first kind. The finite element approach uses basis functions defined by the coefficients of the differential equations. All the discretized boundary value problems are of inverse isotone type and so are its exact discretizations which involve tridiagonal M-matrices in the linear case and M-functions in the nonlinear case.

Introduction In recent years, domain decomposition methods have been used extensively to efficiently solve boundary value problems for partial differential equations in complex{form domains. On the other hand, multilevel techniques on hierarchical data structures also have developed into an effective tool for the construction and analysis of fast solvers. But direct realization of multilevel techniques on a parallel computer system for the global problem in the original domain involves difficult communication problems. I this paper, we present and analyze a combination of these two approaches: domain decomposition and multilevel decomposition on hierarchical structures to design optimal preconditioning operators.

The experimental program ¨SPC-PM Po 3D¨ is part of the ongoing research of the Chemnitz research group Scientific Parallel Computing (SPC) into finite element methods for problems over three dimensional domains. The package in its version 2.0 is documented in two manuals. The User's Manual provides an overview over the program, its capabilities, its installation, and handling. Moreover, test examples are explained. The aim of the Programmer's Manual is to provide a description of the algorithms and their realization. It is written for those who are interested in a deeper insight into the code, for example for improving and extending. In Version 2.0 the program can solve the Poisson equation and the Lam\'e system of linear elasticity with in general mixed boundary conditions of Dirichlet and Neumann type. The domain $\Omega\subset\R^3$ can be an arbitrarily bounded polyhedron. The input is a coarse mesh, a description of the data and some control parameters. The program distributes the elements of the coarse mesh to the processors, refines the elements, generates the system of equations using linear or quadratic shape functions, solves this system and offers graphical tools to display the solution. Further, the behavior of the algorithms can be monitored: arithmetic and communication time is measured, the discretization error is measured, different preconditioners can be compared. We plan to extend the program in the next future by including a multigrid solver, an error estimator and adaptive mesh refinement, as well as the treatment of coupled thermo-elastic problems. The program has been developed for MIMD computers; it has been tested on Parsytec machines (GCPowerPlus-128 with Motorola Power PC601 processors and GCel-192 on transputer basis) and on workstation clusters using PVM. The special case of only one processor is included, that means the package can be compiled for single processor machines without any change in the source files.

Abstract P-version finite element method for the second order elliptic equation in an arbitrary sufficiently smooth domain is studied in the frame of DD method. Two types square reference elements are used with the products of the integrated Legendre's polynomials for the coordinate functions. There are considered the estimates for the condition numbers, preconditioning of the problems arising on subdomains and the Schur complement, the derivation of the DD preconditioner. For the result we obtain the DD preconditioner to which corresponds the generalized condition number of order (logp )2 . The paper consists of two parts. In part I there are given some preliminary re- sults for 1D case, condition number estimates and some inequalities for 2D reference element.

P-version finite element method for the second order elliptic equation in an arbitrary sufficiently smooth domain is studied in the frame of DD method. Two types square reference elements are used with the products of the integrated Legendre's polynomials for the coordinate functions. There are considered the estimates for the condition numbers, preconditioning of the problems arising on subdomains and the Schur complement, the derivation of the DD preconditioner. For the result we obtain the DD preconditioner to which corresponds the generalized condition number of order (logp )2 . The paper consists of two parts. In part I there are given some preliminary results for 1D case, condition number estimates and some inequalities for 2D reference element. Part II is devoted to the derivation of the Schur complement preconditioner and conditionality number estimates for the p-version finite element matrixes. Also DD preconditioning is considered.

Using material models on the basis of the flow theory of plasticity the asymptotic behaviour of solid mechanics solutions in crack tips, interface corners etc. strongly depends on the local realized load trajectory. For incrementally proportional load paths the equations determining the asymptotic fields are very simple ones. The paper considers two-dimensional statements in the neighbourhood of an interface corner consisting of two material ranges. At a distance from the corner the finite element nodes of a regular net are established in a polar co-ordinate system together with the displacement degrees of freedom. The main idea of the presented singular and non-singular stress and deformation field calculation at interface corners characterizes an replacement of the corner neighbourhood effect to the surrounding body by introducing stiffness actions which in usual manner can be assembled together with the other element stiffness matrices to the global stiffness matrix of the body. According to this there exists an in teresting invariant stiffness independence in corner and crack neighbourhoods. The applied technique allows extensions to non-proportional local load increments simplifying the mathematical calculations for the presentation of stress and strain fields in this general case. All computations are made on modern parallel computers. Concrete examples show the advantages of the presented approach.

In this paper we compare the performance, scalability, and robustness of different parallel algorithms for the numerical solution of nonlinear boundary value problems arising in the magnetic field computation and in solid mechanics. These problems are discretized by using the finite element method with triangular meshes and piecewise linear functions. The nonlinearity is handled by a nested Newton solver, and the linear systems of algebraic equations within each Newton step are solved by means of various iterative solvers, namely multigrid methods and conjugate gradient methods with preconditioners based on domain decomposition, multigrid, or BPX techniques, respectively. The basis of the implementation of all solvers is a non-overlapping domain decomposition data structure such that they are well-suited for parallel machines with MIMD architecture.

The realization of the finite element method on parallel computers is usually based on a domain decomposition approach. This paper is concerned with the problem of finding an optimal decomposition and an appropriate mapping of the subdomains to the processors. The quality of this partitioning is measured in several metrics but it is also expressed in the computing time for solving specific systems of finite element equations. The software environment is first described. In particular, the data structure and the accumulation algorithm are introduced. Then several partitioning algorithms are compared. Spectral bisection was used with different modifications including Kernighan-Lin refinement, post-processing techniques and terminal propagation. The final recommendations should give good decompositions for all finite element codes which are based on principles similar to ours. The paper is a shortened English version of Preprint SFB393/96-18 (Uwe Reichel: Partitionierung von Finite-Elemente-Netzen), SFB 393, TU Chemnitz-Zwickau, December 1996. To be selfcontained, some material of Preprint SPC95_5 (see below) is included. The paper appeared as Preprint SFB393/96-18a, SFB 393, TU Chemnitz-Zwickau, January 1997.

Die vorliegende Arbeit befasst sich mit verschiedenen Aspekten der hp-FEM. Insbesondere werden hierbei folgende Punkte genauer untersucht: 1. Das effiziente Aufstellen der Steifigkeitsmatrix auf Referenzelementen, die keine Tensorproduktstruktur besitzen. 2. Eine lokale Konvergenzbetrachtung für die randkonzentrierte hp-FEM. 3. Ein Multilevel-Löser für die randkonzentrierte hp-FEM. 4. Eine Strategie für hp-Adaptivität.