The axisymmetric contact problem of the theory of elasticity for a half-space in the presence of body forces

We consider the axisymmetric problem of determination of the stress-strain state in an elastic half-space in the case of a circular line of separation of the boundary conditions on the boundary plane z = 0. We assume that on the entire boundary z = 0 the tangential stress _rz = 0, while inside the circle r a (z = 0) the normal displacement u_z is known and in its exterior the normal stress _z is given. In addition, we assume that body forces are acting in the half-space. The investigation of problems of similar kind presents interest in connection with the application of A. A. Il'yushin's method of elastic solutions to the problem of the indentation of punches into a nonlinear-elastic, in particular, into an elastoplastic half-space.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 110–115, July–August, 1971.The author wishes to thank M. Ya. Leonov for his valuable suggestions during the preparation of this paper.     

GREEN FUNCTIONS FOR A DECAGONAL QUASICRYSTALLINE MATERIAL WITH A PARABOLIC BOUNDARY

This investigation presents the Green functions for a decagonal quasicrystalline material with a parabolic boundary subject to a line force and a line dislocation by means of the complex variable method. The surface Green functions are treated as a special case, and the explicit expressions of displacements and hoop stress at the parabolic boundary are also given.Finally, the stresses and displacements induced by a phonon line force acting at the origin of the lower half-space are presented.     

Thermoelastostatic equilibrium of a spatial quadrant: Green’s function and solution in integrals

In this study, a new Green??s function and a new Green-type integral formula for a 3D boundary value problem (BVP) in thermoelastostatics for a quarter-space are derived in closed form. On the boundary half-planes, twice mixed homogeneous mechanical boundary conditions are given. One boundary half-plane is free of loadings and the normal displacements and the tangential stresses are zero on the other one. The thermoelastic displacements are subjected by a heat source applied in the inner points of the quarter-space and by mixed non-homogeneous boundary heat conditions. On one of the boundary half-plane, the temperature is prescribed and the heat flux is given on the other one. When the thermoelastic Green??s function is derived, the thermoelastic displacements are generated by an inner unit point heat source, described by ??-Dirac??s function. All results are obtained in elementary functions that are formulated in a special theorem. As a particular case, when one of the boundary half-plane of the quarter-space is placed at infinity, we obtain the respective results for half-space. Exact solutions in elementary functions for two particular BVPs for a thermoelastic quarter-space and their graphical presentations are included. They demonstrate how to apply the obtained Green-type integral formula as well as the derived influence functions of an inner unit point body force on volume dilatation to solve particular BVPs of thermoelasticity. In addition, advantages of the obtained results and possibilities of the proposed method to derive new Green??s functions and new Green-type integral formulae not for quarter-space only, but also for any canonical Cartesian domain are also discussed.     

Determination of all stress components from measurements of the stress invariant by the thermoelastic stress method

A method for the resolution of all stress components from the first invariant J₁ measured by thermoelastic stress analyzer is described. This method may be used to determine, not only surface stress, but also internal stress and stress on the underside.The method is based on the following procedure:

1. (1) Pick an arbitrary domain Ω, within the structure, for which the stresses are required.

2. (2) Measure J₁ on the surface of Ω.

3. (3) Determine the optimum traction along the boundary Γ, which is a part of Ω, by the least squares method such that the difference between the measured J₁ and the calculated J₁ is at a minimum. Either FEM or BEM may be used for this calculation.

Examples of stress resolution for a two-dimensional stress concentration problem and a three-dimensional stress concentration problem are shown. The accuracy of the stress resolution is discussed.     

Three-dimensional solutions for general anisotropy

The Stroh formalism is extended to provide a new class of three-dimensional solutions for the generally anisotropic elastic material that have polynomial dependence on x₃, but which have quite general form in x₁,x₂. The solutions are obtained by a sequence of partial integrations with respect to x₃, starting from Stroh's two-dimensional solution. At each stage, certain special functions have to be introduced in order to satisfy the equilibrium equation. The method provides a general analytical technique for the solution of the problem of the prismatic bar with tractions or displacements prescribed on its lateral surfaces. It also provides a particularly efficient solution for three-dimensional boundary-value problems for the half-space. The method is illustrated by the example of a half-space loaded by a linearly varying line force.     

Mechanics of adhesive contact on a power-law graded elastic half-space

