A continuum method of pile driving analysis: Comparison with the wave equation method

What may be called a ‘Continuum’ method of finite element analysis is used in this paper to predict the behaviour of a pile during driving. In this both the pile and the soil are treated as two distinct parts of the same solid continuum, but with different properties. The behaviour of the soil medium, assumed to be semi-infinite and nonlinear, is represented by the hyperbolic stress-strain relationship. Discretizing the pile-soil system in turn by conventional axi-symmetric and mapped finite elements, the problem is solved in the time-domain using the central difference scheme.

The example considered is that of a fully embedded steel pipe pile for which both field test data and Wave Equation solution are available. Results show that: the Continuum Method is capable of a greater degree of accuracy than the conventional Wave Equation Method, but it is far more expensive than the latter in terms of computational effort needed; the effects of radiation damping and wave dispersion in the soil are found to be small; and the mapped finite elements give significantly better results than conventional elements.     

A comparative study for the estimation of geodetic point velocity by artificial neural networks

Space geodesy era provides velocity information which results in the positioning of geodetic points by considering the time evolution. The geodetic point positions on the Earth’s surface change over time due to plate tectonics, and these changes have to be accounted for geodetic purposes. The velocity field of geodetic network is determined from GPS sessions. Velocities of the new structured geodetic points within the geodetic network are estimated from this velocity field by the interpolation methods. In this study, the utility of Artificial Neural Networks (ANN) widely applied in diverse fields of science is investigated in order to estimate the geodetic point velocities. Back Propagation Artificial Neural Network (BPANN) and Radial Basis Function Neural Network (RBFNN) are used to estimate the geodetic point velocities. In order to evaluate the performance of ANNs, the velocities are also interpolated by Kriging (KRIG) method. The results are compared in terms of the root mean square error (RMSE) over five different geodetic networks. It was concluded that the estimation of geodetic point velocity by BPANN is more effective and accurate than by KRIG when the points to be estimated are more than the points known.     

Modeling the particle breakage of rockfill materials with the cohesive crack model

A practical combined finite–discrete element method was developed to simulate the breakage of irregularly shaped particles in granular geomaterials, e.g., rockfill. Using this method, each particle is discretized into a finite element mesh. The potential fracture paths are represented by pre-inserted cohesive interface elements (CIEs) with a progressive damage model. The Mohr–Coulomb model with a tension cut-off is employed as the damage initiation criterion to rupture the predominant failure mode occurs at the particle scale. Two series of biaxial tests were simulated for both the breakable and unbreakable particle assemblies. The two assemblies have identical configurations, with the exception that the former is inserted with CIEs and is breakable. The simulated stress–strain–dilation responses obtained for both assemblies are in agreement with experimental observations. We present a comprehensive study of the role of particle breakage on the mechanical behavior of rockfill materials at both the macroscopic and microscopic scales. The underlying mechanism of particle breakage can be explained by the force chain in the assemblies.     

Elastically-homogeneous lattice models of damage in geomaterials

This study involves the development of the auxiliary stress approach for producing elastically-homogeneous lattice models of damage in geomaterials. The lattice models are based on random, three-dimensional assemblages of rigid-body-spring elements. Unlike conventional lattice or particle models, the elastic constants of a material (e.g., Young’s modulus and Poisson’s ratio) are represented properly in both global and local senses, without any need for calibration. The proposed approach is demonstrated and validated through analyses of homogeneous and heterogeneous systems under uni- and tri-axial loading conditions. Comparisons are made with analytical solutions and finite element results. Thereafter, the model is used to simulate a series of standard laboratory tests: (a) split-cylinder tests, and (b) uniaxial compressive tests of sedimentary rocks at the Horonobe Underground Research Laboratory in Hokkaido, Japan. Model inputs are based on physical quantities measured in the experiments. The simulation results agree well with the experimental results in terms of pre-peak stress-strain/displacement responses, strength measurements, and failure patterns.     

Performance evaluation of artificial neural networks for planimetric coordinate transformation—a case study,Ghana

