A THERMO-PLASTIC/VISCOPLASTIC DAMAGE MODEL FOR GEOMATERIALS

A thermo-plastic/viscoplastic damage coupled model was formulated to describe the time independent and time dependent behaviors of geomaterials under temperature effect. The plastic strain was divided into instantaneous plastic strain and creep plastic strain. To take temperature effect into acconnt, a temperature variable was introduced into the instantaneous and creep plastic behavior descriptions and damage characterization, and a linear thermal expansion law was used in constitutive equation formulation. According to the mechanical behavior of rock salt, a specific model was proposed based on the previous model and applied to Avery rock salt, in which the numerical results obtained from our model had a good agreement with the data from experiments.

Modelling the viscoplastic behaviour of Callovo-Oxfordian claystone with consideration of damage effect

In order to evaluate the performance of deep geological disposal of radioactive waste,an underground research laboratory(URL)was constructed by Andra in the Callovo-Oxfordian(COx)claystone formation at the Meuse/Haute-Marne(MHM).The construction of URL induced the excavation damage of host formations,and the ventilation in the galleries desaturated the host formation close to the gallery wall.Moreover,it is expected that the mechanical behaviour of COx claystone is time-dependent.This study presents a constitutive model developed to describe the viscoplastic behaviour of unsaturated and damaged COx claystone.In this model,the unsaturation effect is considered by adopting the Bishop effective stress and the van Genuchten(VG)water retention model.In terms of the viscoplastic behaviour,the nonstationary flow surface(NSFS)theory for unsaturated soils is used with consideration of the coupled effects of strain rate and suction on the yield stress.A progressive hardening law is adopted.Meanwhile,a non-associated flow rule is used,which is similar to that in Barcelona basic model(BBM).In addition,to describe the damage effect induced by suction change and viscoplastic loading,a damage function is defined based on the crack volume proportion.This damage function contains two variables:unsaturated effective stress and viscoplastic volumetric strain,with the related parameters determined based on the mercury intrusion porosimetry(MIP)tests.For the model validation,different tests on COx claystone under different loading paths are simulated.Comparisons between experimental and simulated results indicated that the present model is able to well describe the viscoplastic behaviour of damaged COx claystone,including swelling/shrinkage,triaxial extension and compression,and triaxial creep.

Viscoplastic solutions of time-dependent deformation for tunnels in swelling rock mass considering stress release

Excavation and control of tunneling responses in swelling soft-rock tunnels of Sichuan-Tibet railway under seepage conditions were studied.For this,a fractional viscoplastic(FVP)model for swelling soft rocks was established by introducing Abel dashpot and unsteady viscosity coefficient,considering additional swelling deformation and damage of rock caused by humidity effect.In view of the FVP model,the viscoplastic deformation solutions for rock mass surrounding tunnel under seepage conditions were derived and long-term mechanical responses of swelling rocks upon tunnel excavation were analyzed.Next,a stress release coefficient considering seepage and creep was proposed,based on which control responses considering stress release and failure mechanism of stress release measures were analyzed.The results showed that:(i)The one-dimensional(1D)FVP model has a good application for swelling rock and the three-dimensional(3D)FVP model could well describe the whole creep process of rock mass despite a much higher creep attenuation rate in the first stage of creep;and(ii)An appropriate stress release and deformation of surrounding rocks could effectively reduce the supporting resistance.However,upon a large stress release,the radius of plastic region could increase significantly,and the strength of the surrounding rock mass decreases greatly.The proposed solution could provide a theoretical framework for capturing the excavation and support responses for tunneling in swelling rock mass in consideration of time effect.

Viscoplastic friction and microstructural evolution behavior of laser-clad Co-Cr-Ni-Mo coating

The viscoplastic friction and nanostructure formation mechanism of laser-clad Co-based coating were studied by rotary friction between laser-clad Co-Cr-Ni-Mo coating and WC-Co rod.The friction coefficient,friction interface temperature and axial displacement—time curves during rotary friction process were measured.The results showed that all the curves firstly experienced rising stage and then steady stage.The rising stage corresponded to sliding friction while the steady stage corresponded to viscoplastic friction.After viscoplastic friction processing,three typical zones of viscoplastic deformation zone,thermo-mechanically affected zone,and original laser-clad zone can be observed successively from the friction surface to the interior.The viscoplastic deformation significantly crushed the network M23C7 phase in original laser-clad zone and made it dispersively distributed with equiaxial shape and in nano-scale.The viscoplastic zone,in width of 37-131 μm,is mainly characterized by refined M23C7 and α-Co phase with grain size bellow 50 nm,and even a small quantity of amorphous.Thus,the hardness of viscoplastic zone about HV997 was improved compared with the hardness of original laser-clad zone about HV600.