We consider adhesive contact between a rigid sphere of radius R and a graded elastic half-space with Young's modulus varying with depth according to a power law E=E₀(z/c₀)^k (0<k<1) while Poisson's ratio ν remaining a constant. Closed-form analytical solutions are established for the critical force, the critical radius of contact area and the critical interfacial stress at pull-off. We highlight that the pull-off force has a simple solution of P_cr=−(k+3)πRΔγ/2 where Δγ is the work of adhesion and make further discussions with respect to three interesting limits: the classical JKR solution when k=0, the Gibson solid when k→1 and ν=0.5, and the strength limit in which the interfacial stress reaches the theoretical strength of adhesion at pull-off.     

Subinterface crack in an anisotropic piezoelectric bimaterial

In this paper, the problem of a subinterface crack in an anisotropic piezoelectric bimaterial is analyzed. A system of singular integral equations is formulated for general anisotropic piezoelectric bimaterial with kernel functions expressed in complex form. For commonly used transversely isotropic piezoelectric materials, the kernel functions are given in real forms. By considering special properties of one of the bimaterial, various real kernel functions for half-plane problems with mechanical traction-free or displacement-fixed boundary conditions combined with different electric boundary conditions are obtained. Investigations of half-plane piezoelectric solids show that, particularly for the mechanical traction-free problem, the evaluations of the mechanical stress intensity factors (electric displacement intensity factor) under mechanical loadings (electric displacement loading) for coupled mechanical and electric problems may be evaluated directly by considering the corresponding decoupled elastic (electric) problem irrespective of what electric boundary condition is applied on the boundary. However, for the piezoelectric bimaterial problem, purely elastic bimaterial analysis or purely electric bimaterial analysis is inadequate for the determination of the generalized stress intensity factors. Instead, both elastic and electric properties of the bimaterial’s constants should be simultaneously taken into account for better accuracy of the generalized stress intensity factors.     

Small edge crack in a semi-infinite solid subjected to thermal stratification

The mode I stress intensity factor for a small edge crack in an elastic half-space is found when the space is in contact with two stratified fluids of different temperatures, the boundary between the fluids oscillating sinusoidally over the solid surface. The variation in the stress intensity factor, which may lead to thermal fatigue crack growth, is examined as a function of time, crack depth, amplitude and temporal frequency of oscillation, surface heat transfer coefficient and material properties of the half-space. It is shown how this ‘boundary layer’ solution may be applied to problems involving finite geometries.     

Effect of rigid boundary on the propagation of torsional waves in a homogeneous layer over a heterogeneous half-space

In the present paper we study the effect of rigid boundary on the propagation of torsional waves in a homogeneous layer over a semi-infinite heterogeneous half-space, where the heterogeneity is both in rigidity and density. The present study demonstrates that torsional waves can propagate in the layer. The velocities of torsional waves have been calculated numerically as a functions of KH, (where K is the wave number and H is the thickness of the layer) and are presented in a number of graphs. It is also observed that, for a layer over a homogeneous half-space, the velocity of torsional waves does not coincide with that of Love waves in the presence of the rigid boundary whereas it does at the free boundary.     

Influence of stress triaxiality and temperature on void growth and ductile/brittle transition behavior of 40 Cr steel

Examined experimentally are the influence of stress triaxiality and temperature on the growth of microvoids and the ductile/brittle transition (DBT) macrobehavior of 40 Cr steel subjected to two different heat treatments. This is accomplished by testing more than 300 smooth and notched specimens over a temperature range of 20°C to −196°C. Changes in the microstructure morphology are examined by scanning electron microscopy (SEM) and identified with fracture data on a surface constructed from the uniaxial strain ε_c at fracture, the stress triaxiality R_σ and the temperature T. While stress triaxiality has a significant influence on the DBT temperature T_c, it does not affect the ratio of the average radius of voids R_o to that of inclusions R_i. The ratio R_o/R_i is found to increase with temperature and remains constant in specimens with different notch radii regardless of the temperature. Empirical relations between T_c and R_σ⁻ and R_o/R_i and T are proposed to better understand how macrofracture parameters are influenced by microstructure entities.     

DNS of a spatially developing turbulent boundary layer with passive scalar transport

