Motion of a solitary wave over a submerged ridge

The two- and three-dimensional motion of a stationary wave in a layer of liquid with a variable depth was studied earlier by Grimshaw [1, 2], who showed that if the unperturbed state of the liquid were a state of rest, then a certain integrated quantity characterizing the energy of the wave (referred to one element on the leading edge of the wave) was conserved during the motion. In the present investigation (which is based on this property) we shall establish the shape of the wave front for the case of steady motion over an infinite submerged cylindrical ridge; we shall present a model describing a wave of limiting amplitude and shall consider the transient perturbation of the steady-state solution. We should note that one particular case (the motion of a straight front, orthogonal to the axis of the ridge, in a direction parallel to this axis) was considered by another method in [3].     

The moving-load problem in soil dynamics—the vertical displacement approximation

The moving point load problem in soil dynamics is analyzed in the vertical particle displacement approximation. Prior to its motion, the load is stationary. From the instant at which it is set into motion it moves, with constant speed, along a straight path on the (horizontal) planar surface of a semi-infinite elastic medium. The modified Cagniard method for solving transient wave problems is used to determine closed-form expressions for the vertical component of the particle displacement from the elastodynamic wave equation of which only the vertical component is taken into account. The relevant approximation is standard in soil dynamics. Both the cases of “subsonic” and “supersonic” surface load speeds are considered. Methods to include losses in the model are briefly discussed. The study has been initiated with a view to the application of the results to the analysis of the ground motion generated by high-speed trains traveling on a poorly consolidated soil.     

Stochastic seismic response of long-span structures

In this paper, dynamic responses of long-span structures subjected to the action of earthquake with realistic wave speed are analysed. The horizontal or vertical ground motion due to earthquakes is assumed to be a stationary stochastic process, and the seismic waves travel along a horizontal straight line. Expressions for calculating the psd(power spectral density) matrices of structural displacements and internal forces are derived based on three dimensional FEM structural models with the ground motion phase-lags taken into account. A numerical example is given which shows that it is of great importance to consider the effect of the ground motion phase-lags for long-span structures.This work was supported by the National Natural Science Foundation of China.     

Convective Plume Paths in Anisotropic Porous Media

Free convection motions induced by point sources or horizontal line sources of heat are usually assumed to take the form of a vertically orientated plume. In this paper we consider how material anisotropy affects the path of the plume centreline and we show that it is strongly affected by both the anisotropy and the presence of impermeable bounding surfaces. The plume path is a straight line whose angle from the vertical is determined by a balance between the upward buoyancy force, the need for the plume to entrain equal amounts of fluid from the external regions either side of it, and the ease of fluid motion in the direction of the principal axis with the highest permeability.     

Wave resistance of a viscous liquid

In this paper we examine the resistance encountered by a system of normal stresses during its rectilinear motion along the surface of a viscous liquid of infinite depth. The problem is solved in the linear formulation, i.e., it is assumed that amplitudes of the waves which arise are small and the waves are shallow. The solution for the two-and three-dimensional problems is obtained in the general case in closed form. In the two-dimensional case a detailed study is made of the case when a constant pressure p₀, moving with the constant velocity U, is given on a segment of length 2l. In the three-dimen-sional problem the case is studied when the normal stress is concentrated on a segment of a straight line of length 2l, which can replace a ship moving along a straight course with the constant velocity U. The integrals obtained in both cases are studied using the stationary phase method, the application of which for the three-dimensional integrals with respect to a volume with boundaries is justified in §1 of the paper. As a result we obtain equations for the wave resistance in the two- (§2) and three-dimensional (§3) cases.     

On a 3D moving load problem for an elastic half space

The paper deals with 3D dynamic response of an elastic half-space loaded by a point force moving at a constant speed along a straight line on the surface. The problem is formulated within the framework of the asymptotic hyperbolic–elliptic model developed earlier by two of the authors. The validity of the model is restricted to the range of speeds close to the Rayleigh wave speed. Steady-state near-field solutions are derived in terms of elementary functions. Transient analysis of surface motion illustrates peculiarities of the resonance associated with the Rayleigh wave.     

