The extension of Beltrami's ellipsoid to anisotropic bodies

Summary The spectral decomposition of the compliance, stiffness, and failure tensors for transversely isotropic materials was studied and their characteristic values were calculated using the components of these fourth-rank tensors in a Cartesian frame defining the principal material directions. The spectrally decomposed compliance and stiffness or failure tensors for a transversely isotropic body (fiber-reinforced composite), and the eigenvalues derived from them define in a simple and efficient way the respective elastic eigenstates of the loading of the material. It has been shown that, for the general orthotropic or transversely isotropic body, these eigenstates consist of two double components, ₁ and ₂ which are shears (₂ being a simple shear and ₁, a superposition of simple and pure shears), and that they are associated with distortional components of energy. The remaining two eigenstates, with stress components ₃, and ₄, are the orthogonal supplements to the shear subspace of ₁ and ₂ and consist of an equilateral stress in the plane of isotropy, on which is superimposed a prescribed tension or compression along the symmetry axis of the material. The relationship between these superimposed loading modes is governed by another eigenquantity, the eigenangle .The spectral type of decomposition of the elastic stiffness or compliance tensors in elementary fourth-rank tensors thus serves as a means for the energy-orthogonal decomposition of the energy function. The advantage of this type of decomposition is that the elementary idempotent tensors to which the fourth-rank tensors are decomposed have the interesting property of defining energy-orthogonal stress states. That is, the stress-idempotent tensors are mutually orthogonal and at the same time collinear with their respective strain tensors, and therefore correspond to energy-orthogonal stress states, which are therefore independent of each other. Since the failure tensor is the limiting case for the respective _x, which are eigenstates of the compliance tensor S, this tensor also possesses the same remarkable property.An interesting geometric interpretation arises for the energy-orthogonal stress states if we consider the projections of _x in the principal₃D stress space. Then, the characteristic state ₂ vanishes, whereas stress states ₁, ₃ and ₄ are represented by three mutually orthogonal vectors, oriented as follows: The ₃ and ₄ lie on the principal diagonal plane (₃₁₂) with subtending angles equaling (–/2) and (-), respectively. On the positive principal ₃-axis, is the eigenangle of the orthotropic material, whereas the ₁-vector is normal to the (₃₁₂)-plane and lies on the deviatoric -plane. Vector ₂ is equal to zero.It was additionally conclusively proved that the four eigenvalues of the compliance, stiffness, and failure tensors for a transversely isotropic body, together with value of the eigenangle , constitute the five necessary and simplest parameters with which invariantly to describe either the elastic or the failure behavior of the body. The expressions for the _x-vector thus established represent an ellipsoid centered at the origin of the Cartesian frame, whose principal axes are the directions of the ₁-, ₃- and ₄-vectors. This ellipsoid is a generalization of the Beltrami ellipsoid for isotropic materials.Furthermore, in combination with extensive experimental evidence, this theory indicates that the eigenangle alone monoparametrically characterizes the degree of anisotropy for each transversely isotropic material. Thus, while the angle for isotropic materials is always equal to _i = 125.26° and constitutes a minimum, the angle || progressively increases within the interval 90–180° as the anisotropy of the material is increased. The anisotropy of the various materials, exemplified by their ratiosE _L/2G_L of the longitudinal elastic modulus to the double of the longitudinal shear modulus, increases rapidly tending asymptotically to very high values as the angle approaches its limits of 90 or 180°.     

Propagation of a spherical wave in nonlinearly compressible and elastoplastic materials

The problem of spherical wave propagation in soil under the action of an intense uniformly decreasing load ₀(t) applied to the boundary of a cavity with radius r₀ is considered. Soil with a high stress level is modeled either by ideally nonlinearly compressible or elastoplastic material, taking account of linear irreversible unloading for the material. In contrast to [1–7], in order to describe material movement use is made of strain theory [8] with determining functions = (), _i=_i(_i), where , _i, , _i are the first and second invariants of strain and stress tensors. During material loading these functions are presented in the form of polynomials ()=(_i+₂¦¦), _ii)=(_i-_2i)_i, in which constant coefficients _i, _i=1, 2) are determined by experiment, taking account of the triaxial stressed state of soil. Solution of the problem is constructed by an analytically reversible method, with prescribed shape for the shock-wave (SW) surface in the form of a second-degree polynomial relating to time t and a numerical method of characteristics for a prescribed arbitrarily decreasing load _i(t). On the basis of the analytical equations obtained, calculations are carried out for material parameters (including loading profile) in a computer and stresses and mass velocity of plastic and elastoplastic materials are compared.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 4, pp. 95–100, July–August, 1986.The authors express their sincere thanks to Kh. A. Rakhmatulin for discussing the results of this work.     

