Finite-element modeling of a complex deep excavation in Shanghai

The excavation of the north square underground shopping center of Shanghai South Railway Station is a complex deep excavation using the top-down construction method. The excavation has a considerable size and is close to the operating Metro Lines. In order to predict the performance of the excavation more accurately, 3D finite-element analyses are conducted to simulate the construction of this complex excavation. The effects of the anisotropic soil stiffness, the adjacent excavation, and zone excavation on the wall deformation are investigated. It is shown that the numerical simulation with anisotropic soil stiffness yields a more reasonable prediction of the wall deflection than the case with isotropic soil stiffness. The deformation of the shared diaphragm wall between two excavations is influenced by the construction sequence of the two excavations. The zoned excavation can greatly reduce the diaphragm wall deformation. However, only the zoned excavation at the first excavation stage affects the deformation of the walls significantly. When the depth of the excavation increases, the zoned excavation has minor effect on the deformation of diaphragm walls.     

Variation with time of the factor of safety of slopes excavated in unsaturated soils

While the evolution with time of the factor of safety against slope failure is well understood for excavations performed under fully saturated conditions, this evolution has not been extensively studied for excavations carried out in unsaturated soils. The objective of this paper is to investigate the relationship between the factor of safety and time for an excavation performed in an unsaturated silty soil, employing the finite element method. A hypothetical boundary value problem is considered and two types of analyses performed; in the first one, unsaturated soil behaviour is modelled through appropriate constitutive and soil–water retention curve models, whereas in the second one full saturation is assumed. The effect that the saturated soil permeability, its variation with suction, the increase of apparent cohesion due to suction, the depth of the groundwater table and the hydraulic hysteresis have on the results of the unsaturated analysis is examined in a parametric study. The analyses results demonstrate that for unsaturated soils the factor of safety may increase with time, in contrast to what is commonly accepted to be the case in fully saturated soils. Furthermore, it is not possible to know in advance which one of the two types of analysis performed (i.e. the saturated or the unsaturated) will produce conservative results. It is, therefore, advisable to perform unsaturated analysis in geotechnical practice, when dealing with excavations in such soils.     

软土边坡的有限元分析

边坡是否稳定与边坡的应力、应变（位移）状态有十分密切的联系。传统的条分法由于自身的限制,计算结果有很大的限制。本文提出了采用有限元方法分析软土边坡土体的应力－应变状态以确定滑裂面空间位置的新方法。     

Hydrous partial melting in the sheeted dike complex at fast spreading ridges: experimental and natural observations

In ophiolites and in present-day oceanic crust formed at fast spreading ridges, oceanic plagiogranites are commonly observed at, or close to the base of the sheeted dike complex. They can be produced either by differentiation of mafic melts, or by hydrous partial melting of the hydrothermally altered sheeted dikes. In addition, the hydrothermally altered base of the sheeted dike complex, which is often infiltrated by plagiogranitic veins, is usually recrystallized into granoblastic dikes that are commonly interpreted as a result of prograde granulitic metamorphism. To test the anatectic origin of oceanic plagiogranites, we performed melting experiments on a natural hydrothermally altered dike, under conditions that match those prevailing at the base of the sheeted dike complex. All generated melts are water saturated, transitional between tholeiitic and calc-alkaline, and match the compositions of oceanic plagiogranites observed close to the base of the sheeted dike complex. Newly crystallized clinopyroxene and plagioclase have compositions that are characteristic of the same minerals in granoblastic dikes. Published silicic melt compositions obtained in classical MORB fractionation experiments also broadly match the compositions of oceanic plagiogranites; however, the compositions of the coexisting experimental minerals significantly deviate from those of the granoblastic dikes. Our results demonstrate that hydrous partial melting is a likely common process in the root zone of the sheeted dike complex, starting at temperatures exceeding 850°C. The newly formed melt can either crystallize to form oceanic plagiogranites or may be recycled within the melt lens resulting in hybridized and contaminated MORB melts. It represents the main MORB crustal contamination process. The residue after the partial melting event is represented by the granoblastic dikes. Our results support a model with a dynamic melt lens that has the potential to trigger hydrous partial melting reactions in the previously hydrothermally altered sheeted dikes. A new thermometer using the Al content of clinopyroxene is also elaborated.     

