Finite element simulation of wedge indentation

A finite element simulation of frictionless wedge indentation of a copper strip has been carried out under plane strain conditions. The problem was first modelled using an one-pass contact algorithm. The difficulties associated with using this method to model wedge indentation problems are explained. An alternative procedure which alleviates some of the problems associated with the one-pass contact algorithm is proposed for modelling frictionless wedge indentation. Also, a re-meshing procedure which has to be carried out when the distortion of the elements around the indenter becomes significant, is discussed. A sample problem involving indentation of a 4 mm copper strip by a rigid wedge indenter has been modelled and the results are compared with experimental and theoretical results.

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.

Tunable Brownian vortex at the interface

A general kind of Brownian vortices is demonstrated by applying an external nonconservative force field to a colloidal particle bound by a conservative optical trapping force at a liquid-air interface. As the liquid medium is translated at a constant velocity with the bead trapped at the interface, the drag force near the surface provides enough rotational component to bias the particle's thermal fluctuations in a circulatory motion. The interplay between the thermal fluctuations and the advection of the bead in constituting the vortex motions is studied, and we infer that the angular velocity of the circulatory motion offers a comparative measure of the interface fluctuations.

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.

The role of the interface in foam-bearing glass-epoxy composites subjected to impact loading

foam, either stacked together as three layers (MC) or inserted at three different positions (3L) while arranging the stacking sequence during the fabrication of glass fiber-epoxy composites, form the subject of investigation. This stacking variation resulted in a different interfacial area between these foam materials and the glass-epoxy regions in the laminates. This area in designed to be maximum for the 3L variety. The energy of impact being high enough to cause development of the crack in the samples, how the change in interfacial area affects the traverse of the crack front and the failure feature of the laminated composite are reported in the form of photomacrographs in this work. The results point to significant changes for the impact data, like for instance the peak load attained by the different samples, through thickness crack propagation and tensile fracture features on the non-impacted end for the plain variety, separation about the mid-zone for the MC laminates and two or more layer separations for the 3L variety. The separation for the foam-bearing systems occur invariably at the interface and here again one of the (two identical) interfaces only is chosen for the separation.

Influence of foam/glass-epoxy interface on the macroscopic compressive failure features

The effect on the macroscopic compressive failure features of introduction of two flexible foam layers, either together at mid-region or separately at two locations that are away from the midregion, into a glass-epoxy (G-E) system is studied in this work. In this experimental approach an attempt to look at the possible influence the foam/G-E interface region has on the way the materials respond to compressive loading is made by involving an analyses of macrofractographic features. While foam-free samples fail by extensive ear formation and separation nearer to the mid-region, the foam bearing ones display pronounced interface separation. The positioning of the foam sheet(s) has a bearing on the failure features.

Nonlinear analysis of adhesively bonded lap joints considering viscoplasticity in adhesives

This paper presents nonlinear finite element analysis of adhesively bonded joints considering the elastoviscoplastic constitutive model of the adhesive material and the finite rotation of the joint. Though the adherends have been assumed to be linearly elastic, the yielding of the adhesive is represented by a pressure sensitive modified von Mises yield function. The stress-strain relation of the adhesive is represented by the Ramberg-Osgood relation. Geometric nonlinearity due to finite rotation in the joint is accounted for using the Green-Lagrange strain tensor and the second Piola-Kirchhoff stress tensor in a total Lagrangian formulation. Critical time steps have been calculated based on the eigenvalues of the transition matrices of the viscoplastic model of the adhesive. Stability of the viscoplastic solution and time dependent behaviour of the joints are examined. A parametric study has been carried out with particular reference to peel and shear stress along the interface. Critical zones for failure of joints have been identified. The study is of significance in the design of lap joints as well as on the characterization of adhesive strength. (C) 1999 Elsevier Science Ltd. All rights reserved.

Re-examination of crack opening model of interface fracture

The complex singularity associated with a crack at the interface between two dissimilar, isotropic and homogeneous materials leads to mathematical artefacts, such as stress oscillations and crack face interpenetrations in the vicinity of the crack tip. To avoid these unrealistic features, Sinclair (Sinclair GB. On the stress singularity at an interface crack. International Journal of Fracture 1980;16(2):111-9) assumed a finite crack opening angle (COA) such that the singularity lambda became real equal to 1/2. This paper extends the COA model by considering real singularities not necessarily equal to 1/2. When COA is 0 degrees: the interface crack singularity is complex with a real part equal to 1/2. On increasing COA, the imaginary part of the singularity decreases and becomes zero at a threshold value of COA; at this point, the singularity is a real, repeated value. A further increase in COA results in a pair of real singularities. Different crack opening configurations and material combinations are studied, and results presented for threshold COAs and associated values of singularity. Stress analyses for these three regimes: (a) complex, (b) real pair and (c) real repeated singularities, are reported. It is seen that additional complexities are present in the last case. Typical results for stress fields are also included for comparing with standard fields. (C) 1999 Elsevier Science Ltd. All rights reserved.

The influence of temperature gradient and lowering speed on the melt-solid interface shape of GaxIn1-xSb alloy crystals grown by vertical Bridgman technique

Radially homogeneous bulk alloys of GaxIn1-xSb in the range 0.7 < x < 0.8, have been grown by vertical Bridgman technique. The factors affecting the interface shape during the growth were optimised to achieve zero convexity. From a series of experiments, a critical ratio of the temperature gradient (G) of the furnace at the melting point of the melt composition to the ampoule lowering speed (v) was deduced for attaining the planarity of the melt-solid interface. The studies carried out on directional solidification of Ga0.77In0.23Sb mixed crystals employing planar melt-solid interface exhibited superior quality than those with nonplanar interfaces. The solutions to certain problems encountered during the synthesis and growth of the compound were discussed. (C) 1999 Elsevier Science B.V. All rights reserved.

Influence of interface properties on fracture behaviour of concrete

Hardened concrete is a three-phase composite consisting of cement paste, aggregate and interface between cement paste and aggregate. The interface in concrete plays a key role on the overall performance of concrete. The interface properties such as deformation, strength, fracture energy, stress intensity and its influence on stiffness and ductility of concrete have been investigated. The effect of composition of cement, surface characteristics of aggregate and type of loading have been studied. The load-deflection response is linear showing that the linear elastic fracture mechanics (LEFM) is applicable to characterize interface. The crack deformation increases with large rough aggregate surfaces. The strength of interface increases with the richness of concrete mix. The interface fracture energy increases as the roughness of the aggregate surface increases. The interface energy under mode II loading increases with the orientation of aggregate surface with the direction of loading. The chemical reaction between smooth aggregate surface and the cement paste seems to improve the interface energy. The ductility of concrete decreases as the surface area of the strong interface increases. The fracture toughness (stress intensity factor) of the interface seems to be very low, compared with hardened cement paste, mortar and concrete.

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.