Algorithm Development and Application of Double-Yield-Surface Elastoplastic Constitutive Model for Earth-Rock Materials
-
摘要: 为了丰富ADINA、ABAQUS和ANSYS等大型通用有限元软件中岩土类材料本构模型,增强软件对某些土工结构的精细模拟计算能力,以土石料沈珠江双屈服面弹塑性模型为例,详细介绍了模型在商用有限元软件中的二次开发流程,对应力更新算法进行了改进,通过文献中粗粒土三轴试验数据进行对比分析,验证了改进算法的优越性.同时,采用模型面板坝资料进一步验证二次开发流程和算法的准确性,最后进行了某面板堆石坝工程的应用例证.本文所述应力积分算法和开发流程具有通用性,可供大型通用商业有限元软件开发双屈服面和多屈服面弹塑性本构模型时借鉴使用.Abstract: 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.
-
表 1 蓄水期大坝各物理量极值对比
Table 1. Extreme-value comparison at water-impounding stage
项目 水平位移/cm 沉降位移/cm 面板挠度/cm 面板拉应力/MPa 坝体大主应力/MPa 向上游 向下游 本文结果 4.18 12.98 33.50 10.28 1.69 1.826 文献[14] 4.04 10.67 34.65 10.53 2.16 1.826 误差/% 3.47 21.65 -3.32 -1.42 -21.76 0 表 2 坝体及覆盖层沈珠江模型计算参数
Table 2. Constitutive material parameters
材料 ρ/(g·cm-3) φ0/(°) Δφ/(°) KE nk Rf cd nd Rd Kur 覆盖层 2.22 51.8 8.0 600 0.65 0.82 0.007 5 0.40 0.76 900 主堆石 2.15 56.1 8.0 700 0.35 0.77 0.008 0 0.35 0.66 1 050 次堆石 2.23 47.9 5.0 650 0.32 0.76 0.007 2 0.31 0.61 975 垫层料 2.25 57.3 8.0 800 0.50 0.80 0.005 5 0.75 0.68 1 200 过渡层料 2.10 56.1 8.0 750 0.50 0.79 0.006 8 0.50 0.52 1 125 -
熊堃, 岑威钧, 胡清义, 等.多途径综合开发商业软件精细求解土石坝结构静动力反应[J].岩石力学与工程学报, 2013, 32(1):117-125. doi: 10.3969/j.issn.1000-6915.2013.01.017XIONG Kun, CEN Weijun, HU Qingyi, et al. Static and dynamic response analysis of earth rockfill dams with united multipath development of commercial softwares[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(1):117-125. doi: 10.3969/j.issn.1000-6915.2013.01.017 LIU Changhong, LIU Xiaoling, CHEN Qiu. Interval analysis of the finite element method for stochastic structures[J]. Journal of Southwest Jiaotong University, 2001, 12(1):46-48. http://d.wanfangdata.com.cn/Periodical_xnjtdxxb-e200401008.aspx 熊玉春, 房营光. ADINA有限元软件中材料本构的二次开发[J].岩土力学, 2008, 29(8):2221-2225. doi: 10.3969/j.issn.1000-7598.2008.08.037XIONG Yuchun, FANG Yingguang. Secondary development of material constitutive model in ADINA software[J]. Rock and Soil Mechanics, 2008, 29(8):2221-2225. doi: 10.3969/j.issn.1000-7598.2008.08.037 江守燕, 谢庆明, 杜成斌.基于ABAQUS平台的邓肯-张E-B和E-v模型程序开发[J].河海大学学报:自然科学版, 2011, 39(1):61-65. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=SciencePaper201309240000124134JIANG Shouyan, XIE Qingming, DU Chengbin. Development of program of Duncan-Chang E-B and E-v models based on ABAQUS[J]. Journal of Hohai University:Natural Sciences, 2011, 39(1):61-65. http://www.wanfangdata.com.cn/details/detail.do?_type=perio&id=SciencePaper201309240000124134 孙明权, 陈姣姣, 刘运红.邓肯-张E-B模型的ANSYS二次开发及应用[J].华北水利水电学院学报, 2013, 34(2):30-34. doi: 10.3969/j.issn.1002-5634.2013.02.010SUN Mingquan, CHEN Jiaojiao, LIU Yunhong. The second development of Duncan-Chang E-B model in ANSYS and its application[J]. Journal of North China Institute of Water Conservancy and Hydroelectric Power, 2013, 34(2):30-34. doi: 10.3969/j.issn.1002-5634.2013.02.010 岑威钧, 朱岳明.基于ABAQUS的土石料本构模型二次开发及其应用[J].水利水电科技进展, 2005, 25(6):78-81. doi: 10.3880/j.issn.1006-7647.2005.06.022CEN Weijun, ZHU Yueming. ABAQUS-based secondary development of constitutive model for earth rockfill materials and its application[J]. Advances in Science and Technology of Water Resources, 2005, 25(6):78-81. doi: 10.3880/j.issn.1006-7647.2005.06.022 司海宝, 化西婷.南水模型在ABAQUS中的实现及在工程中的应用[J].南水北调与水利科技, 2010, 8(1):52-55. http://d.old.wanfangdata.com.cn/Periodical/nsbdyslkj201001013SI Haibao, HUA Xiting. Development of NHRI constitutive model in ABAQUS and application in engineering[J]. South-to-North Water Transfers and Water Science & Technology, 2010, 8(1):52-55. http://d.old.wanfangdata.com.cn/Periodical/nsbdyslkj201001013 关云飞, 高峰, 赵维炳, 等. ANSYS软件中修正剑桥模型的二次开发[J].岩土力学, 2010, 31(3):976-980. doi: 10.3969/j.issn.1000-7598.2010.03.053GUAN Yunfei, GAO Feng, ZHAO Weibing, et al. Secondary development of modified Cambridge model in ANSYS software[J]. Rock and Soil Mechanics, 2010, 31(3):976-980. doi: 10.3969/j.issn.1000-7598.2010.03.053 WISSMAN J W, HAUCK C. Efficient elasto-plastic finite element analysis with higher order stress point algorithms[J]. Computers & Structures, 1983, 17(1):89-95. SLOAN S W, ABBO A J. Biot consolidation analysis with automatic time stepping and error control. part Ⅰ:theory[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1999, 23:467-492. doi: 10.1002/(ISSN)1096-9853 BORIJA R I, LEE S R. Cam clay plasticity, part Ⅰ:implicit integration of constitutive relations[J]. Computer Methods in Applied Mechanics and Engineering, 1990, 78(1):49-72. doi: 10.1016/0045-7825(90)90152-C POTTS D, GANENDRA D. An evaluation of substepping and implicit stress point algorithms[J]. Computer Methods in Applied Mechanics and Engineering, 1994, 119(3):341-354. http://www.sciencedirect.com/science/article/pii/0045782594900949 沈珠江.土体应力应变分析的一种新模型[C]//第五届全国土力学及基础工程学术会议论文选集.北京: 中国建筑工业出版社, 1990. 冀春楼.深厚覆盖层上高堆石坝静、动力分析方法的研究-冶勒堆石坝静、动工作性态研究[D].南京: 河海大学, 1995. http://www.wanfangdata.com.cn/details/detail.do?_type=degree&id=Y218966 CEN Weijun, WEN Langshen, ZHANG Ziqi, et al. Numerical simulation on seismic damage and cracking of concrete slab for high concrete face rockfill dams[J]. Water Science and Engineering, 2016, 9(3):205-211. doi: 10.1016/j.wse.2016.09.001