Material point method simulation of hydro-mechanical behaviour in twophase porous geomaterials: A state-of-the-art review

The material point method(MPM)has been gaining increasing popularity as an appropriate approach to the solution of coupled hydro-mechanical problems involving large deformation.In this paper,we survey the current state-of-the-art in the MPM simulation of hydro-mechanical behaviour in two-phase porous geomaterials.The review covers the recent advances and developments in the MPM and their extensions to capture the coupled hydro-mechanical problems involving large deformations.The focus of this review is aiming at providing a clear picture of what has or has not been developed or implemented for simulating two-phase coupled large deformation problems,which will provide some direct reference for both practitioners and researchers.

Numerical Simulation of Surrounding Rock Deformation and Grouting Reinforcement of Cross-Fault Tunnel under Different Excavation Methods

Tunnel construction is susceptible to accidents such as loosening, deformation, collapse, and water inrush, especiallyunder complex geological conditions like dense fault areas. These accidents can cause instability and damageto the tunnel. As a result, it is essential to conduct research on tunnel construction and grouting reinforcementtechnology in fault fracture zones to address these issues and ensure the safety of tunnel excavation projects. Thisstudy utilized the Xianglushan cross-fault tunnel to conduct a comprehensive analysis on the construction, support,and reinforcement of a tunnel crossing a fault fracture zone using the three-dimensional finite element numericalmethod. The study yielded the following research conclusions: The excavation conditions of the cross-fault tunnelarray were analyzed to determine the optimal construction method for excavation while controlling deformationand stress in the surrounding rock. The middle partition method (CD method) was found to be the most suitable.Additionally, the effects of advanced reinforcement grouting on the cross-fault fracture zone tunnel were studied,and the optimal combination of grouting reinforcement range (140°) and grouting thickness (1m) was determined.The stress and deformation data obtained fromon-site monitoring of the surrounding rock was slightly lower thanthe numerical simulation results. However, the change trend of both sets of data was found to be consistent. Theseresearch findings provide technical analysis and data support for the construction and design of cross-fault tunnels.

EXTENSION AND APPLICATION OF NEWTON'S METHOD IN NONLINEAR OSCILLATION THEORY

In this paper we suggest and prove that Newton's method may calculate the asymptotic analytic periodic solution of strong and weak nonlinear nonautonomous systems, so that a new analytic method is offered for studying strong and weak nonlinear oscillation systems. On the strength of the need of our method, we discuss the existence and calculation of the periodic solution of the second order nonhomogeneous linear periodic system. Besides, we investigate the application of Newton's method to quasi-linear systems. The periodic solution of Duffing equation is calculated by means of our method.

Modeling pipe-soil interaction under vertical downward relative offset using B-spline material point method

To analyze the pipeline response under permanent ground deformation,the evolution of resistance acting on the pipe during the vertical downward offset is an essential ingredient.However,the efficient simulation of pipe penetration into soil is challenging for the conventional finite element(FE)method due to the large deformation of the surrounding soils.In this study,the B-spline material point method(MPM)is employed to investigate the pipe-soil interaction during the downward movement of rigid pipes buried in medium and dense sand.To describe the density-and stress-dependent behaviors of sand,the J2-deformation type model with state-dependent dilatancy is adopted.The effectiveness of the model is demonstrated by element tests and biaxial compression tests.Afterwards,the pipe penetration process is simulated,and the numerical outcomes are compared with the physical model tests.The effects of pipe size and burial depth are investigated with an emphasis on the mobilization of the soil resistance and the failure mechanisms.The simulation results indicate that the bearing capacity formulas given in the guidelines can provide essentially reasonable estimates for the ultimate force acting on buried pipes,and the recommended value of yield displacement may be underestimated to a certain extent.

Multi-domain equivalent method for prediction of elastic modulus of complex fractured rock mass

