Finite element simulations of ductile rupture in a constrained metal foil

A numerical study of the ductile rupture in a metal foil constrained between two stiff ceramic blocks is performed. The finite element analysis is carried out under the conditions of mode I, plane strain, small-scale yielding. The rate-independent version of the Gurson model that accounts for the ductile failure mechanisms of microvoid nucleation, growth and coalescence is employed to represent the behavior of the metal foil. Different distributions of void nucleating sites in the metal foil are considered for triggering the initiation of discrete voids. The results clearly show that far-field triaxiality-induced cavitation is the dominant failure mode when the spacing of the void nucleating sites is large. On the contrary, void coalescence near the notch tip is found to be the operative failure mechanism when closely spaced void nucleating sites are considered.

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.

Generation of an eight node brick element using Papcovitch-Neuber functions

Some conventional finite elements suffer from drawbacks, such as shear locking, membrane locking, etc. To overcome them researchers have developed various techniques, termed as tricks by some and variational crimes by others. Many attempts have been made, but satisfactory explanations for why some of these techniques work have not been obtained, especially in the case of solid elements. This paper attempts a simple non-conforming solid element using assumed displacement fields which satisfy the Navier equation exactly. Its behaviour under simple loadings like bending, torsion and tension is examined and comparisons are made with existing elements.

Large deformation finite element analysis of laminated circular composite plates

A 48 d.o.f., four-noded quadrilateral laminated composite shell finite element is particularised to a sector finite element and is used for the large deformation analysis of circular composite laminated plates. The strain-displacement relationships for the sector element are obtained by reducing those of the quadrilateral shell finite element by substituting proper values for the geometric parameters. Subsequently, the linear and tangent stiffness matrices are formulated using conventional methods. The Newton-Raphson method is employed as the nonlinear solution technique. The computer code developed is validated by solving an isotropic case for which results are available in the literature. The method is then applied to solve problems of cylindrically orthotropic circular plates. Some of the results of cylindrically orthotropic case are compared with those available in the literature. Subsequently, application is made to the case of laminated composite circular plates having different lay-up schemes. The computer code can handle symmetric/unsymmetric lay-up schemes. The large displacement analysis is useful in estimating the damage in composite plates caused by low-velocity impact.

Violin string shape functions for finite element analysis of rotating Timoshenko beams

Violin strings are relatively short and stiff and are well modeled by Timoshenko beam theory. We use the static part of the homogeneous differential equation of violin strings to obtain new shape functions for the finite element analysis of rotating Timoshenko beams. For deriving the shape functions, the rotating beam is considered as a sequence of violin strings. The violin string shape functions depend on rotation speed and element position along the beam length and account for centrifugal stiffening effects as well as rotary inertia and shear deformation on dynamic characteristics of rotating Timoshenko beams. Numerical results show that the violin string basis functions perform much better than the conventional polynomials at high rotation speeds and are thus useful for turbo machine applications. (C) 2011 Elsevier B.V. All rights reserved.

A finite element method for compressible viscous fluid flows