Response to harmonic wave excitation of finite or infinite elastic plates on a homogeneous or layered half-space

The soil–structure interaction of elastic plates on homogeneous or layered soils excited by horizontally propagating waves is analysed. Large plates are modelled by a combined finite-element boundary-element method (FEBEM), whereas the response of infinitely long plates is calculated by a numerical integration in the frequency–wavenumber domain. The finite-element boundary-element method yields the complete soil–plate transfer function of frequency and distance whereas the frequency–wavenumber solution of the infinitely long plate can serve as an approximation for long distances on a finitely long plate. The soil–plate transfer function starts to decrease strongly at the coincidence frequency, where the bending stiffness equals the plate inertia. A strong decrease follows at mid frequencies and a strong reduction of less than 0.1 of the ground vibration is reached at high frequencies. Rules for the characteristic frequencies are derived from the numerical results clearly indicating the strongest influence of the soil stiffness and the weaker influence of the bending stiffness of the plate. The influence of the mass, length and width of the plate are shown to be limited in case of realistic parameters, but it should be noted that the reduction effects are less effective for layered soils and for nearer observation points.     

Damage Estimation in Masonry Buildings Based on Excavation-Induced Ground Movements

Excavation-induced ground movements and the resulting damages to adjacent structures and facilities is a source of concern for excavation projects in urban areas. The concern will be even higher if the adjacent structure is old or has low strength parameters like masonry building. Frame distortion and crack generation are predictors of building damage resulted from excavation-induced ground movements, which pose challenges to projects involving excavations. This study is aimed to investigate the relation between excavation-induced ground movements and damage probability of buildings in excavation affected distance. The main focus of this paper is on masonry buildings and excavations stabilized using soil nail wall method. To achieve this purpose, 21 masonry buildings adjacent to 12 excavation projects were studied. Parametric studies were performed by developing 3D FE models of brick walls and excavations stabilized using soil nail wall. Finally, probability evaluations were conducted to analyze the outputs obtained from case studies. Based on the obtained results, simple charts were established to estimate the damage of masonry structures in excavation affected distance with two key parameters including “Displacement Ratio” and “Normalized Distance”. The results also highlight the effects of building distance from excavation wall on its damage probability.

     

Groundwater and geotechnical aspects of deep excavations in Hong Kong

Consideration of groundwater is a key element in almost every construction project. The design of deep excavations for basements or underground railway station concourses below the water table require that the water pressures are taken into account. Whilst these can be considered to be hydrostatic in soil, the decreasing permeability of rock with depth and the fact that groundwater flow is invariably along discrete fractures means that the water pressure is unlikely to be hydrostatic at depth.

Groundwater control for deep excavations can be achieved by a number of methods such as grouting, pumping or structural walls or a combination of these. For tunnelling projects grouting is extensively used, but the development of sophisticated tunnelling machines has led in many cases to the demise of compressed air as a means of groundwater control.     

基坑开挖数值分析中土体硬化模型参数的试验研究

土体硬化模型已成为基坑数值分析中最常用的本构模型之一,其应用的关键是计算参数的确定。首先采用薄壁取土器现场取得土样,通过基于GDS的三轴固结排水剪切试验、三轴固结排水卸载-再加载试验及标准固结试验获得了上海典型土层土体的三轴试验应力-应变曲线和固结试验荷载-应变关系曲线。然后根据曲线确定了各土层的HS模型参数,建立了各土层HS模型参数中切线模量 、割线模量 、加载模量 的比例关系,并探讨了这些模量与压缩模量 之间的比例关系。最后将各层的模型参数与模量间的比例关系与国内外其他软土地区的结果进行了比较。试验结果可作为上海地区及其他软土地区基坑工程数值分析中确定HS模型参数时的参考。     

上海地区基坑开挖数值分析中土体 HS-Small模型参数的研究

