Elastoplastical analysis on cylinder of different moduli in tension and compression

The constitutive equation was deduced. The results were obtained from theelastic and elastoplastical analysis of the cylinder sample applied with different ratios ofload on both the inner and outer sides and different moduli ratio. The factors affecting thedevelopment of the plastic field of the cylinder with finite radius were shown. The resultsshow that the different moduli ratio is the most important factor in the development of theplastic zone. The slight fluctuation of the ratio will bring multiplied increment of the displacement, which may result in the final destruction of the engineering material.

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.

Unified description of different soils based on the superloading and subloading concepts

Geotechnical engineering often involves different types of geomaterials,such as sandy soil and clayey soil.Existing studies have confirmed that these soils have some common features,i.e.their mechanical behaviors depend not only on the inherent characteristics but also on their initial states.To describe the main mechanical behaviors of different soils within a simple and reasonable constitutive framework is of great significance for the numerical analysis on geotechnical engineering.This paper first introduces a model based on the concepts of superloading and subloading,which considers the“state dependence”(effects of overconsolidation and structure)of soil and only adds two material parameters compared with the Cam-Clay model.Secondly,conventional triaxial tests are systematically carried out on four types of soils(i.e.sand,silty clay,clay,and intermediate soil)with different initial void ratios,and the mechanical similarities and differences of these soils are discussed uniformly.After that,six material parameters of these soils are uniformly determined based on the concepts of superloading and subloading,and then used in constitutive calculations to verify the feasibility.The calculated results show a good agreement with test data,indicating that the model based on the concepts of superloading and subloading has great potential for describing the general mechanical behaviors of different soils within a unified framework.This work is expected to be applied to constitutive selection and parameter determination in the geotechnical numerical analysis of complex soil profiles.

FFT-Based Numerical Method for Nonlinear Elastic

In theoretical research pertaining to sealing, a contact model must be used to obtain the leakage channel. However, for elastoplastic contact, current numerical methods require a long calculation time. Hyperelastic contact is typically simplifed to a linear elastic contact problem, which must be improved in terms of calculation accuracy. Based on the fast Fourier transform, a numerical method suitable for elastoplastic and hyperelastic frictionless contact that can be used for solving two-dimensional and three-dimensional (3D) contact problems is proposed herein. The nonlinear elastic contact problem is converted into a linear elastic contact problem considering residual deformation (or the equivalent residual deformation). Results from numerical simulations for elastic, elastoplastic, and hyperelastic contact between a hemisphere and a rigid plane are compared with those obtained using the fnite element method to verify the accuracy of the numerical method. Compared with the existing elastoplastic contact numerical methods, the proposed method achieves a higher calculation efciency while ensuring a certain calculation accuracy (i.e., the pressure error does not exceed 15%, whereas the calculation time does not exceed 10 min in a 64 × 64 grid). For hyperelastic contact, the proposed method reduces the dependence of the approximation result on the load, as in a linear elastic approximation. Finally, using the sealing application as an example, the contact and leakage rates between complicated 3D rough surfaces are calculated. Despite a certain error, the simplifed numerical method yields a better approximation result than the linear elastic contact approximation. Additionally, the result can be used as fast solutions in engineering applications.

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.

Numerical Modeling of the Behaviour of a Road Structure on Compressible Soil: Case of the Road Section at the Beau-Rivage-Djassin Intersection

This document presents a study of the behaviour of a pavement structure on compressible soil and the evaluation of its durability. The objective of this study is to highlight the impact of taking into account the non-linear elastic behaviour of soils and granular materials in the design process. To this end, a numerical modelling of the pavement of the beau-rivage-Djassin crossroads section in Porto-Novo was carried out, based on a compressible soil whose behaviour will be considered elastoplastic. The subgrade soil on the section is made up of several sub-layers. The layer of soft, highly plastic clay was modelled according to a modified Cam Clay behaviour, a model of swelling clay soils. The fine sand layer and the granular layers of the structure are modelled according to Mohr-Coulomb behaviour. The loading is considered to be uniformly distributed according to the assumptions of the Burmister model in the French standard. A first verification with ALIZE allowed to validate the structure on the basis of the rutting deformation at the head of the platform εz = 359.6*10-6 which remains lower than the admissible deformation εz,adm = 360*10-6. The numerical calculation was carried out using the finite element method, the code of which is implemented under the PLAXIS v21 software. A comparative study with the results of the ALIZE design revealed that the numerically calculated strains εz = 585*10-6 are higher than those of ALIZE. These numerical strains, which are higher than the elastic strains, do not meet the validation criteria that the strains under loading must remain below the allowable strains. An evaluation of the pavement durability was carried out and it was found that the pavement would only last under traffic for 3 years before the first fatigue deformations appeared.

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.

基于摩擦补偿的伺服转台自抗扰控制策略研究

摩擦非线性扰动是影响伺服跟踪系统控制性能的主要因素之一。为提高转台伺服系统的跟踪性能,提出了一种基于Elastoplastic摩擦模型的改进自抗扰控制方法。首先,建立了转台伺服系统的状态空间模型;其次,采用Elastoplastic摩擦模型描述系统中的非线性摩擦扰动,并用遗传算法辨识了模型参数;最后,基于辨识获得的Elastoplastic摩擦模型,将位置误差和速度误差作为不同的参数分别应用到扩张状态观测器,设计了一种改进型自抗扰控制器。未引入摩擦补偿时的速度跟踪误差平均值约为0.0024 rad/s,而加入补偿后的速度跟踪误差平均值减少为0.00147 rad/s。仿真和实验结果表明,本文提出的控制方案能够提高转台伺服系统的跟踪性能,验证了所提出控制方法的有效性和鲁棒性。