Limit analysis of viscoplastic thick-walled cylinder and spherical shell under internal pressure using a strain gradient plasticity theory

Plastic limit load of viscoplastic thick-walled cylinder and spherical shell subjected to internal pressure is investigated analytically using a strain gradient plasticity theory. As a result, the current solutions can capture the size effect at the micron scale. Numerical results show that the smaller the inner radius of the cylinder or spherical shell, the more significant the scale effects. Results also show that the size effect is more evident with increasing strain or strain-rate sensitivity index. The classical plastic-based solutions of the same problems are shown to be a special case of the present solution.

A STUDY OF B-P UNIFIED VISCOPLASTIC CONSTITUTIVE MODEL

By combining the Bodner-Partom constitutive model and equivalent stressfunction, finite element methods and program on analyzing non-elastic deformation and stress forthermal viscoplastic material are studied in this paper, and it's the first time that this materialmodel is used in a kind of engineering software-MARC. Thermal viscoplastic behavior of hightemperature alloy GH536 specimen with gap is analyzed by this program. The research results show itis feasible to analyze thermal viscoplastic behavior of specimen or structure by applying B-P model.

RIGID-PLASTIC/RIGID-VISCOPLASTIC FEM BASED ON LINEAR PROGRAMMING—THEORETICAL MODELING AND APPLICATION FOR AXISYMMETRICAL PROBLEMS

Compared with the traditional rigid plastic/rigid viscoplastic(RP/RVP) FEM(based on iteration solution),RP/RVP FEM based on linear programming (LP) has some remarkable advantages,such as it's free of convergence problem and its convenience in contact,rigid zone,and friction force treatment.The numerical model of RP/RVP FEM based on LP for axisymmetrical metal forming simulation is studied,and some related key factors and its treatment methods in formulation of constraint condition are proposed.Some solution examples are provided to validate its accuracy and efficiency.

ELASTIC-VISCOPLASTIC FIELDS NEAR THE TIP OF A PROPAGATING CRACK UNDER ANTI-PLANE SHEAR

The elastic-viscoplastic model proposed by Bingham was used to analyse the stress and strain surrounding the tip of a propagating crack under antiplane shear.The proper displacement pattern was given;the asymptotic equations were derived and solved numerically.The analysis and calculation show that for smaller viscosity the crack-tip possesses logarthmic singularity,and for larger viscosity it possesses power-law singularity.In critical case,the two kinds of singularity are consistent with each other.The result revealed the important role of viscosity for crack-tip field.

OPTIMAL DESIGN PROCEDURE BASED ON VISCOPLASTIC MATERIAL BEHAVIOUR

The optimization problems belong to the family of the most important engineering problems. The objective of this paper is to proposed a very efficient numerical algo- rithm which is involved in an optimal design procedure in the field of viscoplasticity phenomena. Because of the complexity of design procedure, particularly in the field of viscoplas- ticity, the finite element method is used. This method provides a good structural discretization and very efficient mathematical describing of time-rate effects. Accord- ing to the design specifications, which will dictate a limiting design strain, using the proposed algorithm, the design life of the considered structure will be covered. For the justification of the proposed method an example is presented. On the basis of specifications data and environment conditions, an operating pressure of a pressure vessel is obtained.

Thermo-plasticity of high-strength and high-ductility 7050 aluminum ingot

Thermo-plasticity of homogenized 7050 aluminum ingot was investigated by instantaneous tensile tests conducted at different temperatures. The results show that, with the increase of testing temperatures, the strength decreases, and the plasticity increases firstly and then decreases in homogenized 7050 ingot. When the studied alloy is deformed between 380℃ and 420℃, the deformation resistance is lower and plasticity is better. And the actual heating temperature for ingot before hot extrusion should be controlled between 360 ~C and 400 ~C. At low tensile temperatures, the deformation structure is mainly composed of dislocation substructure. With the increase of testing temperatures, transgranular fracture transforms into intergranular fracture progressively during deformation. At high tensile temperatures, the grain boundaries are weakened, deformation is concentrated at the grain boundaries and the re-orientation of equilibrium phases at grain boundaries appears.

An elasto-plastic and viscoplastic damage constitutive model for dilatancy and fracturing behavior of soft rock squeezing deformation

