Discrete element and finite element coupling simulation and experiment of hot granule medium pressure forming

The granule medium of discreteness is supposed to be continuous（Drucker-Prager model） in the existing finite element simulation analysis on the hot granule medium pressure forming（HGMF） process, so the granule medium may produce tensile stress in the process of pressure-transferring and flowing, which does not coincide with the reality. The analysis method, discrete element and finite element（DE-FE） coupling simulation, is proposed to solve the problem. The material parameters of simulation model are obtained by the pressure-transfer performance test of granule medium and the hot uniaxial tensile test of sheet metal. The DE-FE coupling simulation platform is established by adopting Visual Basic language. The features in the process that AA7075-T6 conical parts are formed by the HGMF process are analyzed and verified by the process test. The studies show that the results of DE-FE coupling simulation coincide well with the test results, which provides a new analysis method to solve the mechanics problem in the coupling of discrete and continuum.

A New Elasticity and Finite Element Formulation for Different Young's Modulus When Tension and Compression Loadings

This paper presents a new elasticity and finite element formulation for different Young's modulus when tension and compression loadings in anisotropy media. The case studies, such as anisotropy and isotropy, were investigated. A numerical example was shown to find out the changes of neutral axis at the pure bending beams.

Assessing fracturing mechanisms and evolution of excavation damaged zone of tunnels in interlocked rock masses at high stresses using a finitediscrete element approach

Deep underground excavations within hard rocks can result in damage to the surrounding rock mass mostly due to redistribution of stresses.Especially within rock masses with non-persistent joints,the role of the pre-existing joints in the damage evolution around the underground opening is of critical importance as they govern the fracturing mechanisms and influence the brittle responses of these hard rock masses under highly anisotropic in situ stresses.In this study,the main focus is the impact of joint network geometry,joint strength and applied field stresses on the rock mass behaviours and the evolution of excavation induced damage due to the loss of confinement as a tunnel face advances.Analysis of such a phenomenon was conducted using the finite-discrete element method (FDEM).The numerical model is initially calibrated in order to match the behaviour of the fracture-free,massive Lac du Bonnet granite during the excavation of the Underground Research Laboratory (URL) Test Tunnel,Canada.The influence of the pre-existing joints on the rock mass response during excavation is investigated by integrating discrete fracture networks (DFNs) of various characteristics into the numerical models under varying in situ stresses.The numerical results obtained highlight the significance of the pre-existing joints on the reduction of in situ rock mass strength and its capacity for extension with both factors controlling the brittle response of the material.Furthermore,the impact of spatial distribution of natural joints on the stability of an underground excavation is discussed,as well as the potentially minor influence of joint strength on the stress induced damage within joint systems of a non-persistent nature under specific conditions.Additionally,the in situ stress-joint network interaction is examined,revealing the complex fracturing mechanisms that may lead to uncontrolled fracture propagation that compromises the overall stability of an underground excavation.

A Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations

In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.

Multi-scale Simulation Method with Coupled Finite/Discrete Element Model and Its Application

The existing research on continuous structure is usually analyzed with finite element method (FEM) and granular medium with discrete element method (DEM), but there are few researches on the coupling interaction between continuous structure and discrete medium. To the issue of this coupling interaction, a multi-scale simulation method with coupled finite/discrete element model is put forward, in their respective domains of discrete and finite elements, the nodes follow force law and motion law of their own method, and on the their interaction interface, the touch type between discrete and finite elements is distinguished as two types: full touch and partial touch, the interaction force between them is calculated with linear elastic model. For full touch, the contact force is proportional to the overlap distance between discrete element and finite element patch. For partial touch, first the finite element patch is extended on all sides indefinitely to be a complete plane, the full contact force can be obtained with the touch type between discrete element and plane being viewed as full touch, then the full overlap area between them and the actual overlap area between discrete element and finite element patch are computed, the actual contact force is obtained by scaling the full contact force with a factor which is determined by the ratio of the actual overlap area to the full overlap area. The contact force is equivalent to the finite element nodes and the force and displacement on the nodes can be computed, so the ideal simulation results can be got. This method has been used to simulate the cutter disk of the earth pressure balance shield machine (EPBSM) made in North Heavy Industry (NHI) with its excavation diameter of 6.28 m cutting and digging the sandy clay layer. The simulation results show that as the gradual increase of excavating depth of the cutter head, the maximum stress occurs at the roots of cutters on the cutter head, while for the soil, the largest stress is distributed at the region which directly contacted with the cutters. The proposed research can provide good solutions for correct design and installation of cutters, and it is necessary to design mounting bracket to fix cutters on cutter head.