数值分析方法已成为分析敏感环境下基坑开挖对环境影响的最重要手段,其关键是选择合适的土体本构模型和计算参数。HS-Small模型可以同时考虑土体剪切硬化和压缩硬化,而且还可以考虑土体剪切模量在微小应变范围内随应变衰减的行为,因此,HS-Small模型在基坑开挖数值分析中具有很好的适用性。但HS-Small模型的参数较多,不仅包含了HS模型的参数,还包含了小应变参数,因此,获取完整的HS-Small模型参数存在一定困难。通过上海地区5个基坑开挖数值分析的实例对基于室内试验得出的上海典型土层HS模型参数进行了检验。通过参数的敏感性分析,对HS-Small模型的部分参数的确定方法进行了修正。在此基础上,采用反分析法确定了敏感性最强的小应变参数,从而初步完整地获取了上海典型土层土体HS-Small模型参数。采用文中确定的上海典型土层HS-Small模型参数对该地区若干典型基坑工程实例进行计算,结果表明,计算值与实测值吻合得较好,从而验证了确定的上海典型土层HS-Small模型参数的适用性。     

Reliability analysis of soil nail walls

In India, soil nail walls are being extensively used for supporting vertical excavations below ground level to accommodate construction of one-or two-storied basements. Generally, the depth of excavations for basement construction ranges from 10 m to 15 m. For such large depth of excavation, variability of in-situ soil properties has significant influence on the stability of the soil nail walls. In the present study, using reliability analysis, an attempt is made to study the influence of variability of in-situ soil properties on the stability of soil nail walls. For better understanding, a case of 10 m high soil nail wall constructed to support a vertical cut is considered for the study and its stability is evaluated for various failure modes. Additionally, the influence of correlation among soil parameters on soil nail wall stability is assessed. In-situ soil friction angle and correlation between in-situ soil cohesion and angle of friction are found to influence soil nail wall stability significantly. In general, reliability analysis provided a better insight into the assessment of stability of soil nail wall.     

敏感环境下基坑数值分析中土体本构模型的选择

数值分析已成为敏感环境下基坑工程分析的最重要手段,其关键是选择合适的土体本构模型和计算参数。在分析了岩土数值分析中常用土体本构模型特点的基础上,通过算例较系统地对比了各类模型在基坑开挖数值分析中的适用性。敏感环境下的基坑工程需重点关注墙后土体的变形,从满足工程需要和方便实用的角度出发,建议采用能考虑黏土的塑性和应变硬化特征、能区分加荷和卸荷且刚度依赖于应力水平的硬化类弹塑性模型,如MCC模型和HS模型进行分析。具体工程实例的分析,表明了硬化类弹塑性模型在敏感环境下基坑开挖数值分析中的适用性。     

狭窄基坑平动模式刚性挡墙被动土压力分析

对于地铁车站、地下管道沟槽等狭窄基坑,其被动区土体宽度有限,不满足半无限体的假定,采用经典的库仑、朗肯土压力理论计算挡墙被动土压力是不合适的。首先建立了无黏性土中狭窄基坑刚性挡墙的有限元分析模型,研究了挡墙相对平移时不同宽度土体的被动滑裂面的分布规律;借鉴库仑平面土楔假定,建立了狭窄基坑刚性平动挡墙被动土压力的理论计算模型,推导了被动极限状态下滑裂面倾角及被动土压力系数的解析公式;再采用水平薄层单元法,得到了被动土压力分布、土压力合力作用点高度的理论公式。结合算例,深入研究了这种工程背景下挡墙被动滑裂面倾角的影响因素,以及被动土压力合力、土压力分布及合力作用点位置与经典库仑土压力理论的差别,与数值计算结果的对比验证了该理论方法的合理性。研究发现,当被动区土体宽度小于满足半无限体的临界值、且墙土摩擦角大于0时,被动滑裂面倾角大于传统库仑被动滑裂面倾角,被动土压力大于经典库仑解,合力作用点高度则小于库仑解,且基坑越窄,墙土摩擦角越大,其差别越大。     

