A rotating hot-wire technique for spatial sampling of disturbed and manipulated duct flows

The documentation and control of flow disturbances downstream of various open inlet contractions was the primary focus with which to evaluate a spatial sampling technique. An X-wire probe was rotated about the center of a cylindrical test section at a radius equal to one-half that of the test section. This provided quasi-instantaneous multi-point measurements of the streamwise and azimuthal components of the velocity to investigate the temporal and spatial characteristics of the flowfield downstream of various contractions. The extent to which a particular contraction is effective in controlling ingested flow disturbances was investigated by artificially introducing disturbances upstream of the contractions. Spatial as well as temporal mappings of various quantities are presented for the streamwise and azimuthal components of the velocity. It was found that the control of upstream disturbances is highly dependent on the inlet contraction; for example, reduction of blade passing frequency noise in the ground testing of jet engines should be achieved with the proper choice of inlet configurations.List of symbols K _uv correlation coefficient= - P percentage of time that an azimuthal fluctuating velocity derivative dv/d is found - U streamwise velocity component U=U (, t) - V azimuthal or tangential velocity component due to flow and probe rotation V=V (, t) - mean value of streamwise velocity component - U _m resultant velocity from and - mean value of azimuthal velocity component induced by rotation - u fluctuating streamwise component of velocity u=u(, t) - v fluctuating azimuthal component of velocity v = v (, t) - u phase-averaged fluctuating streamwise component of velocity u=u(0) - v phase-averaged fluctuating azimuthal component of velocity v=v() - û average of phase-averaged fluctuating streamwise component of velocity (u()) over cases I-1, II-1 and III-1 û = û() - average of phase-averaged fluctuating azimuthal component of velocity (v()) over cases I-1, II-1 and III-1 - u fluctuating streamwise component of velocity corrected for non-uniformity of probe rotation and/or phase-related vibration u = u(0, t) - v fluctuating azimuthal component of velocity corrected for non-uniformity or probe rotation and/or phase-related vibration v=v (, t) - u ² rms value of corrected fluctuating streamwise component of velocity - rms value of corrected fluctuating azimuthal component of velocity - phase or azimuthal position of X-probe     

Instability of a plane axisymmetric flow with a shear discontinuity in the velocity profile

It is proposed to investigate the stability of a plane axisymmetric flow with an angular velocity profile (r) such that the angular velocity is constant when r < r_O – L and r > r_O + L but varies monotonically from ₁ to ₂ near the point r_O, the thickness of the transition zone being small L r_O, whereas the change in velocity is not small ¦₂ – ₁¦ ₂, ₁. Obviously, as L O short-wave disturbances with respect to the azimuthal coordinate (k=m/r_O 1/r_O) will be unstable with a growth rate-close to the Kelvin—Helmholtz growth rate. In the case L=O (i.e., for a profile with a shear-discontinuity) we find the instability growth rate _O and show that where the thickness of the discontinuity L is finite (but small) the growth rate does not differ from _O up to terms proportional to kL 1 and 1/m 1. Using this example it is possible to investigate the effect of rotation on the flow stability. It is important to note that stabilization (or destabilization) of the flow in question by rotation occurs only for three-dimensional or axisymmetric perturbations.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 111–114, January–February, 1985.     

An application of nonsymmetric Lax-Milgram lemma to nonassociated plasticity

Based on a general assumption for plastic potential and yield surface, some properties of the nonassociated plasticity are studied, and the existence and uniqueness of the distribution of incremental stress and displacement for work-hardening materials are proved by using nonsymmetric Lax-Milgram lemma, when the work-hardening parameter A>F/Q/–F/, Q/.     

Zum Impuls-, Wärme- und Stoffaustausch in turbulenten Strömungen reagierender Binärgemische

Zusammenfassung Die in Teil I vorgestellten Reynolds 'schen Gleichungen und Transportgleichungen werden für Strömungen mit Grenzschichtcharakter angegeben. Weiter werden Integralbedingungen mitgeteilt. Nach einer Diskussion über die Schließung des Gleichungssystems werden Lösungsverfahren besprochen. Dabei wird speziell auf Integralverfahren eingegangen.

---

About the transfer of momentum, heat and mass in turbulent flows of binary mixturesPart II: Thin shear flow layers

