土石料双屈服面弹塑性模型的二次开发算法与应用

为了丰富ADINA、ABAQUS和ANSYS等大型通用有限元软件中岩土类材料本构模型，增强软件对某些土工结构的精细模拟计算能力，以土石料沈珠江双屈服面弹塑性模型为例，详细介绍了模型在商用有限元软件中的二次开发流程，对应力更新算法进行了改进，通过文献中粗粒土三轴试验数据进行对比分析，验证了改进算法的优越性.同时，采用模型面板坝资料进一步验证二次开发流程和算法的准确性，最后进行了某面板堆石坝工程的应用例证.本文所述应力积分算法和开发流程具有通用性，可供大型通用商业有限元软件开发双屈服面和多屈服面弹塑性本构模型时借鉴使用. 

Algorithm Development and Application of Double-Yield-Surface Elastoplastic Constitutive Model for Earth-Rock Materials

To enrich the constitutive models of geotechnical materials in universal commercial finite-element analysis software suites, e.g. ADINA, ABAQUS, and ANSYS, and to improve the precise simulation ability for certain geotechnical structures, Shen's double-yield-surface elastoplastic model was chosen as an example to introduce in detail the secondary development process of a constitutive model in commercial finite-element software. The stress integration algorithms were modified, and the superiority of the algorithms was investigated by comparing the numerical simulation with triaxial test data for coarse-grained soil. Further, an ideal concrete-face rockfill dam was used to validate the accuracy of the secondary development process and the algorithms. Finally, the numerical simulation of a practical concrete-face rockfill dam was conducted. The results indicate that the stress integration algorithms and the secondary development process presented in this study are universal and can be employed for the secondary development of double or multi-yield-surface elastoplastic models in universal finite-element analysis software.