Stress analysis of the cylinder subjected to arbitrarily distributed axisymmetrical loads

Summary Stress analysis has been carried out for a finite cylinder subjected to arbitrarily distributed axisymmetrical surface loads. Direct stress _x in the axial direction is assumed to be of the form _x = ₀+r ₁ +r ₂ where ₀ to ₂ are functions of x. Using the equations of equilibrium and compatibility the other direct stresses and the shearing stress are expressed by ₁ and ₂. Fundamental equations governing ₁ and ₂ are introduced using the variational principle of complementary energy. From the results of the present analysis it is evident that the boundary conditions can be satisfied completely even for the case where the external forces are specified in complicated form, and that more accurate solutions can easily be obtained by introducing additional terms in _x.

---

Spannungsanalyse für den Zylinder unter axialsymmetrischer Last in beliebiger Verteilung

Übersicht Für einen endlichen Zylinder unter axialsymmetrischer Oberflächenlast in beliebiger Verteilung werden die Spannungen ermittelt. Die Normalspannung in Axialrichtung wird in der Form _x = ₀+r ₁ +r ₂ angesetzt mit ₀, ₁, ₂ als Funktionen von x. Mit Hilfe der Gleichgewichtsund Verträglichkeitsbedingungen werden die anderen Normalspannungen und die Schubspannung durch ₁ und ₂ ausgedrückt. Über das Variationsprinzip für die Komplementärenergie werden die grundlegenden Gleichungen für ₁ und ₂ eingeführt. Die Ergebnisse zeigen, daß die Randbedingungen selbst für komplizierte Belastungsarten vollständig erfüllbar sind und mit zusätzlichen Termen in _x mühelos noch genauere Lösungen bestimmt werden können.

     

Theoretical derivation of tangential velocity profiles in a flat vortex chamber-influence of turbulence and wall friction

Summary As part of a study on the hydrodynamics of a cyclone separator, a theoretical investigation of the flow pattern in a flat box cyclone (vortex chamber) has been carried out. Expressions have been derived for the tangential velocity profile as influenced by internal friction (eddy viscosity) and wall friction. The most important parameter controlling the tangential velocity profile is = –u ₀ R/(v+ ), where u ₀ is the radial velocity at the outer radius R of the cyclone, the kinematic liquid viscosity and is the kinematic eddy viscosity. For values of greater than about 10 the tangential velocity profile is nearly hyperbolic, for smaller than 1 the tangential velocity even decreases towards the centre. It is shown how and also the wall friction coefficient may be obtained from experimental velocity profiles with the aid of suitable graphs. Because of the close relation between eddy viscosity and eddy diffusion, measurements of velocity profiles in flat box cyclones will also provide information on the eddy motion of particles in a cyclone, a motion reducing its separation efficiency.List of symbols A cross-sectional area of cyclone inlet - h height of cyclone - p static pressure in cyclone - p static pressure difference in cyclone between two points on different radius - r radius in cyclone - r ₁ radius of cyclone outlet - R radius of cyclone circumference - u radial velocity in cyclone - u ₀ radial velocity at circumference of flat box cyclone - v tangential velocity - v ₀ tangential velocity at circumference of flat box cyclone - w axial velocity - z axial co-ordinate in cyclone - friction coefficient in flat box cyclone (for definition see § 5) - ₁ value of friction coefficient for ₁<< ₂ - ₂ value of friction coefficient for ₂<<1 - = - ₁ value of for ₁<< ₂ - ₂ value of for ₂<<1 - thickness of laminar boundary layer - =/h - turbulent kinematic viscosity - ratio of z to h - k ratio of height of cyclone to radius R of cyclone - parameter describing velocity profile in cyclone =–u ₀ R/(+) - kinematic viscosity of fluid - density of fluid - ratio of r to R - ₁ value of at outlet of cyclone - ₂ value of at inner radius of cyclone inlet - _w shear stress at cyclone wall - angular momentum in cyclone/angular momentum in cyclone inlet - ₁ value of at = ₁ - ₂ value of at = ₂     

