高度非线性弹塑性模型的半隐式算法实现

doi:10.13722/j.cnki.jrme.2017.0043

Abstract
Figure/Table
References (32)
Related Citation (15)

Download: PDF (436 KB) HTML (1 KB)
Export: BibTeX | EndNote (RIS)

Abstract

When implementing a highly nonlinear stress-strain model to solve boundary-valued problems，one of the major challenges is how to reduce the accumulative error and to maintain the effectiveness of the numerical integration. In general，the explicit algorithm tends to have higher accumulative error，and the requirement for accuracy may not be satisfied even when adopting very small increment sizes. When utilizing the Newton-CPPM implicit numerical integration algorithm，the complex Jacobian matrix for a highly nonlinear model may result in non-convergence or even become singular in the modelling process. In order to deal with these challenges，this paper proposes an improved-implicit algorithm，in which a semi-implicit method is used to determine the initial values for the implicit iterations. Using the highly nonlinear SANICLAY model for structured clay as an example，the convergence，the computational efficiency and the accuracy of four algorithms，namely the explicit algorithm，the semi-implicit algorithm，the implicit algorithm and the improved implicit algorithm，are compared via numerical simulations of single element tests. Comparing with the explicit algorithm，the semi-implicit algorithm effectively reduces error accumulation and improves the computation accuracy. The improved implicit algorithm can effectively improve convergence and avoid singularity of Jacobian matrix.

Key words： soil mechanics structured clays nonlinear elasto-plastic constitutive model semi-implicit algorithm numerical implementation

	Service

	E-mail this article
	Add to my bookshelf
	Add to citation manager
	E-mail Alert
	RSS
	Articles by authors
	GENG Dajiang1
	GUO Peijun1
	2
	ZHOU Shunhua1

Cite this article:

GENG Dajiang1,GUO Peijun1,2, et al. Semi-implicit numerical integration of a highly nonlinear elasto-plastic constitutive model[J]. , 2018, 37(S1): 3178-3186.

URL:

https://rockmech.whrsm.ac.cn/EN/10.13722/j.cnki.jrme.2017.0043 OR https://rockmech.whrsm.ac.cn/EN/Y2018/V37/IS1/3178

