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.     

Wettability of Lithium by Liquid Sodium, Potassium and Rubidium

The wettability of the lithium surface by liquid alkali metals Na, K and Rb is investigated for the first time by the sessile drop method in an all-soldered instrument under conditions of a high vacuum. The presence of the temperature threshold of wetting is found in Li–Na and Li–K systems at temperatures of 325°C and 160°C, respectively. A conclusion is drawn that an abrupt decrease in wetting angles in the investigated systems is associated with a marked decrease in the interfacial tension at lithium-lithium, lithium-potassium, and lithium--rubidium interfaces owing to the beginning of a noticeable mutual solubility of the components at relatively high temperatures.     

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.     

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.     

A direct numerical simulation of transition phenomena in a mixed convection channel flow

This study intends to provide an increased understanding of the laminar-turbulent transition phenomena for the buoyancy-assisted heated vertical channel flow during the early transient stage. The spectral method with weak formulation is applied in the direct numerical simulation. Initial disturbances consist of the finite-amplitude two-dimensional TS wave and a pair of three-dimensional oblique waves for the K-type disturbances. The results from the harmonic energy competitions of different wave modes show that for the buoyancy-assisted heated flow, the (k_x=1, k_z=1) or (1,1) and (1,0) modes would gain energy immediately and start to rise at almost the same rate. This phenomenon is different from that of the buoyancy-opposed flow, where the (1,1) mode decays slowly in the beginning until other modes gain enough energy and then it begins to grow quickly and overtakes the (1,0) mode after a short time period. These different transition patterns match with the experimental results that the flow transition is supercritical and subcritical for the buoyancy-assisted and -opposed flows, respectively. Buoyancy-assisted heated flow transition follows the general trend of an isothermal flow in the beginning, but the thermal-buoyant force is crucial in accelerating the instability and also causing notable differences during the subsequent transition process. All of the results for the vortex structures development, kinetic energy budget of the disturbances, flow visualization by tagged fluid particles, and the local temperature fluctuations are consistent in pointing to a clear pattern for the buoyancy-assisted heated flow transition.     

Connecting the second exhaust-heat boiler to the operating first one under the conditions of flow circuits of combined-cycle plants with two gas-turbine units and one steam turbine

Problems arising with connecting the second exhaust-heat boiler to the first exhaust-heat boiler under load in the case of flow circuits of combined-cycle plants of type PGU-450 are considered. Similar problems exist for two-boiler single-turbine units of type PGU-325 and for two-boiler single-turbine units of traditional steampower plants with an output of 200–300 MW upon connection of the second furnace of a two-furnace boiler to the operating first furnace. Analysis of the temperature-induced stressed state and of the cyclic strength of thick-walled elements of the steam-water conduit of PGU-450 and of two-furnace boilers of two-boiler single-turbine units with an output of 200–300 MW (the check valve, the inlet header of the water economizer) is used for determining the steam temperature differences in the high-pressure conduits to the steam turbine and the order of connection of the second exhaust-heat boiler to the operating fist one. Modes of enabling a two-boiler single-turbine unit of a combined-cycle plant, which are most “rigid” from the standpoint of thermal strength, are determined. __________ Translated from élektricheskie Stantsii, No. 2, February 2006, pp. 11–17.