Theoretical Analysis of Bearing Capacity of Shallowly Embeded Rectangular Footing of Marine Structures

In this paper, the finite element analysis software ABAQUS is used to analyze the ultimate bearing capacity of three-dimensional rectangular footing of marine structures. The deformation law and the failure mode of homogeneous seabed soil beneath the rectangular footing are analyzed in detail. According to the equivalent plastic strain of soil under rectangular footing, an allowable velocity field of homogeneous seabed soil is reasonably constructed. Based on the plastic limit analysis theory of soil mass and by using the Mohr-Coulomb yield criterion, an upper bound solution of the ultimate bearing capacity of three-dimensional rectangular footing on general homogeneous seabed soil is derived, and a correction factor of ultimate bearing capacity of three-dimensional rectangular footing is given. To verify the rationality and applicability of this theoretical solution, some numerical solutions are achieved using the general-purpose FEM analysis package ABAQUS, and comparisons are made among the derived upper bound solution, the solution of Vesic, and the solution of Salgado et al. The results indicate that the upper bound solution of the three-dimensional shallowly embedded rectangular footing proposed in this paper is accurate in calculating the bearing capacity of homogeneous seabed soil. For undrained saturated clay foundation and sandy foundation with smaller internal friction angle, this upper bound solution can evaluate the ultimate bearing capacity of rectangular footing; with the gradual increase of the internal friction angle of the soil, the ultimate bearing capacity of the proposed upper bound solution is slightly higher than that of the rectangular footing.     

基于临界状态的边坡渐进破坏力学模型分析及应用

边坡条分法物理意义是以强度折减表征材料峰值强度的下降程度。以不平衡推力部分强度折减法为例,推广传统假设,应用理想弹塑性模型和全过程剪应力-应变模型,在峰值剪应力和峰值应变相等条件下,模拟边坡的渐进破坏过程;提出了临界状态决定法,随着临界状态逐步移动,边坡不平衡力和驱动力越来越大,抗滑力越来越小,直至整体处于破坏状态。以湖北省恩施州稻池村边坡为例,以3种数值(部分强度折减、理想弹塑性和全过程剪应力-应变本构)模型揭示边坡渐进破坏过程中各参量变化特征。结果表明:3种数值条块法均能描述稻池村边坡渐进破坏过程,基于条分法和理想弹塑性模型部分强度折减系数大于全过程剪应力-应变本构模型稳定系数。      

地震前兆特征与岩样剪切应变率异常数值模拟

利用拉格朗日元法(FLAC)模拟了岩样的宏观力学行为、剪切应变率异常及演化。在峰值强度岩石的本构模型分别取线弹性模型及应变软化模型的条件下,得到如下结论：1)在弹性阶段之后,应力-时步曲线具有一定的周期性。与地震复发的周期性类似。在加载过程中,形变场由均匀逐渐向不均匀转变．最终形成交叉的剪切带网络。被该网络所包围的试样中部可视为背景空区．在空区之外,地震群比较活跃。在应变硬化阶段。可以从未来震源观测到前兆一剪切应变率异常条带。2)在中期阶段之后,空区内地震活动的开始增强明显要比空区外地震活动增强晚。地震的迁移是由外向内。3)短期阶段则相反：在硬化阶段．剪切带内外的剪切应变率的差别不很显著。试样处于相对平静状态。在软化阶段,若采用剪切应变率最大的位置作为震源．可确定震源在试样的中部。源兆和场兆都集中在老断层上。     

Evaluation of the possible slip surface of a highly heterogeneous rock slope using dynamic reduction method

A new method, the dynamic reduction method (DRM) combined with the strain-softening method, was applied to evaluate the possible slip surface of a highly heterogeneous rock slope of the Dagangshan hydropower station in Southwest China. In DRM, only the strength of the failure elements is reduced and the softening reduction factor K is adopted to calculate the strength parameters. The simulation results calculated by DRM show that the further slip surface on the right slope of the Dagangshan hydropower station is limited in the middle part of the slope, while both SRM (strength reduction method) and LEM (limit equilibrium method) predict a failure surface which extends upper and longer. The observations and analysis from the three recorded sliding events indicate that the failure mode predicted by DRM is more likely the scenario. The results in this study illustrate that for highly heterogeneous slopes with geological discontinuities in different length scales, the proposed DRM can provide a reliable prediction of the location of the slip surface.     

