Modelling of the elastoplastic behaviour of the bio-cemented soils using an extended Modified Cam Clay model

An elastoplastic constitutive model based on the Modified Cam Clay(MCC)model is developed to describe the mechanical behaviour of soils cemented via microbially induced calcite precipitation(MICP).It considers the increase of the elastic stiffness,the change of the yield surface due to MICP cementation and the degradation of calcium carbonate bonds during shearing.Specifically,to capture the typical contraction-dilation transition in MICP soils,the original volumetric hardening rule in the MCC model is modified to a combined deviatoric and volumetric hardening rule.The model could reproduce a series of drained triaxial tests on MICP-treated soils with different calcium carbonate contents.Further,we carry out a parametric study and observe numerical instability in some cases.In combination with an analytical analysis,our numerical modelling has identified the benefits and limitations of using MCCbased models in the simulation of MICP-cemented soils,leading to suggestions for further model development.

Finite Element Simulations on Failure Behaviors of Granular Materials with Microstructures Using a Micromechanics-Based Cosserat Elastoplastic Model

This paper presents a micromechanics-based Cosserat continuum model for microstructured granular materials.By utilizing this model,the macroscopic constitutive parameters of granular materials with different microstructures are expressed as sums of microstructural information.The microstructures under consideration can be classified into three categories:a medium-dense microstructure,a dense microstructure consisting of one-sized particles,and a dense microstructure consisting of two-sized particles.Subsequently,the Cosserat elastoplastic model,along with its finite element formulation,is derived using the extended Drucker-Prager yield criteria.To investigate failure behaviors,numerical simulations of granular materials with different microstructures are conducted using the ABAQUS User Element(UEL)interface.It demonstrates the capacity of the proposed model to simulate the phenomena of strain-softening and strain localization.The study investigates the influence of microscopic parameters,including contact stiffness parameters and characteristic length,on the failure behaviors of granularmaterials withmicrostructures.Additionally,the study examines themesh independence of the presented model and establishes its relationship with the characteristic length.A comparison is made between finite element simulations and discrete element simulations for a medium-dense microstructure,revealing a good agreement in results during the elastic stage.Somemacroscopic parameters describing plasticity are shown to be partially related to microscopic factors such as confining pressure and size of the representative volume element.

A new elastoplastic model for bolt-grouted fractured rock

Complexities in mechanical behaviours of rock masses mainly stem from inherent discontinuities,which calls for advanced bolt-grouting techniques for stability enhancement.Understanding the mechanical properties of bolt-grouted fractured rock mass(BGFR)and developing accurate prediction methods are crucial to optimize the BGFR support strategies.This paper establishes a new elastoplastic(E-P)model based on the orthotropic and the Mohr-Coulomb(M-C)plastic-yielding criteria.The elastic parameters of the model were derived through a meso-mechanical analysis of composite materials mechanics(CMM).Laboratory BGFR specimens were prepared and uniaxial compression test and variable-angle shear test considering different bolt arrangements were carried out to obtain the mechanical parameters of the specimens.Results showed that the anisotropy of BGFR mainly depends on the relative volume content of each component material in a certain direction.Moreover,the mechanical parameters deduced from the theory of composite materials which consider the short fibre effect are shown to be in good agreement with those determined by laboratory experiments,and the variation rules maintained good consistency.Last,a case study of a real tunnel project is provided to highlight the effectiveness,validity and robustness of the developed E-P model in prediction of stresses and deformations.

Triaxial elastoplastic damage constitutive model of unreinforced clay brick masonry wall

