波浪引起的海底土层中的应力场和位移场

本文分析了波浪压力引起的海底土层中孔隙水流动的状况,给出了简化的力学模型来描述海底土层中孔隙水流动所引起的应力场和位移场.对波浪的分析采用了线性波的理论假设,波浪引起的压力计算采用高斯分布理论.孔隙流体和土骨架都假设为可压缩的,孔隙水在海底土层中的流动满足达西定律.由表示质量守恒的连续性方程与应力平衡方程出发,得到了由位移所表示的高阶微分方程,并给出了海底土层中位移场和应力场的解.其解依赖于孔隙流体和土骨架的相对压缩性,土的泊桑比和剪切模量,以及土介质的渗透率和各向异性等特性.文章还结合实际情况进行了计算,并分析了上述参数的影响.对如何应用上述结果提出了相应的参考意见.     

关于传递矩阵法分析饱和成层介质响应问题的讨论

水平成层土体的地震响应分析(自由场分析)是地震工程领域地震波散射问题的前提基础,由于饱和多孔方程的复杂性,以往的研究大多集中于干土情形,对于饱和土情形的研究相对较少.而实际工程中,地下水位以下,土体孔隙中充满流体,应考虑饱和多孔介质模型.基于Biot多孔介质模型,考虑饱和土中固液相对运动引起的衰减,采用Thomson--Haskell传递矩阵方法得到了饱和成层土体在地震波入射情形时的稳态反应,经傅里叶反变换,可得到时域暂态反应.通过SV波从基岩入射至上覆饱和土层的数值算例,验证了该方法的有效性.发现和初步阐明了计算中出现的两类违背因果律(即响应先于输入)的现象:(1)当SV波入射角度大于导致基岩中反射P波为非均匀波的临界角时,会使得计算结果违背因果律.因此,当入射角超过临界角时,非均匀波的表示尚需进一步完善;(2)由于P2波的衰减,当与稳态波衰减有关的渗透率、土层厚度、入射波频率等参数导致衰减系数超过计算机表示精度时,会出现结果违背因果律现象,并据此得到了满足因果律的参数范围,该范围可作为实际计算时的一个上界.该工作为采用传递矩阵法分析水平饱和土层自由场响应提供了指导依据,且地下水位以上可采用干土模型,水位以下采用饱和土模型,更符合实际情形.     

P1波入射准饱和土的波型转换问题

针对各向同性准饱和土体中体波的传播问题，本文深入研究了P1波(第一纵波)入射土体时所出现的波型转换现象。在准饱和土波动方程的基础上，通过数值算例，找出标准参数土体模型的临界饱和度及波型转换角，并分析土体孔隙率、渗流系数、水压力和入射波频率等参数对该现象的影响。     

弹性波在准饱和土和弹性土界面的反射与透射

考虑土颗粒和孔隙流体的压缩性,根据准饱和土体的变形连续条件和平衡条件,推导了P1波(快压缩波)、P2波(慢压缩波)和S波的计算公式．根据Snell定理,研究了平面P1波从准饱和土体入射到弹性土层时在界面上的反射和透射,并绘制了反射系数、透射系数和界面应力、位移随入射角的变化曲线,结果表明:饱和度是影响界面效应的一个重要因素,当饱和度低于99％时,反射和透射系数以及应力和位移基本相等,变化不大,但当饱和度为100％时,虽然饱和度比99％提高了1％,但反射和透射系数以及应力和位移值比饱和度为99％时明显降低．     

饱和土–隧道动力响应的2.5维有限元–边界元耦合模型

针对饱和土中异形隧道的三维动力响应问题,建立了2.5维有限元与边界元耦合模型.将隧道结构视为弹性体,采用2.5维有限元建立隧道模型;将地基土视为饱和多孔介质,采用2.5维边界元建立饱和土体模型.借助组合辅助问题基本解消除了边界积分方程的奇异性.利用饱和土与隧道接触面的位移、面力连续和完全透水或完全不透水边界条件,实现2.5维有限元和边界元模型的耦合求解.模型具有计算效率高、适用范围广的优点.通过与完全透水和完全不透水边界条件下轴对称问题的半解析解以及单相介质的2.5维有限元与边界元耦合模型对比,验证了本文模型的正确性.最后利用该模型计算了饱和土体中类矩形隧道在移动载荷作用下的三维动力响应,分析了不同土体渗透性下位移及孔隙水压力沿隧道轴向、环向和深度的分布规律.结果表明:(1)孔隙水压力随土体渗透性增大而显著减小,位移受土体渗透性影响小;(2)位移及孔隙水压力在隧道环向主要分布在距载荷作用点两侧约2 m的范围内;(3)孔隙水压力沿深度的衰减比土体位移快,且孔隙水压力和轴向位移沿深度的分布受土体渗透性影响大.     

