Search results

This paper is concerned with the determination of the transient stress and deformational state of plate‐like discontinua subject to flexural cracking. Such a phenomenon can be easily visualized as the type of fragmentation to floating sea ice impacted by an ice‐breaker or offshore platform. The discrete element method is used to solve the dynamic equilibrium equations for each distinct deformable body and the interaction between bodies. Each body may deform elastically and fracture into further pieces if a brittle failure criterion for flexure is exceeded. The discrete plate element is a hybrid thin‐plate (Kirchhoff) mode lumped at element boundaries with transverse shear deformation computed at element centroids. Errors in computed stresses near point loads and cracks by the current element warrant the use of an improved mixed mode plate element. A three‐dimensional application of the discrete element method is presented for the case of fragmentation of floating sea ice impacting an arctic offshore platform. A semi‐implicit solution scheme is introduced to overcome the stringent explicit time step stability conditions due to stiff members in the discrete element formulation.

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 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 study aims to develop an efficient algorithm for generation of conforming mesh for seepage analysis through 3D discrete fracture networks (DFN).

The algorithm is developed based on a refined conforming Delaunay triangulation scheme, which is then validated using analytical solutions. The algorithm is well able to meet the challenge of meshing complex geometry of DFNs.

A series of sensitivity analysis have been performed to evaluate the effect of meshing parameters on steady state solution of Darcy flow using a finite element scheme. The results show that an optimized minimum internal angle of meshing elements should be predetermined to guarantee termination of the algorithm.

The developed algorithm is computationally efficient, fast and is of low cost. Furthermore, it never changes the geometrical structure and connectivity pattern of the DFN.

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.

Hydrofracturing technology has been widely used in tight oil and gas reservoir exploitation, and the fracture network formed by fracturing is crucial to determining the resources recovery rate. Due to the complexity of fracture network induced by the random morphology and type of fluid-driven fractures, controlling and optimising its mechanisms is challenging. This paper aims to study the types of multiscale mode I/II fractures, the fluid-driven propagation of multiscale tensile and shear fractures need to be studied.

A dual bilinear cohesive zone model (CZM) based on energy evolution was introduced to detect the initiation and propagation of fluid-driven tensile and shear fractures. The model overcomes the limitations of classical linear fracture mechanics, such as the stress singularity at the fracture tip, and considers the important role of fracture surface behaviour in the shear activation. The bilinear cohesive criterion based on the energy evolution criterion can reflect the formation mechanism of complex fracture networks objectively and accurately. Considering the hydro-mechanical (HM) coupling and leak-off effects, the combined finite element-discrete element-finite volume approach was introduced and implemented successfully, and the results showed that the models considering HM coupling and leak-off effects could form a more complex fracture network. The multiscale (laboratory- and engineering-scale) Mode I/II fractures can be simulated in hydrofracturing process.

Based on the proposed method, the accuracy and applicability of the algorithm were verified by comparing the analytical solution of KGD and PKN models. The effects of different in situ stresses and flow rates on the dynamic propagation of hydraulic fractures at laboratory and engineering scales were investigated. when the ratio of in situ stress is small, the fracture propagation direction is not affected, and the fracture morphology is a cross-type fracture. When the ratio of in situ stress is relatively large, the propagation direction of the fracture is affected by the maximum in situ stress, and it is more inclined to propagate along the direction of the maximum in situ stress, forming double wing-type fractures. Hydrofracturing tensile and shear fractures were identified, and the distribution and number of each type were obtained. There are fewer hydraulic shear fractures than tensile fractures, and shear fractures appear in the initial stage of fracture propagation and then propagate and distribute around the perforation.

The proposed dual bilinear CZM is effective for simulating the types of Mode I/II fractures and seizing the fluid-driven propagation of multiscale tensile and shear fractures. Practical fracturing process involves the multi-type and multiscale fluid-driven fracture propagation. This study introduces general fluid-driven fracture propagation, which can be extended to the fracture propagation analysis of potential fluid fracturing, such as other liquids or supercritical gases.

Continua and discontinua coexist in natural rock materials. This paper aims to present an improved approach for addressing the mechanical response of rock masses based on the combined finite-discrete element method (FDEM) proposed by Munjiza.

Several algorithms have been programmed in the new approach. The algorithms include (1) a simpler and more efficient algorithm to calculate the contact force; (2) An algorithm for tangential contact force closer to the actual physical process; (3) a plastic yielding criterion (e.g. Mohr-Coulomb) to modify the elastic stress for fitting the mechanical behavior of elastoplastic materials; and (4) a complete code for the mechanical calculation to be implemented in Matrix Laboratory (MATLAB).

Three case studies, including two standard laboratory experiments (uniaxial compression and Brazilian split test) and one engineering-scale anti-dip slop model, are presented to illustrate the feasibility of the Y-Mat code and its ability to deal with multi-scale rock mechanics problems. The results, including the progressive failure process, failure mode and trajectory of each case, are acceptable compared to other corresponding studies. It is shown that, the code is capable of modeling geotechnical and geological engineering problems.

This article gives an improved FDEM-based numerical calculation code. And, feasibility of the code is verified through three cases. It can effectively solve the geotechnical and geological engineering problems.

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.

This paper aims to investigate the responses of laminated glass under soft body impact, including elastic impact and fracture/fragmentation consideration.

The simulation uses the combined finite-discrete element method (FDEM) which combines finite element mesh into discrete elements, enabling the accurate prediction of contact force and deformation. Material rupture is modelled with a cohesive fracture criterion, evaluating the process from continua to discontinua.

Responses of laminated glass under soft impact (both elastic and fracture) agree well with known data. Crack initiation time in laminated glass increases with the increase of the outside glass thickness. With the increase of Eprojectile, failure mode is changing from flexural to shear, and damage tends to propagate longitudinally when the contact surface increases. Results show that the FDEM is capable of modelling soft impact behaviour of laminated glass successfully.

The work is done in 2D, and it will not represent fully the 3D mechanisms.

Elastic and fracture behaviour of laminated glass under soft impact is simulated using the 2D FDEM. Limited work has been done on soft impact of laminated glass with FDEM, and special research endeavours are warranted. Benchmark examples and discussions are provided for future research.

Contact problems are among the most difficult ones in mechanics. Due to its practical importance, the problem has been receiving extensive research work over the years. The finite element method has been widely used to solve contact problems with various grades of complexity. Great progress has been made on both theoretical studies and engineering applications. This paper reviews some of the main developments in contact theories and finite element solution techniques for static contact problems. Classical and variational formulations of the problem are first given and then finite element solution techniques are reviewed. Available constraint methods, friction laws and contact searching algorithms are also briefly described. At the end of the paper, a bibliography is included, listing about seven hundred papers which are related to static contact problems and have been published in various journals and conference proceedings from 1976.