Due to differences in the properties of composition materials and construction techniques,unreinforced masonry is characterized by low strength,anisotropy,nonuniformity,and low ductility.In order to accurately simulate the mechanical behavior of unreinforced brick masonry walls under static and dynamic loads,a new elastoplastic damage constitutive model was proposed and the corresponding subroutine was developed based on the concrete material constitutive model.In the proposed constitutive model,the Rankine strength theory and the Drucker-Prager strength theory were used to define the tensile and compressive yield surface function of materials,respectively.Moreover,the stress updating algorithm was modified to consider the tensile plastic permanent deformation of masonry materials.To verify the accuracy of the proposed constitutive model,numerical simulations of the brick masonry under monotonic and cyclic uniaxial tension and compression loads were carried out.Comparisons among the numerical and theoretical and experimental results show that the proposed model can properly reflect the masonry material mechanical properties.Furthermore,the numerical models of four pieces of masonry walls with different mortar strengths were established.Low cyclic loadings were applied and the results show that the proposed constitutive model can properly simulate the wall shear failure characteristics,and the force-displacement hysteretic curves obtained by numerical simulation are in good agreement with the tests.Overall,the proposed elastic-plastic damage constitutive model can simulate the nonlinear behavior of unreinforced brick masonry walls very well,and can be used to predict the structural response of masonry walls.

A novel elastoplastic model for Yunnan sandstone under poly-axial loading

The lack of understanding of plastic hardening(softening)laws,especially under anisotropic stress conditions,results in inappropriate geotechnical management.Most of the yielding envelopes do not consider the effect of intermediate principal stress and the influence of Lode's angle.In addition,the application of plastic flow rules regarding yielding surfaces compromises the softening of rock internal friction as well as the influence of Lode's angle on the plastic potential.Moreover,the ductility to brittleness transition in the intermediate principal stress direction still requires a theoretical foundation.In this study,based on poly-axial testing results of Yunnan sandstone,we adopted a failure criterion with the intermediate principal stress proposed by Menétrey and Willam.The proposed new failure envelope was applied to capture the plastic evolution of rock samples.A plastic hardening-softening model is constructed,based on the framework of the plastic theory.The softening envelope is modified to better present the stress drop and considers the deterioration of rock internal friction in the post-peak stage of poly-axial loading.The differential of plastic potential according to the principal stresses is also modified,considering the rotation of Lode's angle in the poly-axial loading tests.The model results were compared with laboratory testing results,which showed great consistency across 9 different loading tests(5 under triaxial stress and 4 under poly-axial stress with 22 stress-strain curves in total).The induced brittleness by the intermediate principal stress is also well captured by the proposed model.

Calibration of an elastoplastic model of sand liquefaction using the swarm intelligence with a multi-objective function

According to post-seismic observations,spectacular examples of engineering failures can be ascribed to the occurrence of sand liquefaction,where a sandy soil stratum could undergo a transient loss of shear strength and even behave as a“liquid”.Therefore,correct simulation of liquefaction response has become a challenging issue in geotechnical engineering field.In advanced elastoplastic models of sand liquefaction,certain fitting parameters have a remarkable effect on the computed results.However,the identification of these parameters,based on the experimental data,is usually intractable and sometimes follows a subjective trial-and-error procedure.For this,this paper presented a novel calibration methodology based on an optimization algorithm(particle swarm optimization(PSO))for an advanced elastoplastic constitutive model.A multi-objective function was designed to adjust the global quality for both monotonic and cyclic triaxial simulations.To overcome computational problem probably appearing in simulation of the cyclic triaxial test,two interrupt mechanisms were designed to prevent the particles from wasting time in searching the unreasonable space of candidate solutions.The Dafalias model has been used as an example to demonstrate the main programme.With the calibrated parameters for the HN31 sand,the computed results were highly consistent with the laboratory experiments(including monotonic triaxial tests under different confining pressures and cyclic triaxial tests in two loading modes).Finally,an extension example is given for Ottawa sand F65,suggesting that the proposed platform is versatile and can be easily customized to meet different practical needs.

Elastoplastic response of functionally graded cemented carbides due to thermal loading

Finite dement formulations are used to simulate the evolution of the elastoplastic response of functionally graded cemented carbides （FGCC） due to thermal loading. The geometry of specimens is an axisymmetric solid cylinder with a two-dimensional gradient. The elastoplastic constitutive relationship is developed by constraint factors. Numerical results show that compressive stresses occur in the surface zone and tensile stresses in the cobalt rich zone when the temperature drops from the initial stress-free temperature of 800 to 0℃. The maximum value of the surface compressive stress is 254 MPa and the maximum value of the tensile stress is 252 MPa in the cobalt rich zones. When the cobalt concentration difference in the specimens is equal to or greater than 0.3, there is pronounced plastic flow in cobalt rich zone. When the temperature heats up from 0 to 800 ℃, the total plastic strain reaches 0.001 4. Plastic flow has a significant effect on the reduction of thermal stress concentration.

