Study of Layering Procedures in Finite-Element Analysis of RC Flexural and Torsional Elements

A new postcracking formulation for concrete, along with both implicit and explicit layering procedures, is used in the analysis of reinforced-concrete (RC) flexural and torsional elements. The postcracking formulation accounts for tension stiffening in concrete along the rebar directions, compression softening in cracked concrete based on either stresses or strains, and aggregate interlock based on crack-confining normal stresses. Transverse shear stresses computed using the layering procedures are included in material model considerations that permit the development of inclined cracks through the RC cross section. Examples of a beam analyzed by both the layering techniques, a torsional element, and a column-slab connection region analyzed by the implicit layering procedure are presented here. The study highlights the primary advantages and disadvantages of each layering approach, identifying the class of problems where the application of either procedure is more suitable.

An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

Lateral response of monopiles using centrifuge testing and finite element analysis

Offshore wind capacity is expected to grow exponentially over the next decade resulting in the production of a considerable amount of renewable energy. Monopiles are currently the most popular type of foundation for supporting offshore wind turbines in shallow to medium depth waters. In this paper, the load-deformation response of a 3.8 m diameter monopile installed in soft clays when subjected to axial and lateral loading is investigated using centrifuge testing and soil pore-fluid coupled three-dimensional finite element analysis. Monopile deformation is principally assessed in terms of its lateral displacements and bending moments. Its behaviour as a short rigid pile is discussed using concepts such as its rotation at mudline and the pile depth at which pivoting occurs. © 2014 Taylor & Francis Group.

A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Comparison of finite element and finite volume methods application in geometrically nonlinear stress analysis

A novel three-dimensional finite volume (FV) procedure is described in detail for the analysis of geometrically nonlinear problems. The FV procedure is compared with the conventional finite element (FE) Galerkin approach. FV can be considered to be a particular case of the weighted residual method with a unit weighting function, where in the FE Galerkin method we use the shape function as weighting function. A Fortran code has been developed based on the finite volume cell vertex formulation. The formulation is tested on a number of geometrically nonlinear problems. In comparison with FE, the results reveal that FV can reach the FE results in a higher mesh density.

Dynamic load-balancing of finite element applications with the DRAMA library

The DRAMA library, developed within the European Commission funded (ESPRIT) project DRAMA, supports dynamic load-balancing for parallel (message-passing) mesh-based applications. The target applications are those with dynamic and solution-adaptive features. The focus within the DRAMA project was on finite element simulation codes for structural mechanics. An introduction to the DRAMA library will illustrate that the very general cost model and the interface designed specifically for application requirements provide simplified and effective access to a range of parallel partitioners. The main body of the paper will demonstrate the ability to provide dynamic load-balancing for parallel FEM problems that include: adaptive meshing, re-meshing, the need for multi-phase partitioning.

Dynamic multi-partitioning for parallel finite element applications

The central product of the DRAMA (Dynamic Re-Allocation of Meshes for parallel Finite Element Applications) project is a library comprising a variety of tools for dynamic re-partitioning of unstructured Finite Element (FE) applications. The input to the DRAMA library is the computational mesh, and corresponding costs, partitioned into sub-domains. The core library functions then perform a parallel computation of a mesh re-allocation that will re-balance the costs based on the DRAMA cost model. We discuss the basic features of this cost model, which allows a general approach to load identification, modelling and imbalance minimisation. Results from crash simulations are presented which show the necessity for multi-phase/multi-constraint partitioning components

Finite element analysis of an embankment on a soft estuarine deposit using an elastic–viscoplastic soil model

A new elastic–viscoplastic (EVP) soil model has been used to simulate the measured deformation response of a soft estuarine soil loaded by a stage-constructed embankment. The simulation incorporates prefabricated vertical drains installed in the foundation soils and reinforcement installed at the base of the embankment. The numerical simulations closely matched the temporal changes in surface settlement beneath the centerline and shoulder of the embankment. More importantly, the elastic–viscoplastic model simulated the pattern and magnitudes of the lateral deformations beneath the toe of the embankment — a notoriously difficult aspect of modelling the deformation response of soft soils. Simulation of the excess pore-water pressure proved more difficult because of the heterogeneous nature of the estuarine deposit. Excess pore-water pressures were, however, mapped reasonably well at three of the six monitoring locations. The simulations were achieved using a small set of material constants that can easily be obtained from standard laboratory tests. This study validates the use of the EVP model for problems involving soft soil deposits beneath loading from a geotechnical structure.

Numerical simulation of a patent technology for sealing of deep-sea oil wells using nonlinear finite element method

For over 50 years bridge plugs and cement have been used for well abandonment and work over and are still the material of choice. However the failures of cement abandonments using bridge plugs has been reported on many occasions, some of which have resulted in fatal consequences. A new patented product is designed to address the shortcomings associated with using bridge plugs and cement. The new developed tools use an alloy based on bismuth that is melted in situ using Thermite reaction. The tool uses the expansion properties of bismuth to seal the well. Testing the new technology in real field under more than 2 km deep sea water can be expensive. Virtual simulation of the new device under simulated thermal and mechanical environment can be achieved using nonlinear finite element method to validate the product and reduce cost. Experimental testing in the lab is performed to measure heat generated due to thermite reaction. Then, a sequential thermal mechanical explicit/implicit finite element solver is used to simulate the device under both testing lab and deep water conditions.

The Behavior of compressively loaded stiffener runout specimens - Part II: Finite element analysis

The termination of stiffeners in composite aircraft structures give rise to regions of high interlaminar shear and peel stresses as the load in the stiffener is diffused into the skin. This is of particular concern in co-cured composite stiffened structures where there is a relatively low resistance to through-thickness stress components at the skin-stiffener interface. In Part I, experimental results of tested specimens highlighted the influence of local design parameters on their structural response. Indeed some of the observed behavior was unexpected. There is a need to be able to analyse a range of changes in geometry rapidly to allow the analysis to form an integral part of the structural design process.

This work presents the development of a finite element methodology for modelling the failure process of these critical regions. An efficient thick shell element formulation is presented and this element is used in conjuction with the Virtual Crack Closure Technique (VCCT) to predict the crack growth characteristics of the modelled specimens. Three specimens were modelled and the qualitative aspects of crack growth were captured successfully. The shortcomings in the quantitative correlation between the predicted and observed failure loads are discussed. There was evidence to suggest that high through-thickness compressive stresses enhanced the fracture toughness in these critical regions.

Simulation of three-dimension cutting force and tool deflection in the end milling operation based on finite element method

This paper presents a 3D simulation system which is employed in order to predict cutting forces and tool deflection during end-milling operation. In order to verify the accuracy of 3D simulation, results (cutting forces and tool deflection) were compared with those based on the theoretical relationships, in terms of agreement with experiments. The results obtained indicate that the simulation is capable of predicting the cutting forces and tool deflection.