Historical sense in the historians of the theory of elasticity

This paper studies the historical and epistemological approaches followed by historians of the theory of elasticity. The authors chosen are Todhunter and Pearson, Love, Timoshenko, and Truesdell, who are perhaps the most important historians in English language since the second half of the nineteenth century. It is shown that the historical purpose has changed from one historian to the next, ranging from history as a compilation of what is known to criticism of the evolution of elasticity. The path of progress is investigated and again there are significant differences between historians although overall they have a sense of confluence towards some form of truth. The relation between history and theory is identified as a recurrent theme. And the epistemology in each author, although not explicitly stated, is seen to change following the main trends in their times or before. Some possible avenues of research are mentioned towards the end.     

Stress concentration in the moment theory of elasticity

A multimoment theory of elasticity of inhomogeneous media is used to develop algorithms for determining a sequence of essentially independent fundamental parameters associated with the structure, composition, and characteristics of the medium components. The general form of stress concentration factors is established by analyzing specific models of media. It is found that the mean tangential stresses in a fibrous medium with a transverse stress state are asymmetric. As an example, the methods are used to solve specific problems for plates and shells __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 1, pp. 66–76, January 2007.     

A conservative and shape-preserving semi-Lagrangian method for the solution of the shallow water equations

Semi-Lagrangian methods are now perhaps the most widely researched algorithms in connection with atmospheric flow simulation codes. In order to investigate their applicability to hydraulic problems, cubic Hermite polynomials are used as the interpolant technique. The main advantage of such an approach is the use of information from only two points. The derivatives are calculated and limited so as to produce a shape-preserving solution. The lack of conservation of semi-Lagrangian methods, however, is widely regarded as a serious disadvantage for hydraulic studies, where non-linear problems in which shocks may develop are often encountered. In this work we describe how to make the scheme conservative using an FCT approach. The method proposed does not guarantee an unconditional shock-capturing ability but is able to correctly reproduce the discontinuous flows common in open channel simulation without any shock-fitting algorithm. It is a cheap way to improve existing 1D semi-Lagrangian codes and allows stable calculations beyond the usual CFL limits. A basic semi-Lagrangian method is presented that provides excellent results for a linear problem: the new techniques allow us to tackle non-linear cases without unduly degrading the accuracy for the simpler problems. Two one-dimensional hydraulic problems are used as test cases, water hammer and dam break. In the latter case, because of the non-linearity, special care is needed with the low-order solution and we show the advantages of using Leveque's large-time step version of Roe's scheme for this purpose.     

一种基于能量法则的爆炸成形弹丸速度的工程计算方法

为了准确预报爆炸成形弹丸的速度,提出了以药型罩中所有内部微元都变处于流动状态的时刻作为爆炸成形弹丸成形结束的标志,并因此得出药型罩自身变形时所消耗的变形能应与使整个药型罩都处于流动状态的熔化能相等。在此基础上利用能量守恒原理,建立了爆炸成形弹丸速度的计算公式,并通过试验和数值模拟验证。实验和计算结果表明,本文中所建立的爆炸成形弹丸速度的计算公式具有较好的工程实用价值。     

Stabilization of Unstable Limit Cycles in Systems with Limited Controllability: Expanding the Basin of Convergence of OGY-Type Controllers

The underlying geometric structure of the standard OGY control scheme is analyzed. Some of the main mechanisms that under certain conditions lead to failure of the control algorithm are revealed. The limited controllability available in a system is investigated and it is shown that it may lead to serious problems that will significantly enlarge the state space region of failure of the standard OGY controller. A minimal distance algorithm is analyzed and shown to be, for some problems, more advantageous than the standard OGY technique. Nevertheless, for a broad category of problems, the minimal distance scheme is also shown to fail. As a solution for these problems, two new techniques are proposed: the penalized minimal distance and the multi-step OGY-type scheme. The standard OGY and minimal distance algorithms are particular cases of the new techniques proposed. Finally, we give a necessary condition that estimates the region of controllability under the multi-step OGY-type control. We demonstrate a significantly improved basin of convergence for the new multi-step OGY-type algorithm.     

