Multiscale modeling of large deformation in geomechanics

Large deformation soil behavior underpins the operation and performance for a wide range of key geotechnical structures and needs to be properly considered in their modeling, analysis, and design. The material point method (MPM) has gained increasing popularity recently over conventional numerical methods such as finite element method (FEM) in tackling large deformation problems. In this study, we present a novel hierarchical coupling scheme to integrate MPM with discrete element method (DEM) for multiscale modeling of large deformation in geomechanics. The MPM is employed to treat a typical boundary value problem that may experience large deformation, and the DEM is used to derive the nonlinear material response from small strain to finite strain required by MPM for each of its material points. The proposed coupling framework not only inherits the advantages of MPM in tackling large deformation engineering problems over the use of FEM (eg, no need for remeshing to avoid mesh distortion in FEM), but also helps avoid the need for complicated, phenomenological assumptions on constitutive material models for soil exhibiting high nonlinearity at finite strain. The proposed framework lends great convenience for us to relate rich grain-scale information and key micromechanical mechanisms to macroscopic observations of granular soils over all deformation levels, from initial small-strain stage en route to large deformation regime before failure. Several classic geomechanics examples are used to demonstrate the key features the new MPM/DEM framework can offer on large deformation simulations, including biaxial compression test, rigid footing, soil-pipe interaction, and soil column collapse.     

Modelling of coupled fluid‐mechanical problems in fractured geological media using enriched finite elements

Geological environments, such as petroleum reservoirs, normally exhibit physical discontinuities, for example, fractures and faults. Because of the reduced thickness of these discontinuities, finite element formulations with strong discontinuity have been applied to the numerical modelling of geological environments. Until now, two relevant characteristics of petroleum reservoirs have not been addressed by these formulations. The first is the pore pressure jump in the direction normal to a discontinuity in a fluid‐mechanical coupling condition, which is present primarily in sealing faults owing to the contrast of permeability with the porous medium. The absence of this jump can affect the prediction of the deformability of a physical discontinuity. Furthermore, reservoir models frequently use coarse meshes. Thus, the method used to evaluate the pore pressure in the discontinuity may exhibit a strong dependence relative to the mesh refinement. Based on these characteristics, in this study, a formulation of an enriched finite element for application to coupled fluid‐mechanical problems with pre‐existing physical discontinuities saturated by a single fluid is presented. The formulation employs discontinuous interpolation functions and enables the reproduction of jumps of displacement and pore pressure associated with a discontinuity inside the element without the need to discretise it. An approximation to estimate the pore pressure in the discontinuity was developed, one which seeks to minimise the influence of refinement. The element's response is verified by comparison with a one‐dimensional analytical solution and simple examples that are simulated using commercial software. Copyright © 2015 John Wiley & Sons, Ltd.     

浮式风机一体化模拟的计算参数设置影响研究

随着海上浮式风机的大型化发展，针对漂浮式风机的一体化耦合分析越来越重要。本文利用Simo-Riflex-Aerodyn 仿真工具建立OC4-Deepcwind 漂浮式风机一体化耦合模型，分析计算时间步长、初始截断时间、弹性结构单元离散数等计算参数对模拟结果的影响，包括浮式基础运动、系缆张力、叶片受力等。结果表明：当计算时间步长为0.005 s、0.01 s、0.02 s 时，浮式风机的响应结果差别较小，而计算消耗时间相差较大，0.01 s、0.02 s 的计算时间分别是0.005 s的70%、37%。相同工况下，不同参量响应达到稳定所需时间不同，纵荡需要的时间较长，最长达200 s；不同工况下，同一参量达到稳定所需时间也不相同，切出工况需要的时间最短，较额定工况快约60 s。结构单元离散数对塔柱受力影响较小，对叶片变形影响相对明显，当叶片离散数目减小时，响应值增大12%。实际中应根据具体工况选择合理的计算时间步长、初始截断时间和弹性结构单元离散数量。     

基于路基土蠕变效应的路基预应力损失模型研究

