Object oriented implementation of distributed finite element analysis in .NET

The paper describes a detailed study into the object-oriented implementation of distributed finite element analysis on desktop computers using the .NET framework. The software design aspects are described in some detail for both direct and iterative solution algorithms. The use of interfaces played an important role in the software design. This, together with the .NET framework, enabled remote objects to be implemented in a relatively seamless fashion. The solution routines were “blind” to whether the objects were local or remote. Numerical tests were carried out and reasonable speed-up was achieved, particularly for direct solution methods. It is concluded that .NET provides a viable framework for implementing distributed computing on networks of personal computers.     

Distributed parameter system identification using finite element differential neural networks

Most of the previous work on identification involves systems described by ordinary differential equations (ODEs). Many industrial processes and physical phenomena, however, should be modeled using partial differential equations (PDEs) which offer both spatial and temporal distributions that are simply not available with ODE models. Systems described by a PDE belong to a class of system called distributed parameter system (DPS). This article presents a method for solving the problem of identification of uncertain DPSs using a differential neural network (DNN). The DPS, assumed to be described by a PDE, is approximated using the finite element method (FEM). The FEM discretizes the domain into a set of distributed and connected nodes, thereby, allowing a representation of the DPS in a finite number of ODEs. The proposed DNN follows the same interconnection structure of the FEM, thus allowing the DNN to identify the FEM approximation of the DPS in both 2D and 3D domains. Lyapunov's second method was used to derive adaptive learning laws for the proposed DNN structure. The identification algorithm, here developed in Nvidia's CUDA/C to reduce the execution time, runs mostly on the graphics processing unit (GPU). A physical experiment served to validate the 2D case. In the experiment, the DNN followed the trajectory of 57 markers that were placed on an undulating square piece of silk. The proposed DNN is compared against a method based on principal component analysis and an artificial neural network trained with group search optimization. In addition to the 2D case, a simulation validated the 3D case, where input data for the DNN was generated by solving a PDE with appropriate initial and boundary conditions over an unitary domain. Results show that the proposed FEM-based DNN approximates the dynamic behavior of both a real 2D and a simulated 3D system.     

EasyFEM—An object-oriented graphics interface finite element/finite volume software

An object-oriented finite element/finite volume software, EasyFEM, has been developed. The software, with a fully interactive graphics interface, analyses heat transfer, solid mechanics, and fluids problems by the finite element and the finite volume methods. The Coad and Yourdon methodology is used to describe the general structure of classes and objects implemented by the software. Details of the object-oriented classes for both the graphics pre- and post-processors, as well as for the analysis solutions, are described. Several case studies with graphical representations of numerical solutions are presented to illustrate some functionalities of the software.     

An object-oriented framework for finite element pavement analysis

In this study, we developed an object-oriented (OO) framework with interactive graphics to assist pavement studies using finite element analysis (FEA). FEA has been proven to be effective in studying various pavement failure problems; however, it is time consuming and error prone to manually generate the load sequences where non-regular tire footprints, non-uniform tire-pavement contact stresses, and transverse wheel wander distributions are used. After FEA, extracting the deformations for failure analysis is necessary but tedious. The OO framework developed in this study handles the preprocessing and postprocessing tasks for the FEA of pavements. It has a graphical user interface and is platform independent. It was successfully used in developing a new criterion for characterizing pavement failures that involved approximately four hundred different FEA simulations.     

Tuning the Schur complement computations for finite element partitions

The domain decomposition method (DDM) is an efficient algorithmic tool for the parallelization of finite element computer codes. A variant of the DDM with direct solution algorithm is based on computation of Schur complement matrices for finite element partitions. This paper describes a simple technique that considerably improves execution rate of computationally intensive routines of the Schur complement computations. The technique uses ‘block of columns’ matrix operations and loop unrolling to reduce load instructions from cache memory and to increase instruction-level parallelism. For superscalar RISC processors, experimental results show that it is possible to improve performance of the DDM solution procedure by several times.     

An object-oriented blackboard based approach for automated finite element modeling and analysis of multichip modules

This paper addresses the application of the blackboard system architecture and object-oriented data abstraction techniques to the domain of finite element modeling and analysis. Specifically, a hierarchical object-oriented database was used to represent the physical system at different levels of abstraction including the user-defined physical system level, a computer-generated, simplified physical model level, and the finite element model level. Object link relationships within a given abstraction level and across different abstraction levels resulted in seamless bidirectional information exchange. The blackboard system architecture employed provided a framework for distributed cooperative problem solving, for the application of simplifying domain-specific modeling assumptions, and for integrating the various software modules that are involved in the entire finite element modeling and analysis process. These methodologies were implemented in a prototype computational tool calledIMCMA theIntelligentMultichipModuleAnalyzer. An example illustrates howIMCMA automates finite element thermal analysis of small integrated circuit features in multichip modules through a two-step finite element submodeling process.     