Two national horizontal geodetic datums, namely, the Accra and Leigon datum, have been the only available datum used in Ghana. These two datums are non-geocentric and were established based on astro-geodetic observations. Relating these different geodetic datums mostly involves the use of conformal transformation techniques which could produce results that are not very often satisfactory for certain geodetic, surveying and mapping purposes. This has been ascribed to the incapability of the conformal models to absorb more of the heterogeneous and local character of deformations existing within the local geodetic networks. Presently, application of new approaches such as artificial neural network (ANN) is highly recommended. Whereas the ANN has been gaining much popularity to solving coordinate transformation-related problems in recent times, the existing researches carried out in Ghana have shown that only three-dimensional conformal transformation methods have been utilized. To the best of our knowledge, plane coordinate transformation between the two local geodetic datums in Ghana has not been investigated. In this paper, an attempt has been made to explore the plane coordinate transformation performance of two different ANN approaches (backpropagation neural network (BPNN) and radial basis function neural network (RBFNN)) compared with two different traditional techniques (six- and four-parameter models) in the Ghana national geodetic reference network. The results revealed that transforming plane coordinates from Leigon to Accra datum, the RBFNN was better than the BPNN and traditional techniques. Transforming from Accra to Leigon datum, both the BPNN and RBFNN produced closely related results and were better than the traditional methods. Therefore, this study will create the opportunity for Ghana to recognize the significance and strength of the ANN technology in solving coordinate transformation problems.     

GPS高程拟合方法精度分析

近年来人们往往利用GPS数据来确定大地高,但大地高不同于正常高,为此,利用多项式拟合与地球重力场模型相结合的数学方法,使GPS所测大地高通过这些数学模型直接转换为具有厘米精度的正常高,将该方法得到的正常高与单独利用多项式拟合和地球重力场模型得到的结果进行了比较,其差值的标准差为±38 cm。     

Effect of inherent microcrack populations on rock tensile fracture behaviour: numerical study based on embedded discontinuity finite elements

Inherent microcrack populations have a significant effect on the fracture behaviour of natural rocks. The present study addresses this topic in numerical simulations of uniaxial tension and three-point bending tests. For this end, a rock fracture model based on multiple intersecting embedded discontinuity finite elements is developed. The inherent (pre-existing) microcrack populations are represented by pre-embedded randomly oriented discontinuity populations. Crack shielding (through spurious locking) is prevented by allowing a new crack to be introduced, upon violation of the Rankine criterion, in an element with an initial crack unfavourably oriented to the loading direction. Rock heterogeneity is accounted for by random clusters of triangular finite elements representing different minerals of granitic numerical rock. Numerical simulations demonstrate the strength lowering effect of initial microcrack populations. This effect is substantially stronger under uniaxial tension, due to the uniform stress state, than in semicircular three-point bending having a non-uniform stress state with a clear local maximum of tensile stress.

     

Finite element algorithm for jointed concrete pavements subjected to moving aircraft

An algorithm based on the finite element method is developed to analyze the dynamic response of multiple, jointed concrete pavements to moving aircraft loads. In the finite element idealization, the pavement-subgrade system is idealized by thin plate finite elements resting on Winkler-type viscoelastic foundation represented by a series of distributed springs and dashpots. The dowel bars at the transverse joints are represented by beam elements. It i assumed that the dowel bar is fully embedded into the pavement thus neglecting dowel-pavement interaction effects. The longitudinal keyed or aggregate interlock joints are modeled by vertical spring elements. The dynamic aircraft-pavement interaction effects are considered in the analysis by modeling the aircraft by masses supported by spring-dashpot systems representing the landing gear of the aircraft. It is assumed that the aircraft travels along a straight line with a specific initial velocity and acceleration. The aircraft-pavement interaction takes the form of two sets of coupled equations which result ina non-symmetric stiffness matrix. An approximate mixed iteration-direct elimination scheme is used to solve for the dynamic equations. The accuracy of the computer code is verified by the available experimental and analytical solutions. A parametric study is conducted to investigate the effects of various parameters on the dynamic response of pavements.     