The oseen flow of finite clusters of spheres and the wake effect on the drag factor

In this paper, the wake effect on drag factor in the axisymmetric Oseen flow of the finite clusters of equally spaced spheres with same size is studied. Putting the Oseen lets on the centres of all the spheres, the series solution of the problem is obtained. By truncating the infinite series and applying the collocation method to solve a set of the linear algebraic equations, the approximate solution of the Oseen flow of finite clusters of spheres and the drag factor for each sphere are presented. The effect of the sphere number and spacing on the drag factor of each sphere under different Reynolds numbers are calculated and the wake effect as well as the shielding effect and the end effect are revealed. The influence of various parameters on the effects is considered and compared with the corresponding results of the Stokes flow. The convergence of the method is also studied numerically in this paper.     

Determining the Axisymmetric Geometrically Nonlinear Thermoviscoelastoplastic State of Laminated Shells by the Theory of Deformation Along Paths of Small Curvature

A technique is developed for determining the thermoviscoelastoplastic geometrically nonlinear axisymmetric stress–strain state of laminar shells of revolution under loads that induce meridional stress and torsion. The technique is based on the hypotheses of rectilinear element for the whole stack of layers. The relations of the theory of deformations along paths of small curvature are used as equations of state. The solution is reduced to the numerical integration of a system of ordinary differential equations. The technique is tried out by a test example and illustrated by determining the geometrically nonlinear thermoviscoelastoplastic state of a corrugated shell     

On the Trajectory Method for the Reconstruction of Differential Equations from Time Series

This work investigates the trajectory method [1] for thereconstruction of ordinary differential equations (ODEs) from timeseries. The potentials of the method are analyzed for dynamical systemsdescribed by second- and third-order ODEs, focusing in particular on therole of the parameters of the method and on the influence of the qualityof the time series in terms of noise, length and sampling frequency.Typical models are investigated, such as the van der Pol, the linearmechanical, the Duffing and the Rössler equations, resulting in arobust and versatile method which is capable of allowing interestingapplications to experimental cases. The method is then applied to themeasured time series of a nonlinear mechanical oscillator, a typicalvelocity oscillation of the bursting phenomenon in near-wall turbulenceand the averaged annual evolution of rainfall, temperature andstreamflow over a hydrological basin.     

A formula of solution to the integral of rational functions

ＡＦＯＲＭＵＬＡＯＦＳＯＬＵＴＩＯＮＴＯＴＨＥＩＮＴＥＧＲＡＬＯＦＲＡＴＩＯＮＡＬＦＵＮＣＴＩＯＮＳＧｕｉＺｕ－ｈｕａ（桂祖华）（ＴｈｅＤｅｐａｒｔｍｅｎｔｏｆＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓ,ＳｈａｎｇｈａｉＪｉａｏ－ｔｏｎｇＵｎｉｖｅｒｓｉｔ...     

Mathematical modeling of the motion of a portion of helium under pulsed injection over a fixed bed of cenospheres

A mathematical model is constructed and an analytical solution is obtained for the problem of a one-dimensional steady flow of a mixture of different gases with hollow permeable particles. The case of a one-dimensional unsteady flow of such a mixture is analyzed numerically. The numerical solutions are compared with experimental data on the motion of the peak concentration of helium in a fixed bed filled with cenospheres (solid hollow permeable spherical particles). The permeability of cenosphere walls and the drag coeficient of cenospheres in the gas flow are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 3, pp. 92–102, May–June, 2007.     

Existence theorems in the theory of mixtures

This paper is concerned with boundary-value problems of the linear theory for binary mixtures of elastic bodies. First, a counterpart of the Boussinesq-Somigliana-Galerkin solution in classical elastostatics is established and the fundamental solutions in the equilibrium theory of homogeneous and isotropic mixtures are derived. Then, representations of Somigliana type for the displacement fields are presented. The potentials of single layer and double layer are used to reduce the boundary-value problems to singular integral equations. Existence and uniqueness results are established.     

