Fracture Analysis of High strength and Ultra high strength Concrete beams by using Finite Element Method

This paper presents the details of nonlinear finite element analysis (FEA) of three point bending specimens made up of high strength concrete (HSC, HSC1) and ultra high strength concrete (UHSC). Brief details about characterization and experimentation of HSC, HSC1 and UHSC have been provided. Cracking strength criterion has been used for simulation of crack propagation by conducting nonlinear FEA. The description about FEA using crack strength criterion has been outlined. Bi-linear tension softening relation has been used for modeling the cohesive stresses ahead of the crack tip. Numerical studies have been carried out on fracture analysis of three point bending specimens. It is observed from the studies that the computed values from FEA are in very good agreement with the corresponding experimental values. The computed values of stress vs crack width will be useful for evaluation of fracture energy, crack tip opening displacement and fracture toughness. Further, these values can also be used for crack growth study, remaining life assessment and residual strength evaluation of concrete structural components.

Operator-splitting finite element algorithms for computations of high-dimensional parabolic problems

An operator-splitting finite element method for solving high-dimensional parabolic equations is presented. The stability and the error estimates are derived for the proposed numerical scheme. Furthermore, two variants of fully-practical operator-splitting finite element algorithms based on the quadrature points and the nodal points, respectively, are presented. Both the quadrature and the nodal point based operator-splitting algorithms are validated using a three-dimensional (3D) test problem. The numerical results obtained with the full 3D computations and the operator-split 2D + 1D computations are found to be in a good agreement with the analytical solution. Further, the optimal order of convergence is obtained in both variants of the operator-splitting algorithms. (C) 2012 Elsevier Inc. All rights reserved.

A Finite Element Method Based Approach for Evaluating the Sensitivity of Faraday-Type Electromagnetic Flow Meter for Liquid Sodium Flow

Faraday-type electromagnetic flow meters are employed for measuring the flow rate of liquid sodium in fast breeder reactors. The calibration of such flow meters, owing to the required elaborative arrangements is rather difficult. On the other hand, theoretical approach requires solution of two coupled electromagnetic partial differential equation with profile of the flow and applied magnetic field as the inputs. This is also quite involved due to the 3D nature of the problem. Alternatively, Galerkin finite element method based numerical solution is suggested in the literature as an attractive option for the required calibration. Based on the same, a computer code in Matlab platform has been developed in this work with both 20 and 27 node brick elements. The boundary conditions are correctly defined and several intermediate validation exercises are carried out. Finally it is shown that the sensitivities predicted by the code for flow meters of four different dimensions agrees well with the results given by analytical expression, thereby providing strong validation. Sensitivity for higher flow rates, for which analytical approach does not exist, is shown to decrease with increase in flow velocity.

An Outward-Wave-Favouring Finite Element-Based Strategy for Exterior Acoustical Problems

This work presents a finite element-based strategy for exterior acoustical problems based on an assumed pressure form that favours outgoing waves. The resulting governing equation, weak formulation, and finite element formulation are developed both for coupled and uncoupled problems. The developed elements are very similar to conventional elements in that they are based on the standard Galerkin variational formulation and use standard Lagrange interpolation functions and standard Gaussian quadrature. In addition and in contrast to wave envelope formulations and their extensions, the developed elements can be used in the immediate vicinity of the radiator/scatterer. The method is similar to the perfectly matched layer (PML) method in the sense that each layer of elements added around the radiator absorbs acoustical waves so that no boundary condition needs to be applied at the outermost boundary where the domain is truncated. By comparing against strategies such as the PML and wave-envelope methods, we show that the relative accuracy, both in the near and far-field results, is considerably higher.

Wave propagation in stiffened structures using spectrally formulated finite element

In this work, the wave propagation analysis of built-up composite structures is performed using frequency domain spectral finite elements, to study the high frequency wave responses. The paper discusses basically two methods for modeling stiffened structures. In the first method, the concept of assembly of 2D spectral plate elements is used to model a built-up structure. In the second approach, spectral finite element method (SFEM) model is developed to model skin-stiffener structures, where the skin is considered as plate element and the stiffener as beam element. The SFEM model developed using the plate-beam coupling approach is then used to model wave propagation in a multiple stiffened structure and also extended to model the stiffened structures with different cross sections such as T-section, I-section and hat section. A number of parametric studies are performed to capture the mode coupling, that is, the flexural-axial coupling present in the wave responses.

Partial delamination modeling in composite beams using a finite element method

