Observations on the theoretical and experimental foundations of thermoplasticity

The thermodynamics of materials with internal state variables has been employed to study the properties of a class of thermoplastic materials in which the evolution equation for the internal variables is given by equation k⁽ⁱ⁾ = g⁽ⁱ⁾(σ_kl, '_kl, k⁽ⁱ⁾, θ, · '_kl) where g⁽ⁱ⁾ is homogeneous of order one in · '_kl. The most general form of the Helmholtz potential consistent with the assumption of insensitivity of the elastic relations to inelastic deformation has been derived and a geometric interpretation of the Clausius-Duhem restriction has been made employing the concept of a thermodynamic reference stress. Experimental results of one of the authors have been correlated with the theory.     

A field model interpretation of crack initiation and growth behavior in ferroelectric ceramics: change of poling direction and boundary condition

Strain energy density expressions are obtained from a field model that can qualitatively exhibit how the electrical and mechanical disturbances would affect the crack growth behavior in ferroelectric ceramics. Simplification is achieved by considering only three material constants to account for elastic, piezoelectric and dielectric effects. Cross interaction of electric field (or displacement) with mechanical stress (or strain) is identified with the piezoelectric effect; it occurs only when the pole is aligned normal to the crack. Switching of the pole axis by 90° and 180° is examined for possible connection with domain switching. Opposing crack growth behavior can be obtained when the specification of mechanical stress σ^∞ and electric field E^∞ or (σ^∞,E^∞) is replaced by strain ε^∞ and electric displacement D^∞ or (ε^∞,D^∞). Mixed conditions (σ^∞,D^∞) and (ε^∞,E^∞) are also considered. In general, crack growth is found to be larger when compared to that without the application of electric disturbances. This includes both the electric field and displacement. For the eight possible boundary conditions, crack growth retardation is identified only with (E_y^∞,σ_y^∞) for negative E_y^∞ and (D_y^∞,ε_y^∞) for positive D_y^∞ while the mechanical conditions σ_y^∞ or ε_y^∞ are not changed. Suitable combinations of the elastic, piezoelectric and dielectric material constants could also be made to suppress crack growth.     

Crack size and speed interaction characteristics at micro-, meso- and macro-scale

The motivation to examine physical events at even smaller size scale arises from the development of use-specific materials where information transfer from one micro- or macro-element to another could be pre-assigned. There is the growing belief that the cumulated macroscopic experiences could be related to those at the lower size scales. Otherwise, there serves little purpose to examine material behavior at the different scale levels. Size scale, however, is intimately associated with time, not to mention temperature. As the size and time scales are shifted, different physical events may be identified. Dislocations with the movements of atoms, shear and rotation of clusters of molecules with inhomogeneity of polycrystals; and yielding/fracture with bulk properties of continuum specimens. Piecemeal results at the different scale levels are vulnerable to the possibility that they may be incompatible. The attention should therefore be focused on a single formulation that has the characteristics of multiscaling in size and time. The fact that the task may be overwhelmingly difficult cannot be used as an excuse for ignoring the fundamental aspects of the problem.Local nonlinearity is smeared into a small zone ahead of the crack. A “restrain stress” is introduced to also account for cracking at the meso-scale.The major emphasis is placed on developing a model that could exhibit the evolution characteristics of change in cracking behavior due to size and speed. Material inhomogeneity is assumed to favor self-similar crack growth although this may not always be the case. For relatively high restrain stress, the possible nucleation of micro-, meso- and macro-crack can be distinguished near the crack tip region. This distinction quickly disappears after a small distance after which scaling is no longer possible. This character prevails for Mode I and II cracking at different speeds. Special efforts are made to confine discussions within the framework of assumed conditions. To be kept in mind are the words of Isaac Newton in the Fourth Regula Philosophandi:

“Men are often led into error by the love of simplicity which disposes us to reduce things to few principles, and to conceive a greater simplicity in nature than there really is … We may learn something of the way in which nature operates from fact and observation; but if we conclude that it operates in such a manner, only because to our understanding that operates to be the best and simplest manner, we shall always go wrong.”––Isaac Newton

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.     

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.     

Dynamic interaction of various beams with the underlying soil – finite and infinite, half-space and Winkler models