Flow in porous media III: Deformable media

Stokes flow in a deformable medium is considered in terms of an isotropic, linearly elastic solid matrix. The analysis is restricted to steady forms of the momentum equations and small deformation of the solid phase. Darcy's law can be used to determine the motion of the fluid phase; however, the determination of the Darcy's law permeability tensor represents part of the closure problem in which the position of the fluid-solid interface must be determined.Roman Letters A interfacial area of the- interface contained within the macroscopic system, m² - A interfacial area of the- interface contained within the averaging volume, m² - A _e area of entrances and exits for the-phase contained within the macroscopic system, m² - A ^* interfacial area of the- interface contained within a unit cell, m² - A _e ^* area of entrances and exits for the-phase contained within a unit cell, m² - E Young's modulus for the-phase, N/m² - e i unit base vectors (i = 1, 2, 3) - g gravity vector, m²/s - H height of elastic, porous bed, m - k unit base vector (=e ₃) - characteristic length scale for the-phase, m - L characteristic length scale for volume-averaged quantities, m - n unit normal vector pointing from the-phase toward the-phase (n = -n ) - p pressure in the-phase, N/m² - P p – g·r, N/m² - r ₀ radius of the averaging volume, m - r position vector, m - t time, s - T total stress tensor in the-phase, N/m² - T ⁰ hydrostatic stress tensor for the-phase, N/m² - u displacement vector for the-phase, m - V averaging volume, m₃ - V volume of the-phase contained within the averaging volume, m³ - v velocity vector for the-phase, m/s Greek Letters V /V, volume fraction of the-phase - mass density of the-phase, kg/m³ - shear coefficient of viscosity for the-phase, Nt/m² - first Lamé coefficient for the-phase, N/m² - second Lamé coefficient for the-phase, N/m² - bulk coefficient of viscosity for the-phase, Nt/m² - T –T ⁰ , a deviatoric stress tensor for the-phase, N/m²     

Evaluation of the dynamic-mechanical properties of solids from a cooperative model for stress relaxation

For many solid materials the stress relaxation process obeys the universal relationF = – (d/d lnt)_max = (0.1 ± 0.01) ( ₀ — _i ), regardless of the structure of the material. Here denotes the stress,t the time, ₀ the initial stress of the experiment and _i the internal stress. A cooperative model accounting for the similarity in relaxation behaviour between different materials was developed earlier. Since this model has a spectral character, the concepts of linear viscoelasticity are used here to evaluate the corresponding prediction of the dynamic mechanical properties, i.e. the frequency dependence of the storageE () and lossE () moduli. Useful numerical approximations ofE () andE () are also evaluated. It is noted that the universal relation in stress relaxation had a counterpart in the frequency dependence ofE (). The theoretical prediction of the loss factor for high-density polyethylene is compared with experimental results. The agreement is good.     

Pressure pulsations during interaction between a subsonic jet and the subsonic section of a supersonic jeton a flat obstacle

Pressure pulsations were measured during in-leakage of a subsonic jet and the subsonic section of a heated supersonic jet on a flat obstacle. Data have been obtained on the total and spectrum levels of the pressure pulsations at different spacings X of the obstacle from the nozzle exit. It is shown that when the obstacle is disposed at the section of the jet where the local velocity is subsonic, the pulsation levels outside the dependence on the conditions at the nozzle exit (Mach number Maxa 0 a 3.0; stagnation temperature T₀=280–1200K) vary in direct proportion to the local velocity head q. The ratio between the total level and q is (/g)=0.2–0.3. It is established that for a subsonic velocity ahead of the obstacle, all the spectra obtained for different values of M _a , T₀, d _a and X in the coordinates Sh=f(d/V) and (₁*/q)(V/d) will lie on a single generalized spectrum. Here ₁* is the pulsation level in a 1-Hz band, and d and V are, respectively, the jet diameter and velocity directly in front of the obstacle.Translated from Izvestiya Akademiya Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 172–174, September–October, 1975.     