Soft rock squeezing deformation mainly consists of pre-peak damage-dilatancy and post-peak fracture-bulking at the excavation unloading instant,and creep-dilatancy caused by time-dependent damage and fracturing.Based on the classic elastoplastic and Perzyna over-stress viscoplastic theories,as well as triaxial unloading confining pressure test and triaxial unloading creep test results,an elastoplastic and viscoplastic damage constitutive model is established for the short-and long-term dilatancy and fracturing behavior of soft rock squeezing deformation.Firstly,the criteria for each deformation and failure stage are expressed as a linear function of confining pressure.Secondly,the total damage evolution equation considering time-dependent damage is proposed,including the initial damage produced at the excavation instant,in which the damage variable increases exponentially with the lateral strain,and creep damage.Thirdly,a transient five-stages elasto-plastic constitutive equation for the short-term deformation after excavation that comprised of elasticity,pre-peak damage-dilatancy,post-peak brittle-drop,linear strain-softening,and residual perfectly-plastic regimes is developed based on incremental elasto-plastic theory and the nonassociated flow rule.Fourthly,regarding the timedependent properties of soft rock,based on the Perzyna viscoplastic over-stress theory,a viscoplastic damage model is set up to capture creep damage and dilatancy behavior.Viscoplastic strain is produced when the stress exceeds the initial static yield surface fs;the distance between the static yield surface fs and the dynamic yield surface fd determines the viscoplastic strain rate.Finally,the established constitutive model is numerically implemented and field applied to the-848 m belt conveyer haulage roadway of Huainan Panyidong Coal Mine.Laboratory test results and in-situ monitoring results validate the rationality of the established constitutive model.The presented model takes both the transient and time-dependent damage and fracturing into consideration.

Perfect elastic-viscoplastic field at mode I dynamic propagating crack-tip

The viscosity of material is considered at propagating crack-tip. Under the assumption that the artificial viscosity coefficient is in inverse proportion to power law of the plastic strain rate, an elastic-viscoplastic asymptotic analysis is carried out for moving crack-tip fields in power-hardening materials under plane-strain condition. A continuous solution is obtained containing no discontinuities. The variations of numerical solution are discussed for mode Ⅰ crack according to each parameter. It is shown that stress and strain both possess exponential singularity. The elasticity, plasticity and viscosity of material at crack-tip only can be matched reasonably under linear-hardening condition. And the tip field contains no elastic unloading zone for mode I crack. It approaches the limiting case, crack-tip is under ultra-viscose situation and energy accumulates, crack-tip begins to propagate under different compression situations.

CRACK PROPAGATION IN POLYCRYSTALLINE ELASTIC-VISCOPLASTIC MATERIALS USING COHESIVE ZONE MODELS

Cohesive zone model was used to simulate two-dimensional plane strain crack propagation at the grain level model including grain boundary zones. Simulated results show that the original crack-tip may not be separated firstly in an elastic-viscoplastic polycrystals. The grain interior＇s material properties （e.g. strain rate sensitivity） characterize the competitions between plastic and cohesive energy dissipation mechanisms. The higher the strain rate sensitivity is, the larger amount of the external work is transformed into plastic dissipation energy than into cohesive energy, which delays the cohesive zone rupturing. With the strain rate sensitivity decreased, the material property tends to approach the elastic-plastic responses. In this case, the plastic dissipation energy decreases and the cohesive dissipation energy increases which accelerates the cohesive zones debonding. Increasing the cohesive strength or the critical separation displacement will reduce the stress triaxiality at grain interiors and grain boundaries. Enhancing the cohesive zones ductility can improve the matrix materials resistance to void damage.

Coupled viscoplasticity damage constitutive model for concrete materials

A coupled viscoplasticity damage constitutive model for concrete materials is developed within the framework of irreversible thermodynamics. Simultaneously the Hehnholtz free energy function and a non-associated flow potential function are given, which include the internal variables of kinematic hardening, isotropic hardening and damage. Results from the numerical simulation show that the model presented can describe the deformatioa properties of the concrete without the formal hypotheses of yield criterion and failure criteria, such as the volume dilatancy under the compression, strain-rate sensitivity, stiffness degradation and stress-softening behavior beyond the peak stress which are brought by damages and fractures. Moreover, we could benefit from the application of the finite element method based oi1 this model under complex loading because of not having to choose different constitutive models based on the deformation level.

CYCLIC VISCOPLASTIC DESCRIPTION FOR THE MEAN STRAIN EFFECT