A new method of modeling partial delamination in composite beams is proposed and implemented using the finite element method. Homogenized cross-sectional stiffness of the delaminated beam is obtained by the proposed analytical technique, including extension-bending, extension-twist and torsion-bending coupling terms, and hence can be used with an existing finite element method. A two noded C1 type Timoshenko beam element with 4 degrees of freedom per node for dynamic analysis of beams is implemented. The results for different delamination scenarios and beams subjected to different boundary conditions are validated with available experimental results in the literature and/or with the 3D finite element simulation using COMSOL. Results of the first torsional mode frequency for the partially delaminated beam are validated with the COMSOL results. The key point of the proposed model is that partial delamination in beams can be analyzed using a beam model, rather than using 3D or plate models. (c) 2013 Elsevier B.V. All rights reserved.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, first a Fortran code is developed for three dimensional linear elastostatics using constant boundary elements; the code is based on a MATLAB code developed by the author earlier. Next, the code is parallelized using BLACS, MPI, and ScaLAPACK. Later, the parallelized code is used to demonstrate the usefulness of the Boundary Element Method (BEM) as applied to the realtime computational simulation of biological organs, while focusing on the speed and accuracy offered by BEM. A computer cluster is used in this part of the work. The commercial software package ANSYS is used to obtain the `exact' solution against which the solution from BEM is compared; analytical solutions, wherever available, are also used to establish the accuracy of BEM. A pig liver is the biological organ considered. Next, instead of the computer cluster, a Graphics Processing Unit (GPU) is used as the parallel hardware. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature. Also, a serial MATLAB code, and both serial and parallel versions of a Fortran code, which can solve three dimensional (3D) linear elastostatic problems using constant boundary elements, are provided as supplementary files that can be freely downloaded.

A study of the speed and the accuracy of the Boundary Element Method as applied to the computational simulation of biological organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated, with specific reference to the speed and the accuracy offered by BEM. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. A pig liver is the biological organ considered. Results indicate that BEM is an interesting choice for the simulation of biological organs. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

A Study of Speed of the Boundary Element Method as applied to the Realtime Computational Simulation of Biological Organs

In this work, possibility of simulating biological organs in realtime using the Boundary Element Method (BEM) is investigated. Biological organs are assumed to follow linear elastostatic material behavior, and constant boundary element is the element type used. First, a Graphics Processing Unit (GPU) is used to speed up the BEM computations to achieve the realtime performance. Next, instead of the GPU, a computer cluster is used. Results indicate that BEM is fast enough to provide for realtime graphics if biological organs are assumed to follow linear elastostatic material behavior. Although the present work does not conduct any simulation using nonlinear material models, results from using the linear elastostatic material model imply that it would be difficult to obtain realtime performance if highly nonlinear material models that properly characterize biological organs are used. Although the use of BEM for the simulation of biological organs is not new, the results presented in the present study are not found elsewhere in the literature.

Prediction of the stiffness of nanoclay-polypropylene composites using a monte carlo finite element analysis approach

The present work deals with the prediction of stiffness of an Indian nanoclay-reinforced polypropylene composite (that can be termed as a nanocomposite) using a Monte Carlo finite element analysis (FEA) technique. Nanocomposite samples are at first prepared in the laboratory using a torque rheometer for achieving desirable dispersion of nanoclay during master batch preparation followed up with extrusion for the fabrication of tensile test dog-bone specimens. It has been observed through SEM (scanning electron microscopy) images of the prepared nanocomposite containing a given percentage (3–9% by weight) of the considered nanoclay that nanoclay platelets tend to remain in clusters. By ascertaining the average size of these nanoclay clusters from the images mentioned, a planar finite element model is created in which nanoclay groups and polymer matrix are modeled as separate entities assuming a given homogeneous distribution of the nanoclay clusters. Using a Monte Carlo simulation procedure, the distribution of nanoclay is varied randomly in an automated manner in a commercial FEA code, and virtual tensile tests are performed for computing the linear stiffness for each case. Values of computed stiffness modulus of highest frequency for nanocomposites with different nanoclay contents correspond well with the experimentally obtained measures of stiffness establishing the effectiveness of the present approach for further applications.

A finite element method for the Saint-Venant torsion and bending problems for prismatic beams

In this work, we present a finite element formulation for the Saint-Venant torsion and bending problems for prismatic beams. The torsion problem formulation is based on the warping function, and can handle multiply-connected regions (including thin-walled structures), compound and anisotropic bars. Similarly, the bending formulation, which is based on linearized elasticity theory, can handle multiply-connected domains including thin-walled sections. The torsional rigidity and shear centers can be found as special cases of these formulations. Numerical results are presented to show the good coarse-mesh accuracy of both the formulations for both the displacement and stress fields. The stiffness matrices and load vectors (which are similar to those for a variable body force in a conventional structural mechanics problem) in both formulations involve only domain integrals, which makes them simple to implement and computationally efficient. (C) 2014 Elsevier Ltd. All rights reserved.

Spectral finite element based on an efficient layerwise theory for wave propagation analysis of composite and sandwich beams

In this paper, we present a spectral finite element model (SFEM) using an efficient and accurate layerwise (zigzag) theory, which is applicable for wave propagation analysis of highly inhomogeneous laminated composite and sandwich beams. The theory assumes a layerwise linear variation superimposed with a global third-order variation across the thickness for the axial displacement. The conditions of zero transverse shear stress at the top and bottom and its continuity at the layer interfaces are subsequently enforced to make the number of primary unknowns independent of the number of layers, thereby making the theory as efficient as the first-order shear deformation theory (FSDT). The spectral element developed is validated by comparing the present results with those available in the literature. A comparison of the natural frequencies of simply supported composite and sandwich beams obtained by the present spectral element with the exact two-dimensional elasticity and FSDT solutions reveals that the FSDT yields highly inaccurate results for the inhomogeneous sandwich beams and thick composite beams, whereas the present element based on the zigzag theory agrees very well with the exact elasticity solution for both thick and thin, composite and sandwich beams. A significant deviation in the dispersion relations obtained using the accurate zigzag theory and the FSDT is also observed for composite beams at high frequencies. It is shown that the pure shear rotation mode remains always evanescent, contrary to what has been reported earlier. The SFEM is subsequently used to study wavenumber dispersion, free vibration and wave propagation time history in soft-core sandwich beams with composite faces for the first time in the literature. (C) 2014 Elsevier Ltd. All rights reserved.

An extended finite element model coupled with level set method for stable pitting corrosion

Mass balance between metal and electrolytic solution, separated by a moving interface, in stable pit growth results in a set of governing equations which are solved for concentration field and interface position (pit boundary evolution), which requires only three inputs, namely the solid metal concentration, saturation concentration of the dissolved metal ions and diffusion coefficient. A combined eXtended Finite Element Model (XFEM) and level set method is developed in this paper. The extended finite element model handles the jump discontinuity in the metal concentrations at the interface, by using discontinuous-derivative enrichment formulation for concentration discontinuity at the interface. This eliminates the requirement of using front conforming mesh and re-meshing after each time step as in conventional finite element method. A numerical technique known as level set method tracks the position of the moving interface and updates it over time. Numerical analysis for pitting corrosion of stainless steel 304 is presented. The above proposed method is validated by comparing the numerical results with experimental results, exact solutions and some other approximate solutions.

Automatic finite element formulation and assembly of hyperelastic higher order structural models

The formulation of higher order structural models and their discretization using the finite element method is difficult owing to their complexity, especially in the presence of non-linearities. In this work a new algorithm for automating the formulation and assembly of hyperelastic higher-order structural finite elements is developed. A hierarchic series of kinematic models is proposed for modeling structures with special geometries and the algorithm is formulated to automate the study of this class of higher order structural models. The algorithm developed in this work sidesteps the need for an explicit derivation of the governing equations for the individual kinematic modes. Using a novel procedure involving a nodal degree-of-freedom based automatic assembly algorithm, automatic differentiation and higher dimensional quadrature, the relevant finite element matrices are directly computed from the variational statement of elasticity and the higher order kinematic model. Another significant feature of the proposed algorithm is that natural boundary conditions are implicitly handled for arbitrary higher order kinematic models. The validity algorithm is illustrated with examples involving linear elasticity and hyperelasticity. (C) 2013 Elsevier Inc. All rights reserved.

Wave propagation analysis in laminated composite plates with transverse cracks using the wavelet spectral finite element method

This paper presents a newly developed wavelet spectral finite element (WFSE) model to analyze wave propagation in anisotropic composite laminate with a transverse surface crack penetrating part-through the thickness. The WSFE formulation of the composite laminate, which is based on the first-order shear deformation theory, produces accurate and computationally efficient results for high frequency wave motion. Transverse crack is modeled in wavenumber-frequency domain by introducing bending flexibility of the plate along crack edge. Results for tone burst and impulse excitations show excellent agreement with conventional finite element analysis in Abaqus (R). Problems with multiple cracks are modeled by assembling a number of spectral elements with cracks in frequency-wavenumber domain. Results show partial reflection of the excited wave due to crack at time instances consistent with crack locations. (C) 2014 Elsevier B.V. All rights reserved.