Transport in ordered and disordered porous media V: Geometrical results for two-dimensional systems

In this paper we continue the geometrical studies of computer generated two-phase systems that were presented in Part IV. In order to reduce the computational time associated with the previous three-dimensional studies, the calculations presented in this work are restricted to two dimensions. This allows us to explore more thoroughly the influence of the size of the averaging volume and to learn something about the use of anon-representative region in the determination of averaged quantities.

Nomenclature

Roman Letters A interfacial area of the– interface associated with the local closure problem, m² - a _i i=1, 2, gaussian probability distribution used to locate the position of particles - l unit tensor - characteristic length for the-phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1,2,3 lattice vectors, m - m convolution product weighting function - m _V special convolution product weighting function associated with a unit cell - n _i i=1, 2 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume,V, m³ - x position of the centroid of an averaging area, m - x ₀ reference position of the centroid of an averaging area, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters V /V, volume average porosity - _{a i} standard deviation ofa _i - _r standard deviation ofr - intrinsic phase average of      

Zur Variationsformulierung nichtlinearer Randwertprobleme

Übersicht Es werden verschiedene Bedingungen aufgestellt, die es erlauben, die durch die beiden (Systeme von) nichtlinearen DifferentialgleichungenA (u, ) = q, B (u, ) = und Randbedingungen zusammen mit den nichtlinearen algebraischen Relationenq = C(u, ), = D(u, ) beschriebene Aufgabe durch äquivalente Variationsprobleme zu ersetzen. Dabei zeigt sich ein enger Zusammenhang mit den in der Festkörpermechanik wohlbekannten Prinzipien der virtuellen Verschiebungen und der virtuellen Kräfte. Die auf systematischem Weg konstruierten Variationsfunktionale enthalten viele in der Physik bekannte Funktionale als Sonderfälle, insbesondere jene, die in der Elastomechanik nach Green, Castigliano, Hellinger, Reißner, Hu und Washizu benannt werden.

---

Summary In this paper there are established various conditions which allow a variational formulation of the problem described by the two (systems of) nonlinear differential equationsA(u, ) = q, B(u, ) = and boundary conditions together with the nonlinear algebraic relationsq = C(u, ), = D(u, ). Besides a close relationship is revealed to the principles of virtual displacements and virtual forces which are wellknown in solid mechanics. The systematically constructed variational functional contain many functionals in physics as special cases, mainly those of Green, Castigliano, Hellinger, Reißner, Hu and Washizu in elastomechanics.

     

Die Temperaturabhängigkeit der Oberflächenspannung von reinen Kältemitteln vom Tripelpunkt bis zum kritischen Punkt

Zusammenfassung Die Oberflächenspannung von sechs reinen Substanzen — SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃ und CHClF₂ — wurde mit Hilfe einer modifizierten Kapillarmethode gemessen. Die zur Berechnung der Oberflächenspannung erforderlichen Sättigungsdichten und wurden teils aus vorhandenen Zustandsgleichungen, teils aus ebenfalls gemessenen Brechungsindizes bestimmt. Die Temperaturabhängigkeit der Oberflächenspannung läßt sich durch einen erweiterten Ansatz nach van der Waals =_O (T_c-T)(1+...) darstellen, wobei bei einfachen Stoffen ein eingliedriger, bei polaren und assoziierenden Stoffen ein zweigliedriger Ansatz notwendig und ausreichend ist. Für den kritischen Exponenten der Oberflächenspannung wurde ein von der molekularen Substanz weitgehend unabhängiger Wert von =1.284±0.005 gefunden.

---

Temperature dependence of surface tension of pure refrigerants from triple point up to the critical point