Review of collapse triggering mechanism of collapsible soils due towetting

Loess soil deposits are widely distributed in arid and semi-arid regions and constitute about 10% of land area of the world.These soils typically have a loose honeycomb-type meta-stable structure that is susceptible to a large reduction in total volume or collapse upon wetting.Collapse characteristics contribute to various problems to infrastructures that are constructed on loess soils.For this reason,collapse triggering mechanism for loess soils has been of significant interest for researchers and practitioners all over the world.This paper aims at providing a state-of-the-art review on collapse mechanism with special reference to loess soil deposits.The collapse mechanism studies are summarized under three different categories,i.e.traditional approaches,microstructure approach,and soil mechanics-based approaches.The traditional and microstructure approaches for interpreting the collapse behavior are comprehensively summarized and critically reviewed based on the experimental results from the literature.The soil mechanics-based approaches proposed based on the experimental results of both compacted soils and natural loess soils are reviewed highlighting their strengths and limitations for estimating the collapse behavior.Simpler soil mechanics-based approaches with less parameters or parameters that are easy-to-determine from conventional tests are suggested for future research to better understand the collapse behavior of natural loess soils.Such studies would be more valuable for use in conventional geotechnical engineering practice applications.

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.

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.

Experimental study and modeling of hydromechanical behavior of concrete fracture

In this study, the hydromechanical behavior of a concrete fracture under coupled compressive and shear stresses was investigated. A special experimental device was designed to create a planar fracture in a cylindrical sample and to carry out different kinds of hydromechanical tests on the fracture. Four series of laboratory tests were performed on an ordinary concrete sample. Hydrostatic compression tests were first conducted to characterize the normal compressibility of the fracture. In the second series, direct shear tests were conducted on the fracture under different normal stresses. The maximal shear stress of the fracture was determined as a function of the normal stress. In the third series, fluid flow tests were carried out in view of characterizing the overall hydraulic conductivity of the fracture as a function of its opening and closure. Shear tests with a constant fluid pressure were finally performed to investigate the influence of fluid pressure on the deformation behavior of concrete fractures. Based on the experimental investigation, an elastoplastic model is proposed. This model takes into account the nonlinear elastic behavior of a fracture under normal compression and the plastic deformation and failure due to shear stress. The model was coupled with the classical Darcy's law to describe the fluid flow along the fracture by considering the variation of permeability with fracture aperture. Numerical results agree with experimental data from various laboratory tests.

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.

A practical constitutive theory based on egg-shaped function in elasto-plastic modeling for soft clay

An exploratory discussion is presented on the application of egg-shaped function in elasto-plastic constitutive analysis for soft clay.Two main tasks of the paper are:1)to propose a complete yield criterion based on egg-shaped function and supplement its definition in the deviatoric section,and then a yield criterion suitable for 3D stress conditions is obtained;2)to elaborate its numerical implementation based on the drained triaxial tests.During the above discussion,a non-associated flow rule is proposed,in which the stress-dilatancy relationship in most classical theory is replaced by a linear dependence between the stress state parameterηand the rotation angleγof the plastic potential surface.Thereafter,isotropic and kinematic hardening behavior is considered by employing the hardening parameter H,which can be expressed as the function of plastic work Wp.Finally,comparisons between numerical results and test data on Taizhou soft clay are made to verify the effectiveness of the proposed model.

The Anomalous Response of Elastoplatic Circular Plates Under Impulsive Loading

The anomalous nonlinear elastic, perfectly plastic response behaviors of circular plate subjected to short transverse pulse load is studied. The plate is assumed fixed-pin along the boundary. 'Anomalous' here means that the final deflection may be in the direction opposite that of the load. It has been found by detailed numerical analyses that there exists anomalous response in some narrow loading ranges, so called slots. By further calculations it is shown that this special dynamic behavior is related to coupling affects of internal forces and large plastic deformation after removal loading. Further plastic dissipation will be lead to anomalous dynamic response. This phenomena could be considered as the coupling of the geometry nonlinearity, material nonlinearity,elastic effects and the irrecoverable of the plastic deformation.

Shakedown Criterion Employing Actual Residual Stress Field and Its Application in Numerical Shakedown Analysis

Construction of the static admissible residual stress field and searching the optimal field are key tasks in the shakedown analysis methods applying the static theorem. These methods always meet dimension obstacles when dealing with complex problems. In this paper, a novel shakedown criterion is proposed employing actual residual stress field based on the static shakedown theorem. The actual residual stress field used here is produced under a specified load path, which is a sequence of proportional loading and unloading from zero to all the vertices of the given load domain. This ensures that the shakedown behavior in the whole load domain can be determined based on the theorem proposed by K6nig. The shakedown criterion is then implemented in numerical shakedown analysis, The actual residual stress fields are calculated by incremental finite element elastic-plastic analysis technique for finite deformation under the specified load path with different load levels. The shakedown behavior and the shakedown limit load are determined according to the proposed criterion. The validation of the criterion is performed by a benchmark shakedown example, which is a square plate with a central hole under biaxial loading. The results are consistent with existing results in the literatures and are validated by full cyclic elastic-plastic finite element analysis. The numerical shakedown analysis applying the proposed criterion avoids processing dimension obstacles and performing full cyclic elastic-plastic analysis under arbitrary load paths which should be accounted for appearing. The effect of material model and geometric changes on shakedown behavior can he considered conveniently.