新型预应力路基技术可用于整治路基病害或强化既有线路基,研究其预应力损失规律对保障预应力路基工程的长期安全具有重要意义。基于弹性理论提出了预应力路基侧压力板与路基边坡接触面中心点后方水平路径上附加应力的计算方法,分析表明该特征路径上水平附加应力随距板−土接触面中心点距离的增加呈现良好的指数函数衰减关系。基于水平附加应力的指数函数扩散规律及路基土的蠕变行为,推导了路基土的蠕变变形计算公式。进一步基于蠕变变形与预应力钢筋回缩变形的协调性,建立了路基预应力的损失模型和计算方法,并结合FLAC3D数值仿真开展了对比分析。研究表明：钢筋预拉力损失的理论预测曲线与数值仿真曲线吻合良好,二者偏差不足5%,论证了路基预应力损失模型的有效性;钢筋预拉力可在锚固后的60 d内达到稳定,且预拉力的稳定值可达其初始值的85%～90%,说明新型预应力路基技术可为路基土提供相当可观的稳定附加围压,从而通过改善路基土的应力状态达到强化路基的目的。     

不同渗透水压下完整泥页岩剪切−渗流特性研究

为了探索高压水力浸润带和消落带岩体在压剪荷载与流体入渗双因素耦合作用下剪切渗流特性的衍化规律,利用自主研发的一套新型岩石剪切−渗流耦合试验系统,对饱水泥页岩进行了常法向荷载条件下不同渗透水压的剪切−渗流耦合试验和剪切蠕变−渗流耦合试验。试验结果表明：（1）随着渗透水压力的增加,饱水泥页岩剪切强度峰值、法向变形峰值以及剪切蠕变破坏强度均有所下降,说明渗透水压力对岩样力学强度具有损伤效果。（2）在剪切−渗流耦合试验过程中,岩样剪应力与法向变形在线弹性变化的末期会出现明显的上下波动现象,并伴有水流流出,剪应力和法向变形的突变点在时间上具有一致性。（3）在剪切蠕变−渗流耦合试验过程中,水流渗出速率与渗透水压力大小呈正相关,随着剪切荷载施加等级的提高,岩样蠕变变形量增大,单位时间内累计渗出水量亦有所提升。     

Telecommunication interactions between microplates and basal mantle LLSVPs.SCIEI北大核心CSCD

基于海洋地质地球物理观测建立的板块构造理论意味着板块和浅部地幔共同演化,然而地幔底部尤其是大型横波低速异常区(LLSVP)与板块(尤其微板块)运动和演化之间是否存在关联仍有争议。一些研究认为LLSVP长期保持稳定,而另一些模型则认为它与各级板块存在相互作用。为此,本文通过总结前人成果,并基于近期发表的板块重建和地幔对流模型进行进一步分析,探讨微板块运动和LLSVP的演化关系。模拟结果表明,微板块与大板块类似,俯冲后通常会下沉至核幔边界。微幔块会推动地幔底部热的物质聚集并形成大的热化学结构。该热化学结构与层析成像揭示的LLSVP基本吻合。下地幔径向流速场和温度场的二阶结构与地表速度场散度的二阶结构随时间的移动轨迹相似,表明深浅部圈层的耦合演化,但是下地幔结构演化一般会滞后于浅表。在微幔块推挤之下,地幔柱优先沿着地幔底部热化学结构的边缘形成,且有时会被推至热化学结构的内部。地幔柱上升至浅部后,能够导致岩石圈弱化甚至裂解或板块边界跃迁,形成微板块。因此,地幔底部LLSVP不是稳定或静止的,而是与微板块动态协同演化,并通过地幔柱与浅表板块边界发生遥相关,从而控制微板块生成场所。     

热-水-力作用下圆孔花岗岩的动态损伤特征及结构模型

基于深层地热能开采时储能区井筒围岩所处的工程环境,采用高温加热、不同温度水浸泡、加热-循环次数及径向冲击加载的方法模拟井筒围岩经历的高温、遇水、循环采热及热冲击等造成的动力扰动等物理力学条件。同时,以不同内孔直径的同心圆孔岩样模拟深层地热井,采用分离式霍普金森压杆试验系统开展热-水-力作用下圆孔花岗岩的动态力学试验,并结合VIC-3D非接触应变测量及数值模拟分析技术监测冲击过程中圆孔岩样裂隙萌发、形成的历程和表面应变演化的规律,揭示热-水-力作用下圆孔岩样的动态损伤破坏机制。研究结果表明：径向冲击荷载作用下圆孔花岗岩先后经历弹性变形、塑性变形、结构失稳破坏3个典型阶段;内孔直径、加热温度、浸水温度、加热-浸水循环次数4因素都弱化了圆孔花岗岩抗外界荷载的能力,但未改变其整体的变形演化规律;圆孔花岗岩的破坏模式是动态拉伸破坏,先沿冲击方向由内孔壁向岩样外壁,再垂直冲击方向由岩样外壁向内孔壁萌发、贯通裂纹,形成近垂直的两组破裂面。最后,基于圆孔花岗岩的损伤变形特征及历程,在一定假设基础上,建立动态损伤结构模型,推演了结构方程,并结合试验结果确定了方程参数,通过对比分析发现,理论拟合曲线与试验...     

