Geometry of “developable cones”

When a thin disc is supported on the rim of a bowl, and its centre is pushed down by a finger, it adopts a characteristic conformation, known as a “developable cone”, and sketched in Fig. 1(a): the main, broadly conical, shape can only form if about one-quarter of the disc buckles upwards. There is a curved intersection between the two parts, which takes the form of a crescent-shaped “crease” near its apex, but with the flanking regions less tightly deformed. The “developable cone” is a recurring motif in a wide range of physical situations—crumpling, buckling, draping—and its mechanics provides a key to understand the phenomena, whether the disc deforms in the elastic or the plastic range. The task of this paper is to study only geometrical features of the “developable cone”. The first step is to replace the actual crease (Fig. 1(a)) by an idealised “sharp” crease (Fig. 1(b)). The second step is to study the apparently “large-rotation” problem of kinematics by means of an adaptation of the classical “yield-line” pattern of folding, but with a crucial added constraint that springs from Gauss's analysis of inextensional deformation. We illustrate the method via a graded sequence of examples, and we close with a discussion.     

Some observations on the CLNA model in fretting fatigue

Using the Atzori–Lazzarin criterion, the author has recently proposed a unified model for Fretting Fatigue denominated Crack-Like Notch Analogue—CLNA model, considering only two possible behaviours: either “crack-like” or “large blunt notch”. In a general FF condition, the former condition is treated with a single contact problem corresponding to the MIT Crack Analogue (CA) improved in some details also by the author. The latter, with a simple peak stress condition, i.e. a simple Notch Analogue model, simply stating that below the fatigue limit, infinite life is predicted for any size of contact. In the typical condition of constant normal load and in phase oscillating tangential and bulk loads, both limiting conditions are immediately written, and the CLNA model permits to collapse the effect of the contact loads on a single closed form equation (differently from many other models which do not permit this flexibility). For not too large contact areas (“crack-like” contact) no dependence at all on geometry is predicted, but only on 3 load factors (bulk stress, tangential load and average pressure) and size of the contact. Only in the “large blunt notch” region occurring typically only at very large sizes of contact does size-effect disappear, but the dependence on all other factors including geometry remains. The model compares favourably with some experimental results in the literature. In this paper, some aspects of the CLNA model are further elucidated.     

Closed-form solution for wedge cutting force through thin metal sheets

A simple kinematic model is developed which describes the main features of the process of the cutting of a plate by a rigid wedge. It is assumed in this model that the plate material curls up into two inclined cylinders as the wedge advances into the plate. This results in membrane stretching up to fracture of the material near the wedge tip, while the “flaps” in the wake of the cut undergo cylindrical bending. Self-consistent, single-term formulas for the indentation force and the energy absorption are arrived at by relating the “far-field” and “near-tip” deformation events through a single geometric parameter, the instantaneous rolling radius. Further analysis of this solution reveals a weak dependence on the wedge angle and a strong dependence on friction coefficient. The final equation for the approximate cutting force over a range of wedge semiangles 10° ≤ θ ≤ 30° and friction coefficients 0.1 ≤ μ ≤ 0.4 is: F = 3.28σ₀(δ_t)^0.2l^0.4t^1.6μ^0.4, which is identical in form and characteristics to the empirical results recently reported by Lu and Calladine [Int. J. Mech. Sci.32, 295–313 (1990)].This analysis is believed to resolve a controversy recently developed in the literature over the interpretation of plate cutting experiments.     

Assessments of the kinetic and dynamic loads sustained by stationary tools during high rate plastic forming

A systematic method for evaluating the kinetic and dynamic loads sustained by stationary tools (as opposed to moving tools for which methods already exist) during high rate plastic forming is examined and exemplified by examples. It is essentially based on the momentum theorem for continua for incompressible flow, utilizing kinematically admissible velocity fields. In steady state forming processes (such as rolling, wire drawing, etc.), the difference between the active load (imposed or calculated a priori) and the reactive load, is formulated rigorously, whereas for non-steady processes (forging, impact extrusion, etc.) the formulation gives merely an approximation to the dynamic effects on the tools. The resulting velocity-dependent reactions on the tools are given in terms of two nondimensional numbers, namely, the “kinetic head” (u₀²/σ₀) (called the Euler Number) and the “dynamic head” (ú₀L/σ₀), which includes the machine speed (u₀), machine acceleration ( ), material density , yield strength ₀ and a characteristic dimension of the product, L. The same two non-dimensional heads emerged previously from energy-balance consideration in Ref. [1], while approximating dynamic loads on moving tools, hence a consistency is demonstrated. These heads are unavoidably multiplied by geometrical functions, which typify the specific process under consideration and may amplify (or diminish) the intensity of the dynamic effects. The present work is focussed on quantifying, by the above method, the inherent difference between the reactive load sustained by the non-moving tool (say, a die) and the acting load carried by the moving tool (piston, ram, etc.) In particular cases of very slow processes, these loads are equal by static equilibrium. In some practical processes (like rolling) their difference appears to be relatively small, whereas in others (like impact extrusion) it appears extremely large.     

