Solution of a contact problem for a plate with a deformable insert

Contact problems with friction are solved for a rectangular plate with a circular hole into which a ring plate (insert) is placed with a small clearance. Two versions of contact boundary conditions are formulated. According to the proposed approximate formulation of the problem, the boundary conditions in both versions are satisfied not at the actual contact points but at specified pairs of points. Therefore, it is sufficient to determine attachment, slip, contact, and contact-free regions on just one of the contacting contours. The finite-element method and the Boussinesq principle are used to solve the problem. One of the versions of boundary conditions, compared to the other, gives smaller values for the strain energies of the plate and insert, the stress-concentration coefficient, and the lengths of attachment and contact regions. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 5, pp. 216–226, September–October, 1999.     

Analysis of nonlinear sliding structures by modified stochastic linearization methods

Probabilistic characteristics of a sliding structure is investigated by using new versions of stochastic linearization technique. The structure is composed of base part and upper part, which are connected to each other in a spring-damping system. Coulomb friction between the base structure and earth ground is considered. Two alternative versions of stochastic linearization approach, suggested by X. Zhang and I. Elishakoff, respectively, are applied to such a sliding structure to evaluate its statistical properties. Compared with the results of Monte Carlo simulation, the two new approaches are performing much better than the conventional one in their applications to the sliding structure. Moreover, numerical results indicate that the criterion proposed by Elishakoff turns out to be superior to all other versions in the problem under study. Numerical results also suggest that the entire structure may be replaced by the rigid body in the sliding problem as long as the difference of velocity responses are considered less important than those of displacement responses.     

One-dimensional filtration of immiscible liquids in a nonuniform porous medium

As is known, the differential equation for two-phase filtration with account for capillarity was obtained in [1], and later integrated numerically for the case of a uniform stratum of finite length in [2]. Other versions of the solution of the Rapoport-Leas equation or the system which is equivalent to it are known [3, 4]. This article presents the results of a numerical solution of an analogous problem with account for nonuniform permeability of the stratum.The authors wish to thank T. V. Startsev and L. Kh. Aminov for assistance in performing the calculations on the Ural-3 computer.     

Estimates of the influence of the transverse strain in a contact problem with free boundary

We use the contact problem with free boundary for a “heavy cylindrical shell-over-support reinforcement ring” system as an example to study the influence of transverse shear on the stress state in the shell. To obtain the equilibrium equations for the shell and the reinforcement ring, we apply an express algorithm developed by one of the authors to take into account transverse strains in Kirchhoff versions of the theory of shells. The contact problem with free boundary is solved by the generalized reaction method proposed earlier. We use numerical examples to show that the corrections introduced in the stress state by taking into account the transverse shears are by an order of magnitude larger than the traditional estimate of errors in the Kirchhoff hypotheses.     

Action of an elliptic punch on a transversally isotropic half-space

The 3D contact problem on the action of a punch elliptic in horizontal projection on a transversally isotropic elastic half-space is considered for the case in which the isotropy planes are perpendicular to the boundary of the half-space. The elliptic contact region is assumed to be given (the punch has sharp edges). The integral equation of the contact problem is obtained. The elastic rigidity of the half-space boundary characterized by the normal displacement under the action of a given lumped force significantly depends on the chosen direction on this boundary. In this connection, the following two cases of location of the ellipse of contact are considered: it can be elongated along the first or the second axis of Cartesian coordinate system on the body boundary. Exact solutions are obtained for a punch with base shaped as an elliptic paraboloid, and these solutions are used to carry out the computations for various versions of the five elastic constants. The structure of the exact solution is found for a punch with polynomial base, and a method for determining the solution is proposed.     

Axisymmetric contact between an annular rough punch and a surface nonuniform foundation

The axisymmetric contact problem of interaction between a two-layer foundation and a rigid annular punch is considered under the assumption that the surface nonuniformity of the upper layer and the shape of the punch base are described by rapidly varying functions. The integral equation of the problem containing two rapidly varying functions is derived, and two versions of the problem are considered. Their solutions were first constructed by the generalized projection method. As an illustration, the model problem is analyzed numerically to demonstrate the high efficiency of the method.     

Thermoelastodynamic disturbances in a half-space under the action of a buried thermal/mechanical line source

