In the process of shield tunneling through soft soil layers,the presence of confined water ahead poses a significant threat to the stability of the tunnel face.Therefore,it is crucial to consider the impact of confine...In the process of shield tunneling through soft soil layers,the presence of confined water ahead poses a significant threat to the stability of the tunnel face.Therefore,it is crucial to consider the impact of confined water on the limit support pressure of the tunnel face.This study employed the finite element method(FEM)to analyze the limit support pressure of shield tunnel face instability within a pressurized water-containing layer.Subsequently,a multiple linear regression approach was applied to derive a concise solution formula for the limit support pressure,incorporating various influencing factors.The analysis yields the following conclusions:1)The influence of confined water on the instability mode of the tunnel face in soft soil layers makes the displacement response of the strata not significant when the face is unstable;2)The limit support pressure increases approximately linearly with the pressure head,shield tunnel diameter,and tunnel burial depth.And inversely proportional to the thickness of the impermeable layer,soil cohesion and internal friction angle;3)Through an engineering case study analysis,the results align well with those obtained from traditional theoretical methods,thereby validating the rationality of the equations proposed in this paper.Furthermore,the proposed equations overcome the limitation of traditional theoretical approaches considering the influence of changes in impermeable layer thickness.It can accurately depict the dynamic variation in the required limit support pressure to maintain the stability of the tunnel face during shield tunneling,thus better reflecting engineering reality.展开更多
Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton...Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation(capillary pressure).So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.展开更多
The high temperature split Hopkinson pressure bar (SHPB) compression experiment is conducted to obtain the data relationship among strain, strain rate and flow stress from room temperature to 550 C for aeronautical ...The high temperature split Hopkinson pressure bar (SHPB) compression experiment is conducted to obtain the data relationship among strain, strain rate and flow stress from room temperature to 550 C for aeronautical aluminum alloy 7050-T7451. Combined high-speed orthogonal cutting experiments with the cutting process simulations, the data relationship of high temperature, high strain rate and large strain in high-speed cutting is modified. The Johnson-Cook empirical model considering the effects of strain hardening, strain rate hardening and thermal softening is selected to describe the data relationship in high-speed cutting, and the material constants of flow stress constitutive model for aluminum alloy 7050-T7451 are determined. Finally, the constitutive model of aluminum alloy 7050-T7451 is established through experiment and simulation verification in high-speed cutting. The model is proved to be reasonable by matching the measured values of the cutting force with the estimated results from FEM simulations.展开更多
In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the sof...In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.展开更多
The calculation of earth pressure is a difficult problem in the pit foundation design. Aiming at the problem, the earth pressure calculation formulas considering the displacement are proposed. A method for determining...The calculation of earth pressure is a difficult problem in the pit foundation design. Aiming at the problem, the earth pressure calculation formulas considering the displacement are proposed. A method for determining the limit displacement is given and it is convenient for use. The result indicates that the earth pressure calculated by the formulas is between the earth pressure at rest and the Rankine earth pressure, the formulas can reflect the change rules of active and passive earth pressures along with the displacement. Moreover, the calculation result using the formulas is approximate to the monitoring result. It confirms the validity of the formula deduction and the rationality of the calculation results. As for the passive earth pressure, the calculation method is theoretically feasible.展开更多
Metals heat-treated under high pressure can exhibit different properties. The heat-induced pressure on 2024 aluminum alloy during restricting expansion-deformation heat-treatment was calculated by using the ABAQUS fin...Metals heat-treated under high pressure can exhibit different properties. The heat-induced pressure on 2024 aluminum alloy during restricting expansion-deformation heat-treatment was calculated by using the ABAQUS finite element software, and the effects of the mould material properties, such as coefficient of thermal expansion (CTE), elastic modulus and yield strength, on the pressure were discussed. The simulated results show that the relatively uniform heat-induced pressure, approximately 503 MPa at 500 ℃, appears on 2024 alloy when 42CrMo steel is as the mould material. The heat-induced pressure increases with decreasing the CTE and the increases of elastic modulus and yield strength of the mould material. The influences of the CTE and elastic modulus on the heat-induced pressure are more notable.展开更多
This paper presents an analytical scheme for predicting the collapse strength of a flexible pipe, which considers the structural interaction between relevant layers. The analytical results were compared with a FEA mod...This paper presents an analytical scheme for predicting the collapse strength of a flexible pipe, which considers the structural interaction between relevant layers. The analytical results were compared with a FEA model and a number of test data, and showed reasonably good agreement. The theoretical analysis showed that the pressure armor layer enhanced the strength of the carcass against buckling, though the barrier weakened this effect. The collapse strength of pipe was influenced by many factors such as the inner radius of the pipe, the thickness of the layers and the mechanical properties of the materials. For example, an increase in the thickness of the barrier will increase contact pressure and in turn reduce the critical pressure.展开更多
This paper proposes a theoretical method using finite element analysis(FEA) to calculate the plastic collapse loads of pressure vessels under internal pressure,and compares the analytical methods according to three cr...This paper proposes a theoretical method using finite element analysis(FEA) to calculate the plastic collapse loads of pressure vessels under internal pressure,and compares the analytical methods according to three criteria stated in the ASME Boiler Pressure Vessel Code. First,a finite element technique using the arc-length algorithm and the restart analysis is developed to conduct the plastic collapse analysis of vessels,which includes the material and geometry non-linear properties of vessels. Second,as the mechanical properties of vessels are assumed to be elastic-perfectly plastic,the limit load analysis is performed by em-ploying the Newton-Raphson algorithm,while the limit pressure of vessels is obtained by the twice-elastic-slope method and the tangent intersection method respectively to avoid excessive deformation. Finally,the elastic stress analysis under working pressure is conducted and the stress strength of vessels is checked by sorting the stress results. The results are compared with those obtained by experiments and other existing models. This work provides a reference for the selection of the failure criteria and the calculation of the plastic collapse load.展开更多
An integrated approach was proposed to evaluate the remaining useful life(RUL)of corroded petroleum pipelines.Two types of failure modes(i.e.,leakage and burst failure)were considered,and the corresponding limit state...An integrated approach was proposed to evaluate the remaining useful life(RUL)of corroded petroleum pipelines.Two types of failure modes(i.e.,leakage and burst failure)were considered,and the corresponding limit state functions(LSFs)were established with the structural reliability theory.A power-law function was applied to model the growth of corrosion defects,and the effect of external environmental factors on the growth of the pipeline s defect was considered.Moreover,the result was compared with the commonly used linear growth model.After that,a finite element simulation model was established to calculate the burst pressure of the pipeline with corrosion defects,and its accuracy was verified through hydraulic burst test and by comparison with international criteria.On that basis,the probability that the pipeline may fail was calculated with Monte Carlo simulation(MCS)and by considering the LSFs,and the pipeline s RUL was obtained accordingly.Furthermore,sensitivity analysis was conducted to determine the sensitivity parameters for the corrosion and RUL of the pipeline.The results indicate that the radial corrosion rate,wall thickness and working pressure have a great influence on the failure probability of the pipeline.Thus,corresponding measures should be adopted during the operation process of the pipeline to reduce the corrosion rate and increase the wall thickness,so as to prolong the pipeline s RUL.展开更多
The fracture initiation behavior for hydraulic fracturing treatments highlighted the necessity of proposing fracture criteria that precisely predict the fracture initiation type and location during the hydraulic fract...The fracture initiation behavior for hydraulic fracturing treatments highlighted the necessity of proposing fracture criteria that precisely predict the fracture initiation type and location during the hydraulic fracturing process.In the present study,a Mohr-Coulomb criterion with a tensile cut-off is incorporated into the finite element code to determine the fracture initiation type and location during the hydraulic fracturing process.This fracture criterion considers the effect of fracture inclination angle,the internal friction angle and the loading conditions on the distribution of stress field around the fracture tip.The results indicate that the internal friction angle resists the shear fracture initiation.Moreover,as the internal friction angle increases,greater external loads are required to maintain the hydraulic fracture extension.Due to the increased pressure of the injected water,the tensile fracture ultimately determines the fracture initiation type.However,the shear fracture preferentially occurs as the stress anisotropy coefficient increases.Both the maximum tensile stress and equivalent maximum shear stress decrease as the stress anisotropy coefficient increases,which indicates that the greater the stress anisotropy coefficient,the higher the external loading required to propagate a new fracture.The numerical results obtained in this paper provide theoretical supports for establishing basis on investigating of the hydraulic fracturing characteristics under different conditions.展开更多
Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the invest...Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.展开更多
Based on the construction bridge of Xiamen-Shenzhen high-speed railway(9-32 m simply-supported beam + 6×32 m continuous beam),the pier-beam-track finite element model,where the continuous beam of the ballast trac...Based on the construction bridge of Xiamen-Shenzhen high-speed railway(9-32 m simply-supported beam + 6×32 m continuous beam),the pier-beam-track finite element model,where the continuous beam of the ballast track and simply-supported beam are combined with each other,was established.The laws of the track stress,the pier longitudinal stress and the beam-track relative displacement were analyzed.The results show that reducing the longitudinal resistance can effectively reduce the track stress and the pier stress of the continuous beam,and increase the beam-track relative displacement.Increasing the rigid pier stiffness of continuous beam can reduce the track braking stress,increase the pier longitudinal stress and reduce the beam-track relative displacement,Increasing the rigid pier stiffness of simply-supported beam can reduce the track braking stress,the rigid pier longitudinal stress and the beam-track relative displacement.展开更多
Isotropic consolidation test and consolidated-undrained triaxial test were first undertaken to obtain the parameters of the modified cam-clay(MCC)model and the behavior of natural clayey soil.Then,for the first time,n...Isotropic consolidation test and consolidated-undrained triaxial test were first undertaken to obtain the parameters of the modified cam-clay(MCC)model and the behavior of natural clayey soil.Then,for the first time,numerical simulation of the two tests was performed by three-dimensional finite element method(FEM)using ABAQUS program.The consolidated-drained triaxial test was also simulated by FEM and compared with theoretical results of MCC model.Especially,the behaviors of MCC model during unloading and reloading were analyzed in detail by FEM.The analysis and comparison indicate that the MCC model is able to accurately describe many features of the mechanical behavior of the soil in isotropic consolidation test and consolidated-drained triaxial test.And the MCC model can well describe the variation of excess pore water pressure with the development of axial strain in consolidated-undrained triaxial test,but its ability to predict the relationship between axial strain and shear stress is relatively poor.The comparison also shows that FEM solutions of the MCC model are basically identical to the theoretical ones.In addition,Mandel-Cryer effect unable to be discovered by the conventional triaxial test in laboratories was disclosed by FEM.The analysis of unloading-reloading by FEM demonstrates that the MCC model disobeys the law of energy conservation under the cyclic loading condition if the elastic shear modulus is linearly pressure-dependent.展开更多
基金Project(ZDRW-ZS-2021-3)supported by the Key Deployment Projects of Chinese Academy of SciencesProjects(52179116,51991392)supported by the National Natural Science Foundation of China。
文摘In the process of shield tunneling through soft soil layers,the presence of confined water ahead poses a significant threat to the stability of the tunnel face.Therefore,it is crucial to consider the impact of confined water on the limit support pressure of the tunnel face.This study employed the finite element method(FEM)to analyze the limit support pressure of shield tunnel face instability within a pressurized water-containing layer.Subsequently,a multiple linear regression approach was applied to derive a concise solution formula for the limit support pressure,incorporating various influencing factors.The analysis yields the following conclusions:1)The influence of confined water on the instability mode of the tunnel face in soft soil layers makes the displacement response of the strata not significant when the face is unstable;2)The limit support pressure increases approximately linearly with the pressure head,shield tunnel diameter,and tunnel burial depth.And inversely proportional to the thickness of the impermeable layer,soil cohesion and internal friction angle;3)Through an engineering case study analysis,the results align well with those obtained from traditional theoretical methods,thereby validating the rationality of the equations proposed in this paper.Furthermore,the proposed equations overcome the limitation of traditional theoretical approaches considering the influence of changes in impermeable layer thickness.It can accurately depict the dynamic variation in the required limit support pressure to maintain the stability of the tunnel face during shield tunneling,thus better reflecting engineering reality.
基金supported by the 973 Program (Grant No.2007CB209505)the National Natural Science Foundation of China (Grant No.40674061,40704019)
文摘Wavefields in porous media saturated by two immiscible fluids are simulated in this paper.Based on the sealed system theory,the medium model considers both the relative motion between the fluids and the solid skeleton and the relaxation mechanisms of porosity and saturation(capillary pressure).So it accurately simulates the numerical attenuation property of the wavefields and is much closer to actual earth media in exploration than the equivalent liquid model and the unsaturated porous medium model on the basis of open system theory.The velocity and attenuation for different wave modes in this medium have been discussed in previous literature but studies of the complete wave-field have not been reported.In our work,wave equations with the relaxation mechanisms of capillary pressure and the porosity are derived.Furthermore,the wavefield and its characteristics are studied using the numerical finite element method.The results show that the slow P3-wave in the non-wetting phase can be observed clearly in the seismic band.The relaxation of capillary pressure and the porosity greatly affect the displacement of the non-wetting phase.More specifically,the displacement decreases with increasing relaxation coefficient.
文摘The high temperature split Hopkinson pressure bar (SHPB) compression experiment is conducted to obtain the data relationship among strain, strain rate and flow stress from room temperature to 550 C for aeronautical aluminum alloy 7050-T7451. Combined high-speed orthogonal cutting experiments with the cutting process simulations, the data relationship of high temperature, high strain rate and large strain in high-speed cutting is modified. The Johnson-Cook empirical model considering the effects of strain hardening, strain rate hardening and thermal softening is selected to describe the data relationship in high-speed cutting, and the material constants of flow stress constitutive model for aluminum alloy 7050-T7451 are determined. Finally, the constitutive model of aluminum alloy 7050-T7451 is established through experiment and simulation verification in high-speed cutting. The model is proved to be reasonable by matching the measured values of the cutting force with the estimated results from FEM simulations.
文摘In order to study the laws of the extrusion pressure changing with the extrusion parameters in the process of hydrostatic extrusion for the tungsten alloys, the large deformation elasto plastic theory and the software of ANSYS 5 5 are used to carry out the numerical simulation research. The laws of the extrusion pressure changing with the extrusion parameters, such as the die angle, extrusion ratio, and friction coefficient, are obtained. The simulation results are in good agreement with the experimental ones, and the simulated results are believable.
基金Supported by the Natural Foundation of Jiangsu Province(BK2004427)~~
文摘The calculation of earth pressure is a difficult problem in the pit foundation design. Aiming at the problem, the earth pressure calculation formulas considering the displacement are proposed. A method for determining the limit displacement is given and it is convenient for use. The result indicates that the earth pressure calculated by the formulas is between the earth pressure at rest and the Rankine earth pressure, the formulas can reflect the change rules of active and passive earth pressures along with the displacement. Moreover, the calculation result using the formulas is approximate to the monitoring result. It confirms the validity of the formula deduction and the rationality of the calculation results. As for the passive earth pressure, the calculation method is theoretically feasible.
基金Project (51071125) supported by the National Natural Science Foundation of ChinaProjects (SKLSP201107, SKLSP 201124) supported by the Fund of the State Key Laboratory of Solidification Processing in Northwestern Polytechnical University, China
文摘Metals heat-treated under high pressure can exhibit different properties. The heat-induced pressure on 2024 aluminum alloy during restricting expansion-deformation heat-treatment was calculated by using the ABAQUS finite element software, and the effects of the mould material properties, such as coefficient of thermal expansion (CTE), elastic modulus and yield strength, on the pressure were discussed. The simulated results show that the relatively uniform heat-induced pressure, approximately 503 MPa at 500 ℃, appears on 2024 alloy when 42CrMo steel is as the mould material. The heat-induced pressure increases with decreasing the CTE and the increases of elastic modulus and yield strength of the mould material. The influences of the CTE and elastic modulus on the heat-induced pressure are more notable.
文摘This paper presents an analytical scheme for predicting the collapse strength of a flexible pipe, which considers the structural interaction between relevant layers. The analytical results were compared with a FEA model and a number of test data, and showed reasonably good agreement. The theoretical analysis showed that the pressure armor layer enhanced the strength of the carcass against buckling, though the barrier weakened this effect. The collapse strength of pipe was influenced by many factors such as the inner radius of the pipe, the thickness of the layers and the mechanical properties of the materials. For example, an increase in the thickness of the barrier will increase contact pressure and in turn reduce the critical pressure.
基金Project (Nos. 2006BAK04A02-02 and 2006BAK02B02-08) supported by the National Key Technology R&D Program, China
文摘This paper proposes a theoretical method using finite element analysis(FEA) to calculate the plastic collapse loads of pressure vessels under internal pressure,and compares the analytical methods according to three criteria stated in the ASME Boiler Pressure Vessel Code. First,a finite element technique using the arc-length algorithm and the restart analysis is developed to conduct the plastic collapse analysis of vessels,which includes the material and geometry non-linear properties of vessels. Second,as the mechanical properties of vessels are assumed to be elastic-perfectly plastic,the limit load analysis is performed by em-ploying the Newton-Raphson algorithm,while the limit pressure of vessels is obtained by the twice-elastic-slope method and the tangent intersection method respectively to avoid excessive deformation. Finally,the elastic stress analysis under working pressure is conducted and the stress strength of vessels is checked by sorting the stress results. The results are compared with those obtained by experiments and other existing models. This work provides a reference for the selection of the failure criteria and the calculation of the plastic collapse load.
基金The National Natural Science Foundation of China(No.71671035,72001039)the Open Fund of Hunan Provincial Key Laboratory of Health Maintenance for Mechanical Equipment(No.201901)the Open Fund of Jiangsu Wind Power Engineering Technology Center(No.ZK19-03-03)。
文摘An integrated approach was proposed to evaluate the remaining useful life(RUL)of corroded petroleum pipelines.Two types of failure modes(i.e.,leakage and burst failure)were considered,and the corresponding limit state functions(LSFs)were established with the structural reliability theory.A power-law function was applied to model the growth of corrosion defects,and the effect of external environmental factors on the growth of the pipeline s defect was considered.Moreover,the result was compared with the commonly used linear growth model.After that,a finite element simulation model was established to calculate the burst pressure of the pipeline with corrosion defects,and its accuracy was verified through hydraulic burst test and by comparison with international criteria.On that basis,the probability that the pipeline may fail was calculated with Monte Carlo simulation(MCS)and by considering the LSFs,and the pipeline s RUL was obtained accordingly.Furthermore,sensitivity analysis was conducted to determine the sensitivity parameters for the corrosion and RUL of the pipeline.The results indicate that the radial corrosion rate,wall thickness and working pressure have a great influence on the failure probability of the pipeline.Thus,corresponding measures should be adopted during the operation process of the pipeline to reduce the corrosion rate and increase the wall thickness,so as to prolong the pipeline s RUL.
基金Project(2017YFC1503102)supported by the National Key Research and Development ProgramProjects(51874065,U1903112)supported by the National Natural Science Foundation of China。
文摘The fracture initiation behavior for hydraulic fracturing treatments highlighted the necessity of proposing fracture criteria that precisely predict the fracture initiation type and location during the hydraulic fracturing process.In the present study,a Mohr-Coulomb criterion with a tensile cut-off is incorporated into the finite element code to determine the fracture initiation type and location during the hydraulic fracturing process.This fracture criterion considers the effect of fracture inclination angle,the internal friction angle and the loading conditions on the distribution of stress field around the fracture tip.The results indicate that the internal friction angle resists the shear fracture initiation.Moreover,as the internal friction angle increases,greater external loads are required to maintain the hydraulic fracture extension.Due to the increased pressure of the injected water,the tensile fracture ultimately determines the fracture initiation type.However,the shear fracture preferentially occurs as the stress anisotropy coefficient increases.Both the maximum tensile stress and equivalent maximum shear stress decrease as the stress anisotropy coefficient increases,which indicates that the greater the stress anisotropy coefficient,the higher the external loading required to propagate a new fracture.The numerical results obtained in this paper provide theoretical supports for establishing basis on investigating of the hydraulic fracturing characteristics under different conditions.
基金Project(2009ZX04004-031-04) supported by the National Science and Technology Major Project of the Ministry of Science and Technology of China
文摘Compaction process simulation and residual stress prediction of green PM compact were carried out with elasto-plastic 3D FEA based on the modified Drueker-Prager Cap model in Abaqus. The model parameters of the investigated powder Distaloy AE were determined as functions of relative density through typical mechanical property tests of powder. The model was implemented as a user subroutine USDFLD. Single sided compaction of a d20 ram^5 mm disk green compact of Distaloy AE was simulated, and the residual stress of the disk after ejection was predicted with FEA. The FEA results of the compaction process and the residual stress of the disk show good agreement with compaction experiments and X-ray diffraction measurements, which validates the model and its parameters. The results indicate that the compressive residual stresses exist mainly in a thin layer on the side surface, but the residual stresses are very small on the top and bottom surfaces.
基金Project(50678176) supported by the National Natural Science Foundation of China
文摘Based on the construction bridge of Xiamen-Shenzhen high-speed railway(9-32 m simply-supported beam + 6×32 m continuous beam),the pier-beam-track finite element model,where the continuous beam of the ballast track and simply-supported beam are combined with each other,was established.The laws of the track stress,the pier longitudinal stress and the beam-track relative displacement were analyzed.The results show that reducing the longitudinal resistance can effectively reduce the track stress and the pier stress of the continuous beam,and increase the beam-track relative displacement.Increasing the rigid pier stiffness of continuous beam can reduce the track braking stress,increase the pier longitudinal stress and reduce the beam-track relative displacement,Increasing the rigid pier stiffness of simply-supported beam can reduce the track braking stress,the rigid pier longitudinal stress and the beam-track relative displacement.
基金Project(2011J01308) supported by the Natural Science Foundation of Fujian Province,China
文摘Isotropic consolidation test and consolidated-undrained triaxial test were first undertaken to obtain the parameters of the modified cam-clay(MCC)model and the behavior of natural clayey soil.Then,for the first time,numerical simulation of the two tests was performed by three-dimensional finite element method(FEM)using ABAQUS program.The consolidated-drained triaxial test was also simulated by FEM and compared with theoretical results of MCC model.Especially,the behaviors of MCC model during unloading and reloading were analyzed in detail by FEM.The analysis and comparison indicate that the MCC model is able to accurately describe many features of the mechanical behavior of the soil in isotropic consolidation test and consolidated-drained triaxial test.And the MCC model can well describe the variation of excess pore water pressure with the development of axial strain in consolidated-undrained triaxial test,but its ability to predict the relationship between axial strain and shear stress is relatively poor.The comparison also shows that FEM solutions of the MCC model are basically identical to the theoretical ones.In addition,Mandel-Cryer effect unable to be discovered by the conventional triaxial test in laboratories was disclosed by FEM.The analysis of unloading-reloading by FEM demonstrates that the MCC model disobeys the law of energy conservation under the cyclic loading condition if the elastic shear modulus is linearly pressure-dependent.