磁性流体密封压差的数值计算

分析并简化了含磁性流体的磁场模型,以密封理论为基础,采用有限元法计算出磁性流体密封的磁场和等压线分布,进而计算出密封压差,分析了密封压差与磁性流体量的关系,给出了多级密封压差的计算方法,同时分析了转速对密封能力的影响。结果发现,数值计算结果与实验结果相一致。     

Green's Function for Orthorhombic Aggregates of Weakly Anisotropic Cubic Crystals

Herein a closed but approximate formula of the Green's function is obtained for orthorhombic aggregates of cubic crystallites. This formula, which includes three material constants and three texture coefficients, accounts for the effects of the orientation distribution function (ODF) up to terms linear in the texture coefficients. Thus it is expected that our formula would be applicable to aggregates with weak texture or to materials such as aluminum whose single crystal has weak anisotropy. The approximate formula remains valid and assumes a simpler form when the polycrystal reduces to a weakly anisotropic cubic crystal. Two examples are presented to compare predictions from our formula with those from Nishioka and Lothe's formula and from Synge's contour integral through numerical integration.     

跨音速压气机级的三维周期性非定常流动计算

对跨音速压气机级动静叶排相干形成的三维非定常流场进行了数值研究，利用时间推进ＬＵ－ＳＧＳ稳式迭代法求解三维非定常欧拉方程，对流项采用高分辨率ＮＮＤ格式离散。对某压气机第一级动静叶排相干非定常流场的计算结果表明，本文方法不仅在收敛速度上明显地优于一般显式方法，而且保持了流场中激波的高分辨率，适于推广到计算量巨大的的多级轴流压气机三维非定常流场的数值分析问题     

海堤越浪的数值模拟

基于RANS方程和两方程湍流模型,采用有限体积法,将人射波波场作为人工的分布源项加人动量方程,提出了适用于VOF方法的源造波一消波技术。通过对行波及驻波的计算,分别考察了数值波浪水槽前端及末端消波段的有效性。在本文建立的数值波浪水槽内对规则波在海堤上爬高和越浪过程进行了数值模拟,并将计算结果与现有实验结果进行了比较。验证计算结果表明,数值模拟结果较好地复演了海堤越浪过程。为了研究模型尺度对越浪量的影响,文中设计了一组满足重力相似但具有不同几何比尺的数值实验模型。系列数值实验结果表明,若按重力相似换算越浪量,计算结果与实验预报值间的偏差随模型比尺的增大和堤前波浪破碎强度的增强而增大,建议在进行越浪物理模型实验时需进一步考虑模型比尺对原型预报值的影响。     

Anisotropic Elastic Materials that Uncouple All Three Displacement Components, and Existence of One-Displacement Green's Function

It was shown in an earlier paper that, under a two-dimensional deformation, there are anisotropic elastic materials for which the antiplane displacement u ₃ and the inplane displacements u ₁, u ₂ are uncoupled but the antiplane stresses σ₃₁, σ₃₂ and the inplane stresses σ₁₁, σ₁₂, σ₂₂ remain coupled. The conditions for this to be possible were derived, but they have a complicated expression. In this paper new and simpler conditions are obtained, and a general anisotropic elastic material that satisfies the conditions is presented. For this material, and for certain monoclinic materials with the symmetry plane at x ₃ = 0, we show that the unnormalized Stroh eigenvectors a _k for k = 1, 2, 3 are all real. The matrix A =[a ₁, a ₂, a ₃] is a unit matrix when the material has a symmetry plane at x ₂ = 0. Thus any one of the u ₁, u ₂, u ₃ can be the only nonzero displacement, and the solution is a one-displacement field. Application to the Green's function due to a line of concentrated force f and a line dislocation with Burgers vector v in the infinite space, the half-space with a rigid boundary, and the infinite space with an elliptic rigid inclusion shows that one can indeed have a one-displacement field u ₁, u ₂ or u ₃. One can also have a two-displacement field polarized on a plane other than the (x ₁, x ₂)-plane. The material that uncouples u ₁, u ₂, u ₃ is not as restrictive as one might have thought. It can be triclinic, monoclinic, orthotropic, tetragonal, transversely isotropic, or cubic. However, it cannot be isotropic. This revised version was published online in August 2006 with corrections to the Cover Date.     

孤立波与多个潜水障碍物相互作用的数值研究

本文以完整二维Navier-Stokes方程为控制方程,采用VOF界面跟踪技术和差分方法,数值计算孤立波与多个淹没水下物体相互作用的问题。本文对潜水物体的高度接近水深时,孤立波通过水下孤立直立方柱、两个间距较大的水下直立方柱和两个间距较小的水下直立方柱等三种情况分别进行了计算,给出了波形随时间的演化图,可以看到反射波、前传波和跟随振荡型小波列的生成。对孤立波通过水下孤立直立方柱情形的计算结果,与实验结果和势流理论结果进行了比较。     

Computations of Interaction of an Incident Oblique Shock Wave with a Turbulent Boundary Layer on a Flat Plate

Results of numerical simulation of interaction between an oblique shock wave and a turbulent boundary layer formed in a supersonic (Mach number M =5) flow past a flat plate are presented. The computations are performed for three cases of interaction of different intensity, which result in an attached or detached flow. Numerical results are compared with experimental data. The effect of flow turbulence and shockwave unsteadiness on flow parameters is studied.     

非均匀损伤介质中波传播的数值解

