首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
A computational method, incorporating the finite element model (FEM) into data assimilation using the particle filter, is presented for identifying elasto‐plastic material properties based on sequential measurements under the known changing traction boundary conditions to overcome some difficulties in identifying the parameters for elasto‐plastic problems from which the existing inverse analysis strategies have suffered. A soil–water coupled problem, which uses the elasto‐plastic constitutive model, is dealt with as the geotechnical application. Measured data on the settlement and the pore pressure are obtained from a synthetic FEM computation as the forward problem under the known parameters to be identified for both the element tests and the ground behavior during the embankment construction sequence. Parameter identification for elasto‐plastic problems, such as soil behavior, should be made by considering the measurements of deformation and/or pore pressure step by step from the initial stage of construction and throughout the deformation history under the changing traction boundary conditions because of the embankment or the excavation because the ground behavior is highly dependent on the loading history. Thus, it appears that sequential data assimilation techniques, such as the particle filter, are the preferable tools that can provide estimates of the state variables, that is, deformation, pore pressure, and unknown parameters, for the constitutive model in geotechnical practice. The present paper discusses the priority of the particle filter in its application to initial/boundary value problems for elasto‐plastic materials and demonstrates a couple of numerical examples. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

2.
This study concerns the identification of parameters of soil constitutive models from geotechnical measurements by inverse analysis. To deal with the non‐uniqueness of the solution, the inverse analysis is based on a genetic algorithm (GA) optimization process. For a given uncertainty on the measurements, the GA identifies a set of solutions. A statistical method based on a principal component analysis (PCA) is, then, proposed to evaluate the representativeness of this set. It is shown that this representativeness is controlled by the GA population size for which an optimal value can be defined. The PCA also gives a first‐order approximation of the solution set of the inverse problem as an ellipsoid. These developments are first made on a synthetic excavation problem and on a pressuremeter test. Some experimental applications are, then, studied in a companion paper, to show the reliability of the method. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

3.
A constitutive model that captures the material behavior under a wide range of loading conditions is essential for simulating complex boundary value problems. In recent years, some attempts have been made to develop constitutive models for finite element analysis using self‐learning simulation (SelfSim). Self‐learning simulation is an inverse analysis technique that extracts material behavior from some boundary measurements (eg, load and displacement). In the heart of the self‐learning framework is a neural network which is used to train and develop a constitutive model that represents the material behavior. It is generally known that neural networks suffer from a number of drawbacks. This paper utilizes evolutionary polynomial regression (EPR) in the framework of SelfSim within an automation process which is coded in Matlab environment. EPR is a hybrid data mining technique that uses a combination of a genetic algorithm and the least square method to search for mathematical equations to represent the behavior of a system. Two strategies of material modeling have been considered in the SelfSim‐based finite element analysis. These include a total stress‐strain strategy applied to analysis of a truss structure using synthetic measurement data and an incremental stress‐strain strategy applied to simulation of triaxial tests using experimental data. The results show that effective and accurate constitutive models can be developed from the proposed EPR‐based self‐learning finite element method. The EPR‐based self‐learning FEM can provide accurate predictions to engineering problems. The main advantages of using EPR over neural network are highlighted.  相似文献   

4.
The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

5.
Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

6.
This paper presents novel visualization techniques to simplify representation of the fourth‐order material stiffness tensor as a set of three‐dimensional geometric objects. Stiffness visualization aids in understanding the complex stiffness characteristics of highly non‐linear constitutive models including modelled material anisotropy and loading path dependent stiffness variation. Stiffness visualization is relevant for understanding the relationship of material stiffness to global behaviour in the analysis of a boundary value problem. The spherical pulse stiffness visualization method, developed in the acoustics field, is extended to visualize stiffness of geomaterials using three three‐dimensional objects. This method is limited to relatively simple constitutive models with symmetric stiffness matrices insensitive to loading magnitude and direction. A strain dependent stiffness visualization method is developed that allows the examination of material stiffness for a range of loading directions and is suitable for highly non‐linear and path dependent material models. The proposed stiffness visualization can be represented as 3‐D, 2‐D and 1‐D objects. The visualization technique is used to represent material stiffness and its evolution during simulated soil laboratory tests and deep excavation construction. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

