Three dimensional boundary element analysis of internal cracks under sliding contact load with a spherical indenter

Using boundary element based three dimensional modelling for linear fracture mechanics, we present an analysis of cracking in a homogeneous medium subject to contact load. The proposed iterative solution procedure allows a simultaneous treatment of a reasonable number of partially closed cracks. It is shown that the most probable direction of propagation of a vertical internal crack is strongly dependent on its size compared to the contact radius and its location with respect to the axis of maximum normal load.     

Bifurcation points and bifurcated branches by an asymptotic numerical method and Padé approximants

The aim of this work is to develop a reliable and fast algorithm to compute bifurcation points and bifurcated branches. It is based upon the asymptotic numerical method (ANM) and Padé approximants. The bifurcation point is detected by analysing the poles of Padé approximants or by evaluating, along the computed solution branch, a bifurcation indicator well adapted to ANM. Several examples are presented to assess the effectiveness of the proposed method, that emanate from buckling problems of thin elastic shells. Especially problems involving large rotations are discussed. Copyright © 2004 John Wiley & Sons, Ltd.     

A p‐version,first order shear deformation,finite element for geometrically non‐linear vibration of curved beams

A p‐version, hierarchical finite element for curved, moderately thick, elastic and isotropic beams is introduced. The convergence properties of the element are analysed and some results are compared with results published elsewhere or calculated using a commercial finite element package. It is verified that, with the proposed element, shear locking does not affect the computation of the natural frequencies and that low dimensional, accurate models are obtainable. Geometrically non‐linear vibrations due to finite deformations, which occur for harmonic excitations with frequencies close to the first three natural frequencies of vibration, are investigated using Newmark's method. The influence of the thickness, longitudinal inertia and curvature radius on the dynamic behaviour of curved beams are studied. Copyright © 2004 John Wiley & Sons, Ltd.     

Second-order accurate constraint formulation for subdivision finite element simulation of thin shells

We present a new method for enforcing boundary conditions within subdivision finite element simulations of thin shells. The proposed framework is demonstrated to be second-order accurate with respect to increasing refinement in the displacement and energy norm for simply supported, clamped, free and symmetric boundary conditions. Second-order accuracy on the boundary is consistent with the accuracy of subdivision-based approaches for the interior of a body. Our proposed framework is applicable to both triangular and quadrilateral refinement schemes, and does not impose any topological requirements upon the underlying subdivision control mesh. Several examples from an obstacle course of benchmark problems are used to demonstrate the convergence of the scheme. Copyright © 2004 John Wiley & Sons, Ltd.     

A direct stiffness analysis of a composite beam with partial interaction

The use of the conventional semi-analytical stiffness method in finite element analysis, in which interpolation polynomials are used to develop the stiffness relationships, leads to problems of curvature locking when beam-type elements are developed for composite members with partial interaction between the materials of which it is comprised. The curvature locking phenomenon that occurs for composite steel–concrete members is quite well reported, and the general approach to minimizing the undesirable ramifications of curvature locking has been to use higher-order polynomials with increasing numbers of internal nodes. This paper presents an alternate formulation based on a direct stiffness approach rather than starting from pre-defined interpolation polynomials, and which does not possess the undesirable locking characteristics. The formulation is based on a more general approach for a bi-material composite flexural member, whose constituent materials are joined by elastic shear connection so as to provide partial interaction. The stiffness relationships are derived, and these are applied to a simply supported and a continuous steel–concrete composite beam to demonstrate the efficacy of the method, and in particular its ability to model accurately both very flexible and very stiff shear connection that causes difficulties when implemented in competitive semi-analytical algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

Application of piece‐wise linear weight functions for 2D 8‐node quadrilateral element in contact problems

The present study is a continuation of our previous work with the aim to reduce problems caused by standard higher order elements in contact problems. The difficulties can be attributed to the inherent property of the Galerkin method which gives uneven distributions of nodal forces resulting in oscillating contact pressures. The proposed remedy is use of piece‐wise linear weight functions. The methods to establish stiffness and/or mass matrix for 8‐node quadrilateral element in 2D are presented, i.e. the condensing and direct procedures. The energy and nodal displacement error norms are also checked to establish the convergence ratio. Interpretation of calculated contact pressures is discussed. Two new 2D 8‐node quadrilateral elements, QUAD8C and QUAD8D, are derived and tested in many examples, which show their good performance in contact problems. Copyright © 2004 John Wiley & Sons, Ltd.     

Constitutive modelling and numerical simulation of multivariant phase transformation in superelastic shape‐memory alloys

This work concerns the micromechanical constitutive modelling, algorithmic implementation and numerical simulation of polycrystalline superelastic alloys under multiaxial loading. The model is formulated in finite deformations and incorporates the effect of texture. The numerical implementation is based on the constrained minimization of the Helmholtz free energy with dissipation. Simulations are conducted for thin tubes of Nitinol under tension–torsion, as well as for a simplified model of a biomedical stent. Copyright © 2004 John Wiley & Sons, Ltd.     

Reduced modified quadratures for quadratic membrane finite elements

Reduced integration is frequently used in evaluating the element stiffness matrix of quadratically interpolated finite elements. Typical examples are the serendipity (Q8) and Lagrangian (Q9) membrane finite elements, for which a reduced 2 × 2 Gauss–Legendre integration rule is frequently used, as opposed to full 3 × 3 Gauss–Legendre integration. This ‘softens’ these element, thereby increasing accuracy, albeit at the introduction of spurious zero energy modes on the element level. This is in general not considered problematic for the ‘hourglass’ mode common to Q8 and Q9 elements, since this spurious mode is non‐communicable. The remaining two zero energy modes occurring in the Q9 element are indeed communicable. However, in topology optimization for instance, conditions may arise where the non‐communicable spurious mode associated with the elements becomes activated. To effectively suppress these modes altogether in elements employing quadratic interpolation fields, two modified quadratures are employed herein. For the Q8 and Q9 membrane elements, the respective rules are a five and an eight point rule. As compared to fully integrated elements, the new rules enhance element accuracy due to the introduction of soft, higher‐order deformation modes. A number of standard test problems reveal that element accuracy remains comparable to that of the under‐integrated counterparts. Copyright © 2004 John Wiley & Sons, Ltd.     

Crushing behaviour of the energy absorbing 'tensor skin' panels

ABSTRACT Research results concerning the simulation of the crushing behaviour of composite systems with energy absorption characteristics are presented in the present work. The study is focused on the ‘tensor skin’ concept, an energy absorbing composite system that was originally developed to improve the crashworthiness of helicopters under water impact and which is promising for utilization in the construction of the lower part of composite fuselage aircraft. The ‘tensor skin’ concept comprises a folded or corrugated composite construction, which upon loading unfolds by forming ‘plastic hinges’, leading to an increase in the load bearing capability of the structure. The numerical modelling issues and the critical aspects of the simulation are discussed. Verification of the numerical simulation procedure is performed by experimental work. The experimental results utilized to assess and validate the numerical procedure were derived within the European Research Project ‘Design for Crash Survivability – CRASURV’ (BRITE – Aeronautics Area). The results of the simulations are generally in good agreement with experimental data.