Characterizing the influence of stress-induced microcracks on the laboratory strength and fracture development in brittle rocks using a finite-discrete element method-micro discrete fracture network FDEM-μDFN approach

Heterogeneity is an inherent component of rock and may be present in different forms including mineralheterogeneity, geometrical heterogeneity, weak grain boundaries and micro-defects. Microcracks areusually observed in crystalline rocks in two forms： natural and stress-induced; the amount of stressinducedmicrocracking increases with depth and in-situ stress. Laboratory results indicate that thephysical properties of rocks such as strength, deformability, P-wave velocity and permeability areinfluenced by increase in microcrack intensity. In this study, the finite-discrete element method （FDEM）is used to model microcrack heterogeneity by introducing into a model sample sets of microcracks usingthe proposed micro discrete fracture network （mDFN） approach. The characteristics of the microcracksrequired to create mDFN models are obtained through image analyses of thin sections of Lac du Bonnetgranite adopted from published literature. A suite of two-dimensional laboratory tests including uniaxial,triaxial compression and Brazilian tests is simulated and the results are compared with laboratory data.The FDEM-mDFN models indicate that micro-heterogeneity has a profound influence on both the mechanicalbehavior and resultant fracture pattern. An increase in the microcrack intensity leads to areduction in the strength of the sample and changes the character of the rock strength envelope. Spallingand axial splitting dominate the failure mode at low confinement while shear failure is the dominantfailure mode at high confinement. Numerical results from simulated compression tests show thatmicrocracking reduces the cohesive component of strength alone, and the frictional strength componentremains unaffected. Results from simulated Brazilian tests show that the tensile strength is influenced bythe presence of microcracks, with a reduction in tensile strength as microcrack intensity increases. Theimportance of microcrack heterogeneity in reproducing a bi-linear or S-shape failure envelope and itseffects on the mechanisms leading to spalling damage near an underground opening are also discussed.

Towards a Micromechanical Understanding of Landslides—Aiming at a Combination of Finite and Discrete Elements with Minimal Number of Degrees of Freedom

In this paper, we propose a combination of discrete elements for the soil and finite elements for the fluid flow field inside the pore space to simulate the triggering of landslides. We give the details for the implementation of third order finite elements (“P2 with bubble”) together with polygonal discrete elements, which allows the formulation with a minimal number of degrees of freedom to save computer time and memory. We verify the implementation with several standard problems from computational fluid dynamics, as well as the decay of a granular step in a fluid as test case for complex flow.

Analysis of regular and chaotic dynamics of the Euler-Bernoulli beams using finite difference and finite element methods

Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.

Numerical simulation of hydraulic fracturing and associated microseismicity using finite-discrete element method

Hydraulic fracturing （HF） technique has been extensively used for the exploitation of unconventional oiland gas reservoirs. HF enhances the connectivity of less permeable oil and gas-bearing rock formationsby fluid injection, which creates an interconnected fracture network and increases the hydrocarbonproduction. Meanwhile, microseismic （MS） monitoring is one of the most effective approaches to evaluatesuch stimulation process. In this paper, the combined finite-discrete element method （FDEM） isadopted to numerically simulate HF and associated MS. Several post-processing tools, includingfrequency-magnitude distribution （b-value）, fractal dimension （D-value）, and seismic events clustering,are utilized to interpret numerical results. A non-parametric clustering algorithm designed specificallyfor FDEM is used to reduce the mesh dependency and extract more realistic seismic information.Simulation results indicated that at the local scale, the HF process tends to propagate following the rockmass discontinuities; while at the reservoir scale, it tends to develop in the direction parallel to themaximum in-situ stress. 2014 Institute of Rock and Soil Mechanics, Chinese Academy of Sciences. Production and hosting byElsevier B.V. All rights reserved.

Asymptotic Behavior of the Finite Difference and the Finite Element Methods for Parabolic Equations

The asymptotic convergence of the solution of the parabolic equation is proved. By the eigenvalues estimation, we obtain that the approximate solutions by the finite difference method and the finite element method are asymptotically convergent. Both methods are considered in continnous time.

DIFFERENCE SCHEME AND NUMERICAL SIMULATION BASED ON MIXED FINITE ELEMENT METHOD FOR NATURAL CONVECTION PROBLEM

The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.

A symplectic finite element method based on Galerkin discretization for solving linear systems

We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.

Fracture Response of Reinforced Concrete Deep Beams Finite Element Investigation of Strength and Beam Size