Recent crustal deformation and seismicity studies at Abu-Dabbab area, Red Sea, Egypt

The Abu-Dabbab area is characterized by high seismicity and complex tectonic setting; for these facts, a local geodetic network consisting of 11 geodetic benchmarks has been established. The crustal deformation data in this area are collected using the GPS techniques. Five campaigns of GPS measurements have been collected, processed, and adjusted to get the more accurate positions of the GPS stations. The horizontal velocity vectors, the dilatational, the maximum shear strains, and the principal strain rates were estimated. The horizontal velocity varies in average between 3 and 6 mm per year across the network. The results of the deformation analysis indicate a significant contraction and extension across the southern central part of the study area which is characterized by high seismic activity represented by the clustering shape of the microearthquakes that trending ENE. The north and northeastern parts are characterized by small strain rates. This study is an attempt to provide valuable information about the present state of the crustal deformation and its relationship to seismic activity and tectonic setting at the Abu-Dabbab area. The present study is the first work demonstrating crustal deformation monitoring at the Abu-Dabbab area. The time interval is relatively short. Actually, these results are preliminary results. So, the continuity of GPS measurements is needed for providing more information about the recent crustal deformation in that area.     

Investigation of a hybrid polygonal finite element formulation for confined and unconfined seepage

This study presents a formulation for field problems using hybrid polygonal finite elements, taking steady state seepage through a porous material as the focus. We make comparisons with a conventional finite element formulation based on a single primary variable, focussing on the advantages of the hybrid formulation in terms of flux field accuracy and extension to convex polygonal shaped elements. For the unconfined case, we adopt a head dependent hydraulic conductivity that does not require remeshing. The performance of the hybrid polygonal element formulation is demonstrated through a series of numerical examples. The results show a sensitivity of the location of the free surface in unconfined seepage to mesh configuration for hybrid quadrilateral meshes with various aspect ratios, but not for hybrid polygonal meshes with various orientations and irregularity. Examination of the free surface location results for several conforming shape function options shows an insensitivity to choice of interpolation function, provided that it conforms with the assumptions in the formulation. Copyright © 2016 John Wiley & Sons, Ltd.     

Beam–solid contact formulation for finite element analysis of pile–soil interaction with arbitrary discretization

This paper presents an embedded beam formulation for discretization independent finite element (FE) analyses of interactions between pile foundations or rock anchors and the surrounding soil in geotechnical and tunneling engineering. Piles are represented by means of finite beam elements embedded within FEs for the soil represented by 3D solid elements. The proposed formulation allows consideration of piles and pile groups with arbitrary orientation independently from the FE discretization of the surrounding soil. The interface behavior between piles and the surrounding soil is represented numerically by means of a contact formulation considering skin friction as well as pile tip resistance. The pile–soil interaction along the pile skin is considered by means of a 3D frictional point‐to‐point contact formulation using the integration points of the beam elements and reference points arbitrarily located within the solid elements as control points. The ability of the proposed embedded pile model to represent groups of piles objected to combined axial and shear loading and their interactions with the surrounding soil is demonstrated by selected benchmark examples. The pile model is applied to the numerical simulation of shield driven tunnel construction in the vicinity of an existing building resting upon pile foundation to demonstrate the performance of the proposed model in complex simulation environments. Copyright © 2014 John Wiley & Sons, Ltd.     

Effect of vehicle–pavement interaction on dynamic response of rigid pavements

In the present paper, dynamic response of rigid pavements subjected to moving vehicular loads is presented using an analytical procedure based on the finite element method. The vehicle–pavement interaction effects are taken into account while developing the solution algorithm. The concrete pavement is discretized by finite and infinite beam elements. Infinite elements are helpful in appropriate modeling of end conditions. The main purpose of the development of infinite elements is to model the unbounded domain. The underlying soil medium is modeled by Pasternak model which assumes the existence of shear interaction between the spring elements. The moving vehicle is represented by a mass supported by a spring-dashpot system. The vehicle -pavement interaction force is modeled with Dirac-delta function. Dynamic equilibrium equation is solved with Newmark-Beta integration scheme. It has been observed that the dynamic interaction between the moving load and the pavement has significant effect on pavement response. Parametric study is carried out to investigate the effect of vehicle–pavement interaction (VPI) and soil parameters on the response of pavement.     