On existence and uniqueness of the solution of elastoplastic contact problems

Contact problems and elastoplastic problems are unified and described by the variational inequality formulation, in which the constraints of the constitutional relations for elastoplastic materials and the contact conditions are relaxed totally. First, the coerciveness of the functional is proved. Then the uniqueness of the solution of variational inequality for the elastoplastic contact problems is demonstrated. The existence of the solution is also demonstrated according to the sufficient conditions for the solution of the elliptic variational inequality. A mathematical foundation is developed for the variational extremum principle of elastoplastic contact problems. The developed variational extremum forms can give an effective and strict mathematical modeling to solve contact problems with mathematical programming.

Determination of Elastoplastic Mechanical Properties of the Weld and Heat Affected Zone Metals in Tailor-Welded Blanks by Nanoindentation Test

The elastoplastic mechanical properties of the weld and heat affected zone metals have comparatively major impact on the forming process of tailor-welded blanks. A few scholars investigated the elastoplastic mechanical properties of the weld and heat affected zone, but they only simply assumed that it was a uniform distribution elastoplastic material different from the base materials. Four types of tailor-welded blanks which consist of ST12 and 304 stainless steel plates are selected as the research objects, the elastoplastic mechanical properties of the tailor-welded blanks weld and heat affected zone metals are obtained based on the nanoindentation tests, and the Erichsen cupping tests are conducted by combining numerical simulation with physical experiment. The nanoindentation tests results demonstrate that the elastoplastic mechanical properties of the weld and heat affected zone metals are not only different from the base materials, but also varying between the weld metals and the heat affected zone metals. Comparing the Erichsen cupping test resulted from numerical with that from experimental method, it is found that the numerical value of Erichsen cupping test which consider the elastoplastic mechanical properties of the weld and heat affected zone metals have a good agreement with the experimental result, and the relative error is only 4.8%. The proposed research provides good solutions for the inhomogeneous elastoplastic mechanical properties of the tailor-welded blanks weld and heat affected zone metals, and improves the control performance of tailor-welded blanks forming accuracy.

Elastoplastic 3D Deformation and Stress Analysis of Strip Rolling

An elastoplastic method for analyzing the 3D deformation, stress and transverse distribution of tension stress during cold strip rolling is developed. The analysis is based on the elastoplastic variational principle in which a kinematically admissible velocity field is constructed with the lateral flow function as an unknown function. The stress distribution and volume strain distribution are obtained by solving the simultaneous equations formed by the longitudinal differential equation of equilibrium and constitutive equations. The lateral flow function is determined by minimizing the total energy dissipation rate. Experimental investigation was carried out on a reversible cold mill. The front tension stress distributions in cold rolled strips were measured by a multi roll segmented tension sensing shapemeter. The calculated results are in good agreement with the measured ones.

Finite element simulation of elastoplastic fieldnear crack tips and results for a central crackedplate of le-lhp material under tension

The elastoplastic field near crack tips is investigated through finite element simulation.A refined mesh model near the crack tip is proposed. In the mesh refining area, element size continuously varies from the nanometer scale to themicrometer scale and the millimeter scale. Graphics of the plastic zone, the crack tip blunting, and the deformed crack tip elements are given in the paper.Based on the curves of stress and plastic strain, closely near the crack tip, the stresssingularity index and the stress intensity factor,as well as the plastic strain singularity index and the plastic strain intensity factor are determined.Thestress and plastic strainsingular index vary with the load, while the dimensions of the stress and the plastic strain intensity factorsdependon the stress and the plastic strain singularity index, respectively. The singular field near the elastoplastic crack tip is characterized by the stress singularity index and the stress intensity factor, or alternativelythe plastic strain singularity index and the plastic strain intensityfactor.At the end of the paper, following Irwin’s concept of fracture mechanics,σδKσδKcriterion andεδQεδQcriterion are proposed.Besides, crack tip angle criterion is also presented.

