Propagation of longitudinal and transverse waves in a multimodulus elastic medium

A problem of propagation of longitudinal and transverse waves in a multimodulus elastic isotropic medium is considered. In the model used, the medium is described by a potential depending on three invariants of strains, which allows the influence of preliminary deformation of the medium on the longitudinal and transverse velocities to be taken into account. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 176–182, July–August, 2009.     

Domain decomposition method with hybrid approximations applied to solve problems of elasticity

To solve two-dimensional boundary-value problems of elasticity, two iteration algorithms of the domain decomposition method are proposed: parallel Neumann–Neumann and sequential Dirichlet–Neumann. They are based on the hybrid boundary–finite-element approximations. The algorithms are proved to converge. The optimal parameters are selected using the minimum-residual and steepest-descent methods. Some plane problems of elasticity are solved as examples, and stationary and nonstationary iteration algorithms in these examples are analyzed for efficiency Translated from Prikladnaya Mekhanika, Vol. 44, No. 11, pp. 18–29, November 2008.     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part I: theoretical analysis and typical solution

Boundary film shear elastic modulus effect is analyzed in a hydrodynamic contact. The contact is one-dimensional composed of two parallel plane surfaces, which are, respectively, rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact is consisted of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. Analysis is thus carried out for a micro Raleigh bearing in this contact. The hydrodynamic contact in this micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity large enough at the film–contact interface. In boundary film area, the viscosity, density and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling, which depend on the film thickness. The analytical approach proposed by Zhang (J Mol Liq 128:60–64, 2006) and Zhang et al. (Int J Fluid Mech Res 30:542–557, 2003) is used for boundary film area. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. Conventional hydrodynamic analysis is used for fluid film area. The present paper presents the theoretical analysis and a typical solution. It is found that for the simulated case the boundary film shear elastic modulus effects on the mass flow through the contact, the overall film thickness of the contact and the carried load of the contact are negligible but the boundary film shear elastic modulus effect on the local film thickness of the contact may be significant when the boundary film thickness is on the 1 nm scale and the contact surfaces are elastic. In Part II will be presented detailed results showing boundary film shear elastic modulus effects in different operating conditions.

---

     

Boundary film shear elastic modulus effect in hydrodynamic contact. Part II: influence of operational parameters

The present paper is the subsequent research of the first part (Theor Comput Fluid Dyn, 2009). It investigates the boundary film shear elastic modulus effect in a hydrodynamic contact in different operating conditions. The hydrodynamic contact is one-dimensional, composed of two parallel plane surfaces, which are respectively rough rigid with rectangular micro projections in profile periodically distributed on the surface and ideally smooth rigid. The whole contact consists of cavitated area and hydrodynamic area. The hydrodynamic area consists of many micro Raleigh bearings which are discontinuously and periodically distributed in the contact. The hydrodynamic contact in a micro Raleigh bearing consists of boundary film area and fluid film area which, respectively, occur in the outlet and inlet zones. In boundary film area, the film slips at the upper contact surface due to the limited shear stress capacity of the film–contact interface, while the film does not slip at the lower contact surface due to the shear stress capacity of the film–contact interface large enough. In boundary film area, the viscosity, density, and shear elastic modulus of the film are varied across the film thickness due to the film–contact interactions, and their effective values are used in modeling which depends on the film thickness. In fluid film area, the film does not slip at either of the contact surfaces, and the shear elastic modulus of the film is neglected. It is found from the simulation results that the boundary film shear elastic modulus influences are normally negligible on the mass flow through the contact, the carried load of the contact and the overall film thickness of the contact, and the boundary film shear elastic modulus would normally influence the local film thickness in an elastic contact when the local film thickness is on the film molecule diameter scale. It is also found that the boundary film shear elastic modulus effect has the tendency of being increased with the reduction of the width of a micro contact. It is increased with the reduction of the boundary film–contact interfacial shear strength or with the increase of the critical boundary film thickness, while it is strongest at certain values of the contact surface roughness, the width ratio of fluid film area to boundary film area, and the lubricant film shear elastic modulus.

---

     

Determining the limiting solutions of nonstationary axisymmetric Hele-Shaw problems

