A fully non‐linear model for three‐dimensional overturning waves over an arbitrary bottom

An accurate three‐dimensional numerical model, applicable to strongly non‐linear waves, is proposed. The model solves fully non‐linear potential flow equations with a free surface using a higher‐order three‐dimensional boundary element method (BEM) and a mixed Eulerian–Lagrangian time updating, based on second‐order explicit Taylor series expansions with adaptive time steps. The model is applicable to non‐linear wave transformations from deep to shallow water over complex bottom topography up to overturning and breaking. Arbitrary waves can be generated in the model, and reflective or absorbing boundary conditions specified on lateral boundaries. In the BEM, boundary geometry and field variables are represented by 16‐node cubic ‘sliding’ quadrilateral elements, providing local inter‐element continuity of the first and second derivatives. Accurate and efficient numerical integrations are developed for these elements. Discretized boundary conditions at intersections (corner/edges) between the free surface or the bottom and lateral boundaries are well‐posed in all cases of mixed boundary conditions. Higher‐order tangential derivatives, required for the time updating, are calculated in a local curvilinear co‐ordinate system, using 25‐node ‘sliding’ fourth‐order quadrilateral elements. Very high accuracy is achieved in the model for mass and energy conservation. No smoothing of the solution is required, but regridding to a higher resolution can be specified at any time over selected areas of the free surface. Applications are presented for the propagation of numerically exact solitary waves. Model properties of accuracy and convergence with a refined spatio‐temporal discretization are assessed by propagating such a wave over constant depth. The shoaling of solitary waves up to overturning is then calculated over a 1:15 plane slope, and results show good agreement with a two‐dimensional solution proposed earlier. Finally, three‐dimensional overturning waves are generated over a 1:15 sloping bottom having a ridge in the middle, thus focusing wave energy. The node regridding method is used to refine the discretization around the overturning wave. Convergence of the solution with grid size is also verified for this case. Copyright © 2001 John Wiley & Sons, Ltd.     

Characterization of shock accelerometers using davies bar and strain-gages

This paper proposes a novel method for evaluating the dynamic characteristics of shock accelerometers under high acceleration levels and a wide frequency bandwidth. High accelerations of 10³∼10⁵m/s² can be generated by the reflection of an elastic wave pulse propagating in a metal bar known as the Davies bar. The elastic wave pulse is produced by the collision of a projectile against one end of the bar, and is detected by straingages. The accelerometer to be characterized is attached to the other end of the bar. The one-dimensional theory of elastic waves enables the derivation of an input acceleration to the accelerometer from the measured strain. The dispersion of the elastic waves caused by the lateral inertia of the bar is compensated for by using a two-dimensional analytical solution. This method was validated by an experiment characterizing a piezoelectric-type accelerometer within the frequency band approximately 1 kHz∼70 kHz.     

A numerical method for the determination of weakly non‐hydrostatic non‐linear free surface wave propagation

A numerical method is described that may be used to determine the propagation characteristics of weakly non‐hydrostatic non‐linear free surface waves over a general, bottom topography. In shallow water of constant undisturbed depth, such waves are equivalent to the familiar cnoidal waves characterized by sharp crests and relatively flat troughs. For a certain range of parameters, these propagate without change of form by virtue of the weakly non‐hydrostatic balance in the vertical momentum equation. Effectively, this counters the tendency for the non‐linearity in a purely hydrostatic theory to lead to a continuously deforming surface wave profile. The realistic representation furnished by cnoidal wave theory of free surface waves in the shallow near‐shore zone has led to its utilization in evaluating their propagation characteristics. Nonetheless, the classic analytical theory is inapplicable to the case of wave propagation over a sloping beach or off‐shore sand bar topography. Under these conditions, a local change in form of the surface wave profile is anticipated before the waves break and knowing this is required in order to evaluate fully the propagation process. The efficacy of the numerical method is first demonstrated by comparing the solution for water of constant depth with the evaluation of the analytical solution expressed in terms of the Jacobian elliptic function cn. The general method described in the paper is then illustrated by experiments to determine the change in profile of weakly non‐hydrostatic non‐linear surface waves propagating over bed forms representative of those found in shallow coastal seas. Copyright © 2001 John Wiley & Sons, Ltd.     