Various beams lying on the elastic half-space and subjected to a harmonic load are analyzed by a double numerical integration in wavenumber domain. The compliances of the beam–soil systems are presented for a wide frequency range and for a number of realistic parameter sets. Generally, the soil stiffness G has a strong influence on the low-frequency beam compliance whereas the beam parameters EI and m^′ are more important for the high-frequency compliance. An important parameter is the elastic length l=(EI/G)^1/4 of the beam–soil system. Around the corresponding frequency ω_l=v_S/l, the wave velocity of the combined beam–soil system changes from the Rayleigh wave v_R≈v_S to the bending wave velocity v_B and the combined beam–soil wave has typically a strong damping. The interaction frequency ω_l is found not far from the characteristic frequency ω₀=(G/m^′)^1/2 where an amplification compared to the static compliance is observed for special parameter constellations. In contrast, real foundation beams show no resonance effects as they are highly damped by the radiation into the soil. At medium and high frequencies, asymptotes for the compliance of the beam–soil system are found, u/P(ρv_Paiω)^−3/4 in case of the dominating damping and u/P(−m^′ω²)^−3/4 for high frequencies. The low-frequency compliance of the coupled beam–soil system can be approximated by u/P1/Gl, but it also depends weakly on the width a of the foundation. All numerical results of different beam–soil systems are evaluated to yield a unique relation u/P₀=f(a/l). The integral transform method is also applied to ballasted and slab tracks of railway lines, showing the influence of train speed on the deformation of the track beam. The presented results of infinite beams on half-space are compared with results of finite beams and with infinite beams on a Winkler support. Approximating Winkler parameters are given for realistic foundation-soil systems which are useful when vehicle-track interaction is analyzed for the prediction of railway induced vibration.     

Explicit expressions of S(v), H(v) and L(v) for anisotropic elastic materials

The three matrices L(v), S(v) and H(v), appearing frequently in the investigations of the two-dimensional steady state motions of elastic solids, are expressed explicitly in terms of the elastic stiffness for general anisotropic materials. The special cases of monoclinic materials with a plane of symmetry at x₃ = 0, x₁ = 0, and x₂ = 0 are all deduced. Results for orthotropic materials appearing in the literature may be recovered from the present explicit expressions.     

ANALYSIS OF FINANCIAL DERIVATIVES BY MECHANICAL METHOD (Ⅱ)-BASIC EQUATION OF MARKET PRICE OF OPTION

The basic equation of market price of option is formulated by taking assumptions based on the characteristics of option and similar method for formulating basic equations in solid mechanics: hv ₀(t) = m ₁ v ₀ ^–1(t) – n ₁ v ₀(t) + F, where h, m ₁, n ₁, F are constants. The main assumptions are: the ups and downs of market price v ₀(t) are determined by supply and demand of the market; the factors, such as the strike price, tenor, volatility, etc. that affect on v ₀(t) are demonstrated by using proportion or inverse proportion relation; opposite rules are used for purchasing and selling respectively. The solutions of the basic equation under various conditions are found and are compared with the solution v _f (t) of the basic equation of market price of futures. Furthermore the one-one correspondence between v _f and v ₀(t) is proved by implicit function theorem, which forms the theoretic base for study of v _f affecting on the market price of option v ₀(t).     

Crack repair using an elastic filler

The effect of repairing a crack in an elastic body using an elastic filler is examined in terms of the stress intensity levels generated at the crack tip. The effect of the filler is to change the stress field singularity from order 1/r^1/2 to 1/r^(1-λ) where r is the distance from the crack tip, and λ is the solution to a simple transcendental equation. The singularity power (1-λ) varies from (the unfilled crack limit) to 1 (the fully repaired crack), depending primarily on the scaled shear modulus ratio γ_r defined by G₂/G₁=γ_rε, where 2πε is the (small) crack angle, and the indices (1, 2) refer to base and filler material properties, respectively. The fully repaired limit is effectively reached for γ_r≈10, so that fillers with surprisingly small shear modulus ratios can be effectively used to repair cracks. This fits in with observations in the mining industry, where materials with G₂/G₁ of the order of 10^-3 have been found to be effective for stabilizing the walls of tunnels. The results are also relevant for the repair of cracks in thin elastic sheets.     