An object-oriented framework for finite element analysis based on a compact topological data structure

This paper describes an ongoing work in the development of a finite element analysis system, called TopFEM, based on the compact topological data structure, TopS [1], [2]. This new framework was written to take advantage of the topological data structure together with object-oriented programming concepts to handle a variety of finite element problems, spanning from fracture mechanics to topology optimization, in an efficient, but generic fashion. The class organization of the TopFEM system is described and discussed within the context of other frameworks in the literature that share similar ideas, such as GetFEM++, deal.II, FEMOOP and OpenSees. Numerical examples are given to illustrate the capabilities of TopS attached to a finite element framework in the context of fracture mechanics and to establish a benchmark with other implementations that do not make use of a topological data structure.     

Evaluation of distributed finite element algorithms on a workstation network

This paper discusses the design, implementation and evaluation of linear finite element programs that distribute their computations over a network of workstations. We consider five different algorithms based on direct, iterative and hybrid equation solvers, each of which partitions and maps the model domain across conventional network hardware. A software architecture based on the client-server model distributes the computations and, at the language level, Berkeley sockets enable communication between processes. We evaluate and describe the performance of these algorithms in terms of execution time and speed-up, and we conclude that distributed solvers, particularly those based on substructuring and static condensation, can be effective even on high-latency communication networks.     

Collaborative object-oriented visualization environment

In this paper, we present a Collaborative Object-oriented Visualization Environment (COVE) which provides a flexible and extensible framework for collaborative visualization. COVE integrates collaborative and parallel computing environments based on a distributed object model. It is built as a collection of concurrent objects: collaborative and application objects which interact with one another to construct collaborative parallel computing environments. The former enables COVE to execute various collaborative functions, while the latter allows it to execute fast parallel visualization in various modes. Also, flexibility and extensibility are provided by plugging the proper application objects into COVE at run-time, and making them interact with one another through collaboration objects. For our experiment, three visualization modes for volume rendering are designed and implemented to support the fast and flexible analysis of volume data in a collaborative environment. This work has been supported by KIPA-Information Technology Research Center, University research program by Ministry of Information & Communication, and Brain Korea 21 projects in 2005.     

A parallel preconditioned conjugate gradient solution method for finite element problems with coarse–fine mesh formulation

The object of this paper is a parallel preconditioned conjugate gradient iterative solver for finite element problems with coarse-mesh/fine-mesh formulation. An efficient preconditioner is easily derived from the multigrid stiffness matrix. The method has been implemented, for the sake of comparison, both on a IBM-RISC590 and on a Quadrics-QH1, a massive parallel SIMD machine with 128 processors. Examples of solutions of simple linear elastic problems on rectangular grids are presented and convergence and parallel performance are discussed.     

A parallel finite element procedure for contact-impact problems

An efficient parallel finite element procedure for contact-impact problems is presented within the framework of explicit finite element analysis with thepenalty method. The procedure concerned includes a parallel Belytschko-Lin-Tsay shell element generation algorithm and a parallel contact-impact algorithm based on the master-slave slideline algorithm. An element-wise domain decomposition strategy and a communication minimization strategy are featured to achieve almost perfect load balancing among processors and to show scalability of the parallel performance. Throughout this work, a prototype code, named GT-PARADYN, is developed on the IBM SP2 to implement the procedure presented, under message-passing paradigm. Some examples are provided to demonstrate the timing results of the algorithms, discussing the accuracy and efficiency of the code.     

An implementation of the object-oriented concurrent programming language SINA

SINA is an object-oriented language for distributed and concurrent programming. The primary focus of this paper is on the object-oriented concurrent programming mechanisms of SINA and their implementation. This paper presents the SINA constructs for concurrent programming and inter-object communication, some illustrative examples and a message-based implementation model for SINA that we have used in our current implementation.     

Distributed evolutionary multi-objective mesh-partitioning algorithm for parallel finite element computations

Majority of the mesh-partitioning algorithms attempt to optimise the interprocessor communications, while balancing the computational load among the processors. However, it is desirable to simultaneously optimise the submesh aspect ratios in order to significantly improve the convergence characteristics of the domain decomposition based Preconditioned-conjugate-gradient algorithms, being used extensively in the state-of-the-art parallel finite element codes. Keeping this in view, a new distributed multi-objective mesh-partitioning algorithm using evolutionary computing techniques is proposed in this paper. Effectiveness of the proposed distributed mesh-partitioning algorithm is demonstrated by solving several unstructured meshes of practical-engineering problems and also benchmark problems.     