A direct numerical simulation (DNS) of a spatially developing turbulent boundary layer over a flat plate under zero pressure gradient (ZPG) has been carried out. The evolution of several passive scalars with both isoscalar and isoflux wall boundary condition are computed during the simulation. The Navier–Stokes equations as well as the scalar transport equation are solved using a fully spectral method. The highest Reynolds number based on the free-stream velocity U_∞ and momentum thickness θ is Re_θ=830, and the molecular Prandtl numbers are 0.2, 0.71 and 2. To the authors’ knowledge, this Reynolds number is to date the highest with such a variety of scalars. A large number of turbulence statistics for both flow and scalar fields are obtained and compared when possible to existing experimental and numerical simulations at comparable Reynolds number. The main focus of the present paper is on the statistical behaviour of the scalars in the outer region of the boundary layer, distinctly different from the channel-flow simulations. Agreements as well as discrepancies are discussed while the influence of the molecular Prandtl number and wall boundary conditions is also highlighted. A Pr scaling for various quantities is proposed in outer scalings. In addition, spanwise two-point correlation and instantaneous fields are employed to investigate the near-wall streak spacing and the coherence between the velocity and the scalar fields. Probability density functions (PDF) and joint probability density functions (JPDF) are shown to identify the intermittency both near the wall and in the outer region of the boundary layer. The present simulation data will be available online for the research community.     

Local weight functions for semi-elliptical surface cracks in finite thickness plates

Weight functions for any local point, 0 < Φ < π/2 along a semi-elliptical surface crack in finite thickness plates were derived from an assumed approximate general weight function and two reference stress intensity factors. The resulting weight functions were verified using available finite element results for two nonlinear stress fields and good agreement was achieved. When used together with weight functions for Φ = 0 and Φ = π/2 the weight functions are suitable for the calculation of stress intensity factors anywhere along the crack front for semi-elliptical surface cracks in complex stress fields with aspect ratios in the range 0 ≤ a/c ≤ 1 and relative depths 0 ≤ a/t ≤ 0.8.     

Stability of a melt/solid interface with reference to a plume channel

The structure of free-convection flow in a plume channel formed as a result of melting above a local heat source placed on the basement of a solid mass is experimentally investigated. The channel shape and the flow pattern in it are functions of the relative power Ka = N/N ₁, where N is the plume source power and N ₁ is the heat removed to the surrounding mass. When the heat is withdrawn from the plume channel by heat conduction, the channel represents a system of convective cells on whose boundaries there are channel constrictions. The temperature fields and the cell flow patterns are investigated. For mantle plumes, such as the Hawaiian, Iceland, and Bouvet plumes and extended igneous provinces, the basement diameter and the values of the criterion Ka are determined.     

Magneto-elastic stress induced by an electric current in an infinite thin plate with an elliptical hole

Two-dimensional magnetic field and stress analyses have been presented for soft ferromagnetic, paramagnetic, and diamagnetic materials of an infinite thin plate with an elliptical hole under steady electric current. The magnetic stress has been analyzed in the Maxwell Stress Model. Except for the approximation of the plane stress state since the plate is the thin plate, any assumption is not made for the stress analysis, though the Maxwell stress components are expressed by nonlinear terms. The boundary condition expressed by Maxwell’s stress is completely satisfied without any linear assumptions on the boundary. Two ways for the boundary condition are stated. The analysis of σ _z in the direction of the plate thickness is also carried out. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived, and the magnitude of the stress intensity factor for the magnetic stress and thermal stress due to the Joule heat caused by the electric current is discussed.     

Scattering of guided waves by through-thickness cavities with irregular shapes

This paper investigates the three-dimensional (3D) scattering of guided waves by a through-thickness cavity with irregular shape in an isotropic plate. The scattered field is decomposed on the basis of Lamb and SH waves (propagating and non-propagating), and the amplitude of the modes is calculated by writing the nullity of the total stress at the boundary of the cavity. In the boundary conditions, the functions depend on the through-thickness coordinate, z, but contrary to the case where the cavity has a circular shape, they also depend on the angular coordinate θ. This is dealt with by projecting the z-dependent functions onto a basis of orthogonal functions, and by expanding the θ-dependent functions in Fourier series. Examples include the scattering of the S₀, SH₀ and A₀ modes by elliptical cavities with different values of aspect ratio, and the scattering of the S₀ mode by a cavity with an arbitrary shape. Results obtained with this model are compared with results obtained with the finite element (FE) method, showing very good agreement.     

On thermodynamic potentials in linear thermoelasticity