Nonlinear analysis of NATM tunnel construction with the boundary element method

This paper presents a novel approach to the simulation of NATM tunnel construction using the Boundary Element Method (BEM) as principal numerical method. This new approach has the advantage that only the excavation surface, the possible plastic zones and the tunnel lining have to be discretised. The whole rock mass is represented by the BEM whereas the Finite Element Method (FEM) is used to represent the tunnel lining only. Thus, a general coupling strategy for coupling three-dimensional boundary elements with shell finite elements (shotcrete) and beam finite elements (steel arches) is presented. To achieve realistic results the effect of hydration of the shotcrete and yielding of the steel arches is considered in the excavation process. Furthermore, the nonlinear rock behaviour is modelled more realistically by using a powerful hierarchical constitutive model which considers a large range of rock materials. The combination of these ideas results in higher user-friendliness and efficiency. Some verification tests and practical applications in tunnelling are presented.     

冀东平原土壤中重金属元素的地球化学特征

通过分析冀东平原土壤中有害重金属元素Pb、Cr、Cd、As、Hg的含量变化特征,显示有害重金属元素Hg、Cd、Pb在冀东平原表层土壤的富集已较为明显,Cr、As也形成局部的异常,原有的地球化学特征发生较大变化。初步认为表层土壤中重金属元素的高含量与当地地质背景和人为活动有关,局部异常成因需要进一步研究。研究区地积累指数污染评价结果表明:以全国土壤背景值计算,重金属未造成污染;以研究区背景值计算,Hg表现为轻度—中度污染状况,As、Cd、Cr、Pb未造成污染。     

Rock slope stability assessment using finite element based modelling – examples from the Indian Himalayas

Numerical modelling of rock slides is a versatile approach to understand the failure mechanism and the dynamics of rock slopes. Finite element slope stability analysis of three rock slopes in Garhwal Himalaya, India has been carried out using a two dimensional plane strain approach. Two different modelling techniques have been attempted for this study. Firstly, the slope is represented as a continuum in which the effect of discontinuities is considered by reducing the properties and strength of intact rock to those of rock mass. The equivalent Mohr-Coulomb shear strength parameters of generalised Hoek-Brown (GHB) criterion and modified Mohr-Coulomb (MMC) criterion has been used for this continuum approach. Secondly, a combined continuum-interface numerical method has been attempted in which the discontinuities are represented as interface elements in between the rock walls. Two different joint shear strength models such as Barton-Bandis and Patton’s model are used for the interface elements. Shear strength reduction (SSR) analysis has been carried out using a finite element formulation provided in the PHASE². For blocky or very blocky rock mass structure combined continuum-interface model is found to be the most suitable one, as this model is capable of simulating the actual field scenario.     

Non-linear analysis of geomechanical problems using coupled finite and infinite elements

In modeling of many geomechanics problems such as underground openings, soil-foundation structure interaction problems, and in wave propagation problems through semi-infinite soil medium the soil is represented as a region of either infinite or semi-infinite extent. Numerical modeling of such problems using conventional finite elements involves a truncation of the far field in which the infinite boundary is terminated at a finite distance. In these problems, appropriate boundary conditions are introduced to approximate the solution of the infinite or semi-infinite boundaries as closely as possible. However, the task of positioning the finite boundary in conventional finite element discretization and the definition of the boundary and its conditions is very delicate and depends on the modeller's skill and intuition. Moreover, such a choice is influenced by the size of the domain to be discretized. Consequently, the dimensions of the global matrices and the time required for solution of the problem will increase considerably and also selection of the arbitrary location of truncated boundary may lead to erroneous result. In order to over come these problems, mapped infinite elements have been developed by earlier researchers (Simoni and Schrefier, 1987). In the present work the applicability of infinite element technique is examined for different geomechanics problems. A computer program INFEMEP is developed based on the conventional finite element and mapped infinite element technique. It is then validated using selected problems such as strip footing and circular footing. CPU time taken to obtain solutions using finite element approach and infinite element approach was estimated and presented to show the capability of coupled modeling in improving the computational efficiency. Mesh configurations of different sizes were used to explore the enhancement of both computational economy and solution accuracy achieved by incorporation of infinite elements to solve elastic and elasto-plastic problems in semi-infinite/finite domain as applied to geotechnical engineering. © Rapid Science Ltd. 1998     