New Explicit Expression of Barnett-Lothe Tensors for Anisotropic Linear Elastic Materials

The three Barnett-Lothe tensors H, L, S appear often in the Stroth formalism of two-dimensional deformations of anisotropic elastic materials [1–3]. They also appear in certain three-dimensional problems [4, 5]. The algebraic representation of H, L, S requires computation of the eigenvalues p^v(v=1,2,3) and the normalized eigenvectors (a, b). The integral representation of H, L, S circumvents the need for computing p ^v(v=1,2,3) and (a, b), but it is not simple to integrate the integrals except for special materials. Ting and Lee [6] have recently obtained an explicit expression of H for general anisotropic materials. We present here the remaining tensors L, S using the algebraic representation. They key to our success is the obtaining of the normalization factor for (a, b) in a simple form. The derivation of L and S then makes use of (a, b) but the final result does not require computation of (a, b), which makes the result attractive to numerical computation. Even though the tensor H given in [6] is in terms of the elastic stiffnesses C^μ v while the tensors L, S presented here are in terms of the reduced elastic compliances s′ _μv , the structure of L, S is similar to that of H. Following the derivation of H, we also present alternate expressions of L, S that remain valid for the degenerate cases p ₁ p ₂ and p₁=p₂ = p ₃. One may want to compute H, L, S using either C _μv or s′ _μv v, but not both. We show how an expression in C^μ v can be converted to an expression in s′ _μv v, and vice versa. As an application of the conversion, we present explicit expressions of the extic equation for p in C^μ v and s′ _μv v. This revised version was published online in August 2006 with corrections to the Cover Date.     

Two-phase geothermal flows with conduction and the connection with Buckley-Leverett theory

Two-phase flows of boiling water and steam in geothermal reservoirs satisfy a pair of conservation equations for mass and energy which can be combined to yield a hyperbolic wave equation for liquid saturation changes. Recent work has established that in the absence of conduction, the geothermal saturation equation is, under certain conditions, asymptotically identical with the Buckley-Leverett equation of oil recovery theory. Here we summarise this work and show that it may be extended to include conduction. In addition we show that the geothermal saturation wave speed is under all conditions formally identical with the Buckley-Leverett wave speed when the latter is written as the saturation derivative of a volumetric flow.Roman Letters C(P, S,q) geothermal saturation wave speed [ms^–1] (14) - c _t (P, S) two-phase compressibility [Pa^–1] (10) - D(P, S) diffusivity [m s^–2] (8) - E(P, S) energy density accumulation [J m^–3] (3) - g gravitational acceleration (positive downwards) [ms^–2] - h _w (P),h _w (P) specific enthalpies [J kg^–1] - J _M (P, S,P) mass flow [kg m^–2 s^–1] (5) - J _E (P, S,P) energy flow [J m^–2s^–1] (5) - k absolute permeability (constant) [m²] - k _w (S),k _s (S) relative permeabilities of liquid and vapour phases - K formation thermal conductivity (constant) [Wm^–1 K^–1] - L lower sheetC<0 in flow plane - m, c gradient and intercept - M(P, S) mass density accumulation [kg m^–3] (3) - O flow plane origin - P(x,t) pressure (primary dependent variable) [Pa] - q volume flow [ms^–1] (6) - S(x, t) liquid saturation (primary dependent variable) - S ^*(x,t) normalised saturation (Appendix) - t time (primary independent variable) [s] - T temperature (degrees Kelvin) [K] - T _sat(P) saturation line temperature [K] - TdT _sat/dP saturation line temperature derivative [K Pa^–1] (4) - T _c ,T _D convective and diffusive time constants [s] - u _w (P),u _s (P),u _r (P) specific internal energies [J kg^–1] - U upper sheetC > 0 in flow plane - U(x,t) shock velocity [m s^–1] - x spatial position (primary independent variable) [m] - X representative length - x, y flow plane coordinates - z depth variable (+z vertically downwards) [m] Greek Letters _P , _S remainder terms [Pa s^–1], [s^–1] - double-valued saturation region in the flow plane - h =h _s –h _w latent heat [J kg^–1] - = _w – _s density difference [kg m^–3] - line envelope - =D _K /D ₀ diffusivity ratio - porosity (constant) - _w (P), _s (P), _t (P, S) dynamic viscosities [Pa s] - v _w (P),v _s (P) kinematic viscosities [m²s^–1] - v ₀ =kh/KT kinematic viscosity constant [m² s^–1] - ₀ =v ₀ dynamic viscosity constant [m² s^–1] - _w (P), _s (P) density [kg m^–3] Suffixes r rock matrix - s steam (vapour) - w water (liquid) - t total - av average - 0 without conduction - K with conduction     