ＣＹＣＬＩＣＶＩＳＣＯＰＬＡＳＴＩＣＤＥＳＣＲＩＰＴＩＯＮＦＯＲＴＨＥＭＥＡＮＳＴＲＡＩＮＥＦＦＥＣＴＣＹＣＬＩＣＶＩＳＣＯＰＬＡＳＴＩＣＤＥＳＣＲＩＰＴＩＯＮＦＯＲＴＨＥＭＥＡＮＳＴＲＡＩＮＥＦＦＥＣＴ￥ＹａｎｇＸｉａｎｊｉｅ（ＳｏｕｔｈｗｅｓｔＪ...

Impact Buckling and Post-Buckling of Elasto-Viscoplastic Cylindrical Shell Under Torsion

The dynamic torsional buckling and post-buckling of elasto-plastic circular cylindrical shell is studied. By the usage of the Bodner-Partom constitutive relation, the present visco-plastic problem is treated as an orthotropic elastic problem. Based on the nonlinear large deflection shell therory, the governing dynamic buckling equation is solved by the Runge-Kutta method, and the critical buckling stress is determined by the B-R criterion. The present paper is mainly focused on the imperfection sensitivity, the strain rate sensitivity, and the dynamic post-buckling characters for the shell loaded dynamically with a constant shear strain rate. Numerical results show that the critical shear stress and imperfection sensitivity will increase with the increase of strain rate. It is found that the elastic buckling stress is more sensitive to imperfection than the visco-plastic buckling. The present research also reveals some dynamic post-buckling characters for the cylindrical shell, and it is shown that the shell will vibrate with the amplitude proportional to the strain rate in elastic buckling but improportional to the strain rate in visco-plastic buckling. An interesting phenomenon is found that when the strain rate is sufficiently high, there is no dynamic visco-plastic buckling occurrence and the shell produces a harmonic oscillation with a very small amplitude.

SYMMETRIC FORMULATION OF TANGENTIAL STIFFNESSES FOR NON-ASSOCIATED VISCOPLASTICITY WITH AN IMPLICIT TIME INTEGRATION SCHEME

A numerical scheme is presented which enables the use of symmetric equation solvers in tangential stiffness programs for non-associated viscoplastic materials.

VISCOPLASTIC SOLUTION TO FIELD AT STEADILY PROPAGATING CRACK TIP IN LINEARHARDENING MATERIALS

An elastic-viscoplastic constitutive model was adopted to analyze asymptotically the tip-field of moving crack in linear-hardening materials under plane strain condition. Under the assumption that the artificial viscosity coefficient was in inverse proportion to power law of the rate of effective plastic strain, it is obtained that stress and strain both possess power law singularity and the singularity exponent is uniquely determined by the power law exponent of the rate of effective plastic strain. Variations of zoning structure according to each material parameter were discussed by means of numerical computation for the tip-field of mode II dynamic propagating crack, which show that the structure of crack tip field is dominated by hardening coefficient rather than viscosity coefficient. The secondary plastic zone can be ignored for weak hardening materials while the secondary plastic zone and the secondary elastic zone both have important influence on crack tip field for strong hardening materials. The dynamic solution approaches to the corresponding quasi-static solution when the crack moving speed goes to zero, and further approaches to the HR （Hui-Riedel） solution when the hardening coefficient is equal to zero.

Viscoplastic Torsional Buckling of Cylindrical Shell

The Bodner-Partom constitutive equation is used to study the viscoplastic torsional buckling of perfect cylindrical shell. By treating the viscoplastic shell as an orthotropic shell at each moment and neglecting the inertia term, the critical torque is determined from a set of homogeneous linear equations. The strain rate sensitivity is mainly discussed in the present paper.

An Extended Chaboche's Viscoplastic Law at Finite Strains: Theoretical and Numerical Aspects

This paper presents a newly extended Chaboche’s viscoplastic law at finite strains, so that the classical Chaboche’s theories can be applied to the physical and numerical simulation of metals processing and behavior description of spatial metal structures. The extension is based on a new dissipation inequality at finite strains. The evolution equations are formulated in terms of the corotational rates of the logarithmic elastic strain and the strain-like internal variable conjugate to the back stress as well as the material time derivative of the accumulated plastic strain. The stress equation is expressed on the hyperelastic theory. Therefore, the possible inconsistency with elasticity, caused by the hypoelastic equations, is completely removed. A set of numerical examples with finite deformations are presented to prove the effectivities of the new model and numerical algorithms.