饱和黏弹性多孔介质中的平面波及能量耗散

研究了流体饱和不可压黏弹性多孔介质中的非均匀平面波及其能量流和能量耗散规律. 在流 相和固相物质微观不可压、固相骨架宏观服从积分型本构关系和小变形的假定下，利用 Helmholtz分解，得到了饱和黏弹性多孔介质中非均匀平面波的一般解以及纵波、横波相速 度和衰减率等的解析表达式，分析了平面波传播矢量和衰减矢量之间的关系. 数值结果表明 孔隙流体与固相骨架间的相互作用以及固相骨架的黏性对波的相速度、衰减率等有着显著的 影响. 同时，得到了饱和黏弹性多孔介质的能量方程，给出了能量流矢量和能量耗散率. 对 非均匀平面纵波和横波，推导了平均能量流矢量和平均能量耗散率的解析表达式.     

流体饱和多孔介质黏弹性动力人工边界

基于Biot流体饱和多孔介质本构方程,分别考察具有辐射阻尼性质的外行柱面波和球 面波在圆柱面和球面人工边界上引起的法向、切向应力的表达式. 在应力表达形式上,固相 介质和孔隙流体的法向和切向应力都是由两项组成,它们分别与质点的位移和速度成正比, 因此,可在人工边界的法向和切向设置连续分布的并联弹簧------黏滞阻尼器,用来模拟人工边 界以外的无限域介质对来自有限计算域的外行波动的能量吸收作用,从而形成了流体饱和多 孔介质的黏弹性动力人工边界. 流体饱和多孔介质的黏弹性动力人工边界可方便地与大型通 用软件结合,用于分析饱和土中复杂的结构-地基动力相互作用问题. 算例表明流体饱和多 孔介质黏弹性动力人工边界具有较好的精度和稳定性.     

多孔介质内黏弹性流体的热对流稳定性研究

基于修正的Darcy模型, 介绍了多孔介质内黏弹性流体热对流稳定性研究的现状和主要进展. 通过线性稳定性理论, 分析计算多孔介质几何形状(水平多孔介质层、多孔圆柱以及多孔方腔)、热边界条件(底部等温加热、底部等热流加热、底部对流换热以及顶部自由开口边界)、黏弹性流体的流动模型(Darcy-Jeffrey, Darcy-Brinkman-Oldroyd以及Darcy-Brinkman -Maxwell模型)、局部热非平衡效应以及旋转效应对黏弹性流体热对流失稳的临界Rayleigh数的影响. 利用弱非线性分析方法, 揭示失稳临界点附近热对流流动的分叉情况, 以及失稳临界点附近黏弹性流体换热Nusselt数的解析表达式. 采用数值模拟方法, 研究高Rayleigh数下黏弹性流体换热Nusselt数和流场的演化规律,分析各参数对黏弹性流体热对流失稳和对流换热速率的影响.主要结果:(1)流体的黏弹性能够促进振荡对流的发生;(2)旋转效应、流体与多孔介质间的传热能够抑制黏弹性流体的热对流失稳;(3)在临界Rayleigh数附近,静态对流分叉解是超临界稳定的, 而振荡对流分叉可能是超临界或者亚临界的,主要取决于流体的黏弹性参数、Prandtl数以及Darcy数;(4)随着Rayleigh数的增加,热对流的流场从单个涡胞逐渐演化为多个不规则单元涡胞, 最后发展为混沌状态.      

不可压饱和多孔弹性简支梁的动力响应

在杆件弯曲小变形的假定下,考虑杆件的侧向变形因素,根据多孔介质理论,本文首先建立了不可压饱和多孔弹性梁弯曲变形时动力响应的控制方程。其次,基于所建立的控制微分方程,利用变量分离法,研究了两端可渗透的饱和多孔弹性简支梁在梁中间集中载荷作用下的动力响应,得到了不同物性参数下简支梁动态弯曲时挠度和孔隙流体压力等效力偶等随时间的响应曲线。研究发现由于孔隙流体和固相骨架的相互作用,不可压饱和多孔弹性梁挠度的动力响应具有粘性特征,同时,随着时间的增加,饱和多孔弹性梁的挠度、弯矩等最终趋于经典弹性梁的静挠度、弯矩,此时,孔隙流体压力为零,梁的固相骨架承担所有的外载荷。     

