THE EXACT SOLUTIONS OF ELASTIC－PLASTIC CRACK LINE FIELD FOR MODE Ⅱ PLANE STRESS CRACK

ＴＨＥＥＸＡＣＴＳＯＬＵＴＩＯＮＳＯＦＥＬＡＳＴＩＣ－ＰＬＡＳＴＩＣＣＲＡＣＫＬＩＮＥＦＩＥＬＤＦＯＲＭＯＤＥⅡＰＬＡＮＥＳＴＲＥＳＳＣＲＡＣＫＹｉＺｈｉｊｉａｎ（易志坚）ＷａｎｇＳｈｉｊｉｅ（王士杰）ＷａｎｇＸｉａｎｇｊｉａｎ（王向坚）（Ｒｅｃｅｉ...     

THE EXACT SOLUTIONS OF ELASTIC－PLASTIC CRACK LINE FIELD FOR MODE ⅡPLANE STRESS CRACK

ＴＨＥＥＸＡＣＴＳＯＬＵＴＩＯＮＳＯＦＥＬＡＳＴＩＣ－ＰＬＡＳＴＩＣＣＲＡＣＫＬＩＮＥＦＩＥＬＤＦＯＲＭＯＤＥＩＩＰＬＡＮＥＳＴＲＥＳＳＣＲＡＣＫＹｉＺｈｉｊｉａｎ（易志坚）ＷａｎｇＳｈｉｊｉｅ（王士杰）ＷａｎｇＸｉａｎｇｊｉａｎ（王向坚）（Ｒｅｃｅ...     

STRESS－STBAIN FIELD NEAR CRACK TIP AND CALCULATION OF CRITICAL STRESS OF CRACK PROPAGATION IN A PURE BENDING BEAM OF RECTANGULAR SECTION WITH ONE－SIDED MODEICRACK

ＳＴＲＥＳＳ－ＳＴＢＡＩＮＦＩＥＬＤＮＥＡＲＣＲＡＣＫＴＩＰＡＮＤＣＡＬＣＵＬＡＴＩＯＮＯＦＣＲＩＴＩＣＡＬＳＴＲＥＳＳＯＦＣＲＡＣＫＰＲＯＰＡＧＡＴＩＯＮＩＮＡＰＵＲＥＢＥＮＤＩＮＧＢＥＡＭＯＦＲＥＣＴＡＮＧＵＬＡＲＳＥＣＴＩＯＮＷＩＴＨＯＮＥ－Ｓ...     

THE EXISTENCE OF PERIODIC SOLUTIONS FOR A CLASS OF FUNCTIONAL DIFFERENTIAL EQUATIONS AND THEIR APPLICATION

ＴＨＥＥＸＩＳＴＥＮＣＥＯＦＰＥＲＩＯＤＩＣＳＯＬＵＴＩＯＮＳＦＯＲＡＣＬＡＳＳＯＦＦＵＮＣＴＩＯＮＡＬＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＳＡＮＤＴＨＥＩＲＡＰＰＬＩＣＡＴＩＯＮＺｈａｏＪｉｅ－ｍｉｎ（赵杰民）;ＨｕａｎｇＫｅ－ｌｅｉ（黄克累）...     

CONSTRUCTION OF MODIFIED TAYLOR－GALERKIN FINITE ELEMENTS AND ITS APPLICATION IN COMPRESSIBLE FLOW COMPUTATION

ＣＯＮＳＴＲＵＣＴＩＯＮＯＦＭＯＤＩＦＩＥＤＴＡＹＬＯＲ－ＧＡＬＥＲＫＩＮＦＩＮＩＴＥＥＬＥＭＥＮＴＳＡＮＤＩＴＳＡＰＰＬＩＣＡＴＩＯＮＩＮＣＯＭＰＲＥＳＳＩＢＬＥＦＬＯＷＣＯＭＰＵＴＡＴＩＯＮＣＯＮＳＴＲＵＣＴＩＯＮＯＦＭＯＤＩＦＩＥＤＴＡＹＬＯＲ...     

MOTION EQUATIONS OF MULTILAYERED ELASTIC ELECTROCONDUCTIVE PLATES IN A MAGNETIC FIELD

ＭＯＴＩＯＮＥＱＵＡＴＩＯＮＳＯＦＭＵＬＴＩＬＡＹＥＲＥＤＥＬＡＳＴＩＣＥＬＥＣＴＲＯＣＯＮＤＵＣＴＩＶＥＰＬＡＴＥＳＩＮＡＭＡＧＮＥＴＩＣＦＩＥＬＤＧｏｕＸｉｎｇ－ｈｕａ（苟兴华）（ＳｉｃｈｕａｎＵｎｉｖｅｒｓｉｔｙ,Ｃｈｅｎｇｄｕ）ＺｈａｎｇＦａ...     

