截尾Gauss概率密度函数中待定参数的确定

给出一种计算描述标量湍流脉动的截尾Gauss概率密度函数的待定参数的方法. 通过将关于待定参数的代数方程组表示成适于求解的形式，综合应用牛顿法、牛顿下山法和阻尼牛顿法等迭代算法，并恰当地选取待定参数的迭代初值，获得了在标量平均值及其脉动均方值的各种取值条件下待定参数的相应数值.     

非均质流固耦合介质轴对称动力问题时域解

将地基视为流固两相介质并考虑其非均质成层特性，推导了多层地基动力问题时域解．文中首先建立了一组解耦的两相介质动力控制方程；而后利用Ｌａｐｌａｃｅ－Ｈａｎｋｅｌ变换推导了单层地基象空间初参数解答；再利用初参数法及传递矩阵技术导出了任意多层地基瞬态解的一般解析算式．本文获得的解答可方便地退化为现有理想弹性介质的解答     

流固耦合介质轴对称动力问题解法的改进

用直接求解常微分方程组解文［１］所得的控制方程，减少了传递矩阵计算工作量，避免了子阵求逆，使问题的求解得到了简化     

On Fokker–Planck Equations for Turbulent Reacting Flows. Part 2. Filter Density Function for Large Eddy Simulation

The application of large eddy simulation (LES) to turbulent reacting flow calculations is faced with several closure problems. Suitable parametrizations for filtered reaction rates for instance are hardly available in general. A way to overcome these problems is investigated here. This is done by extending LES equations for filtered velocities and scalars (mass fractions of species and temperature) to equations that involve subgrid scale (SGS) fluctuations. Such equations are called filter density function (FDF) methods because they determine the FDF, which is essentially the probability density function of SGS variables. The FDF model considered involves only three parameters: C ₀ that controls the generation of velocity fluctuations and two parameters which determine the relaxation of velocity and scalar fluctuations. The consideration of this model may be seen as the analysis of a limiting case: the implications of the most simple equations for the dynamics of SGS fluctuations are investigated in this way. These equations were proved recently by various simulations. Here, the FDF model is used analytically to improve simpler methods. Existing models for the SGS stress tensor in velocity LES equations and the diffusion coefficient in scalar FDF equations are generalized in this way. The advantages of these models compared to existing ones are pointed out. These investigations provide further evidence for the suitability of the FDF model considered and they provide its parameters. A theoretical value C ₀ = 19/12 is derived, which agrees very well with the results of direct numerical simulation. This estimate implies the same value for the universal Kolmogorov constant of the energy spectrum, which is consistent with the results of many measurements. The other two model parameters can be obtained then by dynamic procedures. Therefore, the closure problems of LES equations are overcome in this way such that adjustable parameters are not involved.     

Single-mode response and control of a hinged–hinged flexible beam

The problem of suppressing the vibrations of a hinged–hinged flexible beam that is subjected to primary and principal parametric excitations is tackled. Different control laws are proposed, and saturation phenomenon is investigated to suppress the vibrations of the system. The dynamics of the beam are modeled with a second-order nonlinear ordinary-differential equation. The method of multiple scales is used to derive two-first ordinary differential equations that govern the time variation of the amplitude and phase of the response. These equations are used to determine the steady-state responses and their stability. The results of perturbation solution have been verified through numerical simulations, where different effects of the system parameters on the steady-state amplitude and on saturation phenomena at resonance have been reported.     

Exact solutions for axisymmetric flexural free vibrations of inhomogeneous circular Mindlin plates with variable thickness

