A multi-parameter nonlinear elasto-plastic constitutive model which can fully capture the three typical features of stress-strain response, linearity, plasticity-like stress plateau and densification phases was develo...A multi-parameter nonlinear elasto-plastic constitutive model which can fully capture the three typical features of stress-strain response, linearity, plasticity-like stress plateau and densification phases was developed. The functional expression of each parameter was determined using uniaxial compression tests for aluminum alloy foams. The parameters of the model can be systematically varied to describe the effect of relative density which may be responsible for the changes in yield stress and hardening-like or softening-like behavior at various strain rates. A comparison between model predictions and experimental results of the aluminum alloy foams was provided to validate the model. It was proved to be useful in the selection of the optimal-density and energy absorption foam for a specific application at impact events.展开更多
In this paper, a method for the design optimization of elasto-plastic truss structures is proposed based on parametric variational principles (PVPs). The optimization aims to find the minimum weight/volume solution ...In this paper, a method for the design optimization of elasto-plastic truss structures is proposed based on parametric variational principles (PVPs). The optimization aims to find the minimum weight/volume solution under the constraints of allowable node displacements. The design optimization is a formulation of mathematical programming with equilibrium constraints (MPECs). To overcome the numerical difficulties of the complementary constraints in optimization, an iteration process, comprising a quadratic programming (QP) and an updating process, is employed as the optimization method. Furthermore, the elasto-plastic buckling of truss mem- bers is considered as a constraint in design optimization. A combinational optimization strategy is proposed for the displacement constraints and the buckling constraint, which comprises the method mentioned above and an optimal criterion. Three numerical examples are presented to show the validity of the methods proposed.展开更多
Based on the elasto-plastic mechanics and continuum damage theory, a yield criterion related to spherical tensor of stress is proposed to describe the mixed hardening of damaged orthotropic materials. Its dimensionles...Based on the elasto-plastic mechanics and continuum damage theory, a yield criterion related to spherical tensor of stress is proposed to describe the mixed hardening of damaged orthotropic materials. Its dimensionless form is isomorphic with the Mises criterion for isotropic materials. Furthermore, the incremental elasto-plastic damage constitutive equations and damage evolution equations are established. Based on the classical nonlinear plate theory, the incremental nonlinear equilibrium equations of orthotropic thin plates considering damage effect are obtained, and solved with the finite difference and iteration methods. In the numerical examples, the effects of damage evolution and initial deflection on the elasto-plastic postbuckling of orthotropic plates are discussed in detail.展开更多
To gain a better understanding about texture evolution during rolling process of AZ31 alloy, polycrystalline plasticity model was implemented into the explicit FE package, ABAQUS/Explicit by writing a user subroutine ...To gain a better understanding about texture evolution during rolling process of AZ31 alloy, polycrystalline plasticity model was implemented into the explicit FE package, ABAQUS/Explicit by writing a user subroutine VUMAT. For each individual grain in the polycrystalline aggregate, the rate dependent model was adopted to calculate the plastic shear strain increment in combination with the Voce hardening law to describe the hardening response, the lattice reorientation caused by slip and twinning were calculated separately due to their different mechanisms. The elasto-plastic self consistent (EPSC) model was employed to relate the response of individual grain to the response of the polycrystalline aggregate. Rolling processes of AZ31 sheet and as-cast AZ31 alloy were simulated respectively. The predicted texture distributions are in aualitative a^reement with experimental results.展开更多
The dynamic buckling of elasto-plastic cylindrical shells under axial fluid-solid impact is investigated theoretically. A simplified liquid- gas- structure model is given to approximately imitate the problem. The basi...The dynamic buckling of elasto-plastic cylindrical shells under axial fluid-solid impact is investigated theoretically. A simplified liquid- gas- structure model is given to approximately imitate the problem. The basic equation of the structure is derived from a minimum principle in dynamics of elasto-plastic continua at finite deformation, and the flow theory of plasticity is employed. The liquid is incompressible and the gas is compressed adiabatically. A number of numerical results are presented and the characteristics of the buckling behavior under fluid-solid impact are illustrated.展开更多
The present paper aims to establish a versatile strength theory suitable for elasto-plastic analysis of underground tunnel surrounding rock. In order to analyze the effects of intermediate principal stress and the roc...The present paper aims to establish a versatile strength theory suitable for elasto-plastic analysis of underground tunnel surrounding rock. In order to analyze the effects of intermediate principal stress and the rock properties on its deformation and failure of rock mass, the generalized nonlinear unified strength theory and elasto-plastic mechanics are used to deduce analytic solution of the radius and stress of tunnel plastic zone and the periphery displacement of tunnel under uniform ground stress field. The results show that: intermediate principal stress coefficient b has significant effect on the plastic range,the magnitude of stress and surrounding rock pressure. Then, the results are compared with the unified strength criterion solution and Mohr–Coulomb criterion solution, and concluded that the generalized nonlinear unified strength criterion is more applicable to elasto-plastic analysis of underground tunnel surrounding rock.展开更多
The relationship between the Hoek-Brown parameters and the mechanical response of circular tunnels is il-lustrated. Closed-form and approximate solutions are given for the extent of the plastic zone and the stress and...The relationship between the Hoek-Brown parameters and the mechanical response of circular tunnels is il-lustrated. Closed-form and approximate solutions are given for the extent of the plastic zone and the stress and dis-placement fields under axisymmetrical and asymmetric stress conditions. For the same rock masses and under axisym-metrical stress conditions,the radius of the plastic zone in terms of Hoek-Brown criterion is generally an approximation of the radius in terms of the Mohr-Coulomb criterion. The radius in terms of the Hoek-Brown criterion is larger under low stress conditions. For poor quality rock masses (GSI<25),measures (such as grouting,setting rock bolts,etc.) that improve the GSI of rock masses are effective in improving the stability of tunnels. It is not advisable to improve the sta-bility of the tunnels by providing a small support resistance p through shotcrete,except for very poor quality jointed rock masses. Without reference to the quality of the rock mass,the disturbance factor D should not less than 0.5. Meas-ures which disturb rock masses during tunnel construction should be taken carefully when the tunnel depth increases.展开更多
In this work, the finite element analysis of the elasto-plastic plate bending problems is carried out using transition rectangular plate elements. The shape functions of the transition plate elements are derived based...In this work, the finite element analysis of the elasto-plastic plate bending problems is carried out using transition rectangular plate elements. The shape functions of the transition plate elements are derived based on a practical rule. The transition plate elements are all quadrilateral and can be used to obtain efficient finite element models using minimum number of elements. The mesh convergence rates of the models including the transition elements are compared with the regular element models. To verify the developed elements, simple tests are demonstrated and various elasto-plastic problems are solved. Their results are compared with ANSYS results.展开更多
A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed.As compared with the existing discontinuous Galerkin finite element methods,the distinct fea...A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed.As compared with the existing discontinuous Galerkin finite element methods,the distinct feature of the proposed method is that the continuity of the displacement vector at each discrete time instant is automatically ensured,whereas the discontinuity of the velocity vector at the discrete time levels still remains.The computational cost is then obviously reduced, particularly,for material non-linear problems.Both the implicit and explicit algorithms to solve the derived formulations for material non-linear problems are developed.Numerical results show a good performance of the present method in eliminating spurious numerical oscillations and providing with much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain.展开更多
Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).How...Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).However,the deformation field obtained by GSRM could not reflect the real deformation of a slope when the slope became unstable.For most slopes,failure occurs once the strength of some regional soil is sufficiently weakened; thus,the local strength reduction method(LSRM)was proposed to analyze slope stability.In contrast with GSRM,LSRM only reduces the strength of local soil,while the strength of other soil remains unchanged.Therefore,deformation by LSRM is more reasonable than that by GSRM.In addition,the accuracy of the slope's deformation depends on the constitutive model to a large degree,and the variable-modulus elasto-plastic model was thus adopted.This constitutive model was an improvement of the Duncan–Chang model,which modified soil's deformation modulus according to stress level,and it thus better reflected the plastic feature of soil.Most importantly,the parameters of the variable-modulus elasto-plastic model could be determined through in-situ tests,and parameters determination by plate loading test and pressuremeter test were introduced.Therefore,it is easy to put this model into practice.Finally,LSRM and the variable-modulus elasto-plastic model were used to analyze Egongdai ancient landslide.Safety factor,deformation field,and optimal reinforcement measures for Egongdai ancient landslide were obtained based on the proposed method.展开更多
The non-linear constitutive model suggested by the authors and the Alonso's elasto-plasticity model of unsaturated soil modified by the authors are introduced into the consolidation theory of unsaturated soil prop...The non-linear constitutive model suggested by the authors and the Alonso's elasto-plasticity model of unsaturated soil modified by the authors are introduced into the consolidation theory of unsaturated soil proposed by CHEN Zheng-han, and the non-linear and the elasto-plasticity consolidation models of unsaturated soil are obtained. Programs related to the two consolidation models are designed, and a 2-D consolidation problem of unsaturated sail is solved using the programs, the consolidation process and the development of plastic;one under multi-grade bad are studied. The above research develops the consolidation theory of unsaturated soil to a new level.展开更多
Owing to the absence of proper analytical solution of cantilever beams for couple stress/strain gradient elasto-plastic theory, experimental studies of the cantilever beam in the micro-scale are not suitable for the d...Owing to the absence of proper analytical solution of cantilever beams for couple stress/strain gradient elasto-plastic theory, experimental studies of the cantilever beam in the micro-scale are not suitable for the determination of material length-scale. Based on the couple stress elasto-plasticity, an analytical solution of thin cantilever beams is firstly presented, and the solution can be regarded as an extension of the elastic and rigid-plastic solutions of pure bending beam. A comparison with numerical results shows that the current analytical solution is reliable for the case of σ0 〈〈 H 〈〈 E, where σ0 is the initial yield strength, H is the hardening modulus and E is the elastic modulus. Fortunately, the above mentioned condition can be satisfied for many metal materials, and thus the solution can be used to determine the material length-scale of micro-structures in conjunction with the experiment of cantilever beams in the micro-scale.展开更多
A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed ...A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis.展开更多
Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed...Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.展开更多
Determination of collapse load-carrying capacity of elasto-plastic material is very important in designing structure. The problem is commonly solved by elasto-plastic finite element method (FEM). In order to deal wi...Determination of collapse load-carrying capacity of elasto-plastic material is very important in designing structure. The problem is commonly solved by elasto-plastic finite element method (FEM). In order to deal with material nonlinear problem involving strain softening problem effectively, a new numerical method-damped Newton method was proposed. The iterative schemes are discussed in detail for pure equilibrium models. In the equilibrium model, the plasticity criterion and the compatibility of the strains are verified, and the strain increment and plastic factor are treated as independent unknowns. To avoid the stiffness matrix being singularity or condition of matrix being ill, a damping factor a was introduced to adjust the value of plastic consistent parameter automatically during the iterations. According to the algorithm, the nonlinear finite element program was complied and its numerical example was calculated. The numerical results indicate that this method converges very fast for both small load steps and large load steps. Compared with those results obtained by analysis and experiment, the predicted ultimate bearing capacity from the proposed method is identical.展开更多
An in-depth analysis of propagation characteristics ofelasto-plastic combined stress waves in circular thin-walled tubeshas been made. In obtaining the simple-wave solution, however, mostresearches have ignored the in...An in-depth analysis of propagation characteristics ofelasto-plastic combined stress waves in circular thin-walled tubeshas been made. In obtaining the simple-wave solution, however, mostresearches have ignored the influence of the circumferential stressrelated to the radial inertial ef- fect in the tubes. In this paperthe incremental elasto-plastic constitutive relations which areconve- nient for dynamic numerical analysis are adopted, and thefinite-difference method is used to study the evolution adpropagation of elasto-plastic combined stress waves in a thin-walledtube with the radial inertial effect of the tube considered. Thecalculation results are compared with those obtained when the radialinertial effect is not considered. The calculation results show thatthe radial inertial effect of a tube has a fairly great influence onthe propagation of elasto-plastic combined stress waves.展开更多
Application research of neural networks to geotechnical engineering has become a hotspot nowadays.General model may not reach the predicting precision in practical application due to different characteristics in diffe...Application research of neural networks to geotechnical engineering has become a hotspot nowadays.General model may not reach the predicting precision in practical application due to different characteristics in different fields.In allusion to this,an elasto-plastic constitutive model based on clustering radial basis function neural network(BC-RBFNN) was proposed for moderate sandy clay according to its properties.Firstly,knowledge base was established on triaxial compression testing data;then the model was trained,learned and emulated using knowledge base;finally,predicting results of the BC-RBFNN model were compared and analyzed with those of other intelligent model.The results show that the BC-RBFNN model can alter the training and learning velocity and improve the predicting precision,which provides possibility for engineering practice on demanding high precision.展开更多
According to the Mohr-Coulomb yield criterion, the stress field of the infinite slope is derived under a vertical uniform load q on the top of the slope. It is indicated that elastic and elasto-plastic states would oc...According to the Mohr-Coulomb yield criterion, the stress field of the infinite slope is derived under a vertical uniform load q on the top of the slope. It is indicated that elastic and elasto-plastic states would occur in the slope. When q is smaller than the critical load, q(p), the slope is in the elastic state. If q equals q(p), the slope is in the critical state, and the plastic deformation would occur along the critical angle. With the increase of q, the plastic zone would extend, and the slope is in the elasto-plastic State. If q equals limit load, the slope is in the limit equilibrium state. The slope may be divided into three zones. Some charts of the critical angle, the critical and limit load coefficients are presented in this paper.展开更多
Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation(ANCF). The int...Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation(ANCF). The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation.A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.展开更多
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 ...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.展开更多
基金Projects (90716005, 10802055, 10972153) supported by the National Natural Science Foundation of ChinaProject (2007021005) supported by the Natural Science Foundation of Shanxi Province, China+2 种基金Project supported by the Postdoctoral Science Foundation of ChinaProject supported by the Homecomings Foundation, ChinaProject supported by the Top Young Academic Leaders of Higher Learning Institutions of Shanxi, China
文摘A multi-parameter nonlinear elasto-plastic constitutive model which can fully capture the three typical features of stress-strain response, linearity, plasticity-like stress plateau and densification phases was developed. The functional expression of each parameter was determined using uniaxial compression tests for aluminum alloy foams. The parameters of the model can be systematically varied to describe the effect of relative density which may be responsible for the changes in yield stress and hardening-like or softening-like behavior at various strain rates. A comparison between model predictions and experimental results of the aluminum alloy foams was provided to validate the model. It was proved to be useful in the selection of the optimal-density and energy absorption foam for a specific application at impact events.
基金Project supported by the National Natural Sciences Foundation of China (Nos. 10372084 and 10572119)the Program for New Century Excellent Talents in University (No. NCET-04-0958)the Open Foundation of State Key Laboratory of Structural Analysis of Industrial Equipment and the Doctorate Foundation of Northwestern Polytechnical University.
文摘In this paper, a method for the design optimization of elasto-plastic truss structures is proposed based on parametric variational principles (PVPs). The optimization aims to find the minimum weight/volume solution under the constraints of allowable node displacements. The design optimization is a formulation of mathematical programming with equilibrium constraints (MPECs). To overcome the numerical difficulties of the complementary constraints in optimization, an iteration process, comprising a quadratic programming (QP) and an updating process, is employed as the optimization method. Furthermore, the elasto-plastic buckling of truss mem- bers is considered as a constraint in design optimization. A combinational optimization strategy is proposed for the displacement constraints and the buckling constraint, which comprises the method mentioned above and an optimal criterion. Three numerical examples are presented to show the validity of the methods proposed.
基金Project supported by the National Natural Science Foundation of China (No.10572049)
文摘Based on the elasto-plastic mechanics and continuum damage theory, a yield criterion related to spherical tensor of stress is proposed to describe the mixed hardening of damaged orthotropic materials. Its dimensionless form is isomorphic with the Mises criterion for isotropic materials. Furthermore, the incremental elasto-plastic damage constitutive equations and damage evolution equations are established. Based on the classical nonlinear plate theory, the incremental nonlinear equilibrium equations of orthotropic thin plates considering damage effect are obtained, and solved with the finite difference and iteration methods. In the numerical examples, the effects of damage evolution and initial deflection on the elasto-plastic postbuckling of orthotropic plates are discussed in detail.
基金Projects(50821003,50405014)supported by the National Natural Science Foundation of ChinaProjects(10QH1401400,10520705000,10JC1407300)supported by Shanghai Committee of Science and Technology,China+1 种基金Project(NCET-07-0545)supported by Program for New Century Excellent Talents in University,ChinaFord University Research Program,China
文摘To gain a better understanding about texture evolution during rolling process of AZ31 alloy, polycrystalline plasticity model was implemented into the explicit FE package, ABAQUS/Explicit by writing a user subroutine VUMAT. For each individual grain in the polycrystalline aggregate, the rate dependent model was adopted to calculate the plastic shear strain increment in combination with the Voce hardening law to describe the hardening response, the lattice reorientation caused by slip and twinning were calculated separately due to their different mechanisms. The elasto-plastic self consistent (EPSC) model was employed to relate the response of individual grain to the response of the polycrystalline aggregate. Rolling processes of AZ31 sheet and as-cast AZ31 alloy were simulated respectively. The predicted texture distributions are in aualitative a^reement with experimental results.
文摘The dynamic buckling of elasto-plastic cylindrical shells under axial fluid-solid impact is investigated theoretically. A simplified liquid- gas- structure model is given to approximately imitate the problem. The basic equation of the structure is derived from a minimum principle in dynamics of elasto-plastic continua at finite deformation, and the flow theory of plasticity is employed. The liquid is incompressible and the gas is compressed adiabatically. A number of numerical results are presented and the characteristics of the buckling behavior under fluid-solid impact are illustrated.
文摘The present paper aims to establish a versatile strength theory suitable for elasto-plastic analysis of underground tunnel surrounding rock. In order to analyze the effects of intermediate principal stress and the rock properties on its deformation and failure of rock mass, the generalized nonlinear unified strength theory and elasto-plastic mechanics are used to deduce analytic solution of the radius and stress of tunnel plastic zone and the periphery displacement of tunnel under uniform ground stress field. The results show that: intermediate principal stress coefficient b has significant effect on the plastic range,the magnitude of stress and surrounding rock pressure. Then, the results are compared with the unified strength criterion solution and Mohr–Coulomb criterion solution, and concluded that the generalized nonlinear unified strength criterion is more applicable to elasto-plastic analysis of underground tunnel surrounding rock.
基金Project 50639100 supported by the National Natural Science Foundation of China
文摘The relationship between the Hoek-Brown parameters and the mechanical response of circular tunnels is il-lustrated. Closed-form and approximate solutions are given for the extent of the plastic zone and the stress and dis-placement fields under axisymmetrical and asymmetric stress conditions. For the same rock masses and under axisym-metrical stress conditions,the radius of the plastic zone in terms of Hoek-Brown criterion is generally an approximation of the radius in terms of the Mohr-Coulomb criterion. The radius in terms of the Hoek-Brown criterion is larger under low stress conditions. For poor quality rock masses (GSI<25),measures (such as grouting,setting rock bolts,etc.) that improve the GSI of rock masses are effective in improving the stability of tunnels. It is not advisable to improve the sta-bility of the tunnels by providing a small support resistance p through shotcrete,except for very poor quality jointed rock masses. Without reference to the quality of the rock mass,the disturbance factor D should not less than 0.5. Meas-ures which disturb rock masses during tunnel construction should be taken carefully when the tunnel depth increases.
文摘In this work, the finite element analysis of the elasto-plastic plate bending problems is carried out using transition rectangular plate elements. The shape functions of the transition plate elements are derived based on a practical rule. The transition plate elements are all quadrilateral and can be used to obtain efficient finite element models using minimum number of elements. The mesh convergence rates of the models including the transition elements are compared with the regular element models. To verify the developed elements, simple tests are demonstrated and various elasto-plastic problems are solved. Their results are compared with ANSYS results.
基金The project supported by the National Natural Science Foundation of China(19832010,50278012,10272027)the National Key Basic Research and Development Program(973 Program,2002CB412709)
文摘A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed.As compared with the existing discontinuous Galerkin finite element methods,the distinct feature of the proposed method is that the continuity of the displacement vector at each discrete time instant is automatically ensured,whereas the discontinuity of the velocity vector at the discrete time levels still remains.The computational cost is then obviously reduced, particularly,for material non-linear problems.Both the implicit and explicit algorithms to solve the derived formulations for material non-linear problems are developed.Numerical results show a good performance of the present method in eliminating spurious numerical oscillations and providing with much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain.
基金Project([2005]205)supported by the Science and Technology Planning Project of Water Resources Department of Guangdong Province,ChinaProject(2012-7)supported by Guangdong Bureau of Highway Administration,ChinaProject(2012210020203)supported by the Fundamental Research Funds for the Central Universities,China
文摘Employing an ideal elasto-plastic model,the typically used strength reduction method reduced the strength of all soil elements of a slope.Therefore,this method was called the global strength reduction method(GSRM).However,the deformation field obtained by GSRM could not reflect the real deformation of a slope when the slope became unstable.For most slopes,failure occurs once the strength of some regional soil is sufficiently weakened; thus,the local strength reduction method(LSRM)was proposed to analyze slope stability.In contrast with GSRM,LSRM only reduces the strength of local soil,while the strength of other soil remains unchanged.Therefore,deformation by LSRM is more reasonable than that by GSRM.In addition,the accuracy of the slope's deformation depends on the constitutive model to a large degree,and the variable-modulus elasto-plastic model was thus adopted.This constitutive model was an improvement of the Duncan–Chang model,which modified soil's deformation modulus according to stress level,and it thus better reflected the plastic feature of soil.Most importantly,the parameters of the variable-modulus elasto-plastic model could be determined through in-situ tests,and parameters determination by plate loading test and pressuremeter test were introduced.Therefore,it is easy to put this model into practice.Finally,LSRM and the variable-modulus elasto-plastic model were used to analyze Egongdai ancient landslide.Safety factor,deformation field,and optimal reinforcement measures for Egongdai ancient landslide were obtained based on the proposed method.
文摘The non-linear constitutive model suggested by the authors and the Alonso's elasto-plasticity model of unsaturated soil modified by the authors are introduced into the consolidation theory of unsaturated soil proposed by CHEN Zheng-han, and the non-linear and the elasto-plasticity consolidation models of unsaturated soil are obtained. Programs related to the two consolidation models are designed, and a 2-D consolidation problem of unsaturated sail is solved using the programs, the consolidation process and the development of plastic;one under multi-grade bad are studied. The above research develops the consolidation theory of unsaturated soil to a new level.
基金the National Natural Science Foundation of China (50479058, 10672032)
文摘Owing to the absence of proper analytical solution of cantilever beams for couple stress/strain gradient elasto-plastic theory, experimental studies of the cantilever beam in the micro-scale are not suitable for the determination of material length-scale. Based on the couple stress elasto-plasticity, an analytical solution of thin cantilever beams is firstly presented, and the solution can be regarded as an extension of the elastic and rigid-plastic solutions of pure bending beam. A comparison with numerical results shows that the current analytical solution is reliable for the case of σ0 〈〈 H 〈〈 E, where σ0 is the initial yield strength, H is the hardening modulus and E is the elastic modulus. Fortunately, the above mentioned condition can be satisfied for many metal materials, and thus the solution can be used to determine the material length-scale of micro-structures in conjunction with the experiment of cantilever beams in the micro-scale.
基金Project supported by the National Natural Science Foundation of China(No.10172052).
文摘A meshless approach based on the moving least square method is developed for elasto-plasticity analysis,in which the incremental formulation is used.In this approach,the dis- placement shape functions are constructed by using the moving least square approximation,and the discrete governing equations for elasto-plastic material are constructed with the direct collo- cation method.The boundary conditions are also imposed by collocation.The method established is a truly meshless one,as it does not need any mesh,either for the purpose of interpolation of the solution variables,or for the purpose of construction of the discrete equations.It is simply formu- lated and very efficient,and no post-processing procedure is required to compute the derivatives of the unknown variables,since the solution from this method based on the moving least square approximation is already smooth enough.Numerical examples are given to verify the accuracy of the meshless method proposed for elasto-plasticity analysis.
基金The Project Supported by National Natural Science Foundation of China
文摘Several effective numerical methods for solving the elasto-plastic contact problems with friction are pres- ented.First,a direct substitution method is employed to impose the contact constraint conditions on condensed finite ele- ment equations,thus resulting in a reduction by half in the dimension of final governing equations.Second,an algorithm composed of contact condition probes and elasto-plastic iterations is utilized to solve the governing equation,which distinguishes two kinds of nonlinearities,and makes the solution unique.In addition,Positive-Negative Sequence Modifica- tion Method is used to condense the finite element equations of each substructure and an analytical integration is intro- duced to determine the elasto-plastic status after each time step or each iteration,hence the computational efficiency is en- hanced to a great extent.Finally,several test and practical examples are pressented showing the validity and versatility of these methods and algorithms.
基金Project(2012CB026200)supported by the National Basic Research Program of ChinaProjects(50978055,50878048)supported by the National Natural Science Foundation of China
文摘Determination of collapse load-carrying capacity of elasto-plastic material is very important in designing structure. The problem is commonly solved by elasto-plastic finite element method (FEM). In order to deal with material nonlinear problem involving strain softening problem effectively, a new numerical method-damped Newton method was proposed. The iterative schemes are discussed in detail for pure equilibrium models. In the equilibrium model, the plasticity criterion and the compatibility of the strains are verified, and the strain increment and plastic factor are treated as independent unknowns. To avoid the stiffness matrix being singularity or condition of matrix being ill, a damping factor a was introduced to adjust the value of plastic consistent parameter automatically during the iterations. According to the algorithm, the nonlinear finite element program was complied and its numerical example was calculated. The numerical results indicate that this method converges very fast for both small load steps and large load steps. Compared with those results obtained by analysis and experiment, the predicted ultimate bearing capacity from the proposed method is identical.
文摘An in-depth analysis of propagation characteristics ofelasto-plastic combined stress waves in circular thin-walled tubeshas been made. In obtaining the simple-wave solution, however, mostresearches have ignored the influence of the circumferential stressrelated to the radial inertial ef- fect in the tubes. In this paperthe incremental elasto-plastic constitutive relations which areconve- nient for dynamic numerical analysis are adopted, and thefinite-difference method is used to study the evolution adpropagation of elasto-plastic combined stress waves in a thin-walledtube with the radial inertial effect of the tube considered. Thecalculation results are compared with those obtained when the radialinertial effect is not considered. The calculation results show thatthe radial inertial effect of a tube has a fairly great influence onthe propagation of elasto-plastic combined stress waves.
基金Project(07031B) supported by the Scientific Research Fund of Central South University of Forestry and TechnologyProject(06C843) supported by the Scientific Research Fund of Hunan Provincial Education Department
文摘Application research of neural networks to geotechnical engineering has become a hotspot nowadays.General model may not reach the predicting precision in practical application due to different characteristics in different fields.In allusion to this,an elasto-plastic constitutive model based on clustering radial basis function neural network(BC-RBFNN) was proposed for moderate sandy clay according to its properties.Firstly,knowledge base was established on triaxial compression testing data;then the model was trained,learned and emulated using knowledge base;finally,predicting results of the BC-RBFNN model were compared and analyzed with those of other intelligent model.The results show that the BC-RBFNN model can alter the training and learning velocity and improve the predicting precision,which provides possibility for engineering practice on demanding high precision.
文摘According to the Mohr-Coulomb yield criterion, the stress field of the infinite slope is derived under a vertical uniform load q on the top of the slope. It is indicated that elastic and elasto-plastic states would occur in the slope. When q is smaller than the critical load, q(p), the slope is in the elastic state. If q equals q(p), the slope is in the critical state, and the plastic deformation would occur along the critical angle. With the increase of q, the plastic zone would extend, and the slope is in the elasto-plastic State. If q equals limit load, the slope is in the limit equilibrium state. The slope may be divided into three zones. Some charts of the critical angle, the critical and limit load coefficients are presented in this paper.
基金supported in part by the National Natural Science Foundation of China (Grants 11290151 and 11221202)supported in part by the Beijing Higher Education Young Elite Teacher Project (Grant YETP1201)
文摘Under the frame of multibody dynamics, the contact dynamics of elasto-plastic spatial thin beams is numerically studied by using the spatial thin beam elements of absolute nodal coordinate formulation(ANCF). The internal force of the elasto-plastic spatial thin beam element is derived under the assumption that the plastic strain of the beam element depends only on its longitudinal deformation.A new body-fixed local coordinate system is introduced into the spatial thin beam element of ANCF for efficient contact detection in the contact dynamics simulation. The linear isotropic hardening constitutive law is used to describe the elasto-plastic deformation of beam material, and the classical return mapping algorithm is adopted to evaluate the plastic strains. A multi-zone contact approach of thin beams previously proposed by the authors is also introduced to detect the multiple contact zones of beams accurately, and the penalty method is used to compute the normal contact force of thin beams in contact. Four numerical examples are given to demonstrate the applicability and effectiveness of the proposed elasto-plastic spatial thin beam element of ANCF for flexible multibody system dynamics.
基金financially supported by the National Natural Science Foundation of China(Grant No.52074258,Grant No.41941018,Grant No.51974289,and Grant No.51874232)the Natural Science Basic Research Program of Shaanxi Province(Shaanxi Coal and Chemical Industry Group Co.,Ltd.Joint Fund Project,Grant No.2021JLM-06)the open project of State Key Laboratory of Shield Machine and Boring Technology(Grant No.E01Z440101)。
文摘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.