The surface tension of six fluids (SF₆, CCl₃F, CCl₂F₂, CClF₃, CBrF₃, CHClF₂) have been measured by means of a modified capillary rise method. The liquid vapor densities, which are needed to calculate the surface tension, have partly been determined by means of refractive indices simultaneously measured in the same apparatus. The temperature dependence of the surface tension is described by an extended van der Waals power law =_O(T_c-T)(1+...). For simple fluids one term and for polar and associating fluids two terms are necessary and sufficient. The critical exponent is found to be 1.284 ± 0.005 and nearly independent of the molecular structure.

---

Formelzeichen a² Laplace-Koeffizient - a Parameter - B_O, B_on Koeffizient der Koexistenzkurve - g Erdbeschleunigung - H Höhe, kapillare Steighöhe - LL Lorentz-Lorenz-Funktion oder Refraktionskonstante - M molare Masse - M Zahl der Meßwerte - N Zahl der unbekannten Parameter - n Brechungsindex - p Druck - R,r Radius - s Entropie - SD Standardabweichung - T, t Temperatur - u innere Energie Griechische Formelzeichen Exponent des Laplace-Koeffizienten - Exponent der Koexistenzkurve - 2. Exponent der Oberflächenspannung - Wellenlänge des Lichts - Exponent der Oberflächenspannung - _D Dipolmoment - , Dichte der Flüssigkeit bzw. des Dampfes - Oberflächenspannung - reduzierte Temperatur (1-T/T_c) - ² gewichtete Varianz Indizes c kritischer Zustand - D Differenz - m Mittelwert - T Isotherme - t Zustand am Tripelpunkt - S Zustand am Schmelzpunkt - bezogen auf Oberfläche     

On Perturbative Expansions to the Stochastic Flow Problem

When analyzing stochastic steady flow, the hydraulic conductivity naturally appears logarithmically. Often the log conductivity is represented as the sum of an average plus a stochastic fluctuation. To make the problem tractable, the log conductivity fluctuation, f, about the mean log conductivity, lnK _G, is assumed to have finite variance, _f ². Historically, perturbation schemes have involved the assumption that _f ²<1. Here it is shown that _f may not be the most judicious choice of perturbation parameters for steady flow. Instead, we posit that the variance of the gradient of the conductivity fluctuation, _f ², is more appropriate hoice. By solving the problem withthis parameter and studying the solution, this conjecture can be refined and an even more appropriate perturbation parameter, , defined. Since the processes f and f can often be considered independent, further assumptions on f are necessary. In particular, when the two point correlation function for the conductivity is assumed to be exponential or Gaussian, it is possible to estimate the magnitude of f in terms of _f and various length scales. The ratio of the integral scale in the main direction of flow ( _x ) to the total domain length (L^*), _x ²=_x/L^*, plays an important role in the convergence of the perturbation scheme. For _x smaller than a critical value _c, _x < _c, the scheme's perturbation parameter is =_f/_x for one- dimensional flow, and =_f/_x ² for two-dimensional flow with mean flow in the x direction. For _x > _c, the parameter =_f/_x ³ may be thought as the perturbation parameter for two-dimensional flow. The shape of the log conductivity fluctuation two point correlation function, and boundary conditions influence the convergence of the perturbation scheme.     

Effect of finite conductivity in a medium on shock-wave interaction with a magnetic field

We consider the problem of shock-wave incidence on a magnetic wall, which has been studied in [1]. It is shown that the dynamics of the processes which take place in this case depend significantly on the behavior of the conductivity-temperature dependence (T) of the medium and also on the magnitude of the magnetic-field intensity H₀.An exact solution of the problem is constructed for a special form of the law (T). For an arbitrary law (T) the problem is studied numerically by means of digital computer computations; the results are compared with the exact results.Analysis of these solutions shows that the dissipative properties of the medium (electrical conductivity, viscosity), which determine the structure of the refracted wave front, affect the nature of the entire flow as a whole.The formulated problem also makes it possible to clarify the characteristic features of the decay of a discontinuity in a conducting medium.The authors wish to thank A. A. Samarskii, L. A. Zaklyaz'minskii, L. M. Degtyarev, and A. P. Favorskii for discussions of the study, D. A. Gol'dina and A. A. Ivanov for carrying out the numerical calculations, and also G. A. Lyubimov for several helpful comments.     