The Mach reflection triple-point locus for internal and external conical diffraction of a moving shock wave

This paper examines the different behavior that occurs for the Mach reflection triple-point loci between the two fundamental axisymmetric cases, these being the external diffraction by a cone and the internal diffraction within a conically contracting channel. From equations derived in this paper using a shock dynamics approach, it has been shown that, for external diffraction over a cone, a possible solution is that the triple-point locus is a straight line which corresponds to the experimental results available, while for internal diffraction along a conically converging channel, it cannot be straight and is, in fact, a convex curve. In the latter case, a transition point is noted on the triple-point locus before which the locus is nearly straight but after which the curvature becomes marked. The second region diminishes as a proportion of the total locus with decreasing half cone angle.For the external case, a set of simple, axisymmetric equations are derived which allow a rapid estimation of the triple point locus angle and the Mach stem strength for any incident shock Mach number and cone angle combination. The equations for internal diffraction are similar and allow a quick computation of both the curved triple-point locus and the strength of the diffracting front of the shock wave. A comparison with experiment has been carried out and agreement is good.     

非保守功能梯度弹性组合曲梁的精确模型及其数值解

基于直法线假设,采用可伸长梁的几何非线性理论,建立了功能梯度材料弹性组合曲梁受切线均布随从力作用下的静态大变形数学模型。该模型不仅计及了轴线伸长,同时也精确地考虑了梁的初始曲率对变形的影响以及轴向变形与弯曲变形之间的耦合效应。用打靶法数值求解了由金属和陶瓷两相材料所构成的一种FGM组合曲梁在沿轴线均布切向随动载荷作用下的非线性平面弯曲问题,给出了不同梯度指标下FGM弹性曲梁随载荷参数大范围变化的平衡路径,并与金属和陶瓷两种单相材料曲梁的相应特性进行了比较。     

Limiting self-similar motion of a solid particle in a free-molecular gas flow escaping from an orifice

A problem of motion of a single spherical solid particle in the far wake of a free-molecular gas flow escaping from an orifice is considered. It is shown that the spatial distributions of macroscopic parameters of the gas are completely determined by functions of one variable: coordinate along an arbitrary straight line normal to the axis of symmetry. Based on this property, a dimensionless equation of particle motion is derived, which has self-similar solutions: trajectories of motion and particle coordinates (traces) on the target for different initial conditions. Conditions of determining the gas flow behavior on the basis of particle traces on the target are considered.     

Numerical simulation of the flow in straight blade agitator with the MPS method

The three‐dimensional flow in a straight blade agitator with free surface on top is simulated using a grid‐free method named moving particle semi‐implicit method. The agitator has six rotor blades matched with six stationary guide blades. The mechanism and phenomena of the flow are investigated in the area between two adjacent stationary guide blades. Eddies near each tip of the rotational blades are predicted, and they move with the motion of the rotor blades but in opposite rotation direction of the rotor. The rotation axis of the eddies is traced and annular vortices, which are made by the eddies, are studied. The pressure pulsation in the rotation process is also predicted with this particle method. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamics of a satellite carrying a point mass moving about it

We study the dynamics of a complex system consisting of a solid and a mass point moving according to a prescribed law along a curve rigidly fixed to the body. The motion occurs in a central Newtonian gravitational field. It is assumed that the orbit of the system center of mass is an ellipse of arbitrary eccentricity.We obtain equations that describe the motion of the carrier (satellite) about its center of mass. In the case of a circular orbit, we present conditions that should be imposed on the law of the relative motion of the mass point carried by the satellite so that the latter preserves a constant attitude with respect to the orbital coordinate system. In the case of a dynamically symmetric satellite, we consider the problem of existence of stationary and nearly stationary rotations for the case in which the carried point moves along the satellite symmetry axis.We consider several problems of dynamics of the satellite plane motion about its center of mass in an elliptic orbit of arbitrary eccentricity. In particular, we present the law of motion of the carried point in the case without eccentricity oscillations and study the stability of the satellite permanent attitude with respect to the orbital coordinate system.     

Theoretical analysis of crack front instability in mode I+III