Dynamic response of an axially loaded rigid sphere embedded in a saturated poroelastic medium

In this paper, an analytical solution for the response of a rigid sphere embedded in a full space poroelastic medium subjected to a dynamic lateral load is derived. The solution is obtained using Biots theory for acoustic waves. In this solution, the displacements of the solid skeleton and the pore pressure are expressed in terms of three scalar potentials. These potentials correspond to the wave velocities of the slow and fast compressional waves and to the shear wave. The governing equation for the dynamic motion is expressed in the frequency domain using Fourier transformation. Different boundary and load conditions were investigated. Curves showing variation in the fluid pressure and solid displacements with the loads frequency were plotted in non-dimensional forms.     

Stability of surfaces of strong discontinuity in gases

The existence of solutions with surfaces of strong discontinuity is one of the principal features of the continua whose motions are described by systems of differential equations of hyperbolic type. Shock waves in gas dynamics, magnetohydrodynamics and in solids, detonation waves and combustion fronts, contact discontinuities, etc. are well-known examples of these surfaces. The discontinuities are usually investigated in accordance with the following scheme: 1) derivation of the boundary conditions on the discontinuity from the input system of differential equations in integral form; 2) verification of the fulfilment of the evolution conditions; 3) solution of the problem of the discontinuity structure and, when the occasion requires, obtaining supplementary boundary conditions; 4) investigation of the stability of the discontinuity. Only after obtaining positive results in all fours stages can we assert that the existence of the discontinuity is theoretically justified and that it can be used for constructing the solutions of particular boundary value problems. In the present paper attention will be concentrated on the problem of the stability of discontinuities, all the material, with the exception of the general results of Sec.1, being concerned with gas media and relating to discontinuities on whose surface the normal mass flow is nonzero. Having no way of exploring all the aspects of the problem of the stability of discontinuities in the same detail within the limited context of this paper, the authors hope to demonstrate the most general ideas and approaches which could subsequently be used to investigate the stability of discontinuities in various particular models of continua.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–22, March–April, 1996.     

Weak shock waves and shear bands in compressible,inextensible thermoelastic solids

A study of weak shock waves propagating into a solid, which is compressible but temperature-dependent extensible in a specified direction is presented. The inextensible solid is also considered. The constitutive equations of constrained thermoelastic material are written as the summation of constrained and unconstrained counterparts of the relevant quantities. The equation of motion of weak shock waves, which is recovered by the theory of singular surfaces, reduces to an eigenvalue problem. The solution of this eigenvalue problem yields the speeds of propagation of weak shock waves. In the case of an undeformed solid, the speeds of these waves are explicitly expressed. Additionally, a discussion on the ductility limits of constrained thermoelastic material subjected to the uniaxial and biaxial extensions is presented.     

曲面曲率对Rayleigh波传播特性的影响

对任意形状的均匀各向同性线弹性曲面物体,用 WKB～（1）方法求解了沿曲面传播的Rayleigh表面波的运动微分方程,同时考虑了波传播方向及其垂直方向曲面曲率对波的穿透性的影, 所获波动方程的势函数解答表明,在一般情况下垂直波传播方向的曲面曲率对波的穿透深度的影响是不容忽视的.进而以同种介质平面表面情况下的Rayleigh面波的传播特性为基准,给出了曲面曲率引起波数或波速变化的解析表达式.通过理论分析和数值算例,描述了曲面上Rayleigh面波传播行为的一些基本特征.     

Experimental investigation of three-dimensional supersonic flow of a gas in the interference region of intersecting surfaces

A contemporary high-speed aircraft represents a complex three-dimensional configuration, where supersonic gas flow is accompanied by numerous local flow interaction zones, in particular, near the intersection of different surfaces. Such a flow is characterized by three-dimensional systems of shock and expansion waves, and close to the surfaces one finds interaction of boundary layers and, above all, interaction of shock waves with the boundary layer. In general, the angular configurations are formed by intersection or contact of nonplanar surfaces with swept-back or blunted leading edges. This makes it practically impossible to obtain a rigorous theoretiical solution to the problem of gas flow over these surfaces, and presents considerable difficulty in an experimental investigation. It is therefore of interest to study the physical features of gas flow in corner configurations of very simple form [1–3]. The present paper examines the results of an experimental investigation of typical features of symmetric and asymmetric interaction of compressive, expansive, and mixed flows in the interference region of planar surfaces intersecting at an angle of less than 180?.     