The transient dynamic coupled-thermoelasticity problem of a half-space under the action of a buried thermal/mechanical source is analyzed here. This situation aims primarily at modeling underground explosions and impulsively applied heat loadings near a boundary. Also, the present basic analysis may yield the necessary field quantities required to apply the Boundary Element Method in more complicated thermoelastodynamic problems involving half-plane domains. A material response for the half-space predicted by Biots thermoelasticity theory is assumed in an effort to give a formulation of the problem as general as possible (within the confines of a linear theory) . The loading consists of a concentrated thermal source and a concentrated force (mechanical source) having arbitrary direction with respect to the half-plane surface. Both thermal and mechanical line sources are situated at the same location in a fixed distance from the surface. Plane-strain conditions are assumed to prevail. Our problem can be viewed as a generalization of the classical Nakano–Lapwood–Garvin problem and its recent versions due to Payton (1968) and Tsai and Ma (1991) . The initial/boundary value problem is attacked with one- and two-sided Laplace transforms to suppress, respectively, the time variable and the horizontal space variable. A 9×9 system of linear equations arises in the double transformed domain and its exact solution is obtained by employing a program of symbolic manipulations. From this solution the two-sided Laplace transform inversion is then obtained exactly through contour integration. The one-sided Laplace transform inversion for the vertical displacement at the surface is obtained here asymptotically for long times and numerically for short times.     

On Singular Limits of Mean-Field Equations

Mean-field equations arise as steady state versions of convection-diffusion systems where the convective field is determined by solution of a Poisson equation whose right-hand side is affine in the solutions of the convection-diffusion equations. In this paper we consider the repulsive coupling case for a system of two convection-diffusion equations. For general diffusivities we prove the existence of a unique solution of the mean-field equation by a variational analysis of a saddle point problem (usually without coercivity). Also we analyze the small-Debye-length limit and prove convergence to either the so-called charge-neutral case or to a double obstacle problem for the limiting potential depending on the data.     

Spatial Evolution of Nonstationary Perturbations in Hamel Flow

The propagation of small perturbations in longitudinally inhomogeneous flows is investigated. The evolution of the perturbations is studied with reference to the radial flow of a viscous incompressible fluid between plane nonparallel walls, the simplest inhomogeneous flow. Using a generalized method of variation of constants, the corresponding boundary-value problem is reduced to an infinite-dimensional evolutionary system of ordinary differential equations for the complex amplitudes of the eigensolutions of a locally homogeneous problem. Physically, the method can be interpreted as a representation of the perturbation evolution process via two concomitant processes: the independent amplification (attenuation) of normal modes of the locally homogeneous problem and the rescattering of these modes into each other on local inhomogeneities of the base flow. The calculations show that reduced versions of the method are promising for describing the linear stage of laminar-turbulent transition in a boundary layer.     

Equilibrium of an internal transverse crack in a semiinfinite elastic body with thin coating

Static elasticity problems for a half-plane and a strip weakened by a rectilinear transverse crack are studied. In each case, the upper boundary of the body is reinforced by a flexible patch. Various versions of conditions on the lower boundary are considered in the case of the strip. The crack is maintained in the open state by distributed normal forces. The method of generalized integral transforms reduces solving the problem for the equations of equilibriumto solving a singular integral equation of the first kind with the Cauchy kernel with respect to the derivative of the crack opening function. The solutions of the integral equation are constructed by the small parameter and collocation methods for various combinations of the geometric and physical parameters of the problem, and the structure of the solutions is analyzed. The values of the stress intensity factor (SIF) near the crack vertex are obtained.     

Numerical investigation of the effect of winglets on the evolution of an aircraft vortex wake

The effect of winglets on the aerodynamic characteristics of a heavy aircraft and the parameters of the vortex wake behind it is investigated for the landing regime. The solution of the complete problem is obtained by breaking up the wake into three regions, namely, the near, intermediate, and far flowfields, in which the corresponding subproblems are successively solved. The wake flowfields are obtained at different distances from the aircraft, together with the distributions of the mean azimuthal velocity over the vortex radius, and the lifetime of the vortex system is estimated for several one-and two-element wingtip versions. All the results are compared with those for the baseline layout with no winglets.     