The Mammoth Peak sheeted complex,Tuolumne batholith,Sierra Nevada,California: a record of initial growth or late thermal contraction in a magma chamber?

The Mammoth Peak sheeted intrusive complex formed in the interior of a ~7–10 km deep magma chamber, specifically in the Half Dome granodiorite of the Tuolumne batholith, central Sierra Nevada, CA (USA). The sheets consist of fractionated melts with accumulated hornblende, biotite, magnetite, titanite, apatite, and zircon. The accumulation, especially of titanite, had a profound effect on minor and trace elements (Nb, Ta, Ti, REE, U, Th, P, Zr, Hf, etc.), increasing their contents up to five to six times. Our thermal–mechanical modeling using the finite element method shows that cooling-generated tensile stresses resulted in the inward propagation of two perpendicular sets of dilational cracks in the host granodiorite. We interpret the sheeted complex to have formed by a crack-seal mechanism in a high strength, crystal-rich mush, whereby outward younging pulses of fractionated magma were injected into these syn-magmatic cracks at the margin of an active magma chamber. Thermal–mechanical instabilities developed after the assembly of the sheeted complex, which was then overprinted by late ~NW–SE magmatic foliation. This case example provides a cautionary note regarding the interpretation that sheeted zones in large granitoid plutons imply a diking mechanism of growth because the sheeted/dike complexes in plutons (1) may display inverse growth directions from the growth of the overall intrusive sequence; (2) need not record initial chamber construction and instead may reflect late pulsing of magma within an already constructed magma chamber; (3) have an overprinting magmatic fabric indicating the continued presence of melt after construction of sheeted complexes and thus a prolonged thermal history as compared to dikes; and (4) because the scale of the observed sheeted complexes may be small (<1%) in comparison to large homogenous parts of plutons, in which there is no evidence for sheeting or diking. Thus, where extensive dike complexes in plutons are absent, such as in much of the Tuolumne batholith, the application of an incremental diking model to explain chamber construction is at best speculative.     

超大跨度深埋地下结构围岩压力计算研究

围岩压力是地下结构设计与施工的重要依据。新建京张城际铁路八达岭车站拟定备选方案之一的单拱超大跨设计方案中,最大的跨度近45 m。对于类似超大跨度的深埋地下结构,传统的围岩压力理论与经验计算方法都存在着明显的局限性,已经不再适用。而结合现场勘查资料和室内外试验成果的数值模拟方法,尽管可以得到相对精确的围岩压力计算结果,但存在着计算过程复杂、计算效率低下的不足。因此,在计算分析了多种传统围岩压力计算方法与数值模拟方法获得的围岩压力所存在差异的基础上,对Q系统的围岩压力公式引入跨度项进行了修正,从而获得了能够快速估算超大跨度深埋地下结构围岩压力的估算方法。最后对不同跨度工程的围岩压力进行估算,并与监测值进行比较,对比结果良好。该方法不仅可以为工程设计人员及现场工作人员提供设计计算依据,还可为相关研究人员提供理论参考。     

An approximate analytical solution for calculating ground surface settlements due to diaphragm walling

This paper focuses on the derivation of an approximate analytical solution to predict ground surface settlements along the centre-line perpendicular to a slurry-supported diaphragm wall panel. An analytical model of trench excavation is constructed by applying linearly-distributed total earth pressure changes on the sidewalls and the base of the trench in a homogeneous, isotropic, elastic half-space. By applying the principle of superposition appropriately, an explicit solution to calculate ground surface settlements is derived by integrating Mindlin’s displacement formulae. The validity of the solution is verified by three-dimensional finite-element analyses and case studies. Based on the results obtained from a parametric study using the newly derived solution, an algebraic formula is obtained as a first approximation for calculating ground surface settlements.     

临近基坑既有建筑物地基的抗剪强度指标取值问题