高应变率下延性多孔介质中孔洞的动态演化

本文提出了一个新的材料延性动态损伤模型，模型中不但包括了率效应，同时还考虑了惯性效应，孔洞表面能变化和材料硬化对孔洞演化的影响。此外，在模型中同时考虑了体应力和偏应力对孔洞演化的作用，从孔洞演化方程地接到了孔洞增长和压缩应力临界表达式，Ｃａｒｒｏｌｌ和Ｈｏｌｔ结果作为该表达式的一个特例而得出，模型的数值分析得出以下结论：①延性孔洞的动太增长对率效应十分敏感，应变率越高，孔洞增长越快；②惯性效应在主     

不可压饱和多孔弹性梁、杆动力响应的数学模型

基于多孔介质理论,首先建立了饱和多孔弹性杆件弯曲与轴向变形时动力响应的数学模型.其次,基于多孔弹性梁弯曲变形的数学模型,利用Laplace变换,分析了两端可渗透的饱和多孔弹性悬臂梁在自由端受阶梯载荷作用下的动静力响应,给出了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等物理量随时间的响应曲线.发现不可压多孔弹性梁的拟静态响应亦存在Mandel-Cryer现象,多孔弹性梁的挠度具有与粘弹性梁挠度类似的蠕变特征,然而,其应力响应不同于粘弹性梁,随着时间的增加,梁拟静态响应的弯矩逐渐增加,并达到一个稳态值.这些结果有助于揭示植物根茎等力学行为的机理.     

对饱和多孔介质波动问题中本构关系的探讨

依据连续介质力学的混合物理论，讨论了饱和多孔介质的波动一般方程。其中考虑了质量耦合、初始压力及其浮力效应等。对于澄清这一研究课题中的某些基本概念，以及完善其理论体系具有重要意义。）其中ｖ￣（２）为速度涨落，于是单位体积内流相的动能为：因此将流相动能的涨落Ｋ￣（２）归于自由能更自然些。又因为当：ｖ￣（２＊）≡ｖ￣１时，即Δｖ＝ｖ￣１－ｖ￣２＝０时ｖ￣（２）＝０，此时Ｋ￣（２）＝０而且显然Ｋ￣（２）随着Δｖ·Δｖ的增大而增大，一个简单的近似是令Ｋ￣（２）＝ρ·Δｖ＿ｋ·Δｖ＿ｋ并且是常数。由此可导出表现质量力。由自由能的非负性质知，所以有．由（４．１３ａ）很容易推出质量密度矩阵的对称性和正定性。参考文献     

多孔介质的流变模型研究

多孔介质在应力作用下具有弹性变形和黏性变形两种完全不同的变形机制，多孔介质的弹性变形是由介质的本体有效力所致，而黏性变表则是由介质的结构有效应力所致。多孔介质的总变形为弹性变形和黏性变形的叠加，计算多孔介质总应变量的流变模型必须同时采用本体有效应力和结构有效应力（双重有效应力），而传统的流变模型仅采用Terzaghi有效应力是不妥当的，它无法正确描述多孔的应变行为，采用了双重有效应力之后的流变模型，通过调节介质特性参数，可以拟合介质的实际应变行为，并且把多孔介质与普通固体联系了起来。     

简支饱和多孔弹性梁的非线性动力响应