CHARACTERISTICS OF SUBDIFFERENTIALS OF FUNCTIONS

ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ（郭兴明）ＣＨＡＲＡＣＴＥＲＩＳＴＩＣＳＯＦＳＵＢＤＩＦＦＥＲＥＮＴＩＡＬＳＯＦＦＵＮＣＴＩＯＮＳ￥ＧｕｏＸｉｎｇｍｉｎｇ（ＲｅｃｅｉｖｅｄＪｕｎｅ１６，１９９５...     

THE EXISTENCE OF PERIODIC SOLUTION OF THE FOURTH ORDINARY NONLINEAR DIFFERENTIAL EQUATION CAUSED BY FLOW－IN

ＴＨＥＥＸＩＳＴＥＮＣＥＯＦＰＥＲＩＯＤＩＣＳＯＬＵＴＩＯＮＯＦＴＨＥＦＯＵＲＴＨＯＲＤＩＮＡＲＹＮＯＮＬＩＮＥＡＲＤＩＦＦＥＲＥＮＴＩＡＬＥＱＵＡＴＩＯＮＣＡＵＳＥＤＢＹＦＬＯＷ－ＩＮＤＵＣＥＤＶＩＢＲＡＴＩＯＮＧｕＱｉｎｇ－ｆａｎｇ（顾清芳）Ｔａ...     

FINITE ELEMENT ANALYSIS FOR THE UNSTEADY NEARSHORE CIRCULATION DUE TO WAVE－CURRENT INTERACTION（II）──TWO－STEP EXPLICIT FINITE ELEMENT METHOD

ＦＩＮＩＴＥＥＬＥＭＥＮＴＡＮＡＬＹＳＩＳＦＯＲＴＨＥＵＮＳＴＥＡＤＹＮＥＡＲＳＨＯＲＥＣＩＲＣＵＬＡＴＩＯＮＤＵＥＴＯＷＡＶＥ－ＣＵＲＲＥＮＴＩＮＴＥＲＡＣＴＩＯＮ（ＩＩ）──ＴＷＯ－ＳＴＥＰＥＸＰＬＩＣＩＴＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＷ...     

DYNAMIC STRESS INTENSITY FACTORS AROUND TWO CRACKS NEAR AN INTERFACE OF TWO DISSIMILAR ELASTIC HALF－PLANES UNDER IN－PLANE SHEAR IMPACT LOAD

ＤＹＮＡＭＩＣＳＴＲＥＳＳＩＮＴＥＮＳＩＴＹＦＡＣＴＯＲＳＡＲＯＵＮＤＴＷＯＣＲＡＣＫＳＮＥＡＲＡＮＩＮＴＥＲＦＡＣＥＯＦＴＷＯＤＩＳＳＩＭＩＬＡＲＥＬＡＳＴＩＣＨＡＬＦ－ＰＬＡＮＥＳＵＮＤＥＲＩＮ－ＰＬＡＮＥＳＨＥＡＲＩＭＰＡＣＴＬＯＡＤＱｉａｎＲ...     

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 = ₂     

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.     

ANALYSIS OF INTERPHASE MECHANICAL BEHAVIOR WITH INTERFACE ELEMENT IN THE COMPOSITE MATERIALS

ＡＮＡＬＹＳＩＳＯＦＩＮＴＥＲＰＨＡＳＥＭＥＣＨＡＮＩＣＡＬＢＥＨＡＶＩＯＲＷＩＴＨＩＮＴＥＲＦＡＣＥＥＬＥＭＥＮＴＩＮＴＨＥＣＯＭＰＯＳＩＴＥＭＡＴＥＲＩＡＬＳＡＮＡＬＹＳＩＳＯＦＩＮＴＥＲＰＨＡＳＥＭＥＣＨＡＮＩＣＡＬＢＥＨＡＶＩＯＲＷＩＴＨＩＮ...     

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.

     

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.     

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.     

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°.     

Cavity nucleation, growth and stress field in elastic-plastic medium

A finite element method is used to study the cavity nucleation and growth in an elastic-plastic medium. The critical cohesion strength, _c at IPM (interface between second phase particle and matrix) is employed as the criterion of cavity formation. Three different values of _c are taken to examine their influence on the overall mechanical behaviour and process of cavity formation.     

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²