基于一阶线性应变软化理论的边坡稳定性研究

在持续降雨或开挖卸荷作用下, 土体的强度指标会发生劣化, 但现今采用的边坡稳定性计算大多直接将其视为一个常数。为接近真实的边坡失稳破坏模式, 基于瑞典条分法以及一阶线性应变软化机制, 提出了一种边坡渐进破坏分析新方法, 推导出应变软化型边坡极限平衡表达式, 并获得了各破坏进度下边坡的安全系数。通过对模拟算例进行分析, 证明了条分-软化法的可靠性, 且计算结果表明渐进破坏过程中安全系数不仅取决于边坡的破坏方式与强度参数, 还与岩土体的软化模量密切相关。同时, 通过与滑坡实际案例的对比验算, 证实其强度指标存在不同的衰减系数, 即黏聚力的衰减系数大于摩擦角。从理论到应用, 最终获得的条分-软化法, 不仅考虑了岩土体的强度劣化效应以及滑动面的渐进发展, 还能有效地服务于实际工程背景下边坡的稳定性分析, 可以为滑坡的预防与治理提供指导建议。      

成层土中轴横受荷桩水平响应的非线性解

对于承受轴向荷载的水平受荷桩，以往研究大多基于线弹性或弹塑性水平荷载传递模型。为提升轴横受荷桩的计算设计水平，采用轴向荷载传递法计算桩身轴力，考虑桩身轴力引起的P-Δ效应，基于双曲线型水平荷载传递模型考虑桩-土体系变形的非线性特征，对成层土中轴横受荷桩的水平响应进行分析求解，得到了轴横荷载作用下桩身变形和内力的非线性有限差分解，并采用MATLAB语言编制了计算程序。使用模型试验算例与基于现场试验的有限元算例对非线性解的准确性进行对比验证，结果表明:计算结果与算例数据吻合良好，可靠性较高；采用不同荷载传递模型的计算结果在不同荷载水平下有所差异，在较大荷载水平下桩-土变形的非线性特点不容忽视。      

融雪条件下新疆伊犁谷地黄土-卵砾石接触面残余强度

为研究融雪作用下黄土-卵砾石接触面残余强度,以新疆伊犁谷地某黄土-卵砾石接触面滑坡为例,通过自制模具,制作黄土-卵砾石接触面环剪试样,开展不同含水率黄土-卵砾石接触面环剪试验,并通过扫描电镜（SEM）,从微观结构角度探究水对黄土-卵砾石接触面残余强度的影响。试验结果表明:黄土-卵砾石接触面抗剪强度随着含水率的增加而减小,剪切过程中存在应变软化现象,随着正应力及含水量的增加,应变软化现象越不明显;黄土-卵砾石接触面残余强度随着法向应力的增大而增大,且存在较好的线性关系,符合摩尔库伦强度准则;黄土-卵砾石接触面残余强度参数随着含水率的增加而降低,以含水率w=18%为界（塑限含水率19.1%附近）,当含水率为10%~18%时,残余内摩擦角φ_r变化较小（Δφ_r=5.4°）,当含水率为18%~26%时,残余内摩擦角φ_r变化较大（Δφ_r=9.0°）;微观结构方面,随着含水率的增加,黄土-卵砾石接触面形成软化"泥膜",部分填充了卵砾石凹凸部分,剪切面较为光滑,在剪切过程中,黏粒更加分散,附着在黄土颗粒表面,部分填充孔隙,起润滑作用,降低了残余强度。本文的研究成果可以为新疆伊犁谷地黄土-卵砾石接触面型滑坡形成机理研究及工程防治提供科学参考。      

Cohesive strength and seismogenic stress pattern along the active basement faults of the Precordillera-Sierras Pampeanas ranges, western Argentina: An experimental analysis by means of numerical model