基于饱和多孔弹性梁大挠度变形的数学模型,利用Galerkin截断法,本文研究了两端可渗透的简支饱和多孔弹性梁分别在突加横向均布常载荷和简谐载荷作用下的动力响应,得到了梁弯曲时挠度、弯矩以及孔隙流体压力等效力偶等随时间的响应,考察了不同载荷下多孔弹性梁弯曲的响应特征.结果表明:随着载荷的增加,在常载荷作用下多孔弹性梁非线性大挠度响应与线性小挠度的差别愈加明显,而在简谐载荷作用下,多孔弹性梁的动力响应呈现较丰富的性态,相图由最初的单一椭圆曲线不断变形,形状随载荷幅值的增加而逐渐复杂,同时,时程曲线也由简单正弦曲线变为具有多峰值特征的一个周期曲线.     

一维流体饱和粘弹性多孔介质层的动力响应

本文研究了不可压流体饱和粘弹性多孔介质层的一维动力响应问题。基于粘弹性理论和多孔介质理论,在流相和固相微观不可压、固相骨架服从粘弹性积分型本构关系和小变形的假定下,建立了不可压流体饱和粘弹性多孔介质层一维动力响应的数学模型,利用Laplace变换,求得了原初边值问题在变换空间中的解析解,并利用Laplace逆变换的Crump数值反演方法,得到原动力响应问题的数值解。数值研究了饱和标准线性粘弹性多孔介质层的动力响应,分析了固相位移、渗流速度、孔隙压力及固相有效应力等的响应特征。结果表明,与不可压流体饱和弹性多孔介质相同,不可压流体饱和粘弹性多孔介质中亦只存在一个纵波,并且固相骨架的粘性对动力行为有显著的影响。     

Investigation of convection in a horizontal cylindrical layer of a saturated porous medium

The structures of the convective motions and the nature of the heat transfer in a horizontal cylindrical layer are studied numerically for the Forchheimer model of a porous medium in the Boussinesq approximation. New asymmetric solutions of the equations of convection flow through a porous medium are found. Their development, domains of existence, and stability are investigated. One consists of a multivortex structure with asymmetric vortices in the near-polar region. Another asymmetric solution is realized at large Grashof numbers in the form of a convective plume deflected from the vertical. The threshold Grashof number of formation of the asymmetric motions depends on the Prandtl number and the cylindrical layer thickness.     

Some current studies on dynamic response of the circular cylindrical sbhells in elastic medium

This paper summarizes and comments on the current situation of study on dynamic response of the circular cylindrical shells in elastic medium. After mentioning a brief introduction about early researches the authers with emphasis recount the general situation about recent development in this field, especially some important and effective methods for solving problems are pointed out. Finally we try to give a look ahead of the study on this subject.Project supported by the Science Foundation of the Chinese Academy of Sciences.     

Problem of filling of a spherical cavity in a Kelvin-Voigt medium

This paper is concerned with mathematical modeling and solution of the problem of the collapse of a spherical cavity in a viscoelastic medium under the action of constant pressure at infinity. A differential equation of motion for the cavity boundary is constructed and solved numerically. The existence of three modes of motion of the boundary is established, and a map of these modes in the plane of the determining parameters is constructed. Asymptotic forms of the solutions of the problem for all modes are constructed. The problem of cavity collapse with capillary forces taken into account is formulated and solved. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 5, pp. 93–101, September–October, 2008.     

Stress distribution around a circular cylindrical cavity in a prestressed plate

A system of equilibrium equations for nonthin transversely isotropic plates with a uniform prestress field is derived by expanding the unknown functions into Fourier-Legendre series. A method of finding the general solution of this system is expounded and used to determine the stress state of a plate with a circular cylindrical cavity __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 1, pp. 28–39, January 2008.     

Acoustic radiation force of a solid elastic sphere immersed in a cylindrical cavity filled with ideal fluid