The influence of the viscosity and the magnetic field on reflection and transmission of waves at interface between magneto-viscoelastic materials

The influence of the viscosity on reflection and refraction of plane shear elastic waves in two magnetized semi-infinite media is investigated. The numerical results for the absolute values of the reflection and refraction coefficients and their relative changes for a particular choice of the media are presented graphically. The relative changes of these coefficients are calculated for two special orientations of the magnetic field. It is found that the absolute values of these coefficients are not only functions of the angle of incidence but they are also functions of both the large primary magnetic field and the viscosity of the media. They also vary with the orientation of the magnetic field. Finally, we show that the results of earlier works could be obtained here as particular cases. Finally, this study is regarded an attempt to accommodate magnetic field in visco-elastic media where reflection and refraction of plane shear waves are considered. Results which are obtained in this investigation useful for practical applications or for understanding some aspects of physical acoustic.     

Modeling of isothermal processes in porous materials on the basis of the concept of an ensemble of pores

A class of models of porous media based on the concept of an ensemble of pores with a certain distribution of the main geometric parameters (for example, the pore size) is considered. The cases of pores saturated with single-and two-phase multicomponent liquid mixtures are investigated. The properties of equilibrium states of the mixture are derived from the minimum free energy condition and the transfer laws from the decreasing free energy condition. The hydrodynamic connectivity of the pores is described by two types of kernels: one describes the spatial connectivity and the other the connectivity in an elementary macrovolume. Analytically and numerically, the one-dimensional problems of establishment of a steady-state regime, propagation of a passive admixture, and two-phase flow (an analog of the Buckley-Leverett problem) are investigated. A relationship between the models in question and relaxational filtration models is demonstrated. A simple model of capillary hysteresis related with the non-monotonicity of the pore area to volume ratio function is proposed.     

Buckling of a beam subjected to electromagnetic force and its stabilization by controlling the perturbation of the bifurcation

For a beam subjected to electromagnetic force, magnetoelastic buckling due to the increase of such force is theoretically investigated by taking account of the nonlinearity of the electromagnetic force and the elastic force of the beam. Using Liapunov-Schmidt method and center manifold theory, the equilibrium space, the bifurcation set and the bifurcation diagram are theoretically derived. Also, the effect of the higher modes other than the buckling mode on the mode shape of the postbuckling state is discussed. Furthermore, a control method to stabilize the magnetoelastic buckling is proposed, and the unstable equilibrium state of the beam in the postbuckling state, i.e., the straight position of the beam, is stabilized by controlling the perturbation of the bifurcation.     

On the completeness of Lamb modes

A completeness result for Lamb modes in homogeneous waveguides is proved. The problem is formulated as a linear eigenvalue problem in an appropriate Hilbert space of functions. Orthogonality and biorthogonality relations are given. A detailed spectral analysis which is necessary for the application of a general completeness theorem is presented.     

Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging

I compare application of the method of multiple scales with reconstitution and the generalized method of averaging for determining higher-order approximations of three single-degree-of-freedom systems and a two-degree-of-freedom system. Three implementations of the method of multiple scales are considered, namely, application of the method to the system equations expressed as second-order equations, as first-order equations, and in complex-variable form. I show that all of these methods produce the same modulation equations.I address the problem of determining higher-order approximate solutions of the Duffing equation in the case of primary resonance. I show that the conclusions of Rahman and Burton that the method of multiple scales, the generalized method of averaging, and Lie series and transforms might lead to incorrect results, in that spurious solutions occur and the obtained frequency–response curves bear little resemblance to the actual response, is the result of their using parameter values for which the neglected terms are the same order as the retained terms. I show also that spurious solutions cannot be avoided, in general, in any consistent expansion and their presence does not constitute a limitation of the methods. In particular, I show that, for the Duffing equation, the second-order frequency–response equation does not possess spurious solutions for the case of hardening nonlinearity, but possesses spurious solutions for the case of softening nonlinearity. For sufficiently small nonlinearity, the spurious solutions are far removed from the actual response. But as the strength of the nonlinearity increases, these solutions move closer to the backbone and eventually distort it. This is not a drawback of the perturbation methods but an indication of an application of the analysis for parameter values outside the range of validity of the expansion.Also, I address the problem of obtaining non-Hamiltonian modulation equations in the application of the method of multiple scales to multi-degree-of-freedom Hamiltonian systems written as second-order equations in time and how this problem can be overcome by attacking the state-space form of the governing equations. Moreover, I show that application of a variation of the method of Rahman and Burton to multi-degree-of-freedom systems leads to results that do not agree with those obtained with the generalized method of averaging.Contributed by Prof. R.A. Ibrahim.     

The theoretical study of the efficiency of diffusion deposition of nanoaerosols in the extended range of the Peclet numbers

The efficiencies of the diffusion deposition of nanoaerosols for a single fiber for the models of aerosol filter and wire mesh screen are studied numerically in the extended range of the Peclet number Pe. The rectangular periodic cell model for fluid flow and convective-diffusive transport of small aerosol particles is used. Most of the previous theoretical and experimental studies of single fiber diffusion deposition efficiency were for the case of Pe > 1. The array with uniform square or chess grid of fibers and of a row of circular cylindrical fibers are considered as the filter and wire mesh screen models. The flow and particles transport equations are solved numerically using the Boundary Element Method.The obtained numerical data are used to derive the approximate formulas for the deposition efficiency in the entire range of the Peclet number for the various porosities of the filter medium or distances between fibers in a wire mesh screen. The derived dependencies take into account nonlinearity of the deposition efficiency at the low Peclet numbers. The obtained analytical dependencies compare well with the numerical and experimental data.     

Projective Submodel of the Ovsyannikov Vortex

A submodel of the Ovsyannikov vortex with projective symmetry is studied. Integration of the factor system of the submodel reduces to solving a first-order differential equation which is not solved with respect to the derivative. The properties of the solutions of this equation are studied. It is shown that the submodel describes gas ow with a nonstationary source and a nonstationary sink. The problem of the motion of a gas volume between pistons of cylindrical shapes is studied, and its solution with an invariant shock wave is obtained.__________Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 4, pp. 3–16, July–August, 2005.     

应用突变理论分析Euler压杆的稳定性

建立了Euler压杆问题的能量表达式。对此模型，应用突变理论进行稳定分析，得出了系统全部的分叉集与突变流形。突变流形为一族层层嵌套、互不交叉的抛物线。分析的结论要比有限差分法、有限单元法、大范围非线形分析等数值计算的结果简洁、明了，有助于压杆后屈曲问题的研究。     

Acceleration of convergence by shifting the spectrum of implicit finite difference operators associated with the equations of gas dynamics

Eigensystem analysis techniques are applied to finite difference formulations of the Navier-Stokes equations in one dimension. Spectra of the resulting implicit difference operators are computed. The largest eigenvalues are calculated by using a combination of the Frechet derivative of the operators and Arnoldi's method. The accuracy of Arnoldi's method is tested by comparing the rate of convergence of the iterative method with the dominant eigenvalue of the original iteration matrix. On the basis of the pattern of eigenvalue distributions for various flow configurations, a shifting of the implicit operators in question is devised. The idea of shifting is based on the power method of linear algebra and is very simple to implement. This procedure has improved the rates of convergence of CFD codes (developed at NASA Ames Research Center) by 20%–50%. The sensitivity of the computed solution with respect to the shift is also studied. Finally, an adaptive shifting of the spectrum together with Wynn's acceleration algorithm are discussed. It turns out that the shifting process is a preconditioner for Wynn's method.