Reliable estimation of deformation and failure behaviors of fractured rock mass is important for practical engineering design.This study proposes a multi-domain equivalent method for fracture network to estimate the deformation properties of complex fractured rock mass.It comprehends both the advantages of the discrete fracture network model and the equivalent continuum model to capture the features of discontinuities explicitly while reducing computational intensity.The complex fracture network is stochastically split into a number of subfracture networks according to the domain,length or angle.An analytical solution is derived to infer theoretically the relationship between the elastic moduli of the original complex fractured rock mass and the split subfractured rock masses by introducing a correction term based on the deformation superposition principle.Numerical simulations are conducted to determine the elastic moduli of split subfractured rock masses using universal distinct element code(UDEC),while the elastic modulus of the original model is estimated based on the currently proposed analytical relationship.The results show that the estimation accuracy with the current domainbased splitting model is far superior compared to those with the other two splitting models.Thus,the estimation method of elastic modulus of complex fractured rock mass based on domain splitting mode of fracture network is identified as the multi-domain equivalent method proposed in this paper.The reliability of this method is evaluated,and its high computational efficiency is demonstrated through exemplification with regard to different geometric configurations for stochastically artificial discrete fracture network.The proposed multi-domain equivalent method constructs the theoretical framework except for the regression analysis hypothesis compared to the density-reduced model equivalent method.

A dynamic large-deformation particle finite element method for geotechnical applications based on Abaqus

In this paper,the application of Abaqus-based particle finite element method(PFEM)is extended from static to dynamic large deformation.The PFEM is based on periodic mesh regeneration with Delaunay triangulation to avoid mesh distortion.Additional mesh smoothing and boundary node smoothing techniques are incorporated to improve the mesh quality and solution accuracy.The field variables are mapped from the old to the new mesh using the closest point projection method to minimize the mapping error.The procedures of the proposed Abaqus-based dynamic PFEM(Abaqus-DPFEM)analysis and its implementation in Abaqus are detailed.The accuracy and robustness of the proposed approach are examined via four illustrative numerical examples.The numerical results show a satisfactory agreement with published results and further confirm the applicability of the Abaqus-DPFEM to solving dynamic large-deformation problems in geotechnical engineering.

Modal and Thermal Analysis of a Modified Connecting Rod of an Internal Combustion Engine Using Finite Element Method

The connecting rod is one of the most important moving components in an internal combustion engine. The present work determined the possibility of using aluminium alloy 7075 material to design and manufacture a connecting rod for weight optimisation without losing the strength of the connecting rod. It considered modal and thermal analyses to investigate the suitability of the material for connecting rod design. The parameters that were considered under the modal analysis were: total deformation, and natural frequency, while the thermal analysis looked at the temperature distribution, total heat flux and directional heat flux of the four connecting rods made with titanium alloy, grey cast iron, structural steel and aluminium 7075 alloy respectively. The connecting rod was modelled using Autodesk inventor2017 software using the calculated parameters. The steady-state thermal analysis was used to determine the induced heat flux and directional heat flux. The study found that Aluminium 7075 alloy deformed more than the remaining three other materials but has superior qualities in terms of vibrational natural frequency, total heat flux and lightweight compared to structural steel, grey cast iron and titanium alloy.

"Meridian-sinew Acupuncture Method"Improves Deformity Joint Dysfunction in Rheumatoid Arthritis

Joint deformity and dysfunction are common and serious complications in the late stage of rheumatoid arthritis,which seriously affect the quality of life of patients.Traditional Chinese medicine(TCM)believes that joint deformity and dysfunction in some patients with rheumatoid arthritis are closely related to the apraxia of meridians and tendons due to enduring illness.Based on the theory of meridians and tendons circulation,using the local and nearby therapeutic effect of acupoints as the treatment method in clinical practice,we conducted penetration needling of Houxi,Baxie,Wailaogong as well as Ashi points of interphalangeal joints of both hands through bilateral Sanjian,and used the uniform reinforcing-reducing method to soothe tendons and meridians,thus effectively improving the dysfunction of deformed joints.

铝合金塑性成形粗晶现象的研究进展

铝合金经过热塑性变形或者塑性变形+热处理后,可能发生晶粒异常长大而形成粗晶。粗晶会大幅降低构件的力学性能、耐腐蚀性能和抗疲劳性能,因此必须得到有效控制。对铝合金塑性成形粗晶的研究主要集中在2×××、6×××、7×××系热处理可强化型铝合金以及5×××系热处理不可强化型铝合金,研究内容已经覆盖塑性成形工艺、热处理制度、微观组织演变、数值仿真模拟等方面。为了明晰铝合金塑性成形构件粗晶现象的研究现状,为相关研究提供参考,对铝合金塑性成形粗晶现象的研究进展进行了梳理。首先,分析了材料成分、变形温度、变形速度、固溶温度、固溶时间和模具结构对晶粒异常长大的影响;其次,从晶粒取向分布、应变诱导储能和第二相钉扎作用的角度总结了晶粒异常长大的机理;再次,归纳了包括改变坯料初始状态、合理的锻造技术和中间形变热处理等有效抑制粗晶产生的方法;最后对未来的研究方向进行了展望。