An expression for the acoustic radiation force function on a solid elastic spherical particle placed in an infinite rigid cylindrical cavity filled with an ideal fluid is deduced when the incident wave is a plane progressive wave propagated along the cylindrical axis. The acoustic radiation force of the spherical particle with different materials was computed to validate the theory. The simulation results demonstrate that the acoustic radiation force changes demonstrably because of the influence of the reflective acoustic wave from the cylindrical cavity. The sharp resonance peaks, which result from the resonance of the fluid-filled cylindrical cavity, appear at the same positions in the acoustic radiation force curve for the spherical particle with different radii and materials. Relative radius, which is the ratio of the sphere radius and the cylindrical cavity radius, has more influence on acoustic radiation force. Moreover, the negative radiation forces, which are opposite to the progressive directions of the plane wave, are observed at certain frequencies.     

Analyses of stresses in cylindrical shells of medium length with arbitrary open-section

Based on (1), we discuss the method of solution for the stress states in cylindrical shells of medium length with arbitrary open-section. As an example, the stresses in a channel with semi-circle section were calculated, and the results were compared with those from the elementary beam theory.     

Tension of a nonthin transversely isotropic plate with a noncircular cylindrical cavity

The three-dimensional stress state of a transversely isotropic plate with a nearly circular cylindrical cavity is examined. The cavity surface is subject to normal and tangential stresses and the plate is subject at infinity to tensile and shear forces. The problem is solved by expanding unknown functions into Fourier-Legendre series in the thickness coordinate and using the boundary-shape perturbation method. The equations and recurrence formulas needed to solve the problem in an arbitrary approximation are presented __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 11, pp. 101–113, November 2006.     

浅埋的圆柱形孔洞对SH波的散射与地震动

研究了浅埋的圆柱形孔洞对以任意方向入射的平面SH波的散射与地震动问题。利用复变函数和多极坐标方法构造了问题的位移解。当入射波的波长与圆孔的半径相比较小时,地震动将受到较大的影响。影响地震动有三个主要参数:(1)SH波的入射角0;(2)入射波波数,即圆柱形孔洞的半径与入射波半波长之比;(3)h/R,即圆柱形孔洞至表面的距离与圆孔半径之比。当较大时,地震动幅值变化激烈,位移幅值可出现跳动和放大的现象。当h/R增大至10～12时,位移幅值变化恢复至半空间的情况,表明圆柱形孔洞的影响可被忽略。     

Level‐set based numerical simulation of a migrating and dissolving liquid drop in a cylindrical cavity

In the present paper the dissolution of a binary liquid drop having a miscibility gap and migrating due to thermo‐solutal capillary convection in a cylindrical cavity is studied numerically. The interest in studying this problem is twofold. From a side, in the absence of gravity, capillary migration is one of the main physical mechanisms to set into motion dispersed liquid phases and from the other side, phase equilibria of multi‐component liquid systems, ubiquitous in applications, often exhibit a miscibility gap. The drop capillary migration is due to an imposed temperature gradient between the cavity top and bottom walls. The drop dissolution is due to the fact that initial composition and volume values, and thermal boundary conditions are only compatible with a final single phase equilibrium state. In order to study the drop migration along the cavity and the coupling with dissolution, a previously developed planar two‐dimensional code is extended to treat axis‐symmetric geometries. The code is based on a finite volume formulation. A level‐set technique is used for describing the dynamics of the interface separating the different phases and for mollifying the interface discontinuities between them. The level‐set related tools of redistancing and off‐interface extension are used to enhance code resolution in the critical interface region. Migration speeds and volume variations are determined for different drop radii. Copyright © 2004 John Wiley & Sons, Ltd.     

封闭圆柱腔内旋转点声源声压的近场频域解

在圆柱腔壁面为声学硬壁面的假定下，利用圆柱腔体空间内的格林函数导出了圆柱腔内旋转运动点声源空间声压的近场频域解；进而利用该频域解计算分析了旋转简谐单极子点源的声学特性。计算分析表明：由于考虑了壁面的影响，在圆柱腔体内，点声源旋转频率和源频率的变化将改变声压的空间指向性；源频率和旋转频率的增加伴随明显的多普勒效应，同时基频声压沿圆柱腔体轴线的分布从比较平缓变为明显波动；沿半径方向基频声压也波动变化。