Tribological behaviour of duplex treated Ti–6Al–4V: combining nitrogen PSII with a DLC coating

The chemical structure and tribological behaviour of Ti–6Al–4V plasma source ion implanted with nitrogen then DLC-coated in an acetylene plus hydrogen-glow discharge (bias voltage −10 to −30 kV) were investigated. The as-modified samples have a TiN/H:DLC multilayer architecture (coating resistivity 1.6×10⁹ to 2.4×10¹¹ Ω/cm) and exhibit higher hardness, especially at low loads or plastic penetrations in the order of deposition bias voltage −10, −20 and −30 kV. At a lower contact load (1 N) and higher sliding speed (0.05 m/s), frictional properties in most cases improved, as did wear properties. At a higher contact load (5 N) and lower sliding speed (0.04 m/s), friction showed almost no improvement, and wear properties deteriorated. When the material of the counterbody was then changed from AISI 52100 to Ti–6Al–4V modified as the disc (contact load 5 N unchanged, sliding speed decreased), the friction coefficient decreased (but showed no improvement compared with the unmodified sample), while wear properties deteriorated further, and wear was changed from just the disc to both disc and ball, abrasive and adhesive dominated. Transfer films, mainly made up of wear debris transferred from the disc wear surfaces, were formed on the wear scars of the counterbodies. The deterioration of wear properties of the modified samples at the higher contact load is considered to be caused by the “thin ice” effect.     

The real area of contact and compliance of rough cylinders in compression

Bearing area analysis has been used to study the real area of contact and compliance of rough turned steel cylinders in compression. Calculations show that the elastic real area of contact is very small compared to the plastic real area of contact, and that local compliance due to flattening of asperity tips is a small proportion of the total compliance obtained from experiments. The fact that increased load brings more and more new asperities under load rather than enlarging the contact spots leads to a rather simple load-compliance relation for a rough cylinder, viz., W' = N_h · K₁δⁿ, where W₀ = K₁δⁿ defines the load-compliance relation of the individual asperities, and N_h represents the number of asperities bearing the load.     

A universal slip-line model with non-unique solutions for machining with curled chip formation and a restricted contact tool

A universal slip-line model and the corresponding hodograph for two-dimensional machining which can account for chip curl and chip back-flow when machining with a restricted contact tool are presented in this paper. Six major slip-line models previously developed for machining are briefly reviewed. It is shown that all the six models are special cases of the universal slip-line model presented in this paper. Dewhurst and Collins's matrix technique for numerically solving slip-line problems is employed in the mathematical modeling of the universal slip-line field. A key equation is given to determine the shape of the initial slip-line. A non-unique solution for machining processes when using restricted contact tools is obtained. The influence of four major input parameters, i.e. (a) hydrostatic pressure (P_A) at a point on the intersection line of the shear plane and the work surface to be machined; (b) ratio of the frictional shear stress on the tool rake face to the material shear yield stress (τ/k); (c) ratio of the undeformed chip thickness to the length of the tool land (t₁/h); and (d) tool primary rake angle (γ₁), upon five major output parameters, i.e. (a) four slip-line field angles (θ, η₁, η₂, ψ); (b) non-dimensionalized cutting forces (F_c/kt₁w and F_t/kt₁w); (c) chip thickness (t₂); (d) chip up-curl radius (R_u); and (e) chip back-flow angle (η_b), is theoretically established. The issue of the “built-up-edge” produced under certain conditions in machining processes is also studied. It is hoped that the research work of this paper will help in the understanding of the nature and the basic characteristics of machining processes.     

Kinematic hardening analysis of ratchet strain in the “pulley test”