深部采矿岩石力学进展

随着煤炭开采日益向深部发展,深部采矿引发的围岩大变形破坏和强冲击动力灾害日益严峻。在深部高地应力、高地温、高渗透压、强采动、强流变及多场耦合的复杂地质力学环境下,深部采场的应力场特征、煤岩体破碎性质、岩层移动及能量的积聚释放规律等均发生了显著变化。针对深部采矿中的岩石力学问题,论述了笔者及团队在深部采煤方法、深部巷道破坏机理与围岩控制、深井热害与地热利用三大方向取得的进展,主要包括:(1)提出了平衡开采理论和实现平衡开采的110/N00工法,进行了千米深井现场工程应用;(2)构建了深部“非均压建井”模式,研发了实现深井稳定提升的SAP系统,形成了可大幅简化井巷工程量和提高矿井采出率的建井方法;(3)研发了多套适用于研究深部岩体在水、高温、高压、结构效应及多场耦合作用下发生宏观破坏的实验系统和可进行微观层面演算的超算系统,揭示了深部软岩大变形破坏机理及多尺度力学特性;(4)研制了深部岩体冲击型和应变型岩爆实验系统,阐述了深部岩体冲击能量沿开挖临空面瞬间释放的非线性动力学行为;(5)提出了深部巷道开挖补偿支护理论,进一步发展了具有高恒阻、高延伸率、强吸能和耐冲击超常力学特性的NPR支护材料和技术;(6)研发了模拟深部高温、高湿和高压环境下的岩体热力学实验系统,提出了热害治理和热能资源化利用方法,建立了深部热害治理与热能综合利用系统(HEMS)。相关研究成果已在深部开采领域得以应用,可为深部采矿面临的复杂岩石力学问题提供借鉴。

洞桩法地铁车站边桩结构受力机制的模型试验设计

为了方便学生更形象地理解洞桩法地铁车站边桩结构的受力机制,该文以广州地区某洞桩法地铁车站工程为依托,设计了一种简化的洞桩法车站边桩结构模型试验。基于该试验方法,演示分析了边桩结构在洞桩法车站开挖过程中的施工力学行为规律,并对比研究了不同桩型布置参数下的桩后土压力、边桩结构的力学和变形规律。通过该试验,为学生展示了洞桩法边桩结构模型的设计、制作、试验具体操作及数据监测的全过程,从而让学生更好地掌握洞桩法车站边桩结构力学行为演变规律,以及采用不同边桩布置型式的力学行为差异。

覆岩导水裂隙带发育高度动态演化规律研究

为了研究覆岩导水裂隙带发育高度动态演化规律,以小保当矿区2^(-2)煤层为研究对象,通过理论分析、相似模拟实验及实例验证的方法,对导水裂隙带动态发育高度进行研究;借助概率积分法预计其上部岩层移动变形的理论公式,给出了1种基于覆岩曲率变形的导水裂隙带动态发育高度预计方法。研究表明:导水裂隙带上部岩层下沉系数是1个与挠度及下部自由空间高度相关的分段函数,上部岩层曲率变形大小是决定导水裂隙带发育高度的关键因素;导水裂隙带发育高度与工作面推进长度相关;理论预计小保当矿区2^(-2)煤层导水裂隙带发育高度为160.8 m,相似模拟实验发育高度为155 m,现场实测发育高度为152.01~175.57 m。

盾构隧道下穿高铁路基变形影响规律及加固优化研究

