A numerical analysis of partial slip problems under Hertzian contacts

In this study the tangential partial slip problems in Hertzian contact regions are treated by a numerical technique. The tangential loading may include tangential forces in the contact plane and a twisting moment normal to the contact plane. The Coulomb’s law of friction and the property that the direction of friction must oppose the relative motion lead to nonlinear equations. The Newton-Raphson method is utilized to solve these nonlinear equations. Numerical results for tangential tractions and sizes of stick and slip zones may be determined, and they agree with existing analytical results for circular contacts.     

Breakage mechanics—Part I: Theory

Different measures have been suggested for quantifying the amount of fragmentation in randomly compacted crushable aggregates. A most effective and popular measure is to adopt variants of Hardin's [1985. Crushing of soil particles. J. Geotech. Eng. ASCE 111(10), 1177-1192] definition of relative breakage ‘B_r’. In this paper we further develop the concept of breakage to formulate a new continuum mechanics theory for crushable granular materials based on statistical and thermomechanical principles. Analogous to the damage internal variable ‘D’ which is used in continuum damage mechanics (CDM), here the breakage internal variable ‘B’ is adopted. This internal variable represents a particular form of the relative breakage ‘B_r’ and measures the relative distance of the current grain size distribution from the initial and ultimate distributions. Similar to ‘D’, ‘B’ varies from zero to one and describes processes of micro-fractures and the growth of surface area. However, unlike damage that is most suitable to tensioned solid-like materials, the breakage is aimed towards compressed granular matter. While damage effectively represents the opening of micro-cavities and cracks, breakage represents comminution of particles. We term the new theory continuum breakage mechanics (CBM), reflecting the analogy with CDM. A focus is given to developing fundamental concepts and postulates, and identifying the physical meaning of the various variables. In this part of the paper we limit the study to describe an ideal dissipative process that includes breakage without plasticity. Plastic strains are essential, however, in representing aspects that relate to frictional dissipation, and this is covered in Part II of this paper together with model examples.     

Mie-Grneisen状态方程可压缩多流体流动的PPM方法

采用流体体积分数的混合型多流体数值模型，将piecewise parabolic method (PPM)方法应用于可压缩多流体流动的数值模拟，拓展了以前提出的模型和数值方法，使它能够处理一般的Mie-Grneisen状态方程。采用双波近似和两层迭代算法求解一般状态方程的Riemann问题；并根据多流体接触界面无振荡原则设计高精度计算格式，对典型的纯界面平移问题可以从理论上证明本算法在接触间断附近压力和速度没有振荡，而且数值模拟结果表明界面数值耗散也被控制在2~3个网格之内。模拟了多种复杂的可压缩多流体流动，算例结果表明本文方法可以有效地处理接触间断、激波等物理问题，且具有耗散小精度高的特点。     

A Technique for Dynamic Tensile Testing of Human Cervical Spine Ligaments

An experimental study is undertaken to examine the dynamic stress–strain characteristics of ligaments from the human cervical spine (neck). Tests were conducted using a tensile split Hopkinson bar device and the engineering strain rates imposed were of the order of 10²∼10³/s. As ligaments are extremely soft and pliable, specialized test protocols applicable to Hopkinson bar testing were developed to facilitate acquisition of reliable and accurate data. Seven primary ligaments types from the cervical spines of three male cadavers were subjected to mechanical tests. These yielded dynamic stress–strain curves which could be approximated by empirical equations. The dynamic failure stress/load, failure stain/deformation, modulus/stiffness, as well as energy absorption capacity, were obtained for the various ligaments and classified according to their location, the strain rate imposed and the cadaveric source. Compared with static responses, the overall average dynamic stress–strain behavior foreach type of ligament exhibited an elevation in strength but reduced elongation.     

Hertzian fracture at unloading

Hertzian fracture through indentation of flat float glass specimens by steel balls has been examined experimentally. Initiation of cone cracks has been observed and failure loads together with contact and fracture radii determined at monotonically increasing load but also during unloading phases. Contact of dissimilar elastic solids under decreasing load may cause crack inception triggered by finite interface friction and accordingly the coefficient of friction was determined by two different methods. In order to make relevant predictions of experimental findings, a robust computational procedure has been developed to determine global and local field values in particular at unloading at finite friction. It was found that at continued loading it is possible to specify in advance how the contact domain divides into invariant regions of stick and slip. The maximum tensile stress was found to occur at the free surface just outside the contact contour, the relative distance depending on the different elastic compliance properties and the coefficient of friction. In contrast, at unloading invariance properties are lost and stick/slip regions proved to be severely history dependant and in particular with an opposed frictional shear stress at the contact boundary region. This causes an increase of the maximum tensile stress at the contour under progressive unloading. Predictions of loads to cause crack initiation during full cycles were made based on a critical stress fracture criterion and proved to be favourable as compared to the experimental results.     