对弹性波在非均匀损伤介质中的传播理论进行了研究。通过将非均匀损伤区域离散成分层均匀的区域,结合相邻区域交界面处的连续条件,推导出了以右行波、左行波为状态向量的波动方程和传递矩阵。对几种非均匀损伤介质中波的传播进行了实例数值计算,并和其解析解的结果进行了比较,讨论了弹性波在非均匀损伤介质中传播的一般性质。     

Deformation of Euler-Bernoulli Beams by Means of Modified Green's Function: Application of Fredholm Alternative Theorem

This article deals with the determination of the static displacement function of an Euler-Bernoulli beam with two guided supports. To this end, the Green's function method is employed and exact solution is obtained. The Green's function of the problem is constructed, using pertinent boundary conditions of the problem. Nevertheless, the problem does not admit a Green's function due to a mathematical contradiction. In order to eliminate the trouble, the Fredholm Alternative Theorem is utilized and the Green's function is modified. In this case, application of this theorem adds a particular term to the Green's function which gives rise to an arbitrary constant in the Green's function. Moreover, it is shown that the problem may have no solution or an infinite number of solutions. Besides, the necessary condition for having any solution is investigated. This requirement, which states a significant rule in the mechanics of solids, is the static equilibrium of vertical forces acting on the beam. Some examples are presented and results are thoroughly discussed.     

周向导波在空心圆柱体中传播的数值模拟研究

在石油、化工、食品和城市供水等行业,由于腐蚀等原因而引起泄漏事故,造成巨大的经济损失和资源浪费.因此,能够预先检测到管道缺陷避免事故发生显得十分重要.基于弹性动力学理论,利用Matlab编程实现周向导波在空心圆柱体中传播的数值模拟.利用特征函数展开法计算出周向导波的频散曲线及0.6MHz和1MHz各个模态的时域波形图并据此分析激励方式同周向导波传播的关系,频率为1MHz或0.6MHz时,1和2模态在空心圆柱壳中传播的周向导波中处于主导地位:各模态的径向或周向位移幅度会受到激励入射角的影响而变化并且周向导波的位移幅度随着激励入射角的增加而增加:得到特定频率下,周向导波各模态径向、周向位移的变化关系.本文结果为进一步利用实验方法研究周向导波在空心圆柱体中传播特性和缺陷检测奠定了基础.     

Wave Resistance of Amphibian Aircushion Vehicles During Motion on Ice Fields

The steady motion of amphibian aircushion vehicles over water covered with continuous ice is studied. The ice sheet is simulated by a viscoelastic ice plate. An analysis is made of the effect of the aspect ratio of the vehicle, the depth of the water reservoir, and ice characteristics on the wave resistance of the vehicle and the speed of the vehicle at which the wave resistance is maximal.     

Numerical Simulation of the Wave Pattern within a Harbor due to Ship Waves

A numerical simulation methodology to predict the wave pattern within a harbor due to ship waves is presented. It is based on the spectral decomposition of the solitary incident wave, assuming that non-linear effects are negligible, which is the usual case. Each component is transformed using Berkhoff's mild slope model. Superposition of the different solutions allows reconstruction of the wave temporal history, from which the maximum wave height at any point within the harbor can be obtained. An application of the proposed methodology to study the protection structure for the Argentine Yatch Club, located within Buenos Aires harbor, is presented. Details about the schematization and the necessary data are explained. The maximum wave height distribution for the projected structure and the original situation are compared. A variant of the project is proposed and analyzed, showing a significant extra attenuation of the maximum wave heights. This illustrates the utility of the methodology for design purposes.     

Validation of the Slow Compressional Wave in Porous Media: Comparison of Experiments and Numerical Simulations

We perform numerical simulation of ultrasonic experiments on poroelastic samples, in which Biot's slow compressional wave had been observed. The simulation is performed using OASES modeling code, which allows to compute elastic wave fields in layered poroelastic media. Modeled were the experiments of Plona (1980), Rasolofosaon (1988), and our own measurements. In all the three situations, a good agreement between experiment and simulations has been observed. This further confirms the fact that Biot's theory of poroelasticity, on which the simulations were based, adequately describes the behavior of the porous materials under investigations at ultrasonic frequencies.     

旋转壳的数值传递函数方法

应用数值传递函数方法建立一种用于分析旋转壳静力、动力响应的截锥壳单元,在本方法中,单元的位移在环向展开为Fourier级数的形式,应用薄壳理论可以得到解耦的微分方程,通过Laplace变换可以将方程转化为频域内的常微分方程,将其表示为状态空间形式后,可以应用数值传递函数方法求解,对复杂的系统可以应用与有限元类似的方法,划分多个单元组合求解,文中给出了几种旋转壳的动力、静力问题的数值算例,并与其它方法进行了比较,表明本文方法具有精度高,计算方便等特点。     

Analysis and Numerical Evaluation of the Drift Function of Darwin and Lighthill for Axisymmetric Flows

Darwin and Lighthill have studied the deformation or drift of material surfaces produced in the course of steady irrotational flows past cylinders and spheres. This paper is devoted to a reconsideration of the same problem for axisymmetric inviscid flows produced by arbitrary shaped bodies. Such problems cannot be solved by analytical techniques and, therefore, a recently developed algorithm has been brought into consideration to obtain the drift function numerically. First, the case of a sphere is considered to check and compare with the results obtained by Lighthill. Second, the developed algorithm has been used to obtain the drift functions in chosen patches of the stream surfaces produced by a prolate ellipsoid. Received 14 November 1996 and accepted 22 April 1997