[1]	CALLISTO L，RAMPELLO S. Shear strength and small-strain stiffness of a natural clay under general stress conditions[J]. Geotechnique，2002，52(8)：547-560.
[1]	孙德安，陈波. 结构性软土力学特性的试验研究[J]. 土木工程学报，2011，44(增2)：65-68.(SUN Dean，CHEN Bo. Experimental study on the mechanical behavior of structural soft clay[J].China Civil Engineering Journal，2011，44(Supp.2)：65-68.(in Chinese))
[2]	JUNG YH，FINNO RJ，CHO W. Stress-strain responses of reconstituted and natural compressible Chicago glacial clay[J]. Engineering Geology，2012，129：9-19.
[3]	LIU WZ，SHI ML，MIAO LC，et al.Constitutive modeling of the destructuration and anisotropy of natural soft clay[J]. Computers and Geotechnics，2013，51：24-41.
[4]	PANAYIDES S，ROUAINIA M，WOOD DM. Influence of degradation of structure on the behaviour of a full-scale embankment[J]. Canadian Geotechnical Journal，2012，49(3)：344-356.
[1]	ROCCHI G，VACIAGO G，FONTANA M，et al. Understanding sampling disturbance and behaviour of structured clays through constitutive modelling[J]. Soils and Foundations，2013，53(2)：315-334.
[2]	TAIEBAT M，DAFALIAS YF，PEEK R. A destructuration theoryand its application to SANICLAY model[J]. International Journal for Numerical and Analytical Methods in Geomechanics，2010，34(10)：1009-1040.
[3]	SUEBSUK J，HORPIBULSUK S，LIU MD. Modified Structured Cam Clay：A generalised critical state model for destructured，naturally structured and artificially structured clays[J]. Computers and Geotechnics，2010，37：956-968.
[4]	祝恩阳，姚仰平. 结构性土UH模型[J]. 岩土力学，2015，36(11)：3101-3110.(ZHU Enyang，YAO Yangping. A UH constitutive model for structured soils[J]. Rock and Soil Mechanics，2015，36(11)：3101-3110.(in Chinese))
[5]	CALLISTO L，GAJO A，WOOD DM. Simulation of triaxial and true triaxial tests on natural and reconstituted Pisa clay[J]. Geotechnique，2002，52(9)：649-666.
[6]	LIU MD，CARTER JP，AIREY DW. Sydney soil model. I：theoretical formulation[J]. International Journal of Geomechanics，2011，11(3)：211-224.
[7]	FARIAS MM，PEDROSO DM，NAKAI T. Automatic substepping integration of the subloading TIJ model with stress path dependent hardening[J]. Computers and Geotechnics，2009，36(4)：537-548.
[8]	SOLOWSKI WT，GALLIPOLI D. Explicit stress integration with error control for the Barcelona basic model. Part I：algorithms formulations[J]. Computers and Geotechnics，2010，37：59-67.
[9]	SOLOWSKI WT，GALLIPOLI D. Explicit stress integration with error control for the Barcelona basic model. Part II：algorithms efficiencyand accuracy[J]. Computers and Geotechnics，2010，37：68-81.
[10]	CATTANEO F，DELLA VG，JOMMI C. Evaluation of numerical stress-point algorithms on elastic-plastic models for unsaturated soils with hardening dependent on the degree of saturation[J]. Computers and Geotechnics，2014，55：404-415.
[11]	SLOAN SW，ABBO AJ，SHENG DC. Refined explicit integration of elastoplastic models with automatic error control[J]. Engineering Computations，2001，18(1/2)：121-154.
[12]	CABOT ML，SLOAN SW，SHENG DC，et al. Error behaviour in explicit integration algorithms with automatic substepping[J]. International Journal for Numerical Methods in Engineering，2016，108：1030-1053.
[1]	BICANICN，PEARCECJ. Computational aspects of a softening plasticity model for plain concrete[J]. Mechanics of Cohesive-frictional Materials，1996，1(1)：75-94.
[2]	STUPKIEWICZ S，DENZER R，PICCOLROAZ A，et al. Implicit yield function formulation for granular and rock-like materials[J]. Computational Mechanics，2014，54(5)：1163-1173.
[3]	VALENTINI B，HOFSTETTER G. Review and enhancement of 3D concrete models for large-scale numerical simulations of concrete structures[J].International Journal for Numerical and Analytical Methods in Geomechanics，2013，37(3)：221-246.
[1]	JIU JL. Integration algorithm for a modified Yoshida-Uemori model to simulate cyclic plasticity in extremely large plastic strain ranges up to fracture[J].Computers and Structures，2014，145：36-46.
[2]	MANZARI MT，YONTEN K. On implementation and performance of an anisotropic constitutive model for clays[J]. International Journal of Computational Methods，2014，11(2)：1-31.
[3]	PENASA M，PICCOLROAZ A，ARGANI L，et al.Integration algorithms of elastoplasticity for ceramic powder compaction[J].Journal of the European Ceramic Society，2014，34(11)：2775-2788.
[4]	WANG W，DATCHEVA M，SCHANZ T，et al. A sub-stepping approach for elasto-plasticity with rotational hardening[J].Computational Mechanics，2006，37(3)：266-278.
[1]	HERNANDEZ JA，OLIVER J，CANTE JC，et al. A robust approach to model densification and crack formation in powder compaction processes[J].International Journal for Numerical Methods in Engineering，2011，87(8)：735-767.
[2]	BRANNON RM，LEELAVANICHKUL S. A multi-stage return algorithm for solving the classical damage component of constitutive models for rocks，ceramics，and other rock-like media[J].International Journal of Fracture，2010，163：133-149.
[1]	HOMEL MA，BRANNON RM. Relaxing the multi-stage nested return algorithm for curved yield surfaces and nonlinear hardening laws[J]. International Journal of Fracture，2015，194(1)：1-7.
[2]	HOMEL MA，GUILKEY JE，BRANNON RM. Numerical solution forplasticity models using consistency bisection and a transformed-space closest-point return：a nongradient solution method[J].Computational Mechanics，2015，56(4)：565-584.
[3]	BILOTTA A，LEONETTI L，GARCEA G. An algorithm for incremental elastoplastic analysis using equality constrained sequentialquadratic programming[J].Computers and Structures，2012，102：97-107.
[4]	CONTRAFATTO L，CUOMO A. A globally convergent numerical algorithm for damaging elasto-plasticity based on the multiplier method[J]. International Journal for Numerical Methods in Engineering，2005，63(8)：1089-1125.
[1]	PEREZ FA，ARMERO F. On the formulation of closest-point projection algorithms in elastoplasticity-part II：Globally convergent schemes[J].International Journal for Numerical Methods in Engineering，2002，53(2)：331-374.
[2]	DAFALIAS YF，MANZARI MT，AKAISHI M. A simple anisotropic clay plasticity model[J].Mechanics Research Communications，2002，29(4)：241-245.