基于理想点-可拓云模型的隧道围岩稳定性评价

为了对隧道围岩稳定性进行准确评价并解决研究过程中存在的模糊性和随机性以及评价指标不相容的问题,减少单一的主客观赋权法所带来的误差,将云模型引入可拓理论中,利用可拓理论能够实现矛盾问题向相容问题转化的特点和云模型具有处理事物双重不确定性的优势并结合理想点组合赋权法对隧道围岩进行稳定性评价研究。通过文献调研、数据统计的方式,选取具有代表性的6个指标组成隧道围岩评价指标体系,用理想点法赋予评价指标组合权重,并构建可拓云模型对隧道围岩稳定性进行综合评价。通过将此方法应用于工程实例并与其他方法对比,结果表明：基于理想点-可拓云模型的隧道围岩稳定性评价方法能够减少评价过程中存在的不确定性问题,克服单一赋权方法的不足,具有良好的适用性,可以应用于实际工程中。     

高烈度区双洞隧道洞口段地震响应分析与振动台模型试验

强震作用下,隧道地下结构损伤严重,其中洞口段更是抗震的薄弱环节,但其影响规律和特点目前尚缺乏系统的计算分析和研究。以强震区扯羊隧道为例,首先采用FLAC^3D对双洞隧道洞口段地震响应进行分析,再用模型试验的方法进行洞口段模型试验,测量其地震响应应变规律,观测其震害发展情况,并与数值分析结果相互验证。结果表明:地震作用下,隧道仰拱横向位移较大;随着隧道埋深增加,内力逐步增加,其中墙角部位内力较大;洞口处隧道围岩在强震作用下会产生贯通性裂缝,影响隧道洞口稳定;地震作用下的明洞与暗洞交接处内力较大;地震作用下的双洞隧道之间存在较强的动力相互作用。     

分布埋入式光纤与隧道衬砌耦合性能试验及应用

为探究埋入式光纤与隧道衬砌的耦合性能,分别从理论与试验两个方面进行研究,并在实际工程中进行了验证。构建了光纤、中间体和基体结构力学分析模型,进行光纤应变传递机制理论分析,计算了光纤应变传递效率;使用钢筋混凝土梁模拟隧道衬砌,进行了2组不同加载速率的试验。其中,在同一根梁内(同一工况)设计6种光纤的布设方式,以位移控制的方式在梁跨中部位进行单点多级加载,使用BOFDA(布里渊散射光频域分析)技术分别对6条光纤进行监测。试验结果表明:6条光纤均可以有效监测梁从开始加载至钢筋开始屈服阶段,光纤与梁耦合性最好;钢筋开始屈服直至梁破坏阶段,光纤应变不再增加甚至减小或呈现出光纤断裂的状态,此过程光纤与梁耦合性较差;除开槽埋入式光纤的有效监测应变差为3 000×10^-6外,其余布设方式光纤有效监测应变差为2 000×10^-6;光纤在长距离(>>146 mm)埋入式布设情况下可认为其应变传递效率接近100%,2组不同试验结果呈现相似规律。在北京市新机场线地铁暗挖隧道CRD工法区间进行了工程应用研究,监测结果表明分布埋入式光纤布设工艺是可行的,可为...     

基于Web的地下管网溯源分析算法研究

利用.NetCore的并行运算机制,实现基于Web的高性能地下管网的溯源分析功能.实践证明当某点管网数据监测出现异常时,能够快速分析所有可能流经该点源头管网信息,并在地图上动态展示分析结果,达到快速定位查询管网源头的区域管理目的.     

一种利用Nyström离散与FFT快速褶积的散射地震波并行计算方法

利用数值方法解Lippermann-Schwinger (L-S)方程的主要困难在于系数矩阵存储和线性方程组求解.这主要是因为L-S方程的积分部分是一个空间褶积,在离散后将导致一个满秩矩阵,进而形成一个大型或超大型代数方程组.因此,在利用L-S解决地震波散射问题时,一般是利用散射级数法而非数值方法.然而,散射级数法的计算精度和收敛性强烈地依赖于速度扰动的强度,而克服这种依赖性的一个可能的途径就是对现有的数值方法进行改进或是建立新的数值求解方案.在这种思想指导下,首先对L-S方程进行改写,得到一个与原L-S方程等价的积分方程(等价L-S方程).然后,对等价L-S方程进行逐点归一化处理,并利用Nyström法对经归一化处理的等价L-S方程(归一化等价L-S方程)进行离散,并用FFT计算空间褶积.之所以这样选择是由于归一化等价L-S方程经Nyström法离散生成的系数阵为一个Toeplitz阵,可利用其Toeplitz性质降低存储空间;而FFT可以将矩矢空间褶积转化为乘积,且积分核部分只要计算一次即可.进一步,为节约正演计算时间,设计了进程级和线程级相结合的MPI+OpenMP并行模式.数值试验表明,与传统的积分方程数值算法相比,利用等价L-S方程、Nyström离散和FFT快速褶积的计算方案可极大地降低存储需求,进而在保证精度的同时提高计算效率.

     

