基于相位差谱的空间相关非平稳地震动场的模拟

地震动的非平稳特性主要是由其相位差谱决定的,相位差谱与相位导数之间存在线性倍数关系。根据相位导数的显式计算公式,从能量的角度解释了相位导数的均值大致决定了时程峰值的发生时刻,相位导数的方差决定了强震段的持续时间。在相关地震动场模拟方法中,首次将相位差谱的统计模型引入空间相关非平稳地震动场的模拟方法之中,利用快速傅立叶变换技术生成地震动场。表示地震动随机特性的随机相位谱利用相位差谱的统计模型生成,生成的地震动场不仅具有空间相关性,而且在时域、频域内均具有非平稳特性。     

油水两相渗流压力脉动的实验及数值分析

本文用高精度的压差传感器和高速度的数据采集设备同时测量出油水两相流体在多孔介质流动时产生的较大压力降和微小压力脉动,应用自行研制的大型数据处理软件对恒流速水驱油的两相渗流压力脉动实验数据进行了分析,发现不同阶段压力脉动具有明显的频谱特性和时间相关特性的不同,在第三阶段（油为主,水增加阶段）谱能增加最大,时间正相关程度最强。     

Approximate method of solving the problem of characteristic oscillations of a liquid in cavities of revolution

The problem of the characteristic oscillations of a liquid in axisymmetric cavities of rotation has been fairly fully studied [1–5], its solution in the general case being found by the variational method. Analysis of numerical results using the variational method shows that to achieve acceptable accuracy it is necessary to retain an appreciable number of coordinate functions, which entails the solution of a matrix eigenvalue problem of high order, this applying especially to the case when it is necessary to determine several eigenfrequencies and the shapes of the oscillations. In the present paper, a method proposed earlier by Shmakov [6] is developed, the velocity potential being sought in the form of a sun of two potentials. The first (base) potential is a solution to the problem of the characteristic oscillations of a liquid in a cavity whose free surface coincides with the free surface of the original cavity, and the second (correcting) potential is chosen in the form of a system of harmonic functions, this system being complete and orthogonal on the wetted surface of the cavity. Cavities of revolution are analyzed as examples, and a detailed investigation of numerical results is made for a spherical cavity. The numerical analysis shows that a sufficiently accurate result in the determination of a frequency is obtained if one term of the base problem is retained and only the correcting potential is used to make this more accurate. As a result, it is only necessary to solve an algebraic equation of first degree in the square of the frequency.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 3–8, September–October, 1983.     

Force reconstruction based on the union of singular spectral analysis and regularization method

Abstract

To overcome the problem of some ill-posed inverse problem of force reconstruction, which is caused by the noise in the measured responses and small singular values of the structure, a technology of force reconstruction based on a hybrid method of singular spectral analysis (SSA) and the Landweber regularization method is proposed in this study for the first time. The SSA is used to filter the structural response before using Landweber regularization. A new choice method of phase space reconstruction dimension is theoretically proposed, and the minimum embedding dimension is determined by the concept of optimizing difference spectrum theory. The feasibility of this method was demonstrated through three kinds of force reconstructions. The numerical simulation results and an acoustic vibration experiment demonstrated that the proposed method is more effective than the traditional method.

Communicated by Wei-Chau Xie.     

A NOVEL SEMI-ANALYTICAL METHOD FOR SOLVING ACOUSTIC RADIATION FROM LONGITUDINALLY STIFFENED INFINITE NON-CIRCULAR CYLINDRICAL SHELLS IN WATER

Based on the extended homogeneous capacity high precision integration method and the spectrum method of virtual boundary with a complex radius vector, a novel semi-analytical method, which has satisfactory computation effectiveness and precision, is presented for solving the acoustic radiation from a submerged infinite non-circular cylindrical shell stiffened by longitudinal ribs by means of the Fourier integral transformation and stationary phase method. In this work,besides the normal interacting force, which is commonly adopted by some researchers, the other interacting forces and moments between the longitudinal ribs and the non-circular cylindrical shell are considered at the same time. The effects of the number and the size of the cross-section of longitudinal ribs on the characteristics of acoustic radiation are investigated. Numerical results show that the method proposed is more efficient than the existing mixed FE-BE method.     

扩散方程单内点精细积分法与差分法比较研究