This work presents a mixed three-dimensional finite element formulation for analyzing compressible viscous flows. The formulation is based on the primitive variables velocity, density, temperature and pressure. The goal of this work is to present a `stable' numerical formulation, and, thus, the interpolation functions for the field variables are chosen so as to satisfy the inf-sup conditions. An exact tangent stiffness matrix is derived for the formulation, which ensures a quadratic rate of convergence. The good performance of the proposed strategy is shown in a number of steady-state and transient problems where compressibility effects are important such as high Mach number flows, natural convection, Riemann problems, etc., and also on problems where the fluid can be treated as almost incompressible. Copyright (C) 2010 John Wiley & Sons, Ltd.

Stabilization of thermal degradation of poly(methyl methacrylate) by polysulfide polymers

The thermal degradation of poly(methyl methacrylate) (PMMA) in the presence of polysulfide polymers, namely, poly( styrene disulfide) (PSD) and poly(styrene tetrasulfide) (PST) was studied using thermogravimetry (TG) and direct pyrolysis-mass spectrometric (DP-MS) analysis. Both PSD and PST were found to stabilizethe PMMA degradation, which was explained by both radical recombination and a chain-transfer mechanism. (C) 1997 John Wiley & Sons, Inc.

A finite element assessment of flexural strength of prestressed concrete beams with fiber reinforcement

This paper presents an assessment of the flexural behavior of 15 fully/partially prestressed high strength concrete beams containing steel fibers investigated using three-dimensional nonlinear finite elemental analysis. The experimental results consisted of eight fully and seven partially prestressed beams, which were designed to be flexure dominant in the absence of fibers. The main parameters varied in the tests were: the levels of prestressing force (i.e, in partially prestressed beams 50% of the prestress was reduced with the introduction of two high strength deformed bars instead), fiber volume fractions (0%, 0.5%, 1.0% and 1.5%), fiber location (full depth and partial depth over full length and half the depth over the shear span only). A three-dimensional nonlinear finite element analysis was conducted using ANSYS 5.5 [Theory Reference Manual. In: Kohnke P, editor. Elements Reference Manual. 8th ed. September 1998] general purpose finite element software to study the flexural behavior of both fully and partially prestressed fiber reinforced concrete beams. Influence of fibers on the concrete failure surface and stress-strain response of high strength concrete and the nonlinear stress-strain curves of prestressing wire and deformed bar were considered in the present analysis. In the finite element model. tension stiffening and bond slip between concrete and reinforcement (fibers., prestressing wire, and conventional reinforcing steel bar) have also been considered explicitly. The fraction of the entire volume of the fiber present along the longitudinal axis of the prestressed beams alone has been modeled explicitly as it is expected that these fibers would contribute to the mobilization of forces required to sustain the applied loads across the crack interfaces through their bridging action. A comparison of results from both tests and analysis on all 15 specimens confirm that, inclusion of fibers over a partial depth in the tensile side of the prestressed flexural structural members was economical and led to considerable cost saving without sacrificing on the desired performance. However. beams having fibers over half the depth in only the shear span, did not show any increase in the ultimate load or deformational characteristics when compared to plain concrete beams. (C) 2002 Published by Elsevier Science Ltd.

A fixed-grid finite element based enthalpy formulation for generalized phase change problems: role of superficial mushy region

The enthalpy method is primarily developed for studying phase change in a multicomponent material, characterized by a continuous liquid volume fraction (phi(1)) vs temperature (T) relationship. Using the Galerkin finite element method we obtain solutions to the enthalpy formulation for phase change in 1D slabs of pure material, by assuming a superficial phase change region (linear (phi(1) vs T) around the discontinuity at the melting point. Errors between the computed and analytical solutions are evaluated for the fluxes at, and positions of, the freezing front, for different widths of the superficial phase change region and spatial discretizations with linear and quadratic basis functions. For Stefan number (St) varying between 0.1 and 10 the method is relatively insensitive to spatial discretization and widths of the superficial phase change region. Greater sensitivity is observed at St = 0.01, where the variation in the enthalpy is large. In general the width of the superficial phase change region should span at least 2-3 Gauss quadrature points for the enthalpy to be computed accurately. The method is applied to study conventional melting of slabs of frozen brine and ice. Regardless of the forms for the phi(1) vs T relationships, the thawing times were found to scale as the square of the slab thickness. The ability of the method to efficiently capture multiple thawing fronts which may originate at any spatial location within the sample, is illustrated with the microwave thawing of slabs and 2D cylinders. (C) 2002 Elsevier Science Ltd. All rights reserved.

Modeling of concrete cracking-A hybrid technique of using displacement discontinuity element method and direct boundary element method

This paper presents a new approach by making use of a hybrid method of using the displacement discontinuity element method and direct boundary element method to model concrete cracking by incorporating fictitious crack model. Fracture mechanics approach is followed using the Hillerborg's fictitious crack model. A boundary element based substructure method and a hybrid technique of using displacement discontinuity element method and direct boundary element method are compared in this paper. In order to represent the process zone ahead of the crack, closing forces are assumed to act in such a way that they obey a linear normal stress-crack opening displacement law. Plain concrete beams with and without initial crack under three-point loading were analyzed by both the methods. The numerical results obtained were shown to agree well with the results from existing finite element method. The model is capable of reproducing the whole range of load-deflection response including strain-softening and snap-back behavior as illustrated in the numerical examples. (C) 2011 Elsevier Ltd. All rights reserved.

On using peak amplitude and rise time for AE source characterization

Acoustic Emission (AE) signals, which are electrical version of acoustic emissions, are usually analysed using a set of signal parameters. The major objective of signal analysis is to study the characteristics of the sources of emissions. Peak amplitude (P-a) and rise time (R-t) are two such parameters used for source characterization. In this paper, we theoretically investigate the efficiency of P-a and R-t to classify and characterize AE sources by modelling the input stress pulse and transducer. Analytical expressions obtained for P-a and R-t clearly indicate their use and efficiency for source characterization. It is believed that these results may be of use to investigators in areas like control systems and signal processing also.

Systematic trends in structural and thermodynamic properties of ternary oxides in the systems Ln-Pd-O (Ln = lanthanide element)

Isothermal sections of the phase diagrams for the systems Ln-Pd-O (Ln = lanthanide element) at 1223 K indicate the presence of two inter-oxide compounds Ln(4)PdO(7) and Ln(2)Pd(2)O(5) for Ln = La, Pr, Nd, Sm, three compounds Ln(4)PdO(7), Ln(2)PdO(4) and Ln(2)Pd(2)O(5) for Ln = Eu, Gd and only one compound of Ln(2)Pd(2)O(5) for Ln = Tb to Ho. The lattice parameters of the compounds Ln(4)PdO(7), Ln(2)PdO(4) and Ln(2)Pd(2)O(5) show systematic nonlinear variation with atomic number. The unit cell volumes decrease with increasing atomic number. The standard Gibbs energies, enthalpies and entropies of formation of the ternary oxides from their component binary oxides (Ln(2)O(3) and PdO) have been measured recently using an advanced version of the solid-state electrochemical cell. The Gibbs energies and enthalpies of formation become less negative with increasing atomic number of Ln. For all the three compounds, the variation in Gibbs energy and enthalpy of formation with atomic number is markedly non-linear. The decrease in stability with atomic number is most pronounced for Ln(2)Pd(2)O(5), followed by Ln(4)PdO(7) and Ln(2)PdO(4). This is probably related to the repulsion between Pd2+ ions on the opposite phases Of O-8 cubes in Ln(2)Pd(2)O(5), and the presence of Ln-filled O-8 cubes that share three faces with each other in Ln4PdO7. The values for entropy of formation of all the ternary oxides from their component binary oxides are relatively small. Although the entropies of formation show some scatter, the average value for Ln = La, Pr, Nd is more negative than the average value for the other lanthanide elements. From this difference, an average value for the structure transformation entropy of Ln(2)O(3) from C-type to A-type is estimated as 0.87 J.mol(-1).K-1. The standard Gibbs energies of formation of these ternary oxides from elements at 1223 K are presented as a function of lanthanide atomic number. By invoking the Neumann-Kopp rule for heat capacity, thermodynamic properties of the inter-oxide compounds at 298.15 K are estimated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Nonlinear Finite-Element Modeling of Batter Piles under Lateral Load

In this paper, a finite-element model is developed in which the nonlinear soil behavior is represented by a hyperbolic relation for static load condition and modified hyperbolic relation, which includes both degradation and gap for a cyclic load condition. Although batter piles are subjected to lateral load, the soil resistance is also governed by axial load, which is incorporated by considering the P-Δ moment and geometric stiffness matrix. By adopting the developed numerical model, static and cyclic load analyses are performed adopting an incremental-iterative procedure where the pile is idealized as beam elements and the soil as elastoplastic spring elements. The proposed numerical model is validated with published laboratory and field pile test results under both static and cyclic load conditions. This paper highlights the importance of the degradation factor and its influence on the soil resistance-displacement (p-y) curve, number of cycles of loading, and cyclic load response.

Localisation of discrete change in a transformer winding: a network-function-loci approach

An entirely different approach for localisation of winding deformation based on terminal measurements is presented. Within the context of this study, winding deformation means, a discrete and specific change externally imposed at a particular position on the winding. The proposed method is based on pre-computing and plotting the complex network-function loci e.g. driving-point impedance (DPI)] at a selected frequency, for a meaningful range of values for each element (increasing and decreasing) of the ladder network which represents the winding. This loci diagram is called the nomogram. After introducing a discrete change, amplitude and phase of DPI are measured. By plotting this single measurement on the nomogram, it is possible to estimate the location and identify the extent of change. In contrast to the existing approach, the proposed method is fast, non-iterative and yields reasonably good localisation. Experimental results for actual transformer windings (interleaved and continuous disc) are presented.