A two-dimensional finite element method （FEM） model that incorporates faults, elastic rock physical properties, topographical load due to gravity and far-field plate velocity boundary conditions was used to recognize the seismogenic stress state along the fold-and-thrust belt of the Precordillera-Sierras Pampeanas ranges of western Argentina. A plane strain model with nine experiments was presented here to examine the fault strength with two major rock phyical properties： cohesion and angle of internal friction. Mohr-Coulomb failure criterion with bulk rock properties were applied to analyse faults. The stress field at any point of the model was assumed to be comprised of gravitational and tectonic components. The analysis was focused to recognize the seismogenic shear strain concentrated in the internal-cristaline domain of the orogene shown by the modeling. Modeling results are presented in terms of four parameters, i. e., （i） distributions, orientations, and magnitudes of principal stresses （σ1 and σ3）, （ii） displacement vector1 （iii） strain distribution, and （iv） maximum shear stress （τmax） contour line within the model. The simulation results show that the compressive stress is distributed in and around the fault systems. The overall orientation of of σ1 is in horizontal directions, although reorientations do occur within some stress weaker parts, especially subsequent to the faults. A large-scale shear stress is accumulating along the active faults of Tapias-Villicum Fault （TVF）, Salinas-Berros Fault （SBF）, Ampacama-Niquizanga Fault （ANF） and Las Charas Fault （CF）, which could act as local stress and strain modulators to localize the earthquakesoccurrence.     

Modelling the flexural buckling failure of stratified rock slopes based on the multilayer beam model

The buckling failure of stratified rock slopes intersected by a set of steep discontinuities that are approximately parallel to the slope surface is frequently encountered while constructing railways and roadways in mountainous areas. In this study, an analytical approach based on the energy equilibrium principle is presented to solve the flexural buckling stability of stratified rock slopes within the framework of multilayer beam theory. The generalized HoekBrown failure criterion is introduced to reflect the influences of slope size(scale effects) on the buckling stability. Subsequently, numerical and physical modellings from previous literatures are employed to validate the proposed approach. Furthermore, a practical case of Bawang Mountain landslide is also used for the comparative analysis. The study shows that the present analytical approach is capable to provide a more reasonable assessment for the buckling failure of stratified rock slopes, compared with several existing analytical approaches. Finally, a detailed parametric study is implemented, and the results indicate that the effects of rock strength, rock deformation modulus, geological strength index, layer thickness and disturbance degree of rock mass on the buckling failure of stratified rock slopes are more significant than that of rock type and slope angle.     

Hydro-mechanical analysis of rainfall-induced fines migration process within unsaturated soils

Seepage-induced fines migration under rainfall infiltration is a main cause leading to shallow failures in loose colluvial slopes. To describe the full process of fines migration within unsaturated soils during rainfall infiltration and the associated hydro-mechanical behaviors, a seepage-erosion-deformation coupled formulation is proposed in this paper. The governing equations proposed are implemented into a finite element code and used to investigate the influences of skeleton deformation on the rainfall infiltration process through unsaturated soil columns. The numerical results were presented in detail for a better understanding of the rainfall-induced fines migration process within unsaturated soils. Further, the obtained results are integrated into an infinite slope model for slope stability analysis. The results show that, the skeleton deformation will affect the rainfall infiltration rate and hence the timing of slope failures; meanwhile their influences are more evident if the fines deposition process is taken into account. Moreover, the slope stability could be reduced gradually due to the soil strength loss along with loss of fine particles. Therefore, particular attentions should be paid to analyzing the stability of soil slopes susceptible to internal erosion.     

Identifying crack initiation stress threshold in brittle rocks using axial strain stiffness characteristics