一维扩散方程初值问题可以用全域或子域精细积分求解。子域积分可以采用不同数量的内点，单内点是其最简单的情况。当单内点精细积分中的传递函数即指数函数用其泰勒展开式的一阶近似来替代时，精细积分转化为差分方程。本文研究了这一对应关系。各种常见差分格式均找到了对应的单点精细积分格式，并在单点精细积分一般公式中得到了统一表达形式     

Distribution of the velocity profiles of the liquid phase in a gas-liquid flow with small gas contents

Artilce [1] gives the results of measurement of the friction at the wall of a channel under bubble conditions, in a wide range of Reynolds numbers. It is shown that the concept of laminar flow conditions has no meaning when it is applied to the flow of a two-phase mixture, since, even with very small Reynolds numbers, the level of the pulsations of the friction is high, and the spectrum of the pulsations of the friction is continuous. In this case, the mean friction is much greater than the calculated; here the value of the resistance coefficient is not a single-valued function of the Reynolds number. The present article gives the results of measurement of the velocity profiles of the liquid phase, carried out using an electro diffusion method. It is shown that, with Reynolds numbers corresponding to turbulent flow conditions, the profile of the velocity in a two-phase mixture is close to turbulent and does not depend on the gas content.     

Boger流体的松弛谱对非线性粘弹性表征的影响

分析了一种Boger流体的松弛时间谱和该流体的粘弹性表征问题。文中采用带正约束的非线性规划方法得到了包含4个、5个、6个松弛时间的描述此Boger流体线性粘弹特性的时间谱。这些新时间谱比前文报道的含3个松弛时间的谱能更好地描述该流体的储能模量G′。用PSM模型结合这些新松弛时间谱都不能很准确地描述这个Boger流体的第一法向应力差N1的实验数据。作者提出的3r模型描述此Boger流体的N1的结果比PSM模型的计算结果更接近实验值。     

Optimization of the collocation inversion method for the linear viscoelastic homogenization

The effective behaviour of linear viscoelastic heterogeneous material can be derived from the correspondence principle and the inversion of the obtained symbolic homogenized behavior. Various numerical methods were proposed to carry out this inversion. The collocation method, widely used, within this framework rests on a discretization of the characteristic spectrum in a sum of discrete lines for which it is necessary to determine the intensities and the positions by the minimization of the difference between the exact temporal function and its approximation. The classical method is based on a priori choice of the lines positions and on the optimization of their intensities. It is shown here that the combined optimization of the positions and the (positive) intensities lead to a minimization problem under constraints. In the simple case of an incompressible isotropic two-phase material, the assessment of the effective relaxation function with a continuum spectra or made up of discrete lines proves that the proposed method improves the predictions of the classical approach.     

Stochastic modeling for starting-time of phase evolution of random seismic ground motions

In response to the challenge inherent in classical high-dimensional models of random ground motions, a family of simulation methods for non-stationary seismic ground motions was developed previously through employing a wave-group propagation formulation with phase spectrum model built up on the frequency components' starting-time of phase evolution. The present paper aims at extending the formulation to the simulation of non-stationary random seismic ground motions. The ground motion records associated with N—S component of Northridge Earthquake at the type-II site are investigated. The frequency components' starting-time of phase evolution of is identified from the ground motion records, and is proved to admit the Gamma distribution through data fitting. Numerical results indicate that the simulated random ground motion features zero-mean, non-stationary, and non-Gaussian behaviors, and the phase spectrum model with only a few starting-times of phase evolution could come up with a sound contribution to the simulation.     

Coupling Eulerian‐Lagrangian method of air‐particle two‐phase flow with population balance equations to simulate the evolution of vehicle exhaust plume

In this paper, we present a new numerical scheme to describe the dynamic evolution of multiphase polydisperse systems in terms of time, space, and properties by coupling the Eulerian‐Lagrangian method for air‐particle two‐phase flow and population balance equations to describe particle property evolution due to microbehaviors (eg, aggregation, breakage, and growth). This coupling scheme was used to comprehensively simulate the two‐phase flow structure, particle size spectrum, particle number, and volume concentrations. These were characterized by a high‐resolution particle tracking using the Lagrangian approach and the high precision of moments of the particle size spectrum by solving the population balance equation with the quadrature method of moments. The algorithm of the coupling scheme was incorporated into the open source computational fluid dynamics software OpenFOAM to simulate the dynamic evolution of vehicle exhaust plume. The impacts of vehicle velocity, exhaust temperature, and aggregation efficiency on the distribution of auto exhaust particles in space and changes in their properties were analyzed. The results indicate that the particle number concentration, volume concentration, and average diameter of particles in the vehicle exhaust plume could be strongly affected by the plume structure and flow properties.     