为研究盾构隧道下穿高铁路基施工时对路基的影响规律,探究在隧道开挖过程中运用不同加固措施对既有高铁路基的影响程度。以西安市某地铁工程为依托,利用室内模型试验对比分析盾构隧道下穿高铁CFG桩复合路基过程中,无加固、超前管幕工法加固与地表袖阀管注浆加固3种隧道施工条件下,地表沉降值和沉降槽的空间分布规律;采用数值模拟分析验证超前管幕工法加固工况下盾构隧道下穿时,高铁路基、道床的位移变化规律。室内模型试验结果表明:盾构隧道施工中采用管幕工法加固后复合地基正上方的地表沉降最大值减小28.6%,采用袖阀管注浆加固后减小18.0%,并且采用两种加固措施后的CFG桩最大附加轴力均减小20%以上。因此,在隧道掘进过程中采用一定加固措施,可改善围岩稳定性,其中管幕工法加固效果更为显著。数值模拟分析表明:采取管幕工法加固施工与不采取加固措施相比,路基最大沉降量减少33.78%,道床最大沉降量减少45.08%。因此,管幕工法加固能够有效减小对既有高铁路基的影响。

基于均匀化理论与上限分析的膨胀土滑坡稳定性分析

在膨胀土和滑坡共同作用下,隧道洞口段施工更容易引发地表开裂甚至滑坡等工程灾害,在隧道内采用微型桩群防治滑坡比抗滑桩具有优势。本文基于均匀化理论与上限分析对某高速公路隧道洞口段膨胀土滑坡的稳定性进行计算,并评价微型桩群和削方卸载不同组合方式的处置效果,计算时将微型桩群及桩周土通过均匀化理论等效为符合摩尔库仑强度准则的等效加固体来,以此提高计算效率,最后通过对现场削方+微型桩群加固处置后的滑坡变形监测来验证计算的合理性,得出如下结论:相较于土的强度参数,等效加固体内摩擦角保持不变,黏聚力从26kPa提高到85.36kPa。处置前后的滑坡稳定性评价结果表明,不做处理时,滑坡滑动面从滑坡上方岩土交界面延续到隧道洞口前;仅采用微型桩群加固时,滑坡安全系数在1.17左右,滑动面从岩土交界面延续到隧道洞口后;同时采用微型桩群加固和削方卸载时,滑坡安全系数提高到1.26~1.28,滑动面上缘由土石交界面前移。现场变形监测表明地表变形与深层土体变形均不超过3mm,该措施能保障滑坡的稳定性,同时也验证了计算方法的合理性,可为同类工程提供参考。

基于颗粒特征与预设变形的人工砂土变形图像生成方法及应用

图像变形分析方法是当前土力学与岩土工程领域的一种重要变形测试技术,散斑图像作为其发展、应用和可靠性评价的常用工具,是否反映真实土体图像特征和变形分析可靠性尚有待解答。以福建标准砂为样本,基于土体图像特征分析及与散斑图像对比,提出了一种描述砂土图像和变形特征的人工图像生成方法,利用建立的4种变形坐标解析公式和生成序列变形图像,对国际代表性RG-DIC和PIVlab法可靠性进行评价。结果表明:散斑图像与真实砂土图像在纹理、圆度、色值等方面存在显著差异,提出的人工图像生成方法能有效控制土颗粒组分、圆度、色值等分布参数,反映了真实土体图像特征;生成的序列变形图像增加了时间变量和任意变形函数功能,为动态和复杂变形分析可靠性评价建立条件。散斑图像明显低估了RG-DIC和PIVlab法对土体变形的分析误差,原因为黑色背景、白色亮斑构成的纹理特征更具辨识性;RG-DIC法的分析准确性与稳定性明显优于PIVlab法,而不同变形条件下两种方法分析误差呈一致趋势,当剪应变≤10^(-3)时误差快速上升。研究方法与成果,为土体图像变形分析方法发展、应用和可靠性评价提供了重要支撑和参考。

考虑节点刚度随机性的焊接钢筋部品变形计算方法研究

为量化钢筋部品刚度,从受力机理角度出发,基于试验和理论分析,提出考虑刚度随机性的焊接钢筋部品变形计算方法。从受力机理的角度提出三向转动刚度独立力学求解模型,考虑钢筋直径的影响,设计了105组T形焊接钢筋试件及两类加载固定试验装置,开展分级加载试验;在此基础上,考虑焊接节点转动刚度随机性,提出基于蒙特卡洛随机抽样的焊接钢筋部品变形计算方法,并分别开展了6个焊接钢筋网片和6个钢筋块体试验进行验证。结果表明:钢筋焊接节点刚度与钢筋直径有关,当焊接条件相同时,钢筋直径越大,焊接节点刚度越大;同一钢筋直径组合的焊接节点转动刚度基本符合正态分布,3个方向的转动刚度R x>R z>R y;钢筋焊接节点刚度与焊接节点断面面积基本呈线性关系;竖筋和横筋连接节点刚度对钢筋部品变形影响显著,所采用的节点刚度取值方法计算的焊接钢筋部品变形偏差均在10%以内,与传统刚接和铰接方法相比,变形计算精度大幅提高,变形计算方法合理。