The Reynolds equations and transport equations given in part I are presented for thin shear flow layers. Integral relations are given. After a discussion of the closure problem methods of solution are described. Specially integral methods are discussed.

---

Formelzeichen c Massenkonzentration der Komponente - c_t charakteristische Konzentrationsschwankung - c_o Bezugskonzentration - c spezifische Wärme bei konstantem Druck - c_f Reibungsbeiwert - c_D Dissipationsintegral - c_E Entrainment-Funktion - c Schubspannungsintegral - D binsrer Diffusionskoeffizient - H Formparameter - H₁₂ Formparameter - H₃₂ Formparameter - j Kassendiffusionsstrom - L Bezugslänge - p Druck - p_t charakteristische Druckschwankung - p_o Bezugsdruck - Pr Prandtl-Zahl - q Wärmestrom - q²/2 kinetische Energie der Schwankungsbewegung - Re_L mit L gebildete Reynolds-Zahl - Re mit gebildete Reynolds-Zahl - Re₂ mit ₂ gebildete Reynolds-Zahl - Sc Schmidt-Zahl - T absolute Temperatur - T_t charakteristische TemperaturSchwankung - T_o Bezugstemperatur - u,v,w Geschwindigkeitskomponenten - u_t charakteristische Geschwindigkeitsschwankung - u_o Bezugsgeschwindigkeit - U=/ü dimensionslose. x-Komponente der Geschwindigkeit - x,y,z Komponenten des Ortsvektors Griechische Symbole Grenzschichtdicke - ₁ Verdrängungsdicke - ₂ Impulsverlustdicke - ₃ Energieverlustdicke - _T Enthalpieverlustdicke - _c Konzentrationsverlustdicke - =d/dx Parameter für die Grenzschichtabsch:atzung - turbulente Impulsaustauschgröße - _D turbulente Stoffaustauschgröße - _q turbulente Energieaustauschgröße - Dissipationsfunktion - Wärmeleitfähigkeit - dynamische Viskosität - v=/ kinematische Viskosität - Dichte - Produktionsdichte - Schubspannung Indizes mol molekularer Anteil - tur turbulenter Anteil - res resultierender Anteil - Außenrand der Grenzschicht - w Wand     

Heat transfer in magnetohydrodynamic flow near a stagnation point

Summary Steady, axisymmetric, magnetohydrodynamic flow with a stagnation point on an infinite plane wall is considered with a magnetic field applied normal to the wall. Solutions are obtained in the form of series for the velocity, magnetic field and temperature when the magnetic field parameter () and the ratio of viscosity to magnetic diffusivity () are small. The case=O(1) is considered briefly when solutions which Meyer³⁾ obtained by physical order-of-magnitude arguments are derived mathematically as expansions in. Some remarks are made on the consistency of extending the results to flow within the boundary layer near the nose of a bluff body.     

Analysis of slip and temperature jump coefficients in a binary gas mixture

The possibility of simplifying the formulas obtained by the Maxwell-Loyalka method for the velocity _u, temperature _T and diffusion _d slip coefficients and the temperature jump coefficient in a binary gas mixture with frozen internal degrees of freedom of the molecules is considered. Special attention is paid to gases not having sharply different physicochemical properties. The formulas are written in a form convenient for use without linearization in the thermal diffusion coefficient. They are systematically analyzed for mixtures of inert gases, N₂, O₂, CO₂, and H₂ at temperatures extending from room temperature to 2500°K. It is shown that for the molecular weight ratios m^* = m₂/m₁ considered the expressions for _u and can be radically simplified. With an error acceptable for practical purposes (up to 10%) it is possible to employ expressions of the same structural form as for a single-component gas: for _u if 1 m^* 6, and for if 1 m^* 3. When 1 m^* 2 the expression for _T can be simplified with a maximum error of 5%. Within the limits of accuracy of the method the expression for _t can be linearized in the thermal diffusion coefficient. An approximate expression convenient for practical calculations is proposed for _d Finally, the , _u, and _T for a single-component polyatomic gas with easy excitation of the internal degrees of freedom of the molecules are similarly analyzed; it is shown that these expressions can be considerably simplified.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 152–159, November–December, 1990.     