MOTIONS, FORCES AND MODE TRANSITIONS IN VORTEX-INDUCED VIBRATIONS AT LOW MASS-DAMPING

These experiments, involving the transverse oscillations of an elastically mounted rigid cylinder at very low mass and damping, have shown that there exist two distinct types of response in such systems, depending on whether one has a low combined mass-damping parameter (low m^*ζ), or a high mass-damping (highm^*ζ ). For our low m^*ζ, we find three modes of response, which are denoted as an initial amplitude branch, an upper branch and a lower branch. For the classical Feng-type response, at highm^*ζ , there exist only two response branches, namely the initial and lower branches. The peak amplitude of these vibrating systems is principally dependent on the mass-damping (m^*ζ), whereas the regime of synchronization (measured by the range of velocity U^*) is dependent primarily on the mass ratio, m^*ζ. At low (m^*ζ), the transition between initial and upper response branches involves a hysteresis, which contrasts with the intermittent switching of modes found, using the Hilbert transform, for the transition between upper–lower branches. A 180° jump in phase angle φ is found only when the flow jumps between the upper–lower branches of response. The good collapse of peak-amplitude data, over a wide range of mass ratios (m^*=1–20), when plotted against (m^*+C_A) ζ in the “Griffin” plot, demonstrates that the use of a combined parameter is valid down to at least (m^*+C_A)ζ 0·006. This is two orders of magnitude below the “limit” that had previously been stipulated in the literature, (m^*+C_A) ζ>0·4. Using the actual oscillating frequency (f) rather than the still-water natural frequency (f_N), to form a normalized velocity (U^*/f^*), also called “true” reduced velocity in recent studies, we find an excellent collapse of data for a set of response amplitude plots, over a wide range of mass ratiosm^* . Such a collapse of response plots cannot be predicted a priori, and appears to be the first time such a collapse of data sets has been made in free vibration. The response branches match very well the Williamson–Roshko (Williamson & Roshko 1988) map of vortex wake patterns from forced vibration studies. Visualization of the modes indicates that the initial branch is associated with the 2S mode of vortex formation, while the Lower branch corresponds with the 2P mode. Simultaneous measurements of lift and drag have been made with the displacement, and show a large amplification of maximum, mean and fluctuating forces on the body, which is not unexpected. It is possible to simply estimate the lift force and phase using the displacement amplitude and frequency. This approach is reasonable only for very low m^*.     

Fatigue and fretting fatigue of ion-nitrided 34CrNiMo6 steel

This work is concerned with the surface treatment (ion nitriding) of fretting fatigue and fatigue resistance of 34CrNiMo6. Tests are made on a servo-hydraulic machine under tension for both treated and non-treated specimens. The test parameters involve the applied displacements δ±80–±170 μm; fretting pressure σ_n=1000–1400 MPa; fatigue stress amplitude σ_a=380–680 MPa and stress ratio R=−1. The ion nitriding process improves both fatigue and fretting fatigue lives. Subsurface crack initiation from internal discontinuities was found for ion-nitrided specimens.     

Three-dimensional elastic behavior of a nonhomogeneous layer

Summary  This paper deals with the theoretical treatment of a three-dimensional elastic problem governed by a cylindrical coordinate system (r,θ,z) for a medium with nonhomogeneous material property. This property is defined by the relation G(z)=G ₀(1+z/a) ^m where G ₀,a and m are constants, i.e., shear modulus of elasticity G varies arbitrarily with the axial coordinate z by the power product form. We propose a fundamental equation system for such nonhomogeneous medium by using three kinds of displacement functions and, as an illustrative example, we apply them to an nonhomogeneous thick plate (layer) subjected to an arbitrarily distributed load (not necessarily axisymmetric) on its surfaces. Numerical calculations are carried out for several cases, taking into account the variation of the nonhomogeneous parameter m. The numerical results for displacement and stress components are shown graphically. Received 10 May 1999; accepted for publication 15 August 1999     