Circular plates with radially varying thickness, stiffness, and density are widely used for the structural optimization in engineering. The axisymmetric flexural free vibration of such plates, governed by coupled differential equations with variable coefficients by use of the Mindlin plate theory, is very difficult to be studied analytically.In this paper, a novel analytical method is proposed to reduce such governing equations for circular plates to a pair of uncoupled and easily solvable differential equations of the Sturm-Liouville type. There are two important parameters in the reduced equations.One describes the radial variations of the translational inertia and flexural rigidity with the consideration of the effect of Poisson's ratio. The other reflects the comprehensive effect of the rotatory inertia and shear deformation. The Heun-type equations, recently well-known in physics, are introduced here to solve the flexural free vibration of circular plates analytically, and two basic differential formulae for the local Heun-type functions are discovered for the first time, which will be of great value in enriching the theory of Heun-type differential equations.     

Parametric vibrations of flexible hoses excited by a pulsating fluid flow,Part I: Modelling,solution method and simulation

The article presents an analysis of a model describing lateral vibrations of a pipe induced by fluid flow velocity pulsation. The motion has been described with a set of two non-linear partial differential equations with periodically variable coefficients. In the analysis Galerkin method has been applied using orthogonal polynomials as shape function. To determine instability regions Floquet theory has been employed. The effect of selected parameters on parametric resonance ranges and regions of increased vibration level has been investigated. The character and form of vibrations have been investigated indicating the possibility of excitation of sub-harmonic and quasi-periodic vibrations in the combination resonance ranges.     

板的非线性热弹耦合振动（Ⅰ）：近似解析解

本文以文[2,3,4]为基础,导出了板的热弹耦合非线性振动控制方程,在采用Galerkin法离散化以后,按各个变量性质分别用多尺度法或正则摄动法求得近似解析解。籍此可揭示系统各参数对非线性热弹耦合振动影响的机理和作出必要的近似计算,对工程实际具有较大的参考价值。     

Dynamic Response of a Multilayered Poroelastic Half-Space to Harmonic Surface Tractions

The propagator matrix method is developed to study the dynamic response of a multilayered poroelastic half-space to time-harmonic surface tractions. In a cylindrical coordinate system, a method of displacement potentials is applied first to decouple the Biot’s wave equations into four scalar Helmholtz equations, and then, general solutions to those equations are obtained. After that, the propagator matrix method and the vector surface harmonics are employed to derive the solutions for a multilayered poroelastic half-space subjected to surface tractions. It is known that the original propagator algorithm has the loss-of-precision problem when the waves become evanescent. At present, an orthogonalization procedure is inserted into the matrix propagation loop to avoid the numerical difficulty of the original propagator algorithm. Finally, a high-order adaptive integration method with continued fraction expansions for accelerating the convergence of the truncated integral is adopted to numerically evaluate the integral solutions expressed in terms of semi-infinite Hankel-type integrals with respect to horizontal wavenumber. Furthermore, to validate the present approach, the response of a uniform poroelastic half-space is examined using the formulation proposed in this article. It is shown that the numerical results computed with this approach agree well with those computed with the analytical solution of a uniform half-space.     

饱和土中的任意形状孔洞对弹性波的散射

根据Biot波动理论建立了求解饱和土中任意形状孔洞对弹性波散射问题的复变函数方法.首先通过引入位移势函数把稳态条件下的Biot波动方程解耦为势函数所满足的Helmholtz方程.利用分离变量方法即得到Helmholtz方程完备的通解.根据所得位移势函数的通解,可得骨架位移、流体相对骨架的位移、应力和孔压的表达式.通过保角变换方法,把物理平面上的孔洞映射到像平面上单位圆.利用土骨架和流体的边界条件,即可确定波函数展开式中的未知系数.给出了一些数值结果.     

White in Time Scalar Advection Model as a Tool for Solving Joint Composition PDF Equations

A rapidly decorrelating velocity field model is used to derive stochastic partial differential equations (SPDE) allowing one to compute the modeled one-point joint probability density function of turbulent reactive scalars. Those SPDEs are shown to be hyperbolic advection/reaction equations. They are dealt with in a generalized sense, so that discontinuities in the scalar fields can be treated. The Eulerian Monte Carlo (EMC) method thus defined is coupled with a RANS solver and applied to the computation of a turbulent premixed methane flame over a backward facing step.     