Flow and fracture of rocks under general triaxial compression

Recent laboratory studies of the flow and fracture of rocks under general triaxial compression are reviewed. New developments in laboratory techniques have made it possible to measure three principal stresses and strains under general triaxial stress states, in which all three principal stresses are different.Strength and ductility of isotropic rocks are markedly affected not only by the least compression ₃, but also by the intermediate compression ₂, although these two effects are rather additional in strength, but opposite in ductility. The experimental results show that dilatancy is highly anisotropic under the general triaxial stress states.Deformational properties of anisotropic rocks have been also measured under the general triaxial compression. In this case, the effect of the intermediate compression markedly depends on the orientations of the weak planes.     

The determination of a general time creep compliance relation of linear viscoelastic materials under constant load and its extension to nonlinear viscoelastic behavior for the Burger model

Linear viscoelastic materials yield a creep function which only depends on time if creep experiments are performed under constant stress ₀. In practice, this condition is very difficult to realize, and as a consequence, the experiments are performed under constant force. For small strains the difference between the conditions of constant stress and constant force is negligible. Otherwise, the decrease in cross-section has to be taken into account and leads to increasing stress in the course of time for creep experiments under constant load. The Boltzmann superposition principle is solved under the condition of constant load and for strains . The creep complicance C(t; ₀) defined by the ratio becomes, in principle, dependent on the initial stress ₀. As a consequence, a set of creep compliance curves cannot be approximated with a simple parameter fit. Already the application of the solution on the Burger model yields a creep compliance curve with all three creep ranges. Furthermore, the mathematical structure of the time creep compliance relation of the Burger model allows nonlinear viscoelastic extension via the introduction of the yield strength _max and a nonlinearity parameter n _l . The creep behavior of PBT and PC can be described in the range of long times up to initial stresses ₀, being 75% for PBT and 60% for PC of the yield stress _max with only two or one free fit parameter, respectively.     

Analytical evaluation of stress-strain test data

An analytical model for deducing the actual stress-strain properties from laboratory test results is discussed. As an illustration, an elastic bilinear material is used for unconfined cylindrical compression test conditions, as simulated with a finite element analysis. The results obtained are applicable for assisting in evaluating measured strength and stiffness properties of some clay soils, concrete test cylinders, concrete cores, and rock cores.The quantitative results of this study can be used for interpreting measured stress-strain data for unconfined compression test conditions. The error in measured results is shown to be influenced by Poisson's ratio, length-to-diameter ratio of the specimen, end condition, and ratio of inelastic modulus to initial elastic modulus. Curves for adjusting the measured results to the theoretical results are presented.Nomenclature D specimen diameter - E _i initial elastic stiffness modulus - E _y elastic stiffness modulus beyond the yield stress, plastic or inelastic modulus - L specimen length - axial strain - _av average strain - _g gage length strain - _y yield strain - Poisson's ratio - compressive stress - _av average stress - _t theoretical compressive stress - _y yield stress - _ym measured stress at the yield strain     

Transport in ordered and disordered porous media II: Generalized volume averaging