This paper focusses on the theoretical prediction of the widely observed crack front instability in mode I+III, that causes both the crack surface and crack front to deviate from planar and straight shapes, respectively. This problem is addressed within the classical framework of fracture mechanics, where the crack front evolution is governed by conditions of constant energy-release-rate (Griffith criterion) and vanishing stress intensity factor of mode II (principle of local symmetry) along the front. The formulation of the linear stability problem for the evolution of small perturbations of the crack front exploits previous results of Movchan et al. (1998) (suitably extended) and Gao and Rice (1986), which are used to derive expressions for the variations of the stress intensity factors along the front resulting from both in-plane and out-of-plane perturbations. We find exact eigenmode solutions to this problem, which correspond to perturbations of the crack front that are shaped as elliptic helices with their axis coinciding with the unperturbed straight front and an amplitude exponentially growing or decaying along the propagation direction. Exponential growth corresponding to unstable propagation occurs when the ratio of the unperturbed mode III to mode I stress intensity factors exceeds some “threshold” depending on Poisson's ratio. Moreover, the growth rate of helical perturbations is inversely proportional to their wavelength along the front. This growth rate therefore diverges when this wavelength goes to zero, which emphasizes the need for some “regularization” of crack propagation laws at very short scales. This divergence also reveals an interesting similarity between crack front instability in mode I+III and well-known growth front instabilities of interfaces governed by a Laplacian or diffusion field.     

Stability of the rectilinear motion of a railway Wheelset

A railway wheelset rolling on rails without slip is studied with consideration of the creep hypothesis. The wheelset is represented by two cones having a common base; the rails are represented by two circular cylinders with parallel axes. The kinematic characteristics of undisturbed rolling for the wheelset are determined when its center of mass moves along a straight line; these characteristics are also determined in the case of disturbed motion when the mass center of the wheelset moves along a sinusoidal trajectory. For these modes of motion, the constraint reactions are found with an accuracy up to the second order of smallness with respect to the values of disturbed variables. When the absolutely rigid point contact is replaced by the elastic contact, the creep hypothesis is used, the method of averaging with respect to the fast variables is applied, and a critical velocity above which the rectilinear rolling becomes unstable is determined on the basis of the averaged equations.     

The nature of the triple point singularity in the case of stationary reflection of weak shock waves

The structure of flow in the vicinity of a triple point in the problem of stationary irregular reflection of weak shock waves is numerically investigated within the framework of the Euler model, including the von Neumann paradox range. To improve the accuracy of the solution near singular points a new technology including a grid contracted toward the triple point and the discontinuity fitting is applied. It is shown that in the four-wave flow pattern the curvatures of the tangential discontinuity and the Mach front at the triple point are finite. The singularity is concentrated only in a sector between the reflected wave front and the expansion fan. When the three-wave flow pattern is realized, the curvatures of the tangential discontinuity and both wave fronts at the triple point are infinite. On the range of weak and moderate shock waves the logarithmic singularity in subsonic sectors near the triple point conserves up to transition to the regular reflection.     

Propagation of shock waves in free space

The problem of the exit of a shock wave from an axisymmetric channel and its propagation in a free space occupied by an ideal gas is examined. This problem has been studied earlier in [1], in which the shock wave front was considered planar, as well as in [2], in which the wave front was regarded as a surface of an ellipsoid of revolution. The solutions obtained in these studies assumed the presence of two regions in the wave-front surface: the region of the original shock wave and a region stemming from the decomposition of an infinitesimally thin annular discontinuity of the gas parameters, with the wave intensity over the front surface in each region being considered constant, i.e., the wave character of the process over the front was not considered. In this study a solution will be achieved by the method of characteristics [3–5] of the equations of motion of the shock-wave front, as obtained in [6, 7]. Flow fields are determined for the region immediately adjacent to the shock-wave front for a wide range of shock-wave Mach numbers M _a =1.6–20.0 for = 1.4. On the basis of the data obtained, by introduction of variables connected with the length of the undisturbed zone, as calculated from the channel cross-section along the x axis, together with the pressure transition at the wave front, relationships are proposed which approach self-similarity.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 163–166, September–October, 1971.In conclusion, the author thanks S. S. Semenov for his valuable advice on this study.     