Mixing dynamics of a vapor-condensate cloud with the ambient air

The mixing of a wet vapor with a gas is studied using analytical and numerical models. The one-dimensional problem of diffusion mixing accompanied by phase transitions is solved in a self-similar formulation. The versions of mixing of the vapor with a cold and warm gas and with a superheated vapor are analyzed. The atmospheric diffusion of immediate emissions containing water vapor and condensate is modeled numerically in a three-dimensional formulation. A study is made of the evolution of hydrodynamic, concentration, and temperature fields as a function of the initial emission parameters (temperature and humidity) and ambient air parameters. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 3, pp. 114–127, May–June, 2008.     

Effective medium method in the problem of axial elastic shear wave propagation through fiber composites

The effective medium method (EMM) is applied to the solution of the problem of monochromatic elastic shear wave propagation through matrix composite materials reinforced with cylindrical unidirected fibers. The dispersion equations for the wave numbers of the mean wave field in such composites are derived using two different versions of the EMM. Asymptotic solutions of these equations in the long and short wave regions are found in closed analytical forms. Numerical solutions of the dispersion equations are constructed in a wide region of frequencies of the incident field that covers long, middle and short wave regions of the mean wave field. Velocities and attenuation factors of the mean wave fields in the composites obtained by different versions of the EMM are compared for various volume concentrations and properties of the inclusions. The main discrepancies in the predictions of different versions of the EMM are indicated, analyzed and discussed.     

Numerical solution of the problem of an ice sheet under a moving load

A method for analyzing the bending of an ice sheet subjected to a moving load is proposed. The problem is solved in a dynamic formulation. The algorithm of solution is based on the finiteelement method and the finite-difference method. The method proposed allows one to determine the stress-strain state of an ice sheet for any law of motion of a load over ice. Two versions of initial conditions are considered. Examples of calculations are given. Komsomol’sk-on-Amur State Technical University, Komsomol’sk-on-Amur 681013. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 4, pp. 243–248, July–August, 1999.     

An iterative stabilized fractional step algorithm for numerical solution of incompressible N–S equations

Stabilized fractional step algorithm has been widely employed for numerical solution of incompressible Navier–Stokes equations. However, smaller time step sizes are required to use for existing explicit and semi‐implicit versions of the algorithm due to their fully or partially explicit nature particularly for highly viscous flow problems. The purpose of this paper is to present two modified versions of the fractional step algorithm using characteristic based split and Taylor–Galerkin like based split. The proposed modified versions of the algorithm are based on introducing an iterative procedure into the algorithm and allow much larger time step sizes than those required to the preceding ones. A numerical study of stability at acceptable convergence rate and accuracy as well as capability in circumventing the restriction imposed by the LBB condition for the proposed iterative versions of the algorithm is carried out with the plane Poisseuille flow problem under different Reynolds numbers ranging from low to high viscosities. Numerical experiments in the plane Poisseuille flow and the lid‐driven cavity flow problems demonstrate the improved performance of the proposed versions of the algorithm, which are further applied to numerical simulation of the polymer injection moulding process. Copyright © 2005 John Wiley & Sons, Ltd.     

Nonlinear Response of a Parametrically Excited System Using Higher-Order Method of Multiple Scales

Two fundamentally different versions of the method of multiple scales (MMS) are currently in use in the study of nonlinear resonance phenomena. While the first version is the widely used reconstitution method, the second version is proposed by Rahman and Burton [1]. Both versions of the second-order MMS are applied to the differential equation obtained for a parametrically excited cantilever beam with a lumped mass at an arbitrary position. The bifurcation and stability of the obtained response show the difference between the two versions. While the Hopf bifurcation phenomena with no jump is found in the case of second-order MMS version I, both jump-up and jump-down phenomena are observed in second-order MMS version II, which closely agree with the experimental findings. The results are compared with those obtained by numerically integrating the original temporal equation.     

Allowing for Rotations about the Normal in Nonlinear Theories of Shells