Numerical solution of the Hele-Shaw problem reduces to solution of three boundary-value problems of determining analytic functions of a complex variable in each time step: conformal mapping of the range of the parametric variable to the physical plane, the Dirichlet problems for determining the electric-field strength, and the Riemann-Hilbert problem for calculating partial time derivatives of the coordinates of points of the interelectrode space (the images of the points on the boundary of the parametric plane are fixed). Unlike in the two-dimensional problem, the electric-field strength is determined using integral transformations of an analytic function. Approximation by spline function is performed, and more accurate and steady (than the well-known ones) general solution algorithms for the nonstationary axisymmetric problems are described. Results of a numerical study of the formation of stationary and self-similar configurations are presented. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 87–99, July–August, 2009.     

On determining stresses in rigid inclusions. Plane problems

Plane problems of determining the stress-strain state of an isotropic elastic domain with a rigid inclusion are considered. It is shown that the stress field in the inclusion is uniquely determined. This field is uniform for a plane with an elliptic inclusion, and the stresses at infinity and in the inclusion are related by mutually single-valued formulas. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 183–186, July–August, 2009.     

Existence of solutions of Kraiko's problem

Three initial-boundary-value problems for the equations of gas dynamics are formulated. Successive solution of these problems yields a solution of Kraiko's problem of the isentropic transition of an ideal gas from a homogeneous state of rest to another state of rest with higher or lower density. Solutions are constructed for plane, cylindrical, and spherical layers of an ideal gas. The existence of locally analytic solutions is proved. Ural State Academy of Service, Ekaterinburg 620034. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 41, No. 3, pp. 48–55, May–June, 2000.     

Plane waves in a heavy two-layer liquid excited by a cylinder moving at an angle to the horizontal

A solution of an initial-boundary-value problem for a system of integrodifferential equations which describes the plane waves excited in an initially stationary heavy two-layer ideal fluid by a cylinder moving at an angle to the horizontal is investigated. The homogeneous fluid fractions of different densities are assumed to be separated by an evolving fluid interface (horizontal plane, if the liquid is at rest). An approximate solution of two problems for the waves excited by a cylinder moving with a constant acceleration and an oscillating cylinder is constructed analytically. Nizhnii Novgorod. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 137–152, July–August, 1998.     

Subcritical crack growth in an aging plate with variable elastic modulus

The paper addresses subcritical growth of a crack in a thin isotropic plate made of an aging viscoelastic material with time-dependent elastic modulus. The behavior of the material is described by Arutyunyan’s creep theory. To simulate fracture, a modified Leonov–Panasyuk–Dugdale model and a critical crack opening displacement criterion are used. An equation describing the subcritical growth of the crack is derived assuming that Poisson’s ratio is constant. As an example, the critical loads are determined, and curves of subcritical crack growth are plotted for a specific material. The results are compared with the case of constant elastic modulus     

Mesomechanics of inhomogeneous deformation and short-term damage of elastic bodies: Coupled equations

The coupled equations of inhomogeneous deformation and short-time damage are formulated for problems of three-dimensional elasticity, plane strain, and plane stress. Initial representations based on the stochastic inhomogeneity of macrostrength are given. A special case of microstrength distribution which simplifies the equations is examined. The cylindrical bending of an elastic layer is considered as a numerical example __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 9, pp. 28–37, September 2007.     

Nonlinear stochastic creep problem for an inhomogeneous plane with the damage to the material taken into account

An analytical method for the solution of two-dimensional nonlinear creep problems is developed using as an example the biaxial extension of a plane from a stochastically inhomogeneous material with damage accumulation and the third stage of creep taken into account. The governing creep relation is adopted in accordance with the energetic version of the nonlinear theory of viscous flow. The stochasticity of the material is defined by two random functions of coordinates. Formulas for calculating the stress variance are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 2, pp. 140–146, March–April, 2007.     

Local stress singularities in mixed axisymmetric problems of the bending of circular cylinders

A method of analyzing the near-edge stress state in mixed problems of the deformation of an isotropic cylindrical body is proposed. The method is based on the expansion of the solution of three-dimensional problems of elasticity into a series of Lurie–Vorovich homogeneous basis functions. An asymptotic analysis is performed to find the principal part of the solution of the infinite systems of linear algebraic systems to which the problems are reduced. The type of the stress singularity at the edge of the cylinder is the same as in the mixed problems for a quarter plane. Kummer’s convergence acceleration method is used. The obtained results are validated by testing the boundary conditions and by comparing with results obtained by other authors     

Propagation of unsteady bulk perturbations in an elastic half-plane