7.
Borehole closure data have often been used to determine the constitutive behaviour of rock salt. This is an inverse problem and analysis of the data is confused by the fact that borehole closure is a structural response that does not uniquely determine constitutive behaviour since both the stress and deformation depend on the constitutive behaviour. This application brief assesses the conclusion of a recent study that the transient power law, ? = K τmtn, which includes no term to account for steady-state deformation, is a better constitutive model for salt than are models that account for steady-state straining, and that an empirical model developed from the borehole closure data could be extrapolated to give conservative predictions of long-term closure of single or widely spaced openings.  相似文献   

8.
Effective capabilities of combined chemo‐elasto‐plastic and unsaturated soil models to simulate chemo‐hydro‐mechanical (CHM) behaviour of clays are examined in numerical simulations through selected boundary value problems. The objective is to investigate the feasibility of approaching such complex material behaviour numerically by combining two existing models. The chemo‐mechanical effects are described using the concept of chemical softening consisting of reduction of the pre‐consolidation pressure proposed originally by Hueckel (Can. Geotech. J. 1992; 29 :1071–1086; Int. J. Numer. Anal. Methods Geomech. 1997; 21 :43–72). An additional chemical softening mechanism is considered, consisting in a decrease of cohesion with an increase in contaminant concentration. The influence of partial saturation on the constitutive behaviour is modelled following Barcelona basic model (BBM) formulation (Géotech. 1990; 40 (3):405–430; Can. Geotech. J. 1992; 29 :1013–1032). The equilibrium equations combined with the CHM constitutive relations, and the governing equations for flow of fluids and contaminant transport, are solved numerically using finite element. The emphasis is laid on understanding the role that the individual chemical effects such as chemo‐elastic swelling, or chemo‐plastic consolidation, or finally, chemical loss of cohesion have in the overall response of the soil mass. The numerical problems analysed concern the chemical effects in response to wetting of a clay specimen with an organic liquid in rigid wall consolidometer, during biaxial loading up to failure, and in response to fresh water influx during tunnel excavation in swelling clay. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

9.
A challenging computational problem arises when a discrete structure (e.g. foundation) interacts with an unbounded medium (e.g. deep soil deposit), particularly if general loading conditions and non‐linear material behaviour is assumed. In this paper, a novel method for dealing with such a problem is formulated by combining conventional three‐dimensional finite‐elements with the recently developed scaled boundary finite‐element method. The scaled boundary finite‐element method is a semi‐analytical technique based on finite‐elements that obtains a symmetric stiffness matrix with respect to degrees of freedom on a discretized boundary. The method is particularly well suited to modelling unbounded domains as analytical solutions are found in a radial co‐ordinate direction, but, unlike the boundary‐element method, no complex fundamental solution is required. A technique for coupling the stiffness matrix of bounded three‐dimensional finite‐element domain with the stiffness matrix of the unbounded scaled boundary finite‐element domain, which uses a Fourier series to model the variation of displacement in the circumferential direction of the cylindrical co‐ordinate system, is described. The accuracy and computational efficiency of the new formulation is demonstrated through the linear elastic analysis of rigid circular and square footings. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

10.
This paper discusses the quality of the procedure employed in identifying soil parameters by inverse analysis. This procedure includes a FEM‐simulation for which two constitutive models—a linear elastic perfectly plastic Mohr–Coulomb model and a strain‐hardening elasto‐plastic model—are successively considered. Two kinds of optimization algorithms have been used: a deterministic simplex method and a stochastic genetic method. The soil data come from the results of two pressuremeter tests, complemented by triaxial and resonant column testing. First, the inverse analysis has been performed separately on each pressuremeter test. The genetic method presents the advantage of providing a collection of satisfactory solutions, among which a geotechnical engineer has to choose the optimal one based on his scientific background and/or additional analyses based on further experimental test results. This advantage is enhanced when all the constitutive parameters sensitive to the considered problem have to be identified without restrictions in the search space. Second, the experimental values of the two pressuremeter tests have been processed simultaneously, so that the inverse analysis becomes a multi‐objective optimization problem. The genetic method allows the user to choose the most suitable parameter set according to the Pareto frontier and to guarantee the coherence between the tests. The sets of optimized parameters obtained from inverse analyses are then used to calculate the response of a spread footing, which is part of a predictive benchmark. The numerical results with respect to both the constitutive models and the inverse analysis procedure are discussed. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