岩土参数随机场离散的三角形单元局部平均法

将不确定性岩土参数建模为随机场而非传统意义上的随机变量,基于随机场的局部平均理论,提出了用于二维随机场离散的三角形单元局部平均法。通过面积坐标变换和高斯数值积分,给出了三角形单元局部平均随机场协方差矩阵的解析计算方法和数值计算方法。采用算例再现了所提方法的分析过程和有效性,并与传统二维随机场四边形单元离散法进行了对比。结果表明:提出的二维随机场三角形单元离散法能与有限元三角形单元离散法完美结合,随机场单元与有限元单元的对应关系清晰,易于随机有限元程序的编制;对于随机场单元的均值,传统四边形单元离散法与所提方法的计算结果相同;对于随机场单元的方差,传统四边形单元离散法计算结果偏小,所提方法显得更加科学、合理。     

Dynamic response of a thick plate on viscoelastic foundation to moving loads

A finite element algorithm is presented to evaluate the dynamic response of pavements including aircraft-guideway-foundation interaction. The pavement-foundation system is modeled by a series of thick plate elements supported by discrete springs and dashpots at the nodal points representing the viscoelastic foundation. The moving aircraft loads are represented by masses each supported by a spring and dashpot suspension system and having a specified horizontal velocity and acceleration. The accuracy of the algorithm is verified by comparing the finite element solution with available analytical results. A parametric study is conducted to determine the effects of various parameters on the dynamic response of pavements to moving loads.     

Qualitative and quantitative analysis of vorticity,strain and area change in general non-isochoric 2D deformation

A scale free representation of a general non-isochoric 2D deformation is presented which is amenable to mathematical analysis. By describing deformation in 2D in terms of polar coordinates the stretching and rotational histories of linear elements separate and are easily analysed both qualitatively and quantitatively. An analysis of finite strain combined with dynamical considerations allows the derivation of equations which may be used to estimate finite strain, area change and kinematic vorticity number. Numerical investigation of method developed here was carried out and it was found to perform well unless large area changes occur in combination with large components of simple shear. A re-analysis of natural data indicates the method is consistent.     

Resilience at the periphery: Insurgency,agency and social-ecological change under armed conflict

Armed conflict has played an increasingly important role in the transformation of key social and environmental systems at multiple spatial and temporal scales. Accelerated resource flows and environmental change dynamics intersect with conflict processes in ways that are substantial and yet inadequately understood. Drawing on research along the Pakistani border in eastern Afghanistan’s embattled province of Nangarhar, we employ a coupled systems approach for understanding the ways in which social-ecological processes shape and are shaped by armed conflict. Based on field surveys, geospatial analysis of land and forest change, and participatory research among local communities, government agencies and military actors, we identify several causal processes linking conflict and dynamics of social-ecological change in the context of multiscalar geopolitical processes. We focus attention on four inter-related elements: (1) transitional modes of resource governance relating to armed militia groups and state intervention, (2) forest changes related to illegal logging and trade networks, (3) the erosion of upper-montane rangelands through encroachment and changing pastoral responses to conflict, and (4) significant land use changes in the agricultural sector toward the cultivation of opium poppy. Our research highlights the importance of center-periphery relations, the problematic nature of local agency, and the ways in which local social-ecological elements—here, particularly, timber and opium—become political objects within competing narratives of (in)security and ongoing state formation.