At present, the problems of unsteady waves initiated by surface perturbations in an elastic half-space have been studied sufficiently well (see, e.g., [1–5]; a detailed bibliography on this problem can be found in [6]). At the same time, the analytical solutions of the corresponding unsteady problems of bulk perturbations are practically absent. It is these questions as applied to the plane problem that are considered in this paper.     

Seismic Reflection from an Interface Between an Elastic Solid and a Fractured Porous Medium with Partial Saturation

The theory of Tuncay and Corapcioglu (Transp Porous Media 23:237–258, 1996a) has been employed to investigate the possibility of plane wave propagation in a fractured porous medium containing two immiscible fluids. Solid phase of the porous medium is assumed to be linearly elastic, isotropic and the fractures are assumed to be distributed isotropically throughout the medium. It has been shown that there can exist four compressional waves and one rotational wave. The phase speeds of these waves are found to be affected by the presence of fractures, in general. Of the four compressional waves, one arises due to the presence of fractures in the medium and the remaining three are those encountered by Tuncay and Corapcioglu (J Appl Mech 64:313–319, 1997). Reflection and transmission phenomena at a plane interface between a uniform elastic half-space and a fractured porous half-space containing two immiscible fluids, are analyzed due to incidence of plane longitudinal/transverse wave from uniform elastic half-space. Variation of modulus of amplitude and energy ratios with the angle of incidence are computed numerically by taking the elastic half-space as granite and the fractured porous half-space as sandstone material containing non-viscous wetting and non-wetting fluid phases. The results obtained in case of porous half-space with fractures, are compared graphically with those in case of porous half-space without fractures. It is found that the presence of fractures in the porous half-space do affect the reflection/transmission of waves, which is responsible for raising the reflection and lowering the transmission coefficients.     

Two-Parameter Model of a Mode I Crack in an Elastoplastic Body under Plane-Strain conditions

The Dugdale crack model is generalized to the case of plane strain. The governing equations are set up to determine the stresses in the plastic zone. Numerical results from specific problems are analyzed and compared with those for plane stress state and other cases. A relationship between the crack model and K _I-T theory is established in the case of small-scale yielding at the crack tip __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 6, pp. 44–55, June 2005.     

Interaction of a die with a layered elastic foundation

Plane and axisymmetric contact problems for a three-layer elastic half-space are considered. The plane problem is reduced to a singular integral equation of the first kind whose approximate solution is obtained by a modified Multhopp-Kalandiya method of collocation. The axisymmetric problem is reduced to an integral Fredholm equation of the second kind whose approximate solution is obtained by a specially developed method of collocation over the nodes of the Legendre polynomial. An axisymmetric contact problem for an transversely isotropic layer completely adherent to an elastic isotropic half-space is also considered. Examples of calculating the characteristic integral quantities are given. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 165–175, May–June, 2006.     

Numerical study of nonlinear dynamics of hydroelastically connected, plane curvilinear rods

The nonlinear dynamics of hydroelastically connected, plane curvilinear rods is studied. We take into account the reciprocal effect of deformation and hydrodynamics processes, large displacements and strains of the rods, the preliminary static stress-strain state, and the nonstationary fluid flow. A method is proposed for numerical solution of initial boundary-value problems. The effects of the hydroelastic interaction are investigated. The effect of various factors on the dynamics of a damaged pipeline is analyzed. Institute of Mechanics, Nizhnii Novgorod State University, Nizhnii Novgorod 603600. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 1, pp. 212–219, January–February, 1999.     

On the solution of plane boundary-value problems of elasticity in a rectangle

An algorithm for solving plane boundary-value problems of elasticity for a rectangular domain is expounded. The algorithm is based on a complex-valued representation of the general solution to the differential equations of the plane problem and on the use of Lagrange polynomials to satisfy the boundary conditions. The algorithm can quite easily be implemented in a computer program. This is probably the simplest way of solving boundary-value problems of this class __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 1, pp. 97–102, January 2006.     

Plane ovsyannikov vortex: Motion properties and exact solutions

The physical properties of ideal plasma flow described by the Ovsyannikov plane vortex are studied. The particle trajectories and magnetic lines are shown to be plane curves, and an algorithm for describing the motion in three-dimensional space is proposed. Some exact solutions of the submodel are obtained and studied. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 6, pp. 55–68, November–December, 2008.