弹性曲杆的稳定性问题

本文给出空间任一曲杆在弯扭联合作用下的稳定性问题的一般讨论,并且给出了曲杆某一平衡状态的扰动量所满足的方程组(28)—(36),在适当的边界条件下,这些扰动量的非零解对应于临界状态。文末用这组方程具体讨论了五个实际例子,这些例子有些结果是新的,有些是用新的方法去处理老问题。     

Modelling of the moisture concentration field due to cyclical hygrothermal conditions in thick laminated pipes

The paper presents a new method to calculate the moisture concentration field induced by cyclical environmental conditions in thick laminated pipes. The solution which is obtained is composed of a transient solution over the interior of the pipe wall and a fluctuating solution within two thin regions, close to the inner and outer lateral surfaces of the pipe wall. The thickness of these two regions is depending on both materials and frequency conditions. The transient solution is determined by using an analytical method based on the solving in average of the field equation. The fluctuating solution is derived from a finite difference scheme. It is shown that after some period of time the transient solution tends towards a permanent time independent solution. In that case, the fluctuating solution becomes a periodic solution which is conditioned by the cyclical boundary conditions. Finally, the effect of particular cyclical conditions on the moisture concentration in thick wall pipes will be tackled.     

Analytical approximate 3D solution for the longitudinally vibrating cylinder

Summary The paper presents a three-dimensional approximate solution of a longitudinally vibrating cylinder. It is based on the frequency equation for free waves travelling along an infinitely long cylinder which follows from the radial boundary condition on the lateral surface of the cylinder. An equivalent longitudinal traction is calculated by integration of the normal stresses over the end cross section. The full solution of the longitudinally vibrating rod adapted to given axial boundary conditions is derived. The solution is compared to numerical results, and good agreement is obtained.     

对称表面动载下有限块体的破裂过程

实验表明,受对称布置的表面爆炸载荷的脆性块体,会在通过载荷作用点连线的某一平面上破裂。为解释这一现象并进一步探讨利用类似方法破碎岩块的可能性,本文介绍了应用弹性动力学理论建立的计算模型以及所进行的动光弹和混凝土块模型试验。研究结果表明:在所述载荷条件下,应力波特别是自由界面反射的拉伸波对块体的破碎起重要作用;由于多向拉应力集中,在块体加载端面梭边中部最易产生径向裂隙,此后在侧面反射波的作用下裂隙延伸,与中部可能存在的裂纹贯通,从而将块体破开;对称加载使应力波在块体内多次叠加有利于块体充分破碎。     

Three-dimensional solutions for general anisotropy

The Stroh formalism is extended to provide a new class of three-dimensional solutions for the generally anisotropic elastic material that have polynomial dependence on x₃, but which have quite general form in x₁,x₂. The solutions are obtained by a sequence of partial integrations with respect to x₃, starting from Stroh's two-dimensional solution. At each stage, certain special functions have to be introduced in order to satisfy the equilibrium equation. The method provides a general analytical technique for the solution of the problem of the prismatic bar with tractions or displacements prescribed on its lateral surfaces. It also provides a particularly efficient solution for three-dimensional boundary-value problems for the half-space. The method is illustrated by the example of a half-space loaded by a linearly varying line force.     

An optimisation method for separating and rebuilding one-dimensional dispersive waves from multi-point measurements. Application to elastic or viscoelastic bars

When using a classical SHPB (split Hopkinson pressure bar) set-up, the useful measuring time is limited by the length of the bars, so that the maximum strain which can be measured in material testing applications is also limited. In this paper, a new method with no time limits is presented for measuring the force and displacement at any station on a bar from strain or velocity measurements performed at various places on the bar. The method takes the wave dispersion into account, as must inevitably be done when making long time measurements. It can be applied to one-dimensional and single-mode waves of all kinds propagating through a medium (flexural waves in beams, acoustic waves in wave guides, etc.). With bars of usual sizes, the measuring time can be up to 50 times longer than the time available with classical methods. An analysis of the sensitivity of the results to the accuracy of the experimental data and to the quality of the wave propagation modelling was also carried out. Experimental results are given which show the efficiency of the method.     