11.
The dynamic response of an end bearing pile embedded in a linear visco‐elastic soil layer with hysteretic type damping is theoretically investigated when the pile is subjected to a time‐harmonic vertical loading at the pile top. The soil is modeled as a three‐dimensional axisymmetric continuum in which both its radial and vertical displacements are taken into account. The pile is assumed to be vertical, elastic and of uniform circular cross section. By using two potential functions to decompose the displacements of the soil layer and utilizing the separation of variables technique, the dynamic equilibrium equation is uncoupled and solved. At the interface of soil‐pile system, the boundary conditions of displacement continuity and force equilibrium are invoked to derive a closed‐form solution of the vertical dynamic response of the pile in frequency domain. The corresponding inverted solutions in time domain for the velocity response of a pile subjected to a semi‐sine excitation force applied at the pile top are obtained by means of inverse Fourier transform and the convolution theorem. A comparison with two other simplified solutions has been performed to verify the more rigorous solutions presented in this paper. Using the developed solutions, a parametric study has also been conducted to investigate the influence of the major parameters of the soil‐pile system on the vertical vibration characteristics of the pile. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

12.
Observations from earthquakes over the past several decades have highlighted the importance of local site conditions on propagated ground motions. Downhole arrays are deployed to measure motions at the ground surface and within the soil profile, and also to record the pore pressure response within the soft soil profiles during excitation. The degradation of soil stiffness as excess pore pressures are generated during earthquake events has also been observed. An inverse analysis framework is developed and demonstrated to directly extract soil material behavior including pore water pressure (PWP) generation from downhole array measurements that can then be readily used in 1D nonlinear site response analysis. The self‐learning simulations (SelfSim) inverse analysis framework, previously developed for total stress site response analysis, is extended to extract PWP generation behavior of the soil in addition to cyclic response during ground shaking. A Neural Network based constitutive model is introduced to represent PWP generation during cyclic loading. A new analysis scheme is introduced that can use data from co‐located piezometer and accelerometer sensors. The successful performance of the proposed framework is demonstrated using four synthetic vertical array recordings. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

13.
Piles may be subjected to lateral soil pressures as a result of lateral soil movements from nearby construction‐related activities such as embankment construction or excavation operations. Three‐dimensional finite element analyses have been carried out to investigate the response of a single pile when subjected to lateral soil movements. The pile and the soil were modelled using 20‐node quadrilateral brick elements with reduced integration. For compatibility between the soil–pile interface elements, 27‐node quadrilateral brick elements with reduced integration were used to model the soil around the pile adjacent to the soil–pile interface. A Mohr–Coulomb elastic–plastic constitutive model with large‐strain mode was assumed for the soil. The analyses indicate that the behaviour of the pile was significantly influenced by the pile flexibility, the magnitude of soil movement, the pile head boundary conditions, the shape of the soil movement profile and the thickness of the moving soil mass. Reasonable agreement is found between some existing published solutions and those developed herein. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

14.
We pay a revisit to some classical geomechanics problems using a novel computational multiscale modelling approach. The multiscale approach employs a hierarchical coupling of the finite element method (FEM) and the discrete element method. It solves a boundary value problem at the continuum scale by FEM and derives the material point response from the discrete element method simulation attached to each Gauss point of the FEM mesh. The multiscale modelling framework not only helps successfully bypass phenomenological constitutive assumptions as required in conventional modelling approaches but also facilitates effective cross‐scale interpretation and understanding of soil behaviour. We examine the classical retaining wall and footing problems by this method and demonstrate that the simulating results can be well validated and verified by their analytical solutions. Furthermore, the study sheds novel multiscale insights into these classical problems and offers a new tool for geotechnical engineers to design and analyse geotechnical applications based directly upon particle‐level information of soils. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

15.
Parameter identification for lined tunnels in a viscoplastic medium   总被引：2，自引：0，他引：2  
This paper is dedicated to the identification of constitutive parameters of elasto‐viscoplastic constitutive law from measurements performed on deep underground cavities (typically tunnels). This inverse problem is solved by the minimization of a cost functional of least‐squares type. The exact gradient is computed by the direct differentiation method and the descent is done using the Levenberg–Marquardt algorithm. The method is presented for lined or unlined structures and is applied for an elastoviscoplastic constitutive law of the Perzyna class. Several identification problems are presented in one and two dimensions for different tunnel geometries. The used measurements have been obtained by a preliminary numerical simulation and perturbed with a white noise. The identified responses match the measurements. We also discuss the usage of the sensitivity analysis of the system, provided by the direct differentiation method, for the optimization of in situ monitoring. The sensitivity distribution in space and time assess the location of the measurements points as well as the time of observation needed for reliable identification. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

