双向粘弹性节理岩体时域等效分析

以双向粘弹性节理岩体为研究对象,基于一种简单等效假定和自适应算法,并结合有限元方法,得到递推格式的等效本构模型和等效粘弹性场计算模型.通过算例,将等效解与基于ANSYS的非均质解进行比较.从计算精度和效率平衡的角度而言,等效模型的结果是令人满意的.     

基于时域精细算法的EFG方法及其在粘弹性问题中的应用

将时域精细算法与EFG方法相结合,求解粘弹性静、动力问题.通过离散时间段上的变量展开,将时空耦合的初边值问题转化为一系列递推形式的边值问题,然后利用EFG方法进行自适应计算,对非线性问题无需进行迭代求解.此外,通过面力耦合技术,将FE方法和EFG方法简明有效地结合起来,避免了界面上形函数的复杂化.数值算例给出令人满意的结果.     

节理岩体边坡的动力稳定性分析

阐述了动力离散元的原理及计算方法, 并将其用于一个受高烈度地震威胁的滑节理岩体边坡的动力响应分析与动力稳定性评价, 阐明了滑坡动力稳定性分析的必要性。     

节理岩体的分维特征及其工程地质意义

本文用分形理论,对节理网络分布和粗糙度曲线进行分形分析。结果表明,节理的空间分布和粗糙度曲线具自相似性;分维是一个表征岩体强度、岩体损伤程度、岩体质量及节理粗糙程度的几何参数,并可用于评价岩体稳定性、断裂体系活动性及划分风化卸荷带等。     

节理岩体体积节理数Jv的新计算公式

岩体体积节理数Jv的定义,按照《工程岩体分级标准》采用节理间距或频率较为恰当;本文研究了岩体内发育一组、二组、三组贯通性节理的情况下,Jv定义的适宜性,指出节理正交的情况下Jv的适宜性最好,推导出Jv的修正系数为K1=|sin(n1,n2[DD(-1][HT5SS]^[HT][][DD)])||cos(n1×n2,n3[DD(-1][HT5SS]^[HT][][DD)])|,研究表明：节理夹角越小,对Jv的影响越大;通过连通率等效途径,断续节理岩体的Jv计算公式为Jv=K1•∑[DD(][]i[DD)]Si•ni(i3)。运用该公式可实现工程岩体裂隙发育程度的科学评价。     

节理岩体断裂的损伤与分形效应

本文在岩石断裂的基础上,利用损伤的概念,建立起岩体断裂能与初始损伤的关系,并借助于岩石断裂能与断裂面分形维数的关系,给出了节理岩体的断裂能与断裂面分形维数之间的关系.本法为直接从岩石的断裂参数推出节理岩体断裂参数,提供了一条简捷的途径.     

节理岩体的爆破松动机理

利用应力波理论对节理岩体的爆破松动机理进行了探讨,同时利用动力有限元数值模拟方法验证了理论分析结果。研究表明,爆破荷载的加载阶段对岩体的挤压蓄能和快速卸载阶段的弹性回复是岩体产生后续松动位移的根本原因;节理岩体的爆破松动在开挖面的法向表现为岩体的剪切错动,在平行于开挖面方向表现为结构面的张开;拉格朗日坐标下的动力有限元数值方法,如ANSYS/LS-DYNA程序,能够较好地用来模拟节理岩体的爆破松动过程。     

基于有限元-时域分段自适应算法的围岩衬砌耦合渐进分析

以开挖时间、一衬时间、二衬时间划分时间区段,考虑围岩与衬砌为弹性/粘弹性介质,将围岩的“位移/应力释放不完全”、开挖等影响模拟为载荷的渐进释放,提出一个有限元-时域分段自适应算法,求解围岩衬砌耦合问题,并通过数值算例对算法进行了验证.所提方法给出了各时间区段的关联条件;可更准确描述各变量随时间的变化;同时将原时空耦合问题转化为一系列空间问题,并利用有限元方法递推求解;当步长变化时可通过自适应计算保证稳定的计算精度.此外,还建立了一个载荷渐进参数的反问题数值求解模型,并给出了相关算例.     

单节理岩体在真三轴压缩状态下强度特性的分析

本对单节理岩体在真三轴压缩状态下的强度特性进行了分析,重点讨论了中间主应力对单节理岩体强度特性的影响诸问题。研究认为,中间主应力对节理岩体的强度特性是有影响的,在一定条件下,甚至可以成为主要的影响因素。     

时域自适应算法求解弹性地基薄板的动力问题