The classical Bree problem—which represents an uniaxial model of a thin tube subjected to combined internal pressure and cyclic thermal stress across its wall—can be simulated by means of the pulley test in which a wire or strip specimen is subjected to combined steady tensile stress and cyclic bending stress. In this paper, accumulation of ratchet strain in the pulley test is investigated using a linear kinematic hardening material model from which perfect plasticity can be generated as a special case. The results of the investigation show that asymptotic ratchet strains are linearly related to the excess in mean stress σ_D above its value σ*_D at the ratchetting limit regardless of the thermal stress amplitude. Comparisons with test results on copper wire specimens—which exhibit non-linear hardening rate—confirm the qualitative validity of this simple relation. Deviations between theory and experiment are attributed to metallic cyclic creep. Further, perfect plasticity results are shown to be well predicted by a linearized lower bound estimate.     

“Antient” philosophy-science in Magna Graecia: Cassiodorus and Pythagoras and 20th c. literary travellers, Gissing and Douglas

The paper's philosophy-science time-scale stretches over a millenium, from 6th c. B.C. Pythagoras and into 6th c. A.D., Cassiodorus. Interest centres on Greek “science”, especially that of Pythagoras, in the colonies in Greater Greece — the foot of Italy — and on the writings of Capella and Boethius in the 5th c. A.D., as appreciated by Cassiodorus. His regard for the latters’ work and interest in its transmission to the universities in the Middle Ages became the Trivium and Quadrivium, the seven Liberal Arts-Sciences. The case is argued that Cassiodorus was probably the instigator of Benedictine regard for scholarship and was pivotal in facilitating the transmission of the latter “Seven Pillars of Wisdom”. Current general interest in the ancient classical scholarship of this region was stimulated by two early 20th c. literary travellers, Gissing and Douglas, whose writings, after visiting the region, we describe.     

Reference stresses and temperatures for cylinders and spheres under internal pressure with a steady heat flow in the radial direction

The paper examines the creep behavior of thick cylinders and spheres subjected to internal pressure and a negative temperature gradient in the radial direction. It is found that at stationary state the rate of radial displacement of the vessel wall is simply proportional to the material creep behavior associated with a single stress and temperature. Such “reference stresses” and “reference temperatures” are defined for spheres and cylinders of varying wall thicknesses. These reference stresses and reference temperatures are valid for any creep problem where the material behavior may be characterized by a function of the form exp (γT)σ^m. The extension of these results to variable pressure and temperature loading cases is discussed.     

尺度相关的分形结合面法向接触刚度模型

基于分形理论,利用双变量Weierstrass-Mandelbrot函数模拟三维分形结合面,建立尺度相关的三维分形结合面法向接触刚度模型。推导出各等级微凸体发生弹性、弹塑性以及完全塑性变形的存在条件。确定结合面上各等级微凸体的面积分布密度函数,推导出法向接触刚度和法向接触载荷的解析表达式。计算结果表明：当结合面上的微凸体只能发生弹性变形,即自身等级小于弹性临界等级的微凸体,该部分微凸体引起的法向接触刚度和对应法向载荷关系呈非线性。当微凸体的等级大于弹性临界等级,在结合面接触过程中,微凸体弹性变形引起的法向接触刚度与对应的法向载荷关系为线性,非弹性变形引起的法向接触刚度与法向载荷关系为非线性。微凸体的等级范围对结合面的刚度影响较大,在相同的法向载荷作用下,高等级微凸体的结合面产生较高的法向接触刚度,即结合面越平整,结合面的法向刚度越高。     

A new algorithm for computing the indentation of a rigid body of arbitrary shape on a viscoelastic half-space

In this paper the contact problem between a rigid indenter of arbitrary shape and a viscoelastic half-space is considered. Under the action of a normal force the penetration of the indenter and the distribution of contact pressure change. We wish to find the relations which link the pressure distribution, the resultant force on the indenter and the penetration on the assumption that the surfaces are frictionless. For indenters of arbitrary shape the problem may be solved numerically by using the Matrix Inversion Method (MIM), extended to viscoelastic case. In this method the boundary conditions are satisfied exactly at specified “matching points” (the mid-points of the boundary elements). It can be validated by comparing the numerical results to the analytic solutions in cases of a spherical asperity (loading and unloading) and a conical asperity (loading only). Finally, the method was implemented for a finite cylindrical shape with its curved face indenting the surface of the half-space. This last example shows the efficiency of the method in case of a prescribed penetration as well as a given normal load history.     

Pressure between elastic bodies having a slender area of contact and arbitrary profiles