基坑开挖不仅改变了半无限空间体边界条件,卸载时还会引起周围地层移动,导致既有建筑物地基承载力的附加损失。在临近基坑既有建筑物地基承载力计算中,其抗剪强度参数的选取有别于斜坡上地基计算参数的常规取值方法。通过探讨临近基坑既有建筑物地基抗剪强度指标的取值问题,建议以地基土的完全软化强度与残余抗剪强度的平均值作为临近基坑既有建筑物地基抗剪强度指标,并据此估算了临近基坑既有建筑物地基承载力。所得结果与斜坡上地基承载力的试验值和理论解对比分析后显示出其具有一定的合理性。所建议的抗剪强度指标取值方法,还有待于试验的进一步验证。     

开挖卸荷对既有群桩竖向承载性状的影响分析

既有建筑下增层开挖对已有桩基础的影响不同于基坑开挖对坑内桩基的影响。基于工程实例验证的有限元参数,用硬化土弹塑性模型模拟土体,用接触面单元模拟桩土相互作用,建立了桩筏基础-地基-增层开挖三维有限元模型,对增层开挖后群桩基础的竖向承载性状进行研究。分析了桩顶刚度、桩身轴力、桩侧摩阻力、桩端阻比以及土体回弹的变化规律,并研究了不同桩端土体刚度和增层开挖深度对这些参数的影响。结果表明,增层开挖后群桩中不同基桩表现出不同的承载性状,增大桩端土体刚度可明显提高单桩承载力和端阻比临界值;随着增层开挖深度的增加,侧阻和端阻的发挥程度也随之提高。研究结果有望为地下室增层开挖施工中的结构托换变形控制和补桩设计提供依据。     

Drained bearing response of shallow foundations on structured soils

This paper examines the drained bearing response of circular footings resting on structured soil deposits. Numerical simulations have been carried out using a finite element formulation of the Structured Cam Clay model. A parametric study was conducted by varying the parameters that govern the behaviour of structured soils and guidelines are given for designers to identify when effects of the soil structure are important. Under fully drained conditions, deformation within the structured soil supporting the footing usually occurs as a local or punching shear failure due to high compressibility of the structured soil and the mobilised bearing pressure increases with the footing movement, without reaching an ultimate value. A novel approximate method is presented to obtain the load–displacement response of a rigid circular footing resting on the surface of a structured soil deposit. This requires the properties of the soil in the reconstituted state and two additional parameters, which govern the natural structure of the soil. The proposed method has been applied to a published case study, where plate load test results are given for rigid circular steel plates resting on structured soil deposits. Fair agreement is observed between the computed and experimental results, suggesting the approximate method may be useful in design studies of foundations on structured soil deposits.     

Monitored and simulated variations in matric suction during rainfall in a residual soil slope

This research combines field, laboratory and numerical investigations to estimate the development of a wetting front within a 1.2 m residual soil mantle on a steep forested slope during rainfall events. The field-monitored variations in matric suction due to rain-water infiltration during various events revealed that the maximum infiltration rate was much higher when the wetting front resided in the upper 20 cm of soil compared to the case when the wetting front advanced to depths > 20 cm. Laboratory investigations on soil hydraulic properties (i.e., soil water characteristic curve, and hydraulic conductivity) were useful to establish the parameters of a multilayer finite-element model for one-dimensional vertical infiltration. These parameters were subsequently calibrated by matching the predicted and field measured transient pore water pressure responses during actual rainstorms with irregular rainfall patterns. The calibrated simulation model was used to assess the migration of the wetting front under uniform rainfall with different intensities. Based on the numerical results, a hyperbolic equation was developed to predict the duration of uniform rainfall required for the propagation of wetting front to a certain depth for a given rainfall intensity. The proposed equation was subsequently tested against field-monitored advancements of the wetting front during real rainstorms with variable rainfall intensity.     

济南泉域岩溶地下水三维等参有限元数值模拟

基于济南泉群研究区域岩溶水贮存和运动特征,建立了三维有限元数学模型,并利用等参有限元方法求解,再根据2001年、2002年的水文资料,识别模型参数,进一步对济南泉域内岩溶地下水的开采方案进行评价,探讨了济南市“保泉”与“供水”的可行性,为科学管理济南市泉群区域内地下水资源提供依据。