Mechanism of Structure Variations at Rifted Margins in the Central Segment of South Atlantic: Insights from Numerical Modeling

Rifted margins in the central South Atlantic portray spatial variability in terms of preserved width and thickness, which relates to complex rift-related fault activities. However, there is still a lack of systematic and quantitative explanations for the causes of the variations that are observed along the paired rifts. To elucidate this issue, 2D viscous-plastic thermomechanical numerical models are applied to capture the behavior of deformation, in which we investigate the effects of extensional rate, crustal strength and thickness on crust-mantle coupling, and timing of transition from rifting to breakup. Our numerical experiments demonstrate that crust-mantle decoupling accounts for crustal hyperextension, and that incorporating moderate-intensity rheology into lower crust may yield insights into the hyper-extended crust and asymmetric architecture observed in the central South Atlantic. The results also suggest that undulations in lithospheric basement cause asymmetric mantle upwelling. The lower crust of fold belts takes priority to be thermally weakened over craton and induces rift migration simultaneously. A new mechanism for the formation of failed rift is described, where the mechanical decoupling derived from thermally weakened lower crust gives access to dual rift migration. These results reinforce the interpretation on how crustal rheology shapes margins architectures and highlight the first-order effects of crust-mantle coupling.     

应用于L波段风廓线雷达的2 kW限幅器设计

边界层风廓线雷达为气象预报服务提供垂直性、连续性、高时空分辨率的探测数据。针对在全国各雷达站集中式雷达运行过程中限幅器频繁烧毁导致雷达停止工作的问题,设计了一种L波段2kW高功率无源限幅器。此款限幅器在电路前端引入双耦合电缆和功率电阻,增加了散热能力,一分为二限幅,使通过PIN二极管的功率降低一倍,电路后端多级二极管并联,逐级平缓限幅;具有承受功率高,限幅电平低,插入损耗小的特点。利用ADS软件进行仿真,并对限幅器进行了测试,结果表明各项指标满足设计要求。此款限幅器在边界层风廓线雷达上的应用,极大地提高了雷达稳定性。     

青藏高原东缘—扬子特提斯构造域深部结构与地壳形变研究

青藏高原东缘和扬子西缘的构造带是中国特提斯构造域的重要组成部分,该构造域受欧亚板块与印度板块陆—陆碰撞、高原隆升、块体裂解或拼接挤压等强烈构造活动的影响,记录和保存了多期次的特提斯构造演化历史痕迹。同时,该研究区域也是中国西部地区地壳形变最强烈的地区之一,其浅表形变特征与深部构造之间存在怎样的关联和制约机制是目前国际地球科学的一个研究热点。本研究依据作者十多年来持续在该区域开展的地质—地球物理研究,通过深部地球物理多参数结构成像、沉积盆地分析、地壳形变和强震孕育机制等综合对比分析,发现在青藏高原东缘的下地壳存在低速和高泊松比异常带,该异常体与来自青藏高原上涌的软流圈热物质汇聚,导致从扬子西缘到青藏高原的下地壳和上地幔的深部结构发生显著变化。沿着龙门山断裂带,中、下地壳存在交叠相间的低速(高泊松比)和高速(低泊松比)区域,这些深部结构分布特征与地表形变及前陆盆地隆坳格局具有较好的一致性。基于上述认识,提出了青藏高原东缘—扬子板块的深部接触模式及其相应的盆山耦合关系,阐明了板块碰撞—耦合的深部动力学过程对剧烈地壳形变、盆地隆坳格局和强震诱发的制约关系。本研究成果将为深入认识青藏高原东缘高原急剧隆升、盆地基底结构与隆拗格局,以及强烈地壳形变的深部动力学机制提供参考信息。