Heat transfer from laminar flow in ducts with elliptic cross-section

Summary The cooling of a hot fluid in laminar Newtonian flow through cooled elliptic tubes has been calculated theoretically. Numerical data have been computed for the two values 1.25 and 4 of the axial ratio of the elliptic cross-section . For =1.25 the influence of non-zero thermal resistance between outmost fluid layer and isothermal surroundings has also been investigated. Special attention has been given to the distribution of heat flux around the perimeter; when increases the flux varies more with the position at the circumference. This positional dependence becomes less pronounced, however, as the (position-independent) thermal resistance of the wall increases.Flattening of the conduit, while maintaining its cross-sectional area constant, improves the cooling. Comparison with rectangular pipes shows that this improvement is not as marked with elliptic as with rectangular pipes.Nomenclature A _k =A _{m, n} coefficients of expansion (6) - a, b half-axes of ellipse, b<a - a _p =a _{r, s} coefficients of representation (V) - D hydraulic diameter, = 4S/P; S = cross-sectional area, P = perimeter - D _e equivalent diameter, according to (13) - n coordinate (outward) normal to the tube wall - T temperature of fluid - T _i temperature of fluid at the inlet - T _s temperature of surroundings - v ₀ mean velocity of fluid - v _z longitudinal velocity of fluid - x, y carthesian coordinates coinciding with axes of ellipse - z coordinate in flow direction - , dimensionless half-axes of ellipse, =a/D and =b/D - _t heat transfer coefficient from fluid at bulk temperature to surroundings; equation (11) - _w heat transfer coefficient at the wall; equation (3) - axial ratio of ellipse, = a/b = / - , , , dimensionless coordinates; =x/D, =y/D, =z/D, =n/D - dimensionless temperature, = (T–T _s)/(T _i–T _s) - ₀ cup-mixing mean value of ; equation (10) - thermal conductivity of fluid - _m,n = _k eigenvalue - c volumetric heat capacity of fluid - _{m, n} = _k = _k eigenfunction; equations (6) and (I) - Nu total Nusselt number, = _t D/ - Nusselt number at large distance from the inlet - Nu _w wall Nusselt number, = _w D/, based on _w - Pé Péclet number, = ₀ Dc/     

Brittle to Quasi-Brittle Transition

Abstract. The present study focuses on the kinetic and non-deterministic aspects of the brittle to quasi-brittle transition. A solid is approximated by a lattice formed by the interacting continuum particles and the evolution of damage is estimated using particle dynamics. The onset of transition is measured by the rate of the change of correlation length. The proposed method is illustrated on the examples of creep rupture, strain localization and dynamic expansion of a circular void in a brittle plate.Sommario. Viene posta l'attenzione sugli aspetti cinetici e non deterministici della transizione dal comportamento fragile a quello quasi-fragile. Un solido viene approssimato da un reticolo formato da particelle interagenti e l'evoluzione del danno viene stimata tramite la dinamica delle particelle. L'inizio della transizione viene misurato tramite la variazione della lunghezza di correlazione. Il metodo proposto viene illustrato su esempi di rottura per creep, localizzazione della deformazione e l'espansione di un foro circolare in una piastra fragile.     

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     

On surface waves in generalized thermoelasticity

Knowles' representation theorem for harmonically time-dependent free surface waves on a homogeneous, isotropic elastic half-space is extended to include harmonically time-dependent free processes for thermoelastic surface waves in generalized thermoelasticity of Lord and Shulman and of Green and Lindsay.r , , r , , .This work was done when author was unemployed.     

Unsteady convective diffusion in a rotating parallel plate channel

The paper presents an exact analysis of the dispersion of a passive contaminant in a viscous fluid flowing in a parallel plate channel driven by a uniform pressure gradient. The channel rotates about an axis perpendicular to its walls with a uniform angular velocity resulting in a secondary flow. Using a generalized dispersion model which is valid for all time, we evaluate the longitudinal dispersion coefficientsK _i (i=1, 2, ...) as functions of time. It is shown thatK ₁=0 andK ₃,K ₄, ... decay rapidly in comparison withK ₂. ButK ₂ decreases with increasing (the dimensionless rotation parameter) for values of upto approximately =2.2. ThereafterK ₂ increases with further increase in and its value gets saturated for large values of (say, 500) and does not change any further with increase in . A physical explanation of this anomalous behaviour ofK ₂ is given.