As an estimate for the in-situ spalling strength around massive underground excavations to moderately jointed brittle rocks, crack initiation stress marks the initiation of rock micro fracturing. It is crucial to accurately identify crack initiation stress level by proper method. In this study, confined compression tests of sandstone samples are used to examine the validity/applicability of proposed axial strain stiffness method. The results show that by highlighting the minuscule changes in stress-strain curve, the axial strain stiffness curve provided further insight into rock failure process and revealed five stages:(a) irregular fluctuation,(b) nearly horizontal regular fluctuation,(c) irregular fluctuation gradually decreasing to zero,(d) extreme fluctuation, and(e) near zero, which mainly correspond to five stages of stress–strain curve. The ratio of crack-initiation stress to peak strength determined using this approach is 0.44–0.51, similar to the ranges previously reported by other researchers. In this method, the key is to accurately detect the end point of the stage(b), "nearly horizontal regular fluctuation" characterized by a sudden change in axial strain stiffness curve, and the sudden change signifies crack initiation in rock sample. Finally, the research indicates that the axial strain stiffness curve can provide a mean to identify the crack-initiation stress thresholds in brittle rocks.     

隧道围岩变形监测反演分析及稳定性研究

针对有限元方法模拟隧道开挖过程受参数取值影响导致精度降低问题,提出利用BP神经网络和数值模拟方法来对隧道围岩力学参数进行反演分析。对比实测值与模拟值,证明了反演参数的可靠性,并得到隧道的最大最小主应力和位移变化值;利用强度折减法,通过数值模拟计算相应的拱顶沉降值,通过离散点拟合来求得隧道围岩的自稳系数,分析隧道围岩的稳定性。     

Discrete numerical modeling of granular materials considering crushability

The aim of this study is to numerically investigate the influence of particle breakage on the mechanical behavior of granular materials using a discrete element method (DEM). To enable particle crushing, non-crushable elementary particles are boned together to represents the granular aggregates which can be crushed when the external force exceeds its strength. The flaw of the aggregate was also modeled by randomly distributed void. Single particle crushing tests were carried out to determine the distribution of particle strength. The results of single particle crushing tests illustrate that the simulated single particle fracture strength and pattern agree well with the Weibull’s distribution equation. Conventional oedometer tests, drained monotonic and cyclic triaxial tests were also carried out to investigate the crushing of the aggregates and the associated mechanical behaviors. The effect of confining pressure on the crushing of aggregates and the mechanical behavior was also analyzed. It was found that the peak stress and dilation decrease significantly when particle crushing was considered. The deformation behavior of the specimen is essentially controlled by two factors: particle rearrangement - induced dilation and particle crushing - induced contraction. The increase of permanent strain and the reduction of dilation were observed during cyclic loading and they tend to reach a stable state after a certain number of cycles. The crushing of aggregate is most significant in the first two cycles. The results also indicated that for the same axial strain the volumetric strain and the bound breakage in the cyclic loading tests are significantly larger than that in the monotonic loading tests, especially at high cyclic stress ratio.     

基于有限元修正节理岩质边坡稳定性计算的解析解

采用有限元方法探讨在人工开挖或自然侵蚀环境下,岩质边坡体内应力场的变化及节理发育形成机理,并采用有限元强度折减法对后缘具有张节理边坡的稳定性影响因素进行敏感性对比分析,进而得出具有非贯通节理边坡稳定性计算的修正解析解.结果表明:卸荷及风化作用导致边坡体由表及里出现应力重分布及应力集中的现象,使边坡后缘由顶部向下发育一簇...     

河西重力变化的小波分解与地震活动关系的研究

讨论了小波分析在地震重力测量数据分析中应用的可能性，并采用基于样条函数的小波分析及其相应的Β小波分别计算了河西地区的重力变化资料，其结果如下：①小波分解可有效分离重力场时间变化的不同空间波长成分，更清晰地看清重力场变化与地震活动的关系。②河西地区重力资料分解结果表明：反映该地区浅部物质变化的重力变化高频部分和反映该地区深部物质变化的重力变化低频部分对永登Ｍｓ５．８级地震的反应不明显，而反映该地区中层物质变化的重力变化二、三阶小波分解结果对该地震的反应较明显。重力小波分解计算结果反映的密度变化深度与地震资料统计的震源深度基本一致。     

Effect of clay content to the strength of gravel soil in the source region of debris flow