Elastoplastic pipe-soil interaction analyses of partially-supported jointed water mains

Water distribution networks are essential components of water supply systems. The combination of pipe structural deterioration and mechanics leads to the failure of pipelines. A physical model for estimating the pipe failure must include both the pipe deterioration model and mechanics model. Winkler pipe-soil interaction (WPSI), an analytical mechanics model developed by Rajani and Tesfamariam (2004), takes external and internal loads, temperature changes, loss of bedding support, and the elastoplastic effect of soil into consideration. Based on the WPSI model, a method to evaluate the elastic and plastic areas was proposed in the present study. An FEM model based on pipe-soil interaction (PSI) element was used to verify the analytical model. Sensitivity analyses indicate that the soft soil, long pipe and high temperature induced the axial plastic deformation more likely, which, however, may not occur in normal scenarios. The soft soil, pipes in small diameters, long unsupported bedding are prone to form flexural plastic area. The results show that the pipes subjected to the same loads have smaller stresses in the elastoplastic analysis than elastic analysis. The difference, however, is slight.

Two-way coupled analysis of lithium diffusion and diffusion induced finite elastoplastic bending of bilayer electrodes in lithium-ion batteries

A fully coupling model for the diffusion induced finite elastoplastic bending of bilayer electrodes in lithium-ion batteries is proposed. The effect of the mechanical stress on the lithium diffusion is accounted for by the mechanical part of the chemical potential derived from the Gibbs free energy along with the logarithmic stress and strain. Eight dimensionless parameters, governing the stress-assisted diffusion and the diffusion induced elastoplastic bending, are identified. It is found that the finite plasticity starting from the interface of the bilayer increases the chemical potential gradient and thereby facilitates the lithium diffusion. The full plastic flow makes the abnormal lithium concentration distribution possible, i.e., the concentration at the lithium inlet can be lower than the concentration at the interface(downstream). The increase in the thickness of the active layer during charging is much larger than the eigen-stretch due to lithiation, and this excess thickening is found to be caused by the lithiation induced plastic yield.

Analysis of elastoplasticity problems using an improved complex variable element-free Galerkin method

In this paper, based on the conjugate of the complex basis function, a new complex variable moving least-squares approximation is discussed. Then using the new approximation to obtain the shape function, an improved complex variable element-free Galerkin（ICVEFG） method is presented for two-dimensional（2D） elastoplasticity problems. Compared with the previous complex variable moving least-squares approximation, the new approximation has greater computational precision and efficiency. Using the penalty method to apply the essential boundary conditions, and using the constrained Galerkin weak form of 2D elastoplasticity to obtain the system equations, we obtain the corresponding formulae of the ICVEFG method for 2D elastoplasticity. Three selected numerical examples are presented using the ICVEFG method to show that the ICVEFG method has the advantages such as greater precision and computational efficiency over the conventional meshless methods.

New finite strain elastoplastic equations for accurately and explicitly simulating pseudoelastic-to-plastic transition effects of shape memory alloys

A new finite strain elatoplastic J2-flow model with coupling effects of both isotropic and anisotropic hardening is proposed with the co-rotational logarithmic rate.In terms of certain single-variable shape functions representing uniaxial loading and unloading curves,explicit multi-axial expressions for the three hardening quantities incorporated in the new model proposed are derived in unified forms for the purpose of automatically and accurately simulating complex pseudoelastic-to-plastic transition effects of shape memory alloys(SMAs)under multiple loading-unloading cycles.Numerical examples show that with only a single parameter of direct physical meaning for each cycle,accurate and explicit simulations may be achieved for extensive data from multiple cycle tests.

Parameter Characterisation of Asphalt Mixtures for Visco-Elastoplastic Analysis