---

Instationäre konvektive Diffusion in einem rotierenden Parallelplattenkanal

Zusammenfassung In dieser Untersuchung wird eine exakte Analyse der Ausbreitung eines passiven Kontaminierungsstoffes in einer zähen Flüssigkeit gegeben, die, befördert durch einen gleichförmigen Druckgradienten, in einem Parallelplattenkanal strömt. Der Kanal rotiert mit gleichförmiger Winkelgeschwindigkeit um eine zu seinen Wänden senkrechte Achse, wodurch sich eine Sekundärströmung ausbildet. Unter Verwendung eines generalisierten, für alle Zeiten gültigen Dispersionsmodells werden die longitudinalen DispersionskoeffizientenK _i (i=1, 2, ...) als Funktionen der Zeit ermittelt. Es wird gezeigt, daßK ₁=0 gilt und dieK ₃,K ₄, ... gegenüberK ₂ schnell abnehmen.K ₂ nimmt ab, wenn , der dimensionslose Rotationsparameter, bis etwa zum Wert 2,2 ansteigt. Danach wächstK ₂ mit bis auf einem Endwert an, der etwa ab =500 erreicht wird. Dieses anomale Verhalten vonK ₂ findet eine physikalische Erklärung.

---

List of symbols C solute concentration - D molecular diffusivity - K _i longitudinal dispersion coefficients - 2L depth of the channel - P ₀ dimensionless pressure gradient along main flow - Pe Péclet number - q velocity vector - Q _x,Q _y mass flux along the main flow and the secondary flow directions - dimensionless average velocity along the main flow direction - (x, y, z) Cartesian co-ordinates Greek symbols dimensionless rotation parameter - the inclination of side walls withx-axis - kinematic viscosity - fluid density - dimensionless time - angular velocity of the channel - dimensionless distance along the main flow direction - dimensionless distance along the vertical direction - dimensionless solute concentration - integral of the dispersion coefficientK ₂() over a time interval     

Asymptotics of the homogenized moduli for the elastic chess-board composite

We find the asymptotic behavior of the homogenized coefficients of elasticity for the chess-board structure. In the chess board white and black cells are isotropic and have Lamé constants (, ,) and (, ) respectively. We assume that the black cells are soft, so 0. It turns out that the Poisson ratio for this composite tends to zero with .     

Hydrodynamical changes due to polymer migration in very dilute solutions

The slip hypothesis, based on thermodynamical arguments, has been extended to obtain the flow characteristics of polymer solutions flowing in a nonhomogeneous flow field. An asymptotic analysis, valid for both channel and falling film flows, is presented that predicts the flow enhancement due to polymer migration. Concentration-viscosity coupling is shown to be a critical factor in the hydrodynamic analysis. The analysis, which essentially provides an upper bound on flow enhancement, explicitly accounts for the influence of wall shear stress, initial polymer concentration etc. A comparison with the pertinent experimental data shows reasonable agreement. c concentration - c ₀ concentration in shear-free region - c _i initial concentration - d rate of deformation tensor - g acceleration due to gravity - g ₁ function defined in eq. [13] or [15] - g ₂ function defined in eq. [18] or [20] - H half-channel thickness or film thickness - K gas law constant - L length of the channel or film - q flow rate per unit width - q ^* normalized flow rate - T temperature - v velocity - V mean velocity - y transverse distance - y _c location of solvent layer - _w /µ _s - _w /c ₀ KT - /t convected derivative - dimensionless cenentration,c/c ₀ - _c dimensionless interface concentration - _w dimensionless wall concentration - relaxation time - µ _eff effective viscosity - µ _s solvent viscosity - dimensionless transverse distance,y/H - _c dimensionless interface location - density - stress tensor - _w wall shear stress - c _i KT/ _w - ns no slip NCL-Communication No. 3155     

Heteroclinic orbits in singular systems: A unifying approach