Waves in a viscoelastic bar surrounded by soils under smooth loading

The results of numerical solution of the nonstationary wave problem of longitudinal wave propagation in the bar-soil system are given. The barmaterial and the soil are assumed to be linearly viscoelastic (standard linear body). A nonlinear interaction condition is satisfied on the bar-soil contact boundary. The laws of propagation of longitudinal waves and variations in the longitudinal stresses in the bar are determined depending on the mechanical characteristics of the soil and the bar and on the interaction parameters.     

Determination of hourglass coefficients in the theory of a Cosserat point for nonlinear elastic beams

The theory of a Cosserat point has recently been used [Int. J. Solids Struct. 38 (2001) 4395] to formulate the numerical solution of problems of nonlinear elastic beams. In that theory the constitutive equations for inhomogeneous elastic deformations included undetermined constants associated with hourglass modes which can occur due to nonuniform cross-sectional extension and nonuniform torsion. The objective of this paper is to determine these hourglass coefficients by matching exact solutions of pure bending and pure torsion applied in different directions on each of the surfaces of the element. It is shown that the resulting constitutive equations in the Cosserat theory do not exhibit unphysical stiffness increases due to thinness of the beam, mesh refinement or incompressibility that are present in the associated Bubnov–Galerkin formulation. Also, example problems of a bar hanging under its own weight and a bar attached to a spinning rigid hub are analyzed.     

Propagation of nonlinear waves along a viscoelastic bar

Various types of nonlinear waves propagating along a viscoelastic bar are considered. The rheological equation of state has strong physical and geometric nonlinearities, and nonisothermal effects are included. Both weak (isentropic) and shock waves of loading and unloading are investigated. It is shown that, for certain rubber-like materials, stable shock waves of extension can exist along with the shock waves of compression at very large strains. We then consider the strike of a viscoelastic bar of finite length against a rigid obstacle. Numerical solutions to this problem illustrate the influence of stress relaxation on nonlinear wave processes. A model for sticking and bouncing off is formulated and the mass-averaged velocity of the bar at the moment when it bounces off the obstacle is calculated.     

Elastic Waves in Swelling Porous Media

Time harmonic waves in a swelling porous elastic medium of infinite extent and consisting of solid, liquid and gas phases have been studied. Employing Eringen’s theory of swelling porous media, it has been shown that there exist three dilatational and two shear waves propagating with distinct velocities. The velocities of these waves are found to be frequency dependent and complex valued, showing that the waves are attenuating in nature. Here, the appearance of an additional shear wave is new and arises due to swelling phenomena of the medium, which disappears in the absence of swelling. The reflection phenomenon of an incident dilatational wave from a stress-free plane boundary of a porous elastic half-space has been investigated for two types of boundary surfaces: (i) surface having open pores and (ii) surface having sealed pores. Using appropriate boundary conditions for these boundary surfaces, the equations giving the reflection coefficients corresponding to various reflected waves are presented. Numerical computations are performed for a specific model consisting of sandstone, water and carbon dioxide as solid, liquid and gas phases, respectively, of the porous medium. The variations of phase speeds and their corresponding attenuation coefficients are depicted against frequency parameter for all the existing waves. The variations of reflection coefficients and corresponding energy ratios against the angle of incidence are also computed and depicted graphically. It has been shown that in a limiting case, Eringen’s theory of swelling porous media reduces to Tuncay and Corapcioglu theory of porous media containing two immiscible fluids. The various numerical results under these two theories have been compared graphically.     

Stress concentration near a thin elastic inclusion under interaction with harmonic waves in the case of smooth contact

We solve the problem on the interaction of plane harmonic waves with a thin elastic plate-shaped inclusion. The ambient medium is assumed to be in plane strain. The smooth contact conditions are satisfied on both sides of the inclusion. The bending displacements of the inclusion are determined from the corresponding differential equation. In the statement of boundary conditions for this equation, one should take into account the transverse forces and bending moments applied to the lateral edges of the inclusion, while the boundary conditions are posed on the midplane of the inclusion. Using the discontinuous solution method, we reduce the problem to a system of two singular integral equations, which are solved numerically by the mechanical quadrature method. We obtain approximate formulas for the stress intensity coefficients near the ends of the inclusion and for the transverse forces and moments applied to the inclusion.