求解水下纵向加肋无限长非圆柱壳声辐射问题的一种新的半解析方法

基于齐次扩容精细积分法和复数矢径虚拟边界谱方法，利用Fourier积分变换和稳相法，提出了一种具有较高效率和精度的新的求解水下纵向加肋无限长非圆柱壳声辐射问题的半解析方法．考虑了非圆柱壳和肋骨之间同时存在多种相互作用力和力偶矩，较已往很多学者仅计及法向相互作用力更加符合实际．不仅比较了该文方法和精确解计算纵向加肋圆柱壳在集中点力激励下的声辐射计算结果，同时还研究了肋骨数量、大小以及椭圆柱壳横截面椭圆度对声辐射特性的影响．数值计算结果表明该文方法较已有的混合FE-BE法更为有效．     

云纹／光栅测试技术的几点讨论

本文用信号处理的方法系统地阐述了云纹／光栅测试技术。以面内云纹／光栅测试为例，首先说明了光栅是位移和变形信息的空间载体；位移和应变分别是对载体的相位和频率调制；而云纹现象则是试件栅频谱向零频处平移的结果或说云纹是对光栅载体信号解调的过程；进而讨论了云纹技术与光栅技术的本质区别。其次分析了云纹／光栅测试系统的性能，讨论了条纹分析方法及数字图象处理手段对云纹／光栅系统的测试精度、灵敏度、空间分辨率的影响。     

The minimal error method for the Cauchy problem in linear elasticity. Numerical implementation for two-dimensional homogeneous isotropic linear elasticity

In this paper, yet another iterative procedure, namely the minimal error method (MEM), for solving stably the Cauchy problem in linear elasticity is introduced and investigated. Furthermore, this method is compared with another two iterative algorithms, i.e. the conjugate gradient (CGM) and Landweber–Fridman methods (LFM), previously proposed by Marin et al. [Marin, L., Háo, D.N., Lesnic, D., 2002b. Conjugate gradient-boundary element method for the Cauchy problem in elasticity. Quarterly Journal of Mechanics and Applied Mathematics 55, 227–247] and Marin and Lesnic [Marin, L., Lesnic, D., 2005. Boundary element-Landweber method for the Cauchy problem in linear elasticity. IMA Journal of Applied Mathematics 18, 817–825], respectively, in the case of two-dimensional homogeneous isotropic linear elasticity. The inverse problem analysed in this paper is regularized by providing an efficient stopping criterion that ceases the iterative process in order to retrieve stable numerical solutions. The numerical implementation of the aforementioned iterative algorithms is realized by employing the boundary element method (BEM) for two-dimensional homogeneous isotropic linear elastic materials.     

相位谱线性重构法去除Lamb波频散研究

Lamb波广泛应用于板壳结构的损伤检测,而Lamb波的频散效应使得板结构中损伤散射信号的渡越时间信息难以准确提取,因而影响了阵列波束成形损伤成像算法对损伤成像定位的效果和精度.提出的相位谱线性重构法,在中心频率处通过对相位谱展开,进行线性化处理,有效地去除了激励信号中的频散效应.该方法保证了中心频率下的信号分量不变,为研究阵列波束成形损伤成像算法对损伤精确成像奠定了基础.采用数值模拟和实验手段,验证所提出的相位谱线性重构法去除频散的效果,进而结合阵列波束成形损伤成像方法分析损伤成像定位的效果和精度.结果 表明,采用的频散去除方法能够提高损伤散射信号的信噪比,提升阵列波束成形损伤成像方法对损伤成像定位的效果和精度.     

组装式微流控系统制备双重包裹微液滴的方法研究

介绍了一种组装式微流控系统制备单、双重包裹微液滴的方法。微系统中用三通接头构成T型微流体通道,使得分散相在连续相强烈的剪切力和压力差作用下断裂形成单个微米级液滴。在制备单个微液滴基础上,用毛细管将两个三通接头串联,通过调控三相流量,可产生双重包裹液滴。结合实验结果,分析了流体粘度比对液滴大小的影响,并得出液滴的尺寸与流量比之间的关系式,为制备不同尺寸的液滴提供了参考依据。对制备的样品进行统计分析,结果显示,液滴的多分散性指数均小于3.2%,表明微液滴的高度均匀性。此外,通过调节三相液体的流量不仅可以控制内外层液滴的大小,还可以调节内层包裹液滴的个数。本文提出的制备方法,设备组装拆卸简便,不需表面亲疏水性处理,装置利用率高,产生的单、双重包裹微液滴可满足高通量的测量分析要求。     