Implementation of a strain energy-based nonlinear finite element in the object-oriented environment

The objective of the paper is to describe a novel finite element computational method based on a strain energy density function and to implement it in the object-oriented environment. The original energy-based finite element was put into the known standard framework of classes and handled in a different manner. The nonlinear properties of material are defined with a modified strain energy density function. The local relaxation procedure proposed as a method used to resolve a nonlinear problem is implemented in C++ language. The hexahedral element with eight nodes as well as the adaptation of the nonlinear finite element is introduced. The chosen numerical model is made of nearly incompressible hyperelastic material. The application of the proposed element is shown on the example of a rectangular parallelepiped with a hollow port.     

Adaptive full domain covering meshes for parallel finite element computations

Summary  In this work, we present a new parallelization concept for adaptive finite element methods. Compared to classical domain decomposition approaches, the concept of adaptive full domain covering meshes reduces the parallel communication overhead. Furthermore, it provides an easy way to transform sequential codes into parallel software by changing only a few lines of source code.      

Generic programming techniques for parallelizing and extending procedural finite element programs

We outline an approach for extending procedural finite-element software components using generic programming. A layer of generic software components consisting of C++ containers and algorithms is used for parallelization of the finite-element solver and for solver coupling in multi-physics applications. The advantages of generic programming in connection with finite-element codes are discussed and compared with those of object-oriented programming. The use of the proposed generic programming techniques is demonstrated in a tutorial fashion through basic illustrative examples as well as code excerpts from a large-scale finite-element program for serial and parallel computing platforms.

Localized tensor-based solvers for interactive finite element applications using C++ and Java

This paper presents techniques for solving systems of equations arising in finite element applications using a localized, tensor-based approach. The approach is localized in that a major part of the solution responsibility is delegated to vector degree-of-freedom objects, rather than residing solely in a global solver working on a monolithic data structure. The approach is tensor-based in that the fundamental quantities used for computation are considered to be second-order tensors. The localized data structure strategy provides the benefits of an efficient sparse and symmetric storage scheme without requiring complex implementation code. The tensor-based aspect of the approach can bring substantial performance benefits by increasing the granularity at which solution algorithms deal with their data. Java and C++ implementations of interactive finite element programs are used to demonstrate performance that is competitive with other available solvers, especially in the case of problems for which interactive analysis is feasible on commonly available hardware.     

分布式计算环境中面向对象的目录服务

为了对分布式计算环境中资源的统一管理，提出一种基于开放式系统目录推荐标准Ｘ．５００的面向对象的资源描述和管理方法，并对面向对象的分布式数据库中的目录服务的实现作了详细的讨论。     

Large deflection analyses of skew plates using hybrid/mixed finite element method

A new hybrid/mixed shell element is developed using oblique coordinate systems to analyze the large deflection behavior of skew plate with various skew angles, length to width ratios, thicknesses and supported edges under uniformly distributed and concentrated loads. The results obtained from the new element are compared with available theoretical and numerical solutions. An excellent agreement is achieved even for coarse meshes. The accuracy and efficiency of the proposed element are demonstrated.     

Biomechanical analyses of static and dynamic fixation techniques of retrograde interlocking femoral nailing using nonlinear finite element methods

Femoral shaft fractures can be treated using retrograde interlocking nailing systems; however, fracture nonunion still occurs. Dynamic fixation techniques, which remove either the proximal or distal locking screws, have been used to solve the problem of nonunion. In addition, a surgical rule for dynamic fixation techniques has been defined based on past clinical reports. However, the biomechanical performance of the retrograde interlocking nailing systems with either the traditional static fixation technique or the dynamic fixation techniques has not been investigated by using nonlinear numerical modeling. Three-dimensional nonlinear finite element models were developed, and the implant strength, fixation stability, and contact area of the fracture surfaces were evaluated. Three types of femoral shaft fractures (a proximal femoral shaft fracture, a middle femoral shaft fracture, and a distal femoral shaft fracture) fixed by three fixation techniques (insertion of all the locking screws, removal of the proximal locking screws, or removal of the distal locking screws) were analyzed. The results showed that the static fixation technique resulted in sufficient fixation stability and that the dynamic fixation techniques decreased the failure risk of the implant and produced a larger contact area of the fracture surfaces. The outcomes of the current study could assist orthopedic surgeons in comprehending the biomechanical performances of both static and dynamic fixation techniques. In addition, the surgeons could also select a fixation technique based on the specific patient situation using the numerical outcomes of this study.