横观各向同性含液饱和多孔介质黏弹性瞬态反应

在考虑横观各向同性含液饱和多孔介质固体骨架和流体可压缩性以及固体骨架的黏弹性特征下,基于横观各向同性含液饱和多孔介质u-w形式的三维动力控制方程,以固相位移u、液相相对位移w为基本未知量,综合运用Laplace变换、双重Fourier变换等方法,在直角坐标系下通过引入中间变量,将六元2阶动力控制方程组化为两组各含4个未知变量的常微分方程组,给出了直角坐标系下横观各向同性含液饱和多孔介质三维黏弹性动力反应的积分形式一般解;作为理论推导的验证,通过引入初始条件和边界条件,对横观各向同性含液饱和多孔介质半空间黏弹性瞬态反应问题进行了求解。解答的退化验证表明,所推导的理论解是正确的。     

层状半空间中洞室对入射平面SV波的散射

利用间接边界元法,在频域内求解了层状弹性半空间中洞室对入射平面SV波的散射问题。通过自由场反应分析,求得假想洞室边界上各点位移和各单元应力响应。在洞室边界各个单元上施加虚拟分布荷载,求得位移和应力的格林函数。根据应力边界条件确定虚拟分布荷载,将自由场位移响应和虚拟分布荷载产生的位移响应叠加起来,即得到问题的解答。比较了层状半空间和均匀半空间中洞室对入射平面SV波的放大作用。结果表明,层状半空间情况有可能导致较大的地表位移幅值,尤其是对于较低频率入射波。     

Horizontal vibrations of a disk on a poroelastic half-space

This paper analytically examines the horizontal vibration of a rigid disk on a saturated poroelastic half-space. The pressure-solid displacement form of the harmonic equations of motion for asymmetric dynamic problem are developed from the form of the equations originally presented by Biot. Making use of a new method the solution of the above equations is obtained. According to the mixed boundary -value conditions, the dual integral equations of the horizontal vibration of a rigid disk on a saturated poroelastic half-space are established. By appropriate transforms, it is shown that the dual integral equations can be reduced to a pair of Fredholm integral equations of the second kind, whose solutions are then computed. Numerical results for the horizontal dynamic compliance coefficient are given at the end of this paper.     

各向异性双相介质中多极源声波测井理论研究

将各向异性地层模拟为横向各向同性双相介质 ,使其对称主轴与井孔中心轴线相垂直 .本文针对井孔模型下的多极源声波测井问题提出了一种解析的摄动方法 .将各向异性双相介质看作是由对应的各向同性双相介质加上摄动部分而组成 ,并引入 5个关于弹性模量偏离各向同性时的摄动量 ,对对称主轴与井轴垂直的各向异性双相介质的多极源声波测井问题进行了系统的理论研究和分析 ,严格求解了井孔内外的零级和一级近似声波场 .发现在n级多极源激发下 ,井孔中不仅出现n级多极声场 ,而且还能激发出与各向同性地层时不相同的高于和低于n级的多极声波场     

利用灰岩断层面形貌特征识别罗云山山前断裂古地震信息

以山西地堑系罗云山山前断裂的基岩断层面为例,用陆基Li DAR扫描获取了断层面形貌学数据,采用各向同性变差函数法计算了断层面形貌的分维值。结果显示,断层面形貌在高度上具有显著的分带性,每个分带的特征分维值随断层面高度的增加呈阶跃式变化。分析认为这种形貌特征可能反映了断层面与多次地震相关的间歇式出露方式。分带高度指示了约3m和1m的2组同震倾滑位移量,各分带之间的分维值过渡带是断层面缓慢剥露的结果。与片麻岩区的断层面形貌学研究结果相比,灰岩断层面形貌学分维值随着暴露时间的增加而减小,2种岩性区的断层面分维值D与断层面高度H呈现类似镜像的关系,在气候构造条件大致相同的情况下,考虑为岩性差异的影响。     

三维层状黏弹性半空间中球面SH、P和SV波源自由场