Consideration is given to three versions of nonlinear strain–displacement relations in the case of small strains and moderately small angles of rotation: (i) relations that neglect rotations about the normal in conformity with the hypotheses of the Donnel–Mushtary–Vlasov theory; (ii) relations, derived from the elasticity equations using Novozhilov's tensor, that exactly allow for rotations; and (iii) relations, proposed by Sanders, that allow for rotations but neglect shear strains. These versions are compared by comparing the solutions of the stability problem for a corrugated cylindrical shell. It is established that the critical loads are close when rotations are allowed for exactly and when Sanders' technique is used     

New integrable problems in the dynamics of rigid bodies with the kovalevskaya configuration. I - The case of axisymmetric forces

Despite the rarity of integrable problems in rigid body dynamics, amazing relative abundance of these problems is observed when the body has the famous Kovalevskaya configuration A=B=2C. In the present work, the general problem of motion of such a rigid body about a fixed point under the action of axially symmetric conservative forces is considered. The classical problem of motion of a heavy rigid body or a gyrostat and the problem of motion of a body in an ideal fluid are special cases.Three new integrable cases valid for Kovalevskaya's configuration are introduced of which one is general and the other two are restricted to the zero level of the cyclic integral. It is also shown that all the previously known results concerning the preesent problem are special versions of four different cases.     

RANS-based numerical simulation and visualization of the horseshoe vortex system in the leading edge endwall region of a symmetric body

This contribution is aimed at analyzing the capabilities of popular two-equation turbulence models to predict features of 3D flow fields and endwall heat transfer near the blunt edge of a symmetric body mounted on a plate. The configuration studied experimentally by Praisner and Smith is considered. Results obtained with the in-house CFD code SINF and the commercial package ANSYS–CFX are presented and compared. Prediction capabilities of the low-Re Wilcox turbulence model and two versions of the Menter SST model, the original and the modified one, are analyzed in comparison with the experimental data. Special attention is paid to grid sensitivity of the numerical solutions. Advanced visualization of the vortex structures computed is performed with author’s visualization tool HDVIS. It has been established that the Wilcox model is not capable of predicting the development of a multiple-vortex system observed in the experiment upstream of the body leading edge. Both versions of the MSST model produce qualitatively correct results, with a considerable superiority of the modified version when compared with the quantitative data.     

Multi-step transversal and tangential linearization methods applied to a class of nonlinear beam equations

A novel and continuously parameterized form of multi-step transversal linearization (MTrL) method is developed and numerically explored for solving nonlinear ordinary differential equations governing a class of boundary value problems (BVPs) of relevance in structural mechanics. A similar family of multi-step tangential linearization (MTnL) methods is also developed and applied to such BVP-s. Within the framework of MTrL and MTnL, a BVP is treated as a constrained dynamical system, i.e. a constrained initial value problem (IVP). While the MTrL requires the linearized solution manifold to transversally intersect the nonlinear solution manifold at a chosen set of points across the axis of the independent variable, the essential difference of the present MTrL method from its previous version [Roy, D., Kumar, R., 2005. A multistep transversal linearization (MTL) method in nonlinear structural dynamics. J. Sound Vib. 17, 829–852.] is that it has the flexibility of treating nonlinear damping and stiffness terms as time-variant damping and stiffness terms in the linearized system. The resulting time-variant linearized system is then solved using Magnus’ characterization [Magnus, W., 1954. On the exponential solution of differential equations for a linear operator. Commun. Pure Appl. Math., 7, 649–673.]. Towards numerical illustrations, response of a tip loaded cantilever beam (Elastica) is first obtained. Next, the response of a simply supported nonlinear Timoshenko beam is obtained using a variationally correct (VC) model for the beam [Marur, S., Prathap, G., 2005. Nonlinear beam vibration problems and simplification in finite element model. Comput. Mech. 35(5), 352–360.]. The new model does not involve any simplifications commonly employed in the finite element formulations in order to ease the computation of nonlinear stiffness terms from nonlinear strain energy terms. A comparison of results through MTrL and MTnL techniques consistently indicate a superior quality of approximations via the transversal linearization technique. While the usage of tangential system matrices is common in nonlinear finite element practices, it is demonstrated that the transversal version of linearization offers an easier and more general implementation, requires no computations of directional derivatives and leads to a consistently higher level of numerical accuracy. It is also observed that higher order versions of MTrL/MTnL with Lagrangian interpolations may not work satisfactorily and hence spline interpolations are suggested to overcome this problem.