A numerical method is presented for calculating the pressure distribution and contact area shape between two elastic bodies of arbitrary profile which make contact over a slender contact area, i.e. where the relative curvature of the two profiles is much smaller in the longitudinal direction than in the transverse. The pressure distribution is assumed to be piecewise-linear in the longitudinal direction and semi-elliptical in the transverse. No a priori relationship is assumed between the shape of the contact area and the longitudinal variation in pressure; they are found simultaneously from dual integral equations for the compatibility of (a) the normal displacement and (b) the transverse curvature along the longitudinal axis of the contact zone.In cases where the profiles of the contacting bodies are smooth and continuous up to, and beyond, the ends of the contact area, the method gives a very reliable measure of the contact pressure distribution. Where discontinuities in profile are present, at roller ends for example, stress singularities are to be expected and like any numerical method, only approximate values of the stress concentration can be found. In the cases studied, the concentration of pressure associated with a “sharp” edge of contact is found to be very local.The method has been applied to both cylindrical and variously “crowned” rollers, also to a ball “over-riding” the edge of a closely conforming groove.     

Tension transmission via an elastic rod gripped by two circular-edged plates

The equilibrium of static friction in tension transmission by an inextensible elastic rod gripped by two plates with round edges was modeled and solved analytically. To proceed, we established equilibrium equations in the deflected and circular edged regions, respectively, and then combined the results together. In both deflected and circular edged region, the condition of cantilever beam with clamped end maintaining a continuous contact was assumed to avoid corresponding mathematical complexity. As a result, we expressed the gripping force in terms of the contact angle, the inclined angle of load, the radius ratio and the initial tension in our results.More specifically, a more precise but still analytical relationship between the incoming and outgoing tensions including the effect of the rod bending rigidity was derived through the analysis of the circular edged region; it turned out that the radius ratio is the only parameter representing the bending rigidity in the tension ratio. To analyze the deflected region, the classical elastica analysis in terms of an elliptic integral was adopted. Based on the results, the effect of radius ratio and frictional coefficient on the derived force ratio was investigated. As a result, the bending rigidity can be ignored above the range of ρ≈100. But the tension transmission decreases by 29.1% relative to the classical case at the range of ρ≈5. In addition, the effect of inclined load on the derived gripping force was also investigated.     

New cylinder liner surfaces for low oil consumption

Anticipated emission legislation and reduced fuel consumption are the main driving forces when developing new engines. Optimization of the active surfaces in the piston system is one possible way to meet the above demands. In this study the effects of surface topography and texture direction of the ring/liner contact on oil film thickness and friction were simulated and experimentally tested. “Low wear” results from the experimental wear tests with “glide honed” smooth liner surfaces supported the “low friction” simulation results. In addition a new wear volume sensitive surface roughness parameter, R_ktot, based on the Abbot–Firestone bearing area curve was introduced.     

Buckling of thin cylindrical shells: an attempt to resolve a paradox

The classical theory of buckling of axially loaded thin cylindrical shells predicts that the buckling stress is directly proportional to the thickness t, other things being equal. But empirical data show clearly that the buckling stress is actually proportional to t^1.5, other things being equal. As is well known, there is wide scatter in the buckling-stress data, going from one half to twice the mean value for a given ratio R/t. Current theories of shell buckling explain the low buckling stress—in comparison with the classical—and the experimental scatter in terms of “imperfection-sensitive”, non-linear behaviour. But those theories always take the classical analysis of an ideal, perfect shell as their point of reference.Our present principal aim is to explain the observed t^1.5 law. So far as we know, no previous attack has been made on this particular aspect of thin-shell buckling. Our work is thus breaking new ground, and we shall deliberately avoid taking the classical analysis as our starting point.We first point out that experiments on self-weight buckling of open-topped cylindrical shells agree well with the mean experimental data mentioned above; and then we associate those results with a well-defined post-buckling “plateau” in load/deflection space, that is revealed by finite-element studies. This plateau is linked with the appearance of a characteristic “dimple” of a mainly inextensional character in the deformed shell wall. A somewhat similar post-buckling dimple is also found by quite separate finite-element studies when a thin cylindrical shell is loaded axially at an edge by a localised force; and it turns out that such a dimple grows under a more-or-less constant force that is proportional to t^2.5, other things being equal.This 2.5-power law can be explained by analogy with the inversion of a thin spherical shell by an inward-directed force. Thus, the deformation of such a shell is generally inextensional except for a narrow “knuckle” or boundary layer in which the combined local elastic energy of bending and stretching is proportional to t^2.5, other things being equal. Similarly, the modes of deformation in the post-buckling dimples in a cylindrical shell are practically independent of thickness, except in the highly deformed boundary-layer regions which separate the inextensionally distorted portions of the shell. These ideas lead in turn to an explanation of the t^1.5 law for the post-buckling stress of open-topped cylindrical shells loaded by their own weight.We attribute the absence of experimental scatter in the self-weight buckling of open-topped cylindrical shells to the statical determinacy of the situation, which allows a post-buckling dimple to grow at a well-defined “plateau load”. Conversely, the large experimental scatter in tests on cylinders with closed ends may be attributed to the lack of statical determinacy there.Our paper contains several arguments that are not mathematically water-tight, in contrast to many reports in the field of mechanics of structures. We plead that the problem which we have tackled is so difficult that the only way forward is one of “over-simplification”. We hope that our work will be judged not with respect to its absence of mathematical precision, but by the light which it sheds upon the problem under investigation.     