为更精确地求解弹性地基薄板的动力响应,发展了一种分段时域自适应算法,通过变量在离散时段内的展开,将时空耦合的初边值问题转化为一系列递推的基于有限元(FEM)的空间问题求解,通过自适应计算保持稳定的计算精度。数值算例表明:本文解与解析解相比最大相对误差不超过3.59%;当步长较大时四阶Runge-Kutta法和Newmark法均失效,本文所提算法仍可得到满意的计算结果。     

可储能原则在一类岩爆问题中的研究和应用

从能量的角度对一类岩爆问题--应变型岩爆进行一些探讨.对此类岩爆能量的来源、能量的去处及能量释放的原因和大小进行了完整的论述.提出了岩体的可储能原则--岩体可存储的能量并非定值,而是随应力状态的变化而改变,当应力发生突变,岩体内实际存储的能量大于可储能时,多余的能量必将以各种形式瞬间释放,其中包括动能.根据此准则对简化了的一维、三维模型进行了详细的分析和计算,并结合离散元软件3DEC再现了三维卸载岩爆的破坏过程.     

岩石爆破开挖诱发振动的等效模拟方法

岩石爆破开挖过程中爆炸荷载的复杂性、岩石介质模型的多变性、本构方程的多样性、以及群孔爆破模拟计算工作量巨大等诸多因素使得准确模拟爆破诱发的围岩振动存在较大的困难,鉴于此提出一种等效数值模拟方法。根据爆破过程炮孔周围岩石破坏范围的空间分布特征,确定群孔起爆条件下爆炸荷载作用的等效弹性边界;通过分析炸药起爆后炮孔空腔动力膨胀、岩体裂纹扩展、炮孔堵塞物运动以及爆生气体的逸出,计算了等效边界上沿炮孔轴向变化的爆炸荷载。结合瀑布沟水电站1# 尾水隧洞爆破开挖监测,基于动力有限元法运用该等效技术模拟了围岩的质点振动速度。数值模拟与现场实测数据的对比表明该等效模拟技术用于岩石爆破开挖围岩响应的计算是合适的。计算结果同时还显示,等效边界邻近岩体力学参数的选取对计算结果有重要影响。     

粘弹性岩体中支护圆形洞室施工过程的解析分析

地下洞室的开挖与支护是逐步的连续过程。对具有流变效应的粘弹性岩体,流变时效与施工效应发生耦合,变形与时间相关。针对深埋圆形洞室的施工,用半径时变函数模拟断面开挖过程。当岩体模拟为任一粘弹性材料时,将方程进行拉普拉斯变换求得位移通解,逆变换后代入边界条件确定待定函数,最终得到用洞周面力表达的围岩应力、位移统一解。区分开挖与支护时段,将半径时变函数、洞周面力不同表达式代入,利用支护后围岩与弹性支护接触条件建立关于支护力的Volterra积分方程。当岩石模拟为Boltzmann粘弹模型时,代入材料参数可求解积分方程得到支护力的确切表达,并进一步求得开挖过程及任意时刻支护后应力、位移分段解析表达式。表达式和算例分析表明:加支护后的径向位移增长呈指数形式变化且最终稳定于某一数值。最终洞型相同时,采用不同断面开挖速度且挖完立即支护时,开挖较快的情况位移变化较剧烈,而支护后最终稳定位移较小;但是,相应支护阶段产生的位移较大,支护力也较大。文中给出的方法可用于计算圆形洞室半径任意开挖并加支护后的应力、位移,适用于任一粘弹模型岩体的施工分析。     

Combined adaptive meshing technique and characteristic-based split algorithm for viscous incompressible flow analysis

A combined characteristic-based split algorithm and an adaptive meshing technique for analyzing two-dimensional viscous incompressible flow are presented.The method uses the three-node triangular element with equal-order interpolation functions for all variables of the velocity components and pressure.The main advantage of the combined method is that it improves the solution accuracy by coupling an error estima- tion procedure to an adaptive meshing technique that generates small elements in regions with a large change in solution gradients,and at the same time,larger elements in the other regions.The performance of the combined procedure is evaluated by analyzing one test case of the flow past a cylinder,for their transient and steady-state flow behaviors.     

基于自适应遗传算法的三轴磁强计校准