PLS法与隧道围岩稳定性分类

围岩稳定性是隧道工程支护设计的基本参数,采用偏最小二乘回归法,综合考虑多因变量影响因素,对隧道围岩稳定性分类进行回归建模,并取得较好效果.     

Stability analysis of tunnels driven in jointed rocks by means of a homogenized limit analysis approach

Conceived as a potential alternative to the classical design methods employed for analyzing the stability of underground structures driven in jointed rocks, the homogenization approach stems from the heuristic idea that, from a macroscopic point of view, a rock mass cut by a network of joints may be perceived as a homogenized continuum. The strength properties of the latter can be theoretically obtained from the failure conditions of its individual constituents: rock matrix and joint interfaces. At the material level, the limit analysis reasoning is used in the context of homogenization to formulate the homogenized strength criterion of a jointed rock mass in the particular situation of a single set of parallel joints. As it could be expected, the obtained closed‐form expressions show the strength anisotropy induced by joint preferential orientation. The support functions (π functions) associated with the homogenized strength criterion are also determined in both plane strain and three‐dimensional cases. This criterion is then applied to the investigation of stability analysis of a tunnel excavated in a jointed rock mass. Upper bounds estimated of the stability factor are derived from the implementation of the kinematic approach directly on the homogenized underground structure. Finally, the approach is applied to analyze and discuss the collapse of the Pinheiros subway station (São Paulo, Brazil). Copyright © 2014 John Wiley & Sons, Ltd.     

A new numerical approach to the hyperstatic reaction method for segmental tunnel linings

This paper proposes a numerical approach to the hyperstatic reaction method (HRM) for the analysis of segmental tunnel linings. The influence of segmental joints has been considered directly using a fixity ratio that is determined on the basis of the rotational stiffness. The parameters necessary for the calculation are presented. A specific implementation has been developed using a FEM framework. This code is able to consider the three‐dimensional (3D) effect of segment joints in successive rings on the tunnel lining behaviour. The present HRM allows one to take an arbitrary distribution of segment joints along the tunnel boundary into consideration. In addition, the rotational stiffness of segment joints has been simulated using nonlinear behaviour, as it is closer to the true behaviour of a joint than linear or bilinear behaviour. The numerical results of three hypotheses on ring interaction, which allow the 3D effect of a segmental tunnel lining to be taken into account, have been compared with data obtained from the shield‐driven tunnel of the Bologna–Florence high‐speed railway line project. The numerical results presented in the paper show that the proposed HRM can be used to effectively estimate the behaviour of a segmental tunnel lining. Copyright © 2014 John Wiley & Sons, Ltd.     

Parallelized combined finite-discrete element (FDEM) procedure using multi-GPU with CUDA

This paper focuses on the efficiency of finite discrete element method (FDEM) algorithmic procedures in massive computers and analyzes the time-consuming part of contact detection and interaction computations in the numerical solution. A detailed operable GPU parallel procedure was designed for the element node force calculation, contact detection, and contact interaction with thread allocation and data access based on the CUDA computing. The emphasis is on the parallel optimization of time-consuming contact detection based on load balance and GPU architecture. A CUDA FDEM parallel program was developed with the overall speedup ratio over 53 times after the fracture from the efficiency and fidelity performance test of models of in situ stress, UCS, and BD simulations in Intel i7-7700K CPU and the NVIDIA TITAN Z GPU. The CUDA FDEM parallel computing improves the computational efficiency significantly compared with the CPU-based ones with the same reliability, providing conditions for achieving larger-scale simulations of fracture.     

Cracking risk of partially saturated porous media—Part I: Microporoelasticity model

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. This is the motivation to develop a thermodynamics‐based microporoelasticity model for the assessment of cracking risk in partially saturated porous geomaterials. The study refers to 3D representative volume elements of porous media, including a two‐scale double‐porosity material with a pore network comprising (at the mesoscale) 3D mesocracks in the form of oblate spheroids, and (at the microscale) spherical micropores of different sizes. Surface tensions prevailing in all interfaces between solid, liquid, and gaseous matters are taken into account. To establish a thermodynamics‐based crack propagation criterion for a two‐scale double‐porosity material, the potential energy of the solid is derived, accounting—in particular—for mesocrack geometry changes (main original contribution) and for effective micropore pressures, which depend (due to surface tensions) on the pore radius. Differentiating the potential energy with respect to crack density parameter yields the thermodynamical driving force for crack propagation, which is shown to be governed by an effective macrostrain. It is found that drying‐related stresses in partially saturated mesocracks reduce the cracking risk. The drying‐related effective underpressures in spherical micropores, in turn, result in a tensile eigenstress of the matrix in which the mesocracks are embedded. This way, micropores increase the mesocracking risk. Model application to the assessment of cracking risk during drying of argillite is the topic of the companion paper (Part II). Copyright © 2009 John Wiley & Sons, Ltd.