结构动力弹塑性与倒塌分析(I)——滞回曲线改进、ABAQUS子程序开发与验证*
柳国环1,2,练继建1,2,国 巍3

(I)——滞回曲线改进、ABAQUS子程序开发与验证*] 柳国环1,2,练继建1,2,国 巍3 (1.天津大学 建筑工程学院,天津300072; 2.天津大学 水利工程仿真与安全国家重点实验室,天津300072; 3.中南大学 土木工程学院,湖南 长沙410075)

滞回模型; 弹塑性; ABAQUS; UMAT; VUMAT

Elasto-plasticity and Collapse Analysis for Large-span and High-rise Structure(I)——Hysteresis Relationship Improvement,ABAQUS Subroutine Development and Verification
LIU Guo-huan1,2,LIAN Ji-jian1,2,GUO Wei3

(1. School of Civil Engineering,Tianjin University,Tianjin 300072,China)(2. State Key Laboratory Hydraulic Engineering Simulation and Safety,Tianjin University,Tianjin 300072,China)(3. School of Civil Engineering and Architecture,Central South University,

hysteristic model; elasto-plasticity; ABAQUS; UMAT; VUMAT

备注

详述并明确了钢材静态损伤滞回模型中3种类型加卸载路径,重点改进了混凝土滞回规则中再加载路径拐点前的二次曲线与拐点后的直线规则,定义了再加载至骨架曲线的刚度退化指数型表达式。进而,开发了直接能够与ABAQUS主程序保持数据相互调用的UMAT隐式和VUMAT显式算法子程序,并给出子程序与ABAQUS主程序链接的框架示意图。最后,通过5种具有试验结果的构件对所开发程序加以检验。结果 表明:(1)改进的混凝土再加载滞回规则可较好描述混凝土加载曲线,克服了二次型规则的再加载曲线与卸载曲线在拐点前相交的可能性,体现了拐点后再加载至骨架曲线的刚度退化;(2)开发实现的隐式算法UMAT与显式算法双精度VUMAT子程序不仅可行,而且具有足够的计算精度。

Firstly,three types of unloading and reloading rules for steel static elasto-plastic damage hysteristic models are discussed in detail. The quadratic curve and linear rules respectively before and after inflection point of reloading curves are improved emphatically in the concrete hysteretic rules,and a new experiential expression for the stiffness degrading from reloading to skeleton curves is proposed as well. Secondly,the subroutines are developed respectively to remain the direct call of data between the ABAQUS main program and implicit UMAT and explicit VUMAT algorithms,the schematic framework of relationship between the subroutines and ABAQUS main program is given in detail. At last,the developed program was tested by five kind of typical frames(steel rebar,plain concrete,reinforced-concrete,steel-pipe concrete and steel-reinforced concrete),which have the test results. The analysis results show that:(1)The improved concrete reloading hysteretic and concrete loading curve consistent well,and can avoid the intersection phenomenon between unloading and reloading curves of quadratic form rule before inflection point and reflects the stiffness degradation from reloading to skeleton curves after inflection point.(2)The developed subroutines which are implicit UMAT and explicit double pricision VUMAT algorithms are feasible and possesses high accuracy.