针对MEMS磁强计的精度无法满足姿态测量系统的航向角测量要求,对磁强计的误差来源进行了模型化分析,设计了一种基于自适应遗传算法的空间椭球磁强计校准方法.首先,采取自适应遗传算法,对磁强计测量的原始数据进行空间椭球的拟合,用估计的参数进行刻度系数、软磁干扰、硬磁干扰与零位偏置的综合误差补偿.其次,利用最小二乘法求解出非正交轴的方向余弦,进行非正交误差和安装误差的补偿.最后,将该方法应用到某姿态测量系统中,分别用未补偿和补偿后的数据进行姿态测量实验,实验结果表明该方法准确计算出磁强计的误差参数,使补偿后的航向角精度提高了6.8倍.     

Sensitivity analysis based on Lanczos algorithm in structural dynamics

IntroductionTheproblemforcomputingeigenpairderivativeshasoccupiedmanyresearchersinthepastseveraldecades.SurveypaperbyAdelmanandHaftka[1]isanexcellentintroductiontothetopicandalsoindicatevariousdisciplinesofapplications.Thereasonwhysomanypeopleareinterestedinthisproblemisthatsensitivitiesplayaveryimportantroleinstructuraldynamics,optimalanalysisandcontrolsystemdesign ,andmayalsobeanexpensiveandtime_consumingtask.Researchonsensitivityofmultipleeigenvalueshasbeenanareaofsignificantinterestinrece…     

基于三维接触有限元分析的管片等效抗弯刚度和错台研究

针对盾构隧道管片接头等效抗弯刚度预测研究中,对梁-弹簧与三维模型整体等效性考虑的不充分,以三维接触有限元计算结果作为测量信息,借助基于挠度等效的反问题求解,提出了一种确定管片接头等效抗弯刚度的新方法;并利用不同轴力-偏心距组合下的反演结果,建立了基于Kriging代理模型的轴力-弯矩-等效抗弯刚度的非线性关系,提出了由此关联的管片结构非线性问题的数值求解方法。与三维有限元结果相比,所提方法可较为准确地预测管片结构的内力和变形,表现出良好的整体等效性。此外,借助三维接触有限元分析,进一步深入探讨了螺栓孔间隙与衬垫本构关系对管片纵向错台量的影响。     

一种基于自适应波峰检测的MEMS计步算法

针对微机电测量系统(MEMS)波峰检测计步算法和自相关分析计步算法仅利用单轴加速度和固定阈值对传感器姿态和运动状态变化适应性较差的问题,提出了一种自适应波峰检测算法。该算法将行人运动状态分为正常状态与非正常状态,根据行人每一步的最大整体加速度与运动状态的内在相关性,获取不同运动状态的波峰检测经验阈值,实现不同运动状态下的自适应计步。通过实验对比分析,自适应波峰检测算法在传感器不同姿态和行人不同运动状态下的计步正确率均可达到99%以上,而常规波峰检测算法和自相关分析算法对正常态的计步精度虽然达到97%和99%以上,但对非正常状态下的计步精度仅有70%和50%,无法适应行人运动状态的变化。结果表明:自适应波峰检测算法对MEMS传感器姿态和运动状态的变化适应性较强,能够实现传感器不同姿态和不同运动状态下的可靠性计步。另外,自适应波峰检测、常规波峰检测、自相关分析算法的时间运算效率分别为0.036 s、0.046 s、0.131 s,自适应波峰检测算法时间效率明显优于其他两种算法。     

An efficient goal‐based reduced order model approach for targeted adaptive observations

An efficient adjoint sensitivity technique for optimally collecting targeted observations is presented. The targeting technique incorporates dynamical information from the numerical model predictions to identify when, where and what types of observations would provide the greatest improvement to specific model forecasts at a future time. A functional (goal) is defined to measure what is considered important in modelling problems. The adjoint sensitivity technique is used to identify the impact of observations on the predictive accuracy of the functional, then placing the sensors at the locations with high impacts. The adaptive (goal) observation technique developed here has the following features: (i) over existing targeted observation techniques, its novelty lies in that the interpolation error of numerical results is introduced to the functional (goal), which ensures the measurements are a distance apart; (ii) the use of proper orthogonal decomposition (POD) and reduced order modelling for both the forward and backward simulations, thus reducing the computational cost; and (iii) the use of unstructured meshes. The targeted adaptive observation technique is developed here within an unstructured mesh finite element model (Fluidity). In this work, a POD reduced order modelling is used to form the reduced order forward model by projecting the original complex model from a high dimensional space onto a reduced order space. The reduced order adjoint model is then constructed directly from the reduced order forward model. This efficient adaptive observation technique has been validated with two test cases: a model of an ocean gyre and a model of 2D urban street canyon flows. Copyright © 2016 John Wiley & Sons, Ltd.