This article presents a finite element analysis of reinforced concrete deep beams using nonlinear fracture mechanics. The article describes the development of a numerical model that includes several nonlinear processes such as compression and tension softening of concrete, bond slip between concrete and reinforcement, and the yielding of the longitudinal steel reinforcement. The development also incorporates the Delaunay refinement algorithm to create a triangular topology that is then transformed into a quadrilateral mesh by the quad-morphing algorithm. These two techniques allow automatic remeshing using the discrete crack approach. Nonlinear fracture mechanics is incorporated using the fictitious crack model and the principal tensile strength for crack initiation and propagation. The model has been successful in reproducing the load deflections, cracking patterns and size effects observed in experiments of normal and high-strength concrete deep beams with and without stirrup reinforcement.

TIME-DISCRETIZATION PROCEDURE FOR FINITE ELEMENT APPROXIMATION OF COMPRESSIBLE MISCIBLE DISPLACEMENT IN POROUS MEDIA

We’ll consider the model of two-phase compressible miscible displacement in porous media which includes molecular diffusion and dispersion in one dimensional space. Time-discretization procedure is established and analysed. The optimal error estimate in L2 norm is proved by introducing a new interpolation operator R.

Rock Cutting Analysis Employing Finite and Discrete Element Methods

The petroleum industry has shown great interest in the study of drilling optimization on pre-salt formations given the low rates of penetration observed so far. Rate of penetration is the key to economically drill the pre-salt carbonate rock. This work presents the results of numerical modeling through finite element method and discrete element method for single cutter drilling in carbonate samples. The work is relevant to understand the mechanics of drill bit-rock interaction while drilling deep wells and the results were validated with experimental data raised under simulated downhole conditions. The numerical models were carried out under different geometrical configurations, varying the cutter chamfer size and back-rake angles. The forces generated on the cutter are translated into mechanical specific energy as this parameter is often used to measure drilling efficiency. Results indicate that the chamfer size does not change significantly the mechanical specific energy values, characteristic. Results also show there is a significant increase although the cutter aggressiveness is influenced by this geometrical in drilling resistance for larger values of back-rake angle.

Three-dimensional finite element analysls for external midface distraction after different osteotomy in patients with cleft lip and palate

Objective To study three - dimensional finite element analysis for external midface distraction after different osteotomy in patients with cleft lip and palate ( Clp) . Methods Three - dimensional Fem models of Le Fort Ⅰ,Ⅱand Ⅲ,osteotomy in Clp patients were estabolished. External midface distraction were simulated. An anteriorly and inferiorly directed 900 g force was

Finite element response sensitivity analysis of three-dimensional soil-foundation-structure interaction (SFSI) systems

The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method（DDM） has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction（SFSI） models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method（FFD）.The relative importance（RI） of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM（e.g.,accurate,efficient,insensitive to numerical noise）.Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.

PARALLEL ANALYSIS OF COMBINED FINITE/DISCRETE ELEMENT SYSTEMS ON PC CLUSTER

A computational strategy is presented for the nonlinear dynamic analysis of large- scale combined finite/discrete element systems on a PC cluster.In this strategy,a dual-level domain decomposition scheme is adopted to implement the dynamic domain decomposition.The domain decomposition approach perfectly matches the requirement of reducing the memory size per processor of the calculation.To treat the contact between boundary elements in neighbouring subdomains,the elements in a subdomain are classified into internal,interfacial and external elements.In this way,all the contact detect algorithms developed for a sequential computation could be adopted directly in the parallel computation.Numerical examples show that this implementation is suitable for simulating large-scale problems.Two typical numerical examples are given to demonstrate the parallel efficiency and scalability on a PC cluster.

An explicit finite element－finite difference method for analyzing the effect of visco－elastic local topography on the earthquake motion

An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.

Transient analysis of 1D inhomogeneous media by dynamic inhomogeneous finite element method

The dynamic inhomogeneous finite element method is studied for use in the transient analysis of one dimensional inhomogeneous media. The general formula of the inhomogeneous consistent mass matrix is established based on the shape function. In order to research the advantages of this method, it is compared with the general finite element method. A linear bar element is chosen for the discretization tests of material parameters with two fictitious distributions. And, a numerical example is solved to observe the differences in the results between these two methods. Some characteristics of the dynamic inhomogeneous finite element method that demonstrate its advantages are obtained through comparison with the general finite element method. It is found that the method can be used to solve elastic wave motion problems with a large element scale and a large number of iteration steps.