采用刚度矩阵方法结合Hankel积分变换,求解了层状黏弹性半空间中球面SH、P和SV波的自由波场.首先,在柱坐标系下建立层状黏弹性半空间的反轴对称(柱面SH波)和轴对称(柱面P-SV波)情况精确动力刚度矩阵.进而由Hankel变换将空间域内的球面波展开为波数域内柱面波的叠加,然后将球面波源所在层的上下端面固定,求得固定层内的动力响应和固定端面反力,将固端反力反向施加到层状黏弹性半空间上,采用直接刚度法求得固端反力的动力响应,叠加固定层内和固端反力动力响应,求得波数域内球面波源动力响应.最后由Hankel积分逆变换求得频率-空间域内球面波源自由场,时域结果由傅里叶逆变换求得.文中验证了方法的正确性,并以均匀半空间和基岩上单一土层中球面SH、P和SV波为例分别在频域和时域内进行了数值计算分析.研究表明基岩上单一土层中球面波自由场与均匀半空间情况有着本质差异;基岩上单一土层中球面波位移频谱峰值频率与场地固有频率相对应,基岩面的存在使得基岩上单一土层地表点的位移时程非常复杂,振动持续时间明显增长;阻尼的增大显著降低了动力响应的峰值,同时也显著减少了波在土层的往复次数.     

The effects of intermediate principle stress on the mechanical behavior of transversely isotropic rocks: Insights from DEM simulations

The discrete element method (DEM) is used to study the response of anisotropic rocks under true triaxial testing. Numerical samples of seven different bedding orientations (β = 0^o, 15^o, 30^o, 45^o, 60^o, 75^o, and 90^o) are created. Six series of test simulations (σ₃ = 0, 10, 30, 50, 70, and 100 MPa) are conducted on each sample, with five different σ₂ values, varied from σ₃ to σ₁. The effects of anisotropy and intermediate stress on the peak strength, brittle-ductile transition, and degree of anisotropy are subsequently explained through underlying micromechanics. Results show a “fan-shaped” variation of the peak strength with σ₂, displaying an ascending-then-descending trend. An increasing brittleness with σ₂ is observed at lower confining pressures for all, but medium anisotropy angles. For higher confining pressures, increasing ductility with σ₂ is seen for every anisotropy angle. A U-shaped variation of peak strength with anisotropy angles is noted that flattens under high intermediate stress. Hence, for numerical models of Posidonia shale under normalized σ₂ higher than 0.76, the anisotropy effect is found to be negligible. Micromechanical analyses reveal that the stress asymmetry, suppression of weak plane action as well as the localization and coalescence of microcracks in the intact rock matrix, due to σ₂, are the contributors towards the obtained trends. Since existing failure criteria do not weigh in these features in geotechnical assessments, this paper helps future studies by providing a deeper understanding of these effects and a comprehensive data set for the analyses of anisotropic rocks under polyaxial stress conditions.     

Displacement ring load Green's functions for saturated porous transversely isotropic tri‐material full‐space

Based on the Biot's poroelastic theory and using scalar potential functions both the ring load and point load displacement Green's functions for a transversely isotropic saturated porous full‐space composed of an upper half‐space, a finite thickness middle layer and a lower half‐space is analytically presented for the first time. It is assumed that each region consists of a different transversely isotropic material. The equations of poroelastodymanics in terms of the solid displacements and the pore fluid pressure are uncoupled with the help of two scalar potential functions, so that the governing equations for the potential functions are either a second order wave equation or a repeated wave‐heat transfer equation of sixth order. With the aid of Fourier expansion with respect to circumferential direction and Hankel integral transforms with respect to the radial direction in cylindrical coordinate system, the response is determined in the form of line integrals in the real space, followed by theorem of inverse Hankel integral transforms. The solutions degenerate to a single phase elastic material, and the results are compared with previous studies, where an excellent agreement may be observed with the results provided in the literature. Some examples of displacement Green's functions are finally given to illustrate the solution. Copyright © 2016 John Wiley & Sons, Ltd.     

Forced vertical and horizontal movements of a rectangular rigid foundation on a transversely isotropic half‐space

An analytical investigation of a half‐space containing transversely isotropic material under forced vertical and horizontal displacements applied on a rectangular rigid foundation is presented in this paper. With the goal of a rigorous solution to the shape‐ and rigidity‐ induced singular mixed boundary value problem, the formulation employs scalar potential representation, the Fourier expansion and the Hankel integral transforms method to obtain the surface arbitrary point‐load solution in cylindrical coordinate system. The obtained Green's functions are rewritten in rectangular coordinate system, allowing the response of the half‐space because of an arbitrary distributed load on a rectangular surface area be given in terms of a double integral. The numerical evaluations of stresses are done with the use of an element, which is singular at the edge and the corner of the rectangle. Upon the imposition of the rigidity displacement boundary condition for a rigid foundation and the use of a set of two‐dimensional adaptive‐gradient elements, which can capture the singular behavior in the contact stress effectively, a set of new numerical results are presented to illustrate the effect of transverse isotropy on the foundation response. Copyright © 2012 John Wiley & Sons, Ltd.