Stress field of a coated arbitrary shape inclusion

This paper presents an effective method for the plane problem of a coated inclusion of arbitrary shape embedded in an isotropic matrix subjected to uniform stresses at infinity. Based on the complex variable method combined with the expansion of Faber series and Laurent series, the complex potentials in the matrix, the coating and the arbitrary shape inclusion are given in the form of series with unknown coefficients. The stress and displacement continuous conditions on the interfaces are then used to produce a set of linear equations containing all the coefficients. Through solving these linear equations, the complex potentials are finally obtained in the three phases. Additionally, numerical results are presented and graphically shown to investigate the influence of inclusion geometry and coating on the stress distribution along the interfaces for the cases of a coated elliptic, square and triangle inclusions, respectively. It is found that the coating has little effects on the interface stress for a hard inclusion, while it impacts greatly for a soft inclusion. Especially, it is also found that the stresses show the nature of intense fluctuations near the corner of the triangle inclusion, since the inclusion in this case is similar to a wedge.     

Dynamics of a rigid cylinder near a plane boundary in the radiation field of an acoustic wave

An approach is described for investigation of the interaction between a rigid body and a viscous fluid boundary under acoustic wave propagation. The influence of the liquid on the rigid body is determined as a mean force, which is a constant in the time component of the hydrodynamic force. This enables the use of a previously developed technique for calculation of pressure in a compressible viscous liquid. The technique takes into account the second-order terms with respect to the wave field parameters and is based on investigation of a system of initially nonlinear hydromechanics equations that can be simplified with respect to the wave motion parameters of the liquid. It has proven possible to retain the second-order terms for determination of stresses in the liquid without having to solve the system of nonlinear equations. The stresses can be expressed in terms of parameters found in the solution of the linearized equations of the compressible viscous liquid. In this way, the solution of linearized equations is expressed in terms of a scalar and vector potentials. The problem statement is derived for a rigid cylinder located near a rigid flat wall under the effects of a wave propagating perpendicular to the wall. The solution for this particular example is obtained.     

Interacting elliptic inclusions by the method of complex potentials

An accurate series solution has been obtained for a piece-homogeneous elastic plane containing a finite array of non-overlapping elliptic inclusions of arbitrary size, aspect ratio, location and elastic properties. The method combines standard Muskhelishvili’s representation of general solution in terms of complex potentials with the superposition principle and newly derived re-expansion formulae to obtain a complete solution of the many-inclusion problem. By exact satisfaction of all the interface conditions, a primary boundary-value problem stated on a complicated heterogeneous domain has been reduced to an ordinary well-posed set of linear algebraic equations. A properly chosen form of potentials provides a remarkably simple form of solution and thus an efficient computational algorithm. The theory developed is rather general and can be applied to solve a variety of composite mechanics problems. The advanced models of composite involving up to several hundred inclusions and providing an accurate account for the microstructure statistics and fiber–fiber interactions can be considered in this way. The numerical examples are given showing high accuracy and numerical efficiency of the method developed and disclosing the way and extent to which the selected structural parameters influence the stress concentration at the matrix–inclusion interface.     

Interaction of a plane harmonic wave with a thin rigid inclusion of the shape of a cylindrical shell

The paper presents the solution of the problem of determining the stress state in an elastic matrix containing a rigid inclusion of the shape of a thin cylindrical shell. It is assumed that harmonic vibrations occur in the matrix under the conditions of axial symmetry (the symmetry axis is the inclusion axis) and the conditions of full adhesion between the inclusion and the matrix are satisfied. The vibrations are caused by the propagation of a plane wave whose front is perpendicular to the inclusion axis. The solution method is based on representing the displacements in the matrix as discontinuous solutions of the equations of axisymmetric oscillations of an elastic medium with unknown stress jumps on the inclusion surface. The realization of the boundary conditions for these jumps leads to a system of integral equations. Its solution is constructed numerically by the mechanical quadrature method with the use of special quadrature formulas for specific integrals. It is numerically investigated how the ratio of the inclusion geometric dimensions and the propagating wave frequency affect the stress concentration near the inclusion.     