Vortex pairing in an axisymmetric jet using two-frequency acoustic forcing at low to moderate strouhal numbers

 Hot-wire measurement and multi-smoke wire flow visualization method are employed to study vortex pairing in the jet column mode under two-frequency forcing with controlled initial phase differences. For the range of 0.3<St _D <0.6, vortex pairing can be easily controlled by means of the fundamental and its subharmonic forcing with varying initial phase differences. As stable vortex pairing dominates, the variation of the subharmonic component with the initial phase difference changes from a sine shape to a cusp-like shape. The harmonics of the subharmonic also show similar trends. The detuning induces the amplitude and phase modulations of the u-signal in the time trace and the sideband growth in the spectra. The u-signal reflects the subharmonic variation with the initial phase difference in its envelope. For 0.6<St _D <0.9, non-pairing advection of vortices due to improper phase difference is sometimes observed under single-frequency forcing. In this case, vortex pairing can be made to occur by the addition of a subharmonic with very small amplitude. As the initial level of this subharmonic is increased, the onset position of vortex pairing moves upstream. In this range, the initial phase difference is not an effective parameter in controlling vortex pairing. Received: 22 May 1997 / Accepted: 16 October 1997     

Non-invasive Mechanical Measurement for Transparent Objects by Digital Holographic Interferometry Based on Iterative Least-Squares Phase Unwrapping

A non-invasive measurement method for transparent objects using digital holographic interferometry based on iterative least-squares phase unwrapping is presented. Holograms are captured and processed by digital method. The phase reconstructed from holograms is wrapped. Continual and real phase is obtained by phase unwrapping. In this paper iterative unwrapping of phase difference is combined with least-squares unwrapping to eliminate errors. The relations between the phase variation of object wave and the thickness deformation and stress of transparent object are given. A polymethyl methacrylate (PMMA) specimen with a hole under uniform tensile force is tested by this method. The experimental results are in accordance with the theoretical ones. That means the method is correct.     

精细积分法应用于地震碰撞力反应谱计算研究

相邻结构的碰撞时程分析是计算地震碰撞力反应谱的基础。结构碰撞时程分析要求采用稳定性好、精度高及计算效率高的数值分析方法。精细积分法将二阶动力微分方程通过增元降阶的方式转换成Hamilton对偶变量体系,得到了动力微分方程的精确解。基于此,本文将精细积分法引入结构碰撞时程分析及地震碰撞力反应谱计算中。在推导精细积分法公式的基础上,在MATLAB环境下编制了结构碰撞时程分析程序和碰撞力反应谱计算程序,并实现了碰撞力反应谱程序的并行化。经算例验证,精细积分法应用于结构碰撞时程分析及地震碰撞力反应谱计算是可行的,程序计算结果准确。     

Investigation of two-dimensional transonic flow of a gas through a turbine cascade by a combined method of integral equations and stabilization

Two numerical methods are basically used to calculate two-dimensional separationless steady flow of an inviscid compressible fluid through a turbine cascade. For the subsonic flow over the profile, the method of integral equations [1] is used. It has a high accuracy in the region of the entry edge because of the concentration on it of the computational points, and for an incompressible fluid it requires comparatively little computer time compared with other methods. For all velocities, the stabilization method [2] is used successfully. Application of a finite-difference through-computation scheme makes it possible to carry out calculations without explicit separation of singularities (discontinuities) in the flow. Shock waves are obtained in this case as narrow regions (a few cells of the difference mesh) with large gradients of the parameters. In the existing variant of the method [3], the computational mesh is constructed manually, which entails much time and requires considerable experience of such work. It appears expedient to combine these methods and use the method of integral equations to specify the initial flow field and an automated construction of the difference mesh on which the calculation is then continued by the stabilization method. In the greater part of the flow, the difference mesh is constructed from the equipotentials and streamlines of the flow of an incompressible fluid. This ensures a cell shape which is nearly orthogonal, and this shortens the computing time and raises the accuracy of the results in the stabilization method.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 4, pp. 105–111, July–August, 1980.We are grateful to A. N. Kraiko for suggesting the subject and for constant interest and to G. Yu. Stepanov for interest and advice.