饱和土地基在简谐荷载作用下的解析分析

基于多孔介质混合物理论,用解析的方法研究了不可压饱和土地基受到简谐荷载作用下的动力响应问题。利用Fourier积分变换求解耦合方程组,得到了二维饱和土介质在简谐荷载作用下的通解。针对表面透水的具有下卧基岩的饱和土层以及半无限饱和土地基的边界条件,获得了固体骨架位移、孔隙流体位移、固体骨架有效应力以及孔隙流体压力的积分形式解答,并通过数值算例分析了饱和土地基在简谐荷载作用下的响应。     

多孔饱和半空间上弹性圆板垂直振动的积分方程

应用新的方法求解多孔饱和固体的动力基本方程－Ｂｉｏｔ波动方程,首先把Ｂｉｏｔ波动方程化为仅有土骨架位移和孔隙水压力的偏微分方程组,并且逐次解耦方法（不引入位移势函数）求解此偏微分方程组,然后按混合边值条件建立多孔饱和半空间上弹性圆板垂直振动的对偶积分方程,用Ａｂｅｌ变换化对偶积分方程为第二类Ｆｒｅｄｈｏｌｍ积分方程。文中考虑两种孔隙流体的表面边界条件：（ａ）半空间表面（包括圆板与半空间的接触面）是     

一维渗透力与浮力

为简化分析,针对一维渗流问题,研究提出了土力学中渗透力和浮力的两种推导方法,作为传统的宏观尺度孔隙水隔离体法的有益补充.弹性力学平衡微分方程、土力学Terzaghi有效应力方程和流体力学简化Bernoulli方程构成本文分析渗透力和浮力的3个基本方程.在基本方程基础上,容易推导出相应的骨架和孔隙水两种隔离体的平衡微分方程,从而在静力平衡范畴内揭示渗透力和浮力的内涵.单位体积饱和土体的渗透力,源于总水头压力的梯度,而浮力则源于位置水头压力在竖向的梯度,这两者统一于骨架或孔隙水的平衡微分方程.实际工程关注的有效应力计算问题,一般可以直接应用3个基本方程来确定;只有在简化条件下可按渗透力和浮力计算土体中有效应力分布规律.还讨论了若干研究热点问题,重点探讨了当前一种渗透力新定义j=nγ_wi在形式上的合理性以及在实际应用中可能存在的风险,并验证了一维渗透力的一种经典精细化表述结果中考虑渗流速度时间导数的严谨性.在土力学渗透力和浮力问题研究中应重视和正确应用Terzaghi有效应力方程.     

饱和土和衬砌相互作用的圆形隧洞动力特性

利用 Darcy 渗透定律,通过引入衬砌和土体的相对渗透系数,建立了隧洞边界部分透水条件。将土骨架视为具有分数导数粘弹性本构关系的粘弹性体,基于 Biot 理论,通过界面连续性条件在频率域内给出了简谐轴对称荷载或流体压力作用下饱和粘弹性土-弹性衬砌系统耦合振动时饱和土、衬砌的位移、应力、孔隙水压力表达式,并通过算例分别考察了简谐轴对称荷载、流体压力作用下的分数导数阶数、材料参数比、渗透系数对系统的径向位移幅值 U 、孔隙水压力幅值P 的影响,结果表明：随着分数导数阶数和材料参数比的增加,系统的响应幅值逐渐减小;随着渗透系数κ的增加,轴对称荷载作用下的土体位移幅值和孔隙水压力幅值逐渐减小,当κ大于100时,U 和 P 值无明显变化,流体压力作用下的土体位移幅值和孔隙水压力幅值逐渐增大,当κ大于1时, U 和 P 值无明显变化。     

水饱和岩石中爆炸应力波传播的数值模拟

基于连续介质力学和不相融混合物理论,假定组分间无相对运动,采用B-W-N-B有效应力准则,在屈服面中引入孔隙影响因子,提出了一种多孔含水介质流固耦合的本构模型,并给出了孔隙的演化方程,对水饱和凝灰岩介质中的爆炸应力波传播作了数值模拟。数值计算给出的应力波形与实测结果有良好的一致性。采用本文中所述流固耦合本构模型,可很好地预测水饱和岩石中爆炸波的演化规律。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).