Predictive model to estimate the stress–strain curves of bulk metals using nanoindentation

Many studies have shown that finite element modeling (FEM) can be used to fit experimental load–displacement data from nanoindentation tests. Most of the experimental data are obtained with sharp indenters. Compared to the spherical case, sharp tips do not directly allow the behavior of tested materials to be deduced because these produce a nominally-constant plastic strain impression. The aim of this work is to construct with FEM an equivalent stress–strain response of a material from a nanoindentation test, done with a pyramidal indenter. The procedure is based on two equations which link the parameters extracted from the experimental load–displacement curve with material parameters, such as Young's modulus E, yield stress Y₀ and tangent modulus E_T. We have already tested successfully the relations on well-known pure metallic surfaces. However, the load–displacement curve obtained using conical or pyramidal indenters cannot uniquely determine the stress–strain relationship of the indented material. The non-uniqueness of the solution is due to the existence of a characteristic point (ε_c, σ_c); for a given elastic modulus, all bilinear stress–strain curves that exhibit the same true stress σ_c at the specific true strain ε_C lead to the same loading and unloading indentation curve. We show that the true strain ε_c is constant for all tested materials (Fe, Zn, Cu, Ni), with an average value of 4.7% for a conical indenter with a half-included angle θ=70.3°. The ratio σ_c/ε_c is directly related to the elastic modulus of the indented material and the tip geometry.     

Asymptotic analysis of extrusion of porous metals

Plane strain extrusion of fully dense and porous metals is analysed using asymptotic techniques. The extrusion die is assumed to taper gradually down the extrusion axis. The asymptotic expansions are based on a small parameter ε which is defined as the ratio of the total reduction of the original cross-section to the length of the reduction region. Coulomb's law is used to model the frictional forces that develop along the metal-die interface and the coefficient of friction is assumed to be of order ε. Analytical solutions for the first two terms in the expansions are obtained. In the case of the fully dense metals, it is shown that the leading order [O(1)] solution involves “slab flow.” It is also shown that the next term in the expansion of the solution is O(ε²), and this provides a theoretical justification for the use of the so-called “slab methods” of analysis for dies of moderate slope. An asymptotic analysis of the extrusion of porous metals with dilute concentration of voids is also carried out. Gurson's plasticity model is used to describe the constitutive behavior of the material. The leading order solution is the same as that of the fully dense material and the effects of porosity enter as an O(ε) correction. In order to verify the asymptotic solutions developed, detailed finite element calculations are carried out for both the fully dense and the porous material. The asymptotic solutions agree well with the results of the finite element calculations.     

A simple relation in elastic contact problems

In this paper, a new relation is introduced that simplifies the determination of the Muskhelishvili's potential function in plane contact problems. The relation is Φ(z)=1/2[p(z)-iq(z)], which is correct for all uncoupled contact problems. The new relation is proved in a mathematical way and utilized to obtain the potential function in several contact problems. A complete agreement has been observed between our results and the potential functions that have been obtained from complicated methods in the past. Utilization of the new relation simplifies the solution of contact problems and analytical calculation of the stress and displacement fields, which may lead to the design of better cutting tools or fretting fatigue test pads.     

Kinetics of adhesion on a viscoelastic sample by force microscopy

Indentation of the tip of an atomic force microscope is carried out under purely adhesive force, on a sample that is in a viscoelastic state. A characteristic size of the contact area is also measured by recording the tangential force while a small lateral modulation of the sample position is applied. Both the contact size and the penetration depth of the tip follow the compliance function of the material. It is shown that the kinetics of the adhesion can be described as the adhesive equilibrium of the tip on an elastic material whose compliance is the instantaneous value of the compliance function of the viscoelastic material. Implications for the analysis of the force–distance curves are discussed.