16.
In this paper, the seismic response of ‘infinitely’ long slopes is numerically analysed via the formulation of a 1D analytical/numerical model, in which the soil mechanical behaviour is assumed to be elasto‐perfectly viscoplastic and simple shear (SS) kinematical constraints are imposed. In order to simplify the problem, a theoretically based procedure to set up a fully 1D shear constitutive model is defined, within which the mechanical response of a multiaxial relationship is condensed. The use of a 1D shear constitutive model is aimed at reducing the number of unknowns and, therefore, the computational costs. In particular, the case of the Mohr–Coulomb yield criterion is considered, while an enhanced Taylor–Galerkin finite element algorithm is employed to simulate the seismic wave propagation within the soil stratum. The proposed ‘condensation’/calibration procedure captures both the ‘pseudo’‐hardening pre‐failure behaviour and the influence of dilation on the occurrence of strain‐localization, which characterize, under SS conditions, the static response of virgin perfectly plastic soils. The effectiveness of the conceived method is shown with reference to freshly deposited deposits, while, in the case of highly overconsolidated strata, some difficulties arise because of the brittle behaviour induced both by unloading and non‐associativeness. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

17.
The dynamic response of a soft soil layer of finite thickness under the mutual effects of flow and periodical disturbance at the free surface is discussed in this work. The homogeneous water is governed by potential theory and the soil layer obeys Biot's theory of poroelasticity. The boundary‐value problem is solved by an analytical algorithm, in which the wave number is found first. Secondly, the closed form solutions are found by a two‐parameter perturbation method with the boundary‐layer correction. The results are also compared with those of the poroelastic soil layer of infinite thickness to show the impermeable rigid boundary effect. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献   

18.
Data assimilation, using the particle filter and incorporating the soil‐water coupled finite element method, is applied to identify the yield function of the elastoplastic constitutive model and corresponding parameters based on the sequential measurements of hypothetical soil tests and an actual construction sequence. In the proposed framework of the inverse analysis, the unknowns are both the particular parameter within the exponential contractancy model, nE, which parameterizes various shapes for the yield function of the competing constitutive models, including the original/the modified Cam‐Clay models and in‐between models and the parameters of the corresponding constitutive model. An appropriate set, consisting of the yield function of the constitutive model and the parameters of the constitutive model, can be simultaneously identified by the particle filter to describe the most suitable soil behavior. To examine the validity of the proposed procedure, hypothetical and actual measurements for the displacements of a soil specimen were obtained for consolidated and undrained tests through a synthetic FEM computation and for consolidated and drained tests, respectively. After examining the applicability of the proposed procedure to these test results, the present paper then focuses on the actual measured data, ie, the settlement behavior including the lateral deformation of the Kobe Airport Island constructed on reclaimed land.  相似文献   

19.
A hybrid finite element method and differential quadrature method (DQM) is developed to estimate the dynamic response of two‐dimensional multilayered half‐spaces subjected to impulsive point loading. Nonreflecting absorbing boundary conditions consist of appropriate springs, and dampers are considered. The capabilities of the finite element method for solving boundary value problems with general domain, loading and systematic boundary treatment are combined with accurate and stable time marching capabilities of the DQM to develop an accurate and efficient numerical technique. The capability, efficiency, robustness and convergence of the DQM for solving the dynamic problem are demonstrated through numerical simulations of various half‐spaces with different time increments and layer arrangement. Also, comparison study when using Newmark's time integration scheme for the same problem is done. It can be concluded that the DQM as an unconditionally stable method is suitable for solving such a problem. Also, parametric study is performed to show the effect of the absorbing boundary conditions on the dynamic response. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

20.
An objective of this paper is to demonstrate that the small strain model developed by the authors can be incorporated into the conventional kinematic hardening plasticity framework to predict pre‐failure defor mations. The constitutive model described in this paper is constituted by three elliptical yield surfaces in triaxial stress space. Two inner surfaces are rotated ellipses of the same shape, representing the boundaries of the linear elastic and small strain regions, while the third surface is the modified Cam clay large‐scale yield surface. Within the linear elastic region, the soil behaviour is elastic with cross‐coupling between the shear and volumetric stress–strain components. Within the small strain region, the soil behaviour is elasto‐plastic, described by the kinematic hardening rule with an infinite number of loading surfaces defined by the incremental energy criterion. Within the large‐scale yield surface, the soil behaviour is elasto‐plastic, described by kinematic and isotropic hardening of the small strain region boundary. Since the yield surfaces have different shapes, the uniqueness of the plastic loading condition imposes a restriction on the ratio between their semi‐diameters. The model requires 12 parameters, which can be determined from a single consolidated undrained triaxial compression test. Copyright © 2000 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号