On Asymptotics of Solutions of Nonlinear Differential Equations of the Second Order

We study the problem of asymptotics of unbounded solutions of differential equations of the form y″ = α₀ p(t)ϕ(y), where α₀ ∈ {−1, 1}, p: [a, ω[→]0, +∞[, −∞ < a < ω ≤ +∞, is a continuous function, and ϕ: [y ₀, +∞[→]0, +∞[ is a twice continuously differentiable function close to a power function in a certain sense.__________Translated from Neliniini Kolyvannya, Vol. 8, No. 1, pp. 18–28, January–March, 2005.     

Symmetry and secondary bifurcation of elastic structures

We consider non-linear bifurcation problems for elastic structures modeled by the operator equation F[w;α]=0 where F:X×R^k→Y,X,Y are Banach spaces and XY. We focus attention on problems whose bifurcation equations are of the form

f_i(α₁,α₂;λ,μ)=(a_iμ+b_iλ)α_i+p_iα_i³+q_iα_i∑_j=1,j≠i^kα_j+1²+α_ih_i(λ,μ;α₁,α₂,…α_k) i=1,2,…k

which emanates from bifurcation problems for which the linearization of F is Fredholm operators of index 0. Under the assumption of F being odd we prove an important theorem of existence of secondary bifurcation. Under this same assumption we prove a symmetry condition for the reduced equations and consequently we got an existence result for secondary bifurcation. We also include a stability analysis of the bifurcating solutions.     

Induced parallel vortex shedding from a circular cylinder at Re of O(10) by using the cylinder end suction technique

In the present work, the objective is to attempt to induce parallel vortex shedding at a moderately high Reynolds number (=1.578 × 10⁴) by using the cylinder end suction method, and measure the associated aerodynamic parameters.We first measured the aerodynamic parameters of a single circular cylinder without end suction, and showed that the quantities measured are in good agreement with equivalent data in the published literature. Next, by using different amount of end suction which resulted in increasing the cylinder end velocity by 1%, 2% and 2.5%, we were able to show that the above corresponded to the situation of under suction, optimal suction and over suction, respectively. With optimal suction, we demonstrated that the end suction method works at Re = 1.578 × 10⁴. The shape of the primary vortex shed became straighter than when there is no end suction, and parameters like cylinder surface pressure distribution, drag force per unit span, as well as vortex shedding frequency all showed negligible spanwise variation. Further careful analyses showed that when compared to the naturally existing curved vortex shedding, with parallel vortex shedding the mid-span drag per unit span became slightly smaller, but the drag averaged over the cylinder span became slightly larger. For cylinder surface pressure, it was found that cylinder end effects mainly influenced the surface pressure in the angular ranges −180°  β < −60° and 60° < β  180°. Without end suction, the cylinder surface pressure in the above ranges was found to increase (become less negative) slightly with |z/d|, but such increase disappeared when optimal end suction was applied, and the cylinder surface pressure distribution became spanwise location independent. As for the vortex shedding frequency (Strouhal number), although the Strouhal number showed spanwise variation when there is no end suction and negligible spanwise variation when optimal suction was applied, the difference between the spanwise averaged Strouhal number was quite negligible. With under suction, the spanwise dependence of various aerodynamic parameters existed, but was found to be not as significant as when no end suction was applied at all. With over suction, the flow situation was found to be practically no change from the optimal suction situation.     

On the number of droplets in aerosols

The number of droplets which may be formed with a supersaturated vapor in presence of a gas cannot exceed a number proportional to (p_v−p_v0)⁴ where p_v and p_v0 denote at the same temperature the pressure of the supersaturated vapor–gas mixture and the pressure of the saturated vapor–gas mixture. The energy necessary to the droplet formation is also bounded by a number proportional to (p_v−p_v0)².