Dynamic dip testing as a method to assess the condensate drainage behavior from the air-side surface of compact heat exchangers

A new method to assess the condensate drainage behavior of the air-side surface of compact heat exchangers—dynamic dip testing—is introduced. The new method is shown to provide highly repeatable data for real-time drainage. Results from experiments with more than 20 flat-tube and round-tube-and-fin heat exchangers are presented, and the data clearly show geometrical effects such as the impact of the tube type on condensate drainage. By comparing the results from dip testing to wind-tunnel experiments for the same heat exchangers, we find dip testing can serve as a powerful tool for assessing the condensate retention behavior. The coils retaining the most and the least condensate in a steady-state wind-tunnel test, likewise held the most and the least in a dip test. However, different amounts of water are retained on the air-side surface during dip tests and wind-tunnel tests. A model based on gravity, surface tension and drag effects is developed to help understand and predict the drainage behavior of heat exchangers. The new model and experimental approach are useful in screening heat exchangers for condensate retention and for assessing off-cycle drainage behavior.     

Characteristic-value analysis for plastic dynamic buckling of columns under elastoplastic compression waves

The characteristic-value analysis of plastic dynamic buckling is presented for columns under the action of elastoplastic compression wave caused by an axial-step load. Two critical conditions constituting a dynamic instability criterion are derived on the basis of transformation and conservation of energy. The governing equations, the boundary conditions and the continuity conditions derived by the use of the first critical condition are the same as those given by the adjacent-equilibrium criterion and are insufficient for determining two characteristic parameters involved in the governing equations. A supplementary restraint equation for buckling deformations at the plastic-wave front and the elastic-wave front is derived by the use of the second critical condition. Then, a couple of characteristic equations for two characteristic parameters are derived on the condition that the governing equations have non-trivial solutions satisfying the boundary conditions, the continuity conditions and the supplementary restraint equation. The critical-load parameters, dynamic characteristic parameter (exponent parameter of inertia term) and dynamic buckling modes are calculated from the solutions of the characteristic equations.     

On the Lagrange problem about the optimal form for circular hollow columnsSur le problème de Lagrange de la forme optimale d'une colonne circulaire aux parois minces

A new approach to the problem of finding the form of the strongest clamped circular column with thin walls of fixed volume and height is proposed. The same model describes also the form of the horizontal beam with rectangular vertical sections of a fixed height and variable widths. The existence and the uniqueness of the optimal column are proved here for the first time and the optimal form is found. To cite this article: Yu.V. Egorov, C. R. Mecanique 331 (2003).     

Recovery of small electromagnetic inhomogeneities from partial boundary measurementsIdentification de petites inhomogénéités diélectriques à partir de mesures partielles

We consider for the inverse problem of identifying locations and certain properties of the shapes of small dielectric inhomogeneities in a homogeneous background medium from boundary measurements on part of the boundary or dynamic boundary measurements for a finite time interval. Using as weights particular background solutions we develop asymptotic methods based on appropriate averaging of the data. To cite this article: H. Ammari, A.G. Ramm, C. R. Mecanique 330 (2002) 199–205.     

Flow past an axisymmetric body embedded in a saturated porous mediumÉcoulement derrière un corps à symétrie axiale plongé dans un milieu poreux saturé

The problem of viscous fluid past an axisymmetric body embedded in a fluid saturated porous medium is studied using the Brinkman's extension. A general formula for the drag on the body is derived in the form of a limit of an expression involving the stream function characterizing the flow. The flow past an axisymmetric approximate sphere is also considered. The stream function in this case is obtained in terms of Bessel functions and Gegenbauer's functions. The drag acting on the body is evaluated by using the formula derived. Its variation is studied with respect to geometric and permeability parameters. The special cases of flow past a sphere and a spheroid are obtained from the present analysis. To cite this article: D. Srinivasa Charya, J.V. Ramana Murthy, C. R. Mecanique 330 (2002) 417–423.