In this paper we develop the averaged form of the Stokes equations in terms of weighting functions. The analysis clearly indicates at what point one must choose a media-specific weighting function in order to achieve spatially smoothed transport equations. The form of the weighting function that produces the cellular average is derived, and some important geometrical theorems are presented.Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - A _p surface area of a particle, m² - d _p 6V _p/A_p, effective particle diameter, m - g gravity vector, m/s² - I unit tensor - K _m permeability tensor for the weighted average form of Darcy's law, m² - L general characteristic length for volume averaged quantities, m - L _p general characteristic length for volume averaged pressure, m - L characteristic length for the porosity, m - L _v characteristic length for the volume averaged velocity, m - l characteristic length (pore scale) for the-phase - l _i i=1, 2, 3 lattice vectors, m - (y) weighting function - m(–y) (y), convolution product weighting function - _v special weighting function associated with the traditional averaging volume - m _v special convolution product weighting function associated with the traditional averaging volume - m _g general convolution product weighting function - m _V unit cell convolution product weighting function - m _C special convolution product weighting function for ordered media which produces the cellular average - m _D special convolution product weighting function for disordered media - m _M master convolution product weighting function for ordered and disordered media - n unit normal vector pointing from the-phase toward the-phase - p pressure in the-phase, N/m² - p_m superficial weighted average pressure, N/m² - p _m intrinsic weighted average pressure, N/m² - p traditional intrinsic volume averaged pressure, N/m² - p p p _m , spatial deviation pressure, N/m² - r ₀ radius of a spherical averaging volume, m - r _m support of the convolution product weighting function, m - r position vector, m - r position vector locating points in the-phase, m - V averaging volume, m³ - V volume of the-phase contained in the averaging volume, m³ - V _cell volume of a unit cell, m³ - V velocity vector in the-phase, m/s - vm superficial weighted average velocity, m/s - v _m intrinsic weighted average velocity, m/s - V volume of the-phase contained in the averaging volume, m³ - V _p volume of a particle, m³ - v traditional superficial volume averaged velocity, m/s - v v p _m spatial deviation velocity, m/s - x position vector locating the centroid of the averaging volume or the convolution product weighting function, m - y position vector relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - _m m * . weighted average porosity - mass density of the-phase, kg/m³ - viscosity of the-phase, Ns/m² - V /V, volume fraction of the-phase     

Transport in ordered and disordered porous media IV: Computer generated porous media for three-dimensional systems

In the method of volume averaging, the difference between ordered and disordered porous media appears at two distinct points in the analysis, i.e. in the process of spatial smoothing and in the closure problem. In theclosure problem, the use of spatially periodic boundary conditions isconsistent with ordered porous media and the fields under consideration when the length-scale constraint,r ₀L is satisfied. For disordered porous media, spatially periodic boundary conditions are an approximation in need of further study.In theprocess of spatial smoothing, average quantities must be removed from area and volume integrals in order to extractlocal transport equations fromnonlocal equations. This leads to a series of geometrical integrals that need to be evaluated. In Part II we indicated that these integrals were constants for ordered porous media provided that the weighting function used in the averaging process contained thecellular average. We also indicated that these integrals were constrained by certain order of magnitude estimates for disordered porous media. In this paper we verify these characteristics of the geometrical integrals, and we examine their values for pseudo-periodic and uniformly random systems through the use of computer generated porous media.

Nomenclature

Roman Letters A interfacial area of the- interface associated with the local closure problem, m² - A _e area of entrances and exits for the-phase contained within the averaging system, m² - a _i i=1, 2, 3 gaussian probability distribution used to locate the position of particles - I unit tensor - L general characteristic length for volume averaged quantities, m - L characteristic length for , m - L characteristic length for , m - characteristic length for the -phase particles, m - ⁰ reference characteristic length for the-phase particles, m - characteristic length for the-phase, m - _i i=1, 2, 3 lattice vectors, m - m convolution product weighting function - m _v special convolution product weighting function associated with the traditional volume average - n _i i=1, 2, 3 integers used to locate the position of particles - n unit normal vector pointing from the-phase toward the-phase - n _e outwardly directed unit normal vector at the entrances and exits of the-phase - r _p position vector locating the centroid of a particle, m - r gaussian probability distribution used to determine the size of a particle, m - r ₀ characteristic length of an averaging region, m - r position vector, m - r _m support of the weighting functionm, m - averaging volume, m³ - V volume of the-phase contained in the averaging volume,, m³ - x positional vector locating the centroid of an averaging volume, m - x ₀ reference position vector associated with the centroid of an averaging volume, m - y position vector locating points relative to the centroid, m - y position vector locating points in the-phase relative to the centroid, m Greek Letters indicator function for the-phase - Dirac distribution associated with the- interface - V /V, volume average porosity - /L, small parameter in the method of spatial homogenization - standard deviation ofa _i - _r standard deviation ofr - _r intrinsic phase average of      

Viscous properties of soil