The four thermodynamic potentials, the internal energy u=u(ε_ij,s), the Helmholtz free energy f=f(ε_ij,T), the Gibbs energy g=g(σ_ij,T) and the enthalpy h=h(σ_ij,s) are derived, independently of each other, by using the Duhamel–Neumann extension of Hooke's law and an assumed linear dependence of the specific heat on temperature. A systematic procedure is then presented to express all thermodynamic potentials in terms of four possible pairs of independent state variables. This procedure circumvents a tedious transition from one potential to another, based on the formal change of variables, and inversions of the stress–strain and entropy–temperature relations. The general results are applied to uniaxial loading paths under isothermal, adiabatic, constant stress, and constant strain conditions. An interplay of adiabatic and isothermal elastic constants in the expressions for exchanged heat along certain thermodynamic paths is indicated.     

Influence of roughness on attachment line boundary layer transition in hypersonic flow

 Attachment line boundary layer transition on swept cylinders is studied in a low enthalpy hypersonic wind tunnel at M _∞=7.14. Sweep angles of 60° and 70° are used and transition is detected by means of heat flux measurements. The influence on attachment line transition of single 2D-roughness elements, in the form of tripwires or slots, as well as 3D obstacles is determined and the results are analyzed with respect to Poll’s criterion. Received: 16 January 1996 / Accepted: 12 July 1996     

Effect of surface layers on the constriction resistance of an isothermal spot. Part I: Reduction to an integral equation and numerical results

The problem of the constriction resistance of a circular spot on a half-space covered with a uniform layer of different material is considered. For the general case of any specified axisymmetric distributions of temperature over the spot and heat flux over the rest of the surface, the mixed boundary value problem governing the heat flow from the spot to the underlying layer-substrate composite is converted to a non-homogeneous Fredholm integral equation of the second kind. For the particular case of isothermal spot on otherwise insulated surface, the evaluation of the constriction resistance is reduced to the reciprocal of a simple integral with the solution of the relevant integral equation as integrand. The integral equation is solved numerically and very accurate results are obtained for the constriction resistance over four orders of magnitude variation of the ratio, , of layer thickness to spot radius and the ratio, k_r, of layer to substrate conductivities for both conducting (k_r > 1) and insulating (k_r < 1) layers. An extensive discussion of the numerical results is presented with particular emphasis on their implications for the contact resistance of practical joints in the presence of interfacial layers. Further, in the light of the numerical results, two widely used analytical approximations for the constriction resistance – the first of which results from replacing the isothermal condition over the spot by a special flux (herein called the Equivalent Isothermal Heat Flux, EIHF) condition which is believed to render the spot nearly isothermal and the second is a consequence of the assumption (herein termed the Thin Insulating Layer Approximation, TILA) that, for thin insulating layers like oxide films, the heat flow in the layer region right beneath the spot is purely axial – are assessed as to their levels of accuracy and ranges of applicability with respect to both and k_r.     

Radiative and Thermal Dispersion Effects on Non-Darcy Natural Convection with Lateral Mass Flux for Non-Newtonian Fluid from a Vertical Flat Plate in a Saturated Porous Medium

The problem of natural convective heat transfer for a non-Newtonian fluid from an impermeable vertical plate embedded in a fluid-saturated porous medium has been analyzed. Non-Darcian, radiative and thermal dispersion effects have been considered in the present analysis. The governing boundary layer equations and boundary conditions are cast into a dimensionless form and simplified by using a similarity transformation. The resulting system of equations is solved by using a double shooting Runge–Kutta method. The effect of viscosity index n, the conduction–radiation parameter R, the non-Darcy parameter Gr*, the thermal dispersion parameter Ds and the suction/injection parameter f_w on the fluid velocities, temperatures and the local Nusselt number are discussed.     

Computation of conjugate heat transfer in the turbulent mixed convection regime in a vertical channel with multiple heat sources

This paper presents the results of a comprehensive numerical study to analyze conjugate, turbulent mixed convection heat transfer from a vertical channel with four heat sources, uniformly flush-mounted to one of the channel walls. The results are presented to study the effect of various parameters like thermal conductivity of wall material (k _s), thermal conductivity of flush-mounted discrete heat source (k _c), Reynolds number of fluid flow (Re _s), modified Richardson number (Ri ⁺) and aspect ratio (AR) of the channel. The standard k-ε turbulence model, modified by including buoyancy effects with physical boundary conditions, i.e. without wall functions, has been used for the analysis. Semi-staggered, non-uniform grids are used to discretise the two dimensional governing equations, using finite volume method. A correlation, encompassing a wide range of parameters, is developed for the non-dimensional maximum temperature (T ^*) using the asymptotic computational fluid dynamics (ACFD) technique.