Finite amplitude vibrations of an orthotropic circular plate with an isotropic core

The free and forced non-linear vibrations of a fixed orthotropic circular plate, with a concentric core of isotropic material, are studied. Existence of harmonic vibrations is assumed and thus the time variable is eliminated by a Ritz-Kantorovich method. Hence, the governing non-linear partial equations for the axisymmetric vibration of the composite circular plate are reduced to a set of ordinary differential equations which form a non-linear eigen-value problem. Solutions are obtained by utilizing the related initial-value problems in conjunction with Newton's integration method. The results reveal the effects of finite amplitude and anisotropy of materials upon the dynamic responses. Further, the method developed in this paper, which is used to solve the title problem, is one of some generality. It can be applied to many differential eigenvalue problems with piecewise continuous functions.     

A numerical treatment of the periodic solutions of non-linear vibration systems

Direct numerical integration can be used to find the periodicsolutions for the equations of motion of nonlinear vibrationsystems.The initial conditions are iterated so that theycoincide With the terminal conditions.The time interval ofthe integration(i.e.,the period)and certain parameters ofthe equations of motion can be included in the iterations.Theintegration method has a variable stoplength.This Sbooting method can produce periodic solutions witha shorter computex time.The only error occurs in the numeri-cal integration and it can therefore be estimated and madesmall enough.Using this method one can treat a variety ofvibration problems.such as free conservative.forced.para-meter-excited and self-sustained vibrations with one or se-veral degrees-of-freedom.Unstable solutions and those Whichare sensitive to parameter Changes can also be calculated.Thestability of the solutions is investigated based on the thecryof differential equations with periodic coefficients.The ex-trapolation method and the proc     

Homogenization of scalar wave equations with hysteresis

The paper deals with a scalar wave equation of the form where is a Prandtl–Ishlinskii operator and are given functions. This equation describes longitudinal vibrations of an elastoplastic rod. The mass density and the Prandtl–Ishlinskii distribution function are allowed to depend on the space variable x. We prove existence, uniqueness and regularity of solution to a corresponding initial-boundary value problem. The system is then homogenized by considering a sequence of equations of the above type with spatially periodic data and , where the spatial period tends to 0. We identify the homogenized limits and and prove the convergence of solutions to the solution of the homogenized equation. Received June 17, 1999     

Modeling of properties and motions of an inhomogeneous one-dimensional continuum of a complicated Cosserat-type microstructure

We consider an approach to modeling the properties of the one-dimensional Cosserat continuum [1] by using the mechanical modeling method proposed by Il’yushin in [2] and applied in [3]. In this method, elements (blocks, cells) of special form are used to develop a discrete model of the structure so that the average properties of the model reproduced the properties of the continuum under study. The rigged rod model, which is an elastic structure in the form of a thin rod with massive inclusions (pulleys) fixed by elastic hinges on its elastic line and connected by elastic belt transmissions, is taken to be the original discrete model of the Cosserat continuum. The complete system of equations describing the mechanical properties and the dynamical equilibrium of the rigged rod in arbitrary plane motions is derived. These equations are averaged in the case of a sufficiently smooth variation in the parameters of motion along the rod (the long-wave approximation). It was found that the average equations exactly coincide with the equations for the one-dimensional Cosserat medium [1] and, in some specific cases, with the classical equations of motion of an elastic rod [4–6]. We study the plane motions of the one-dimensional continuum model thus constructed. The equations characterizing the continuum properties and motions are linearized by using several assumptions that the kinematic parameters are small. We solve the problem of natural vibrations with homogeneous boundary conditions and establish that each value of the parameter distinguishing the natural vibration modes is associated with exactly two distinct vibration mode shapes (in the same mode), each of which has its own frequency value.