Search results

In recent years, developments in the field of lightweight armour have been of primary importance to the defence industry. This necessity has led to many organisations adopting composite armours comprising both the traditional heavy armours and new lighter weight ceramic armours. The numerical modelling of metal based armour systems has been well documented over the years using purely continuum based methods; and also the modelling of brittle systems using discrete element methods, therefore it is the objective of this paper to demonstrate how a coupled finite and discrete element approach, can be used in the further understanding of the quantitative response of ceramic systems when subjected to dynamic loadings using a combination of adaptive continuum techniques and discrete element methods. For the class of problems encountered within the defence industry, numerical modelling has suffered from one principal weakness; for many applications the associated deformed finite element mesh can no longer provide an accurate description of the deformed material, whether this is due to large ductile deformation, or for the case of brittle materials, degradation into multiple bodies. Subsequently, two very different approaches have been developed to combat such deficiencies, namely the use of adaptive remeshing for the ductile type materials and a discrete fracture insertion scheme for the modelling of material degradation. Therefore, one of the primary objectives of this paper is to present examples demonstrating the potential benefits of explicitly coupling adaptive remeshing methods to the technique of discrete fracture insertion in order to provide an adaptive discontinuous solution strategy, which is computationally robust and efficient.

This paper gives a bibliographical review of the finite element and boundary element parallel processing techniques from the theoretical and application points of view. Topics include: theory – domain decomposition/partitioning, load balancing, parallel solvers/algorithms, parallel mesh generation, adaptive methods, and visualization/graphics; applications – structural mechanics problems, dynamic problems, material/geometrical non‐linear problems, contact problems, fracture mechanics, field problems, coupled problems, sensitivity and optimization, and other problems; hardware and software environments – hardware environments, programming techniques, and software development and presentations. The bibliography at the end of this paper contains 850 references to papers, conference proceedings and theses/dissertations dealing with presented subjects that were published between 1996 and 2002.

The purpose of this paper is to develop a method to model entire structures on a large scale, at the same time taking into account localized non-linear phenomena of the discrete microstructure of cohesive-frictional materials.

Finite element (FEM) based continuum methods are generally considered appropriate as long as solutions are smooth. However, when discontinuities like cracks and fragmentation appear and evolve, application of models that take into account (evolving) microstructures may be advantageous. One popular model to simulate behavior of cohesive-frictional materials is the discrete element method (DEM). However, even if the microscale is close to the macroscale, DEMs are computationally expensive and can only be applied to relatively small specimen sizes and time intervals. Hence, a method is desirable that combines efficiency of FEM with accuracy of DEM by adaptively switching from the continuous to the discrete model where necessary.

An existing method which allows smooth transition between discrete and continuous models is the quasicontinuum method, developed in the field of atomistic simulations. It is taken as a starting point and its concepts are extended to applications in structural mechanics in this paper. The kinematics in the method presented herein is obtained from FEM whereas DEM yields the constitutive behavior. With respect to the constitutive law, three levels of resolution – continuous, intermediate and discrete – are introduced.

The overall concept combines model adaptation with adaptive mesh refinement with the aim to obtain a most efficient and accurate solution.

Gives a bibliographical review of the finite element methods (FEMs) applied for the linear and nonlinear, static and dynamic analyses of basic structural elements from the theoretical as well as practical points of view. The range of applications of FEMs in this area is wide and cannot be presented in a single paper; therefore aims to give the reader an encyclopaedic view on the subject. The bibliography at the end of the paper contains 2,025 references to papers, conference proceedings and theses/dissertations dealing with the analysis of beams, columns, rods, bars, cables, discs, blades, shafts, membranes, plates and shells that were published in 1992‐1995.

The main aim of the current paper is to find a numerical plan for hydraulic fracturing problem with application in extracting natural gases and oil.

First, time discretization is accomplished via Crank-Nicolson and semi-implicit techniques. At the second step, a high-order finite element method using quadratic triangular elements is proposed to derive the spatial discretization. The efficiency and time consuming of both obtained schemes will be investigated. In addition to the popular uniform mesh refinement strategy, an adaptive mesh refinement strategy will be employed to reduce computational costs.

Numerical results show a good agreement between the two schemes as well as the efficiency of the employed techniques to capture acceptable patterns of the model. In central single-crack mode, the experimental results demonstrate that maximal values of displacements in x- and y- directions are 0.1 and 0.08, respectively. They occur around both ends of the line and sides directly next to the line where pressure takes impact. Moreover, the pressure of injected fluid almost gained its initial value, i.e. 3,000 inside and close to the notch. Further, the results for non-central single-crack mode and bifurcated crack mode are depicted. In central single-crack mode and square computational area with a uniform mesh, computational times corresponding to the numerical schemes based on the high order finite element method for spatial discretization and Crank-Nicolson as well as semi-implicit techniques for temporal discretizations are 207.19s and 97.47s, respectively, with 2,048 elements, final time T = 0.2 and time step size τ = 0.01. Also, the simulations effectively illustrate a further decrease in computational time when the method is equipped with an adaptive mesh refinement strategy. The computational cost is reduced to 4.23s when the governed model is solved with the numerical scheme based on the adaptive high order finite element method and semi-implicit technique for spatial and temporal discretizations, respectively. Similarly, in other samples, the reduction of computational cost has been shown.

This is the first time that the high-order finite element method is employed to solve the model investigated in the current paper.

Gives a bibliographical review of the error estimates and adaptive finite element methods from the theoretical as well as the application point of view. The bibliography at the end contains 2,177 references to papers, conference proceedings and theses/dissertations dealing with the subjects that were published in 1990‐2000.

The simulation of fracture processes for discrete crack propagation is well established for linear‐elastic cracking problems. Applying finite element techniques for the numerical formulation, at every incremental macro‐crack step the element mesh has to be adapted such that the crack path remains independent of the initial mesh. The accuracy of the obtained results has to be controlled by suitable error estimators and error indicators. Considering the dependence of the predicted crack path on the precision of the displacement and stress computation, quality measures for the computed results are recommended. In this research the use of the Babuska/Rheinboldt error indicator in combination with linear‐elastic crack propagation problems is demonstrated. Based on this error measure an adaptive mesh refinement technique is developed. In comparison with classical discrete crack propagation simulations the advantages of the new concept can be clearly observed.

Gives introductory remarks about chapter 1 of this group of 31 papers, from ISEF 1999 Proceedings, in the methodologies for field analysis, in the electromagnetic community. Observes that computer package implementation theory contributes to clarification. Discusses the areas covered by some of the papers ‐ such as artificial intelligence using fuzzy logic. Includes applications such as permanent magnets and looks at eddy current problems. States the finite element method is currently the most popular method used for field computation. Closes by pointing out the amalgam of topics.

This paper discusses the issues involved in the development of combined finite/discrete element methods; both from a fundamental theoretical viewpoint and some related algorithmic considerations essential for the efficient numerical solution of large scale industrial problems. The finite element representation of the solid region is combined with progressive fracturing, which leads to the formation of discrete elements, which may be composed of one or more deformable finite elements. The applicability of the approach is demonstrated by the solution of a range of examples relevant to various industrial sections.

This paper outlines a dynamic domain decomposition‐based parallel strategy for combined finite/discrete element analysis of multi‐fracturing solids and discrete systems. Attention is focused on the parallelised interaction detection between discrete objects. Two graph representation models for discrete objects in contact are proposed which lay the foundation of the current development. In addition, a load imbalance detection and re‐balancing scheme is also suggested to enhance the parallel performance. Finally, numerical examples are provided to illustrate the parallel performance achieved with the current implementation.