We consider singularly perturbed systems , such that=f(, _o, 0). _o ^m , has a heteroclinic orbitu(t). We construct a bifurcation functionG(, ) such that the singular system has a heteroclinic orbit if and only ifG(, )=0 has a solution=(). We also apply this result to recover some theorems that have been proved using different approaches.     

On an inversion theorem for conical regions in elasticity theory

Summary Two-dimensional stress singularities in wedges have already drawn attention since a long time. An inverse square-root stress singularity (in a 360° wedge) plays an important role in fracture mechanics.Recently some similar three-dimensional singularities in conical regions have been investigated, from which one may be also important in fracture mechanics.Spherical coordinates are r, , . The conical region occupied by the elastic homogeneous body (and possible anisotropic) has its vertex at r=0. The mantle of the cone is described by an arbitrary function f(, )=0. The displacement components be u. For special values of (eigenvalues) there exist states of displacements (eigenstates) % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],which may satisfy rather arbitrary homogeneous boundary conditions along the generators.The paper brings a theorem which expresses that if is an eigenvalue, then also-1- is an eigenvalue. Though the theorem is related to a known theorem in Potential Theory (Kelvin's theorem), the proof has to be given along quite another line.

---

Zusammenfassung Zwei-dimensionale Spannungssingularitäten in keilförmigen Gebieten sind schon längere Zeit untersucht worden und neuerdings auch ähnliche drei-dimensionale Singularitäten in konischen Gebieten.Kugelkoordinaten sind r, , . Das konische Gebiet hat seine Spitze in r=0. Der Mantel des Kegels lässt sich beschreiben mittels einer willkürlichen Funktion f(, )=0. Die Verschiebungskomponenten seien u. Für spezielle Werte von (Eigenwerte) bestehen Verschiebunszustände % MathType!MTEF!2!1!+-% feaafiart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqef0uAJj3BZ9Mz0bYu% H52CGmvzYLMzaerbd9wDYLwzYbItLDharqqr1ngBPrgifHhDYfgasa% acOqpw0xe9v8qqaqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbba9q8Wq% Ffea0-yr0RYxir-Jbba9q8aq0-yq-He9q8qqQ8frFve9Fve9Ff0dme% GabaqaaiGacaGaamqadaabaeaafiaakabbaaa6daaahjxzL5gapeqa% aiaadwhadaWgaaWcbaGaeqOVdGhabeaakiabg2da9iaadkhadaahaa% WcbeqaaiabeU7aSbaakiaadAgadaWgaaWcbaGaeqOVdGhabeaakiaa% cIcacqaH7oaBcaGGSaGaeqiUdeNaaiilaiabfA6agjaacMcaaaa!582B!\[u_\xi = r^\lambda f_\xi (\lambda ,\theta ,\Phi )\],welche homogene Randwerte der Beschreibenden des Kegels entlang genügen.Das Bericht bringt ein Theorem, welches aussagt, das und =–1– beide Eigenwerte sind.

     

Über den Einfluß von Querkrümmung und Dichte bei inhomogenen turbulenten Freistrahlen

Zusammenfassung Zur Berechnung turbulenter Strömungen wird das k--Modell im Ansatz für die turbulente Scheinzähigkeit erweitert, so daß es den Querkrümmungs- und Dichteeinfluß auf den turbulenten Transportaustausch erfaßt. Die dabei zu bestimmenden Konstanten werden derart festgelegt, daß die bestmögliche Übereinstimmung zwischen Berechnung und Messung erzielt wird. Die numerische Integration der Grenzschichtgleichungen erfolgt unter Verwendung einer Transformation mit dem Differenzenverfahren vom Hermiteschen Typ. Das erweiterte Modell wird auf rotationssymmetrische Freistrahlen veränderlicher Dichte angewendet und zeigt Übereinstimmung zwischen Rechnung und Experiment.

---

On the influence of transvers-curvature and density in inhomogeneous turbulent free jets

The prediction of turbulent flows based on the k- model is extended to include the influence of transverse-curvature and density on the turbulent transport mechanisms. The empirical constants involved are adjusted such that the best agreement between predictions and experimental results is obtained. Using a transformation the boundary layer equations are solved numerically by means of a finite difference method of Hermitian type. The extended model is applied to predict the axisymmetric jet with variable density. The results of the calculations are in agreement with measurements.

---

Bezeichnungen Wirbelabsorptionskoeffizient - c_i Massenkonzentration der Komponente i - c_D, c_L, c, c₁, c₂ Konstanten des Turbulenzmodells - d Düsendurchmesser - E bezogene Dissipationsrate - f bezogene Stromfunktion - f Korrekturfunktion für die turbulente Scheinzähigkeit - j turbulenter Diffusionsstrom - k Turbulenzenergie - k_i Schrittweite in -Richtung - K dimensionslose Turbulenzenergie - L turbulentes Längenmaß - M_i Molmasse der Komponente i - p Druck - allgemeine Gaskonstante - r Querkoordinate - r_0,5 Halbwertsbreite der Geschwindigkeit - r_0,5c Halbwertsbreite der Konzentration - T Temperatur - u Geschwindigkeitskomponente in x-Richtung - v Geschwindigkeitskomponente in r-Richtung - x Längskoordinate - y allgemeine Funktion - Y_i diskreter Wert der Funktion y - Relaxationsfaktor für Iteration - turbulente Dissipationsrate - transformierte r-Koordinate - kinematische Zähigkeit - Exponent - transformierte x-Koordinate - Dichte - _k, Konstanten des Turbulenzmodells - Schubspannung - allgemeine Variable - Stromfunktion - Turbulente Transportgröße Indizes 0 Strahlanfang - m auf der Achse - r mit Berücksichtigung der Krümmung - t turbulent - mit Berücksichtigung der Dichte - im Unendlichen - Schwankungswert oder Ableitung einer Funktion - – Mittelwert Herrn Professor Dr.-Ing. R. Günther zum 70. Geburtstag gewidmet     

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.     

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.     

On the Ellipticity of the Middle Surface of a Shell and its Application to the Asymptotic Analysis of ‘Membrane’ Shells

We consider a surface S = (), where ² is a bounded, connected, open set with a smooth boundary and : ³ is a smooth map; let () denote the components of the two-dimensional linearized strain tensor of S and let 0 with length 0 > 0. We assume the the norm ,|| ()||0, in the space V0() = { H¹() × H¹() × L²(); = 0 on 0 } is equivalent to the usual product norm on this space. We then establish that this assumption implies that the surface S is uniformly elliptic and that we necessarily have 0 = .     

On the continuum modelling of porous media containing fluid: a molecular viewpoint with particular attention to scale

Mass conservation and linear momentum balance relations for a porous body and any fluid therein, valid at any given length scale in excess of nearest-neighbour molecular separations, are established in terms of local weighted averages of molecular quantities. The mass density field for the porous body at a given scale is used to identify its boundary at this scale, and a porosity field is defined for any pair of distinct length scales. Specific care is paid to the interpretation of the stress tensor associated with each of the body and fluid at macroscopic scales, and of the force per unit volume each exerts on the other. Consequences for the usual microscopic and macroscopic viewpoints are explored.Nomenclature material system; Section 2.1. - porous body (example of a material system); Sections 2.1, 3.1, 4.1 - fluid body (example of a material system); Sections 2.1, 3.1, 4.1 - weighting function; Sections 2.1, 2.3 - ,h weighting function corresponding to spherical averaging regions of radius and boundary mollifying layer of thicknessh; Section 3.2 - Euclidean space; Section 2.1 - V space of all displacements between pairs of points in; Section 2.1 - mass density field corresponding to; (2.3)₁ - ^P , ^f mass density fields for , ; (4.1) - P momentum density field corresponding to; (2.3)₂ - v velocity field corresponding to; (2.4) - S _r (X) interior of sphere of radiusr with centre at pointx; (3.3) - boundary ofany region - region in which ^p > 0 with = ,h; (3.1) - subset of whose points lie at least+h from boundary of ; (3.4) - abbreviated versions of ; Section 3.2, Remark 4 - strict interior of ; (3.7) - analogues of for fluid system ; Section 3.2 - general version of corresponding to any choice of weighting function; (4.6) - interfacial region at scale; (3.8) - ₀ scale of nearest-neighbour separations in ; Section 3.2. Remark 1 - porosity field at scales ( ₁; ₂); (3.9) - pore space at scales ( ₁; ₂); (3.12)