Moving singularities in elasto-diffusive solids with applications to fracture propagation

Solutions are derived for steady-state motion of a singularity class, which includes point sources and dislocations, through a medium in which the elastic stress-field can evolve with time due to the diffusion of an internal second-phase species, such as a pore-fluids and lattice impurity concentrations in a crystalline solid, or the transfer of heat. The technique is to integrate the known influence functions for a stationary singularity. Attention is focused on the most tractable aspect, namely the stress field on the trajectory of motion: this suffices for simulation of growing shear and tensile fractures (e.g. in a porous fluid-saturated solid). Continuous densities of fluid sources and point discontinuities (dislocations) are suitably distributed (as determined by solving the resulting singular integral equations) to satisfy solid stress and fluid pressure or flow conditions on the fracture surfaces. Alternative methods for finding the complete dislocation influence function are discussed and comparisons with existing source solutions are made. Substantial stabilization effects are found in fracture propagation.     

Superposition of finite-amplitude waves propagating in a principal plane of a deformed Mooney–Rivlin material

Here we consider finite-amplitude wave motions in Mooney–Rivlin elastic materials which are first subjected to a static homogeneous deformation (prestrain). We assume that the time-dependent displacement superimposed on the prestrain is along a principal axis of the prestrain and depends on two spatial variables in the principal plane orthogonal to this axis. Thus all waves considered here are linearly polarized along this axis. After retrieving known results for a single homogeneous plane wave propagating in a principal plane, a superposition of an arbitrary number of sinusoidal homogeneous plane waves is shown to be a solution of the equations of motion. Also, inhomogeneous plane wave solutions with complex wave vector in a principal plane and complex frequency are obtained. Moreover, appropriate superpositions of such inhomogeneous waves are also shown to be solutions. In each case, expressions are obtained for the energy density and energy flux associated with the wave motion.     

Analytical solution of a problem of a collision of two hypersonic gas flows from symmetric sources

An asymptotic solution of the Euler equations that describe stationary interaction of two hypersonic gas flows from two identical spherically symmetric sources and an integral equation determining the shock wave shape are obtained with the use of a modified method of expansion of the sought functions with respect to a small parameter, which is the ratio of gas densities in the incoming flow and behind the shock wave. The solution of this equation near the axis of symmetry allows the shock wave stand-off distance from the contact plane and the radius of its curvature to be found. It is shown that the solution obtained agrees well with the known numerical solutions.     

Periodic oscillation and fine structure of wedge-induced oblique detonation waves

An oblique detonation wave for a Mach 7 inlet flow over a long enough wedge of 30 turning angle is simulated numerically using Euler equation and one-step rection model.The fifth-order WENO scheme is adopted to capture the shock wave.The numerical results show that with the compression of the wedge wall the detonation wave front structure is divided into three sections:the ZND model-like strcuture,single-sided triple point structure and dual-headed triple point strucuture.The first structure is the smooth straight,and the second has the characteristic of the triple points propagating dowanstream only with the same velocity,while the dual-headed triple point structure is very complicated.The detonation waves facing upstream and downstream propagate with different velocities,in which the periodic collisions of the triple points cause the oscillation of the detonation wave front.This oscillation process has temporal and spatial periodicity.In addition,the triple point trace are recorded to obtain different cell structures in three sections.     

Stress waves in pyramids by photoelasticity

The propagation of stress waves in pyramids was studied photoelastically with the application of a laser-photomultiplier tube system and an internal polariscope for recording moving fringes. Dispersion and attenuation of stress waves were considered in a straight bar, a 5-deg pyramid, and a 20-deg pyramid made of Hysol 4290 epoxy plastic. In the straight bar and 5-deg pyramid, longitudinal waves propagate without any dispersion even though the waves attenuate as they progress down the models; in the 20-deg pyramid, however, the dispersion of the stress waves is quite significant. The distributions of the axial and radial stresses and the photoelastic fringe patterns obtained on the 20-deg pyramid show that the stress wave front is spherical with the maximum stress along the central axis of the pyramid. A one-dimensional theory of wave propagation without correction factors in a small-angle infinite cone compares well with the experimental results.