基于物质点和深度积分耦合模型的滑坡数值分析

深度积分算法可将滑坡沿地表滑动的三维模型化简为二维模型进行求解,通过减少控制方程未知量的个数以提升求解效率。物质点法(material point method,MPM)具有无网格法和有网格法的双重优势,模拟滑坡大变形问题时可避免网格畸变现象。采用深度积分耦合物质点法建立滑坡数值模型,给出算法实现具体流程,基于影响域改进的物质点法(influence domain material point method,IDMPM),针对两个典型无倾角底面光滑和有倾角底面不光滑滑坡算例进行基准测试。在计算精度方面,深度积分耦合物质点法模型能较好地预测远端距离、流速、深度等滑移特征参数;在计算效率方面,与常规物质点法求解格式相比,深度积分耦合物质点法模型可大幅度提高运行效率。该研究成果可为滑坡地质灾害破坏范围的分析预测、危害评估、应急抢险提供有效理论依据和时间保障。

三维可变形圆化多面体离散单元法

为真实模拟岩体的变形特性与运动形态,将圆化多面体离散元法与有限单元法结合,提出一种三维可变形圆化多面体离散单元法。该方法既能真实表征块体的不规则特征,又降低了接触判断的难度,同时能准确反映块体的变形特性。在求解切向接触力时,将接触判断对象从接触对简化为单元整体,显著提高了计算效率。为分析块体变形特性,在块体离散单元内部划分有限元网格,将最外层有限元网格作为最小接触单元。采用直接平均法,将接触力转化为等效节点力,并采用非线性有限单元法实现对单元变形特性的精确模拟,克服了圆化多面体无法反映单元变形的缺陷。通过5个算例论证了新方法在捕捉单元变形、运动形态以及其力学特征等方面的准确性和高效性。

个性化头部生物力学模型的开发及验证

本文提出了一种个性化头部生物力学模型建模方法,并对所生成模型进行了有效性验证。该方法以图斯特50th百分位头部生物力学模型为基础,基于目标模型头部CT数据,采用三维点云配准和自由网格变形法,实现具有详实脑组织结构的个性化头部生物力学模型的快速建模。通过重构经典尸体试验,发现该方法构建的个性化头部生物力学模型在运动学和生物力学响应方面与尸体试验结果高度一致,且与通过逆向工程方法开发的头部生物力学模型并无明显差异,验证了本文方法所开发模型的有效性。因此,本文方法可用于快速构建具有详细解剖学结构的个性化头部生物力学模型,为损伤生物力学、临床医学和法医鉴定等领域的数字化需求提供基础计算分析工具。

库车坳陷西部阿瓦特—博孜地区构造变形差异及数值模拟分析

库车坳陷西部的阿瓦特-博孜地区位于构造挤压转换区,沿走向构造变形差异显著。本文基于最新的三维地震资料和离散元数值模拟方法,对该区域进行了精细的地质建模,模拟了其形成演化过程,讨论了差异变形的控制因素。阿瓦特地区盐下沿克深断裂远距离逆冲推覆,形成上下叠置的两套沿侏罗系煤层滑脱的逆冲体系;盐上发育厚度较小的阿瓦特凹陷;盐构造以不整合型为主。博孜地区盐下在克深断裂以北为高角度基底卷入变形,以南为沿煤层的滑脱,整体为楔形叠瓦冲断构造;盐上发育厚度巨大的拜城凹陷;盐构造以整合型为主。模拟结果表明,温宿古隆起限制了中生界至古近系的沉积范围,阻碍了变形的向前传递,是造成阿瓦特-博孜地区变形差异性的最主要原因。先存盐底辟构造和南天山的差异挤压作用对该区域的差异变形也具有一定的控制作用。综合地震剖面解释和数值模拟结果的演化过程表明,该区域构造活动具有由弱到强的变化趋势,构造变形开始于早中新世吉迪克组沉积时期,在上新世早期康村组沉积之后发生显著的构造加速事件。