The production of runoff in the source area of a debris flow is the consequence of a reduction in soil strength. Gravel soil is widely distributed in the source region, and the influence of its clay content on soil strength is one of the important questions regarding the formation mechanism of debris flows. In this paper, the clay content in gravel soil is divided into groups of low clay content(1%, 2, 5%), moderate clay content(3.75%, 5.00%, 6.25%, 7.5%) and high clay content(10.0%, 12.5%, 15%). Tests of the unconsolidated undrained shear strength and consolidated drained shear strength were performed. The unconsolidated undrained shearing(UU) experiment simulates the rapid shear failure of loose gravel soil under the conditions of brief heavy rainfall. The consolidated drained shearing(CD) experiment simulates creep failure of consolidated sediment during extended rainfall. The pore water pressure first increased and then decreased as the clay content increased, and the increase in pore pressure was relatively high in the gravel soil sample when the clay content is in the range of 3.25-7.50%, and stress in the gravel soil is relatively low for a moderate clay content. Gravelly soils with a moderate clay content are moreprone to debris-flow initiation. This paper presents a mathematical formula for the maximum shear stress and clay content of gravel soil under two conditions. The key processes whereby the soil fails and triggers a debris flow—volume contraction of soil, expansion of clay soil, and rise of pore pressure―cause reductions in the soil friction force and enhancement of the water content in the clay particles, and subsurface erosion of soil reduces the soil viscosity, which eventually reduces the soil strength so that the soil loses its stability, liquefies and generates a debris flow.     

Model refinement for offshore platforms: Experimental study

Offshore jacket platforms are widely used in offshore oil and gas exploitation. Finite element models of such structures need to have many degrees of freedom (DOFs) to represent the geometrical detail of complex structures, thereby leading to incompatibility in the number of DOFs of experimental models. To bring them both to the same order while ensuring that the essential eigen- properties of the refined model match those of experimental models, an extended model refinement procedure is presented in this paper. Vibration testing of an offshore jacket platform model is performed to validate the applicability of the proposed approach. A full-order finite element model of the platform is established and then tuned to meet the measured modal properties identified from the acceleration signals. Both model reduction and modal expansion methods are investigated, as well as various scenarios of sensor arrangements. Upon completion of the refinement, the updated jacket platform model matches the natural frequencies of the measured model well.     

对流边界层理论中不稳定非自治系统的数值解（英文）

为了求解对流边界层理论中一个非自治微分方程系统,作者采用伽略金有限元方法,此方法是通过将无限区间上的三阶非线性微分方程转化成有限区间上的二阶微分形式,并构造出相应的伽略金有限元方程来求得数值解,该数值解与先前一些作者的结果一致,并且计算效率显高于其它数值方法.     

Adaptive finite element modeling of marine controlled-source electromagnetic fields in two-dimensional general anisotropic media

In this paper, we extend the scope of numerical simulations of marine controlled-source electromagnetic (CSEM) fields in a particular case of anisotropy (dipping anisotropy) to the general case of anisotropy by using an adaptive finite element approach. In comparison to a dipping anisotropy case, the first order spatial derivatives of the strike-parallel components arise in the partial differential equations for generally anisotropic media, which cause a non-symmetric linear system of equations for finite element modeling. The adaptive finite element method is employed to obtain numerical solutions on a sequence of refined unstructured triangular meshes, which allows for arbitrary model geometries including bathymetry and dipping layers. Numerical results of a 2D anisotropic model show both anisotropy strike and dipping angles have great influence on the marine CSEM responses.     

Elastic-plastic analysis for pile-anchor supporting system of deep foundation pit

By using numerical analysis methods to simulate the deep excavation, a lot of analyses are established on the basis of two-dimensional plane strain, ignoring the fact that foundation pit possesses three dimensions. For soil constitutive relation, people always take linear and nonlinear model, without considering the plastic behavior of soil. Using plastic-elastic hardening model to simulate constitutive relation of soil characteristics, the authors carried out mechanical analysis for pit excavation and support. The results show that the analysis for the stress state of pile anchor system is an effective way which provides theoretical basis for calculation of soil displacement.