Dynamic problems connected with the wave propagation in soils not saturated with water and with wave interaction with obstacles and structural elements at the present time are solved on the basis of models in which plastic but not viscous soil properties are taken into account [1–5]. An analysis of experimental data and their comparison with the calculated results [4, 5] confirms that it is permissible to apply the model of an elasticplastic medium to soils in problems concerning the interaction of waves and structures. At the same time plane-wave damping in soils takes place more intensively than would follow from calculations carried out on the basis of models of an elastic-plastic medium. For example, if in a section of a poured sandy soil, taken as the initial section, the maximum stress in the wave is _m=ll kgf/cm² and its duration is 6=8 msec, then at a distance of 25 cm the calculations give _m=9.5 kgf/cm², while the experiment gives _m= 5 kgf/cm². If in the initial section _m= 20 kgf/cm² and =6 msec, then at a distance of 35 cm the calculation gives _m= l7 kgf/cm², while the experiment gives _m= 9 kgf/cm². In the calculations it was assumed that unloading takes place with a constant strain. This deviation of the calculated results from the experiment can be explained, in the first place, by the dependence of the () on the strain rate , which is not taken into account in the model of an elastic-plastic medium. The viscous properties cause additional energy losses and a more intensive damping of the waves. Experimentally the dependence of the () curves on the strain rate has been investigated for many soils [5–8]. The dynamic load on the test sample was produced by a body falling from a height or being accelerated by some method. Below we present test results of viscous soil properties when the test sample is compressed by an air shock wave. Compression curves and approximate numerical values of the coefficient of viscosity are obtained.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 9, No. 4, pp. 68–71, July–August, 1968.The author thanks A. I. Shishikin for his participation in the experiments.     

Relation between viscoelasticity and shear-thinning behaviour in liquids

Summary The shear-thinning behaviour of a liquid is represented in terms of a relaxation time, defined by the ratio ₀/G₀ of initial viscous and elastic constants. The relationship provides a very simple basis for the evaluation of andG ₀ from viscosity/shear data. Results are compared with relaxation times and moduli from primary normal-stress measurement, from stress relaxation and from direct measurement of recoverable shear strain. Good agreement is found but there is experimental evidence the recoverable shear strain _e is related to normal stressN ₁ and shear stress by _e = N₁/3, which does not agree with the theoretical prediction of eitherWeissenberg orLodge.

---

Zusammenfassung Das Scherentzähungsverhalten einer Flüssigkeit wird mittels einer Relaxationszeit beschrieben, die durch das Verhältnis der Anfangswerte von Viskosität und Elastizitätsmodul ₀/G₀ definiert ist. Diese Beziehung eröffnet eine einfache Methode zur Bestimmung von undG ₀ aus Scherviskositätsmessungen. Die damit erhaltenen Ergebnisse werden mit Relaxationszeiten und Moduln verglichen, die durch Messung der ersten Normalspannungsdifferenz, der Spannungsrelaxation und der Scherdehnungsrückstellung (recoverable shear strain) gewonnen worden sind. Es wird eine gute Übereinstimmung gefunden, zugleich aber wird der experimentelle Nachweis geführt, daß die Scherdehnungsrückstellung _e mit der ersten NormalspannungsdifferenzN ₁ und der Schubspannung durch die Beziehung _e = N₁/3 verknüpft ist, was sowohl zu der theoretischen Voraussage vonWeissenberg als auch zu derjenigen vonLodge im Widerspruch steht.

---

---

With 10 figures and 1 table     

Three-dimensional effects in beams: Isodyne assessment of a plane solution

The stress distribution in a homogeneous beam subjected to three-point bending is investigated using the method of optical isodynes. The three stress components _xx,_yy and _xy acting in the planes formed by the longitudinal and vertical axes of the beam are determined in three planes situated at different through the thickness locations with respect to the beam's midplane. The experimental results are subsequently correlated with the two-dimensional elasticity solution. It is illustrated that at locations sufficiently removed from the centrally applied concentrated load, good correlation between theory and experiment is obtained. In the regions where high stress gradients exist however, differences are observed in the in-plane stress distributions in the different planes. These differences are explained by the presence of the out of plane normal stress _zz using the relations of optical isodynes. Greatest differences between theory and experiment are obtained for the in-plane shear stress component _xy.