Asphalt mixtures exhibit strong viscous properties under repetitive loads. This phenomenon can be simulated in creep and recovery tests. By applying the visco elastoplastic model proposed, data recorded in the tests are interpreted. It is emphasised that applicability of the visco elastoplastic model depends on the proper method of characterising parameters involved in the constitutive equations. Since two or more strain components of elasticity and viscosity coexist during the loading or unloading, a measuring system of two loggers is introduced to decouple these components. Test data are analysed in three steps as elastic, visco elastic and visco plastic evaluations, using the records collected at the moment of unloading, during the periods of recovery and creep respectively. Factors that may influence the accuracy of data analysis are also discussed.

Unified elastoplastic finite difference and its application

Two elastoplastic constitutive models based on the unified strength the- ory （UST） are established and implemented in an explicit finite difference code, fast Lagrangian analysis of continua （FLAC/FLAC3D）, which includes an associated/non- associated flow rule, strain-hardening/softening, and solutions of singularities. Those two constitutive models are appropriate for metallic and strength-different （SD） materials, respectively. Two verification examples are used to compare the computation results and test data using the two-dimensional finite difference code FLAC and the finite element code ANSYS, and the two constitutive models proposed in this paper are verified. Two application examples, the large deformation of a prismatic bar and the strain-softening be- havior of soft rock under a complex stress state, are analyzed using the three-dimensional code FLAC3D. The two new elastoplastic constitutive models proposed in this paper can be used in bearing capacity evaluation or stability analysis of structures built of metallic or SD materials. The effect of the intermediate principal stress on metallic or SD mate- rial structures under complex stress states, including large deformation, three-dimensional and non-association problems, can be analyzed easily using the two constitutive models proposed in this paper.

Multi-scale elastoplastic mechanical model and microstructure damage analysis of solid expandable tubular

We present an in-depth study of the failure phenomenon of solid expandable tubular (SET) due to large expansion ratio in open holes of deep and ultra-deep wells. By examining the post-expansion SET, lots of microcracks are found on the inner surface of SET. Their morphology and parameters such as length and depth are investigated by use of metallographic microscope and scanning electron microscope (SEM). In addition, the Voronoi cell technique is adopted to characterize the multi-phase material microstructure of the SET. By using the anisotropic elastoplastic material constitutive model and macro/microscopic multi-dimensional cross-scale coupled boundary conditions, a sophisticated and multi-scale finite element model (FEM) of the SET is built successfully to simulate the material microstructure damage for different expansion ratios. The microcrack initiation and growth is simulated, and the structural integrity of the SET is discussed. It is concluded that this multi-scale finite element modeling method could effectively predict the elastoplastic deformation and the microscopic damage initiation and evolution of the SET. It is of great significance as a theoretical analysis tool to optimize the selection of appropriate tubular materials and it could be also used to substantially reduce costly failures of expandable tubulars in the field. This numerical analysis is not only beneficial for understanding the damage process of tubular materials but also effectively guides the engineering application of the SET technology.

Elastoplastic analysis of knee bracing frame

The knee bracing steel frame (KBF) is a new kind of energy dissipating frame, which combines excellent ductility and lateral stiffness. As the structural fuse of the frame, the knee element will yield first during a severe earthquake so that no damage occurs to the major structural members and the rehabilitation is easy and economical. To help fully understand the relations be- tween its seismic performance and the structural parameters, systematic elastoplastic analysis of the KBF structure with finite element method was conducted in this work. Finally, general design recommendations were made according to the results of the analysis.

Stability Loss of Rotating Elastoplastic Discs of the Specific Form

A method of calculating a possible stability loss by a rotating circular annular disc of variable thickness is suggested within the theory of perfect plasticity with the help of small parameter method. A characteristic equation for a critical radius of a plastic zone is obtained as a first approximation. The formula for the critical angular velocity, determining the stability loss of the disc according to the self-balanced form, is derived. The method using which we can take into account the disc’s geometry and loading parameters is also specified. The efficiency of the proposed method is shown in Section 5 while considering an illustrative example. The values of critical angular velocity of rotating are found numerically for different parameters of the disc.