Page 73 - 2024年第55卷第10期
P. 73
顶部施加，单元底部 节点 采用 法 向 固 定 约 束，底 部 一 角 点 采 用 三 向 链 杆 的 固 定 约 束，以 商 业 软 件
ＡＢＡＱＵＳ进行数值计算，所用的单元类型为 Ｃ３Ｄ８Ｒ八节点六面体单元。混凝土材料采用 ＣＤＰ模型，
根据实际试件的混凝土材料特性，选用的计算参数见表 １。
表 １ 混凝土数值模拟计算的力学参数 ［３０ － ３１］
弹性模量 泊松比 剪胀角 Ψ 偏心率 应力比 抗拉强度 ｆ 抗压强度 ｆ 断裂能 Ｇ ｆ
ｃ
ｔ
材料参数 Ｋ ｃ
?ＧＰａ ｖ ?（°） η σ σ ?ＭＰａ ?ＭＰａ ?（Ｎ?ｍ）
珚 ｃｃ ? 珚 ｃ
单向拉伸循环 ３０．５ ０．２ ３０ ０．１ ０．６６７ １．１６ ３．５ ４０ ５５
单向压缩循环 ３０ ０．２ ３０ ０．１ ０．６６７ １．１６ ２．８ ２８ ２００
图 ４和图 ５分别展示了采用 Ｎａｊａｒ定义的损伤变量式（８），以及本文定义的式（１２），包括 λ ＝ １和
根据试验曲线由式（ １４）计算的 λ值等三种情况下，损伤分析的结果，并与试验曲线 ［３０］ 进行比较。其
－ ６
ｐｔ
ε
中，根据受拉试验曲线实测的等效塑性应变 珔 ＝ ６２．４ × １０ ，由式（１４）计算得到的 λ ＝ ２．６７７３ ；根据受
ｐｔ
－ ６
压试验曲线实测的等效塑性应变 珔 ＝ ５６４．２ × １０ ，由式（１４）得到 λ ＝ １．４５８０。表 ２列出了图 ４（ａ）卸载
ε
过程中 ｂ的变化和图 ５（ａ）卸载过程中 ｂ的变化。
ｃ
ｔ
图 ４ 不同定义损伤因子时的受拉试验曲线损伤分析
图 ５ 不同定义损伤因子时的受压试验曲线损伤分析
表 ２ 循环荷载作用下卸载过程中 ｂ、ｂ数值
ｔ ｃ
受拉时（ λ ＝ ２．６７７３）ｂ ｔ 受压时（ λ ＝ １．４５８０）ｂ ｃ
卸载次数 １ ２ ３ ４ ５ １ ２ ３ ４
试验值 ０．８７５３ ０．７９１３ ０．７８４５ ０．８０３４ ０．７８５１ ０．６５０１ ０．５６８９ ０．５７８７ ０．５９１８
数值模拟 ０．７９３３ ０．６７１１ ０．７０１０ ０．７２８０ ０．７６２８ ０．６３５８ ０．５８７４ ０．５８２７ ０．５９１６
误差?％ － ９．４ － １５．２ － １０．６ － ９．４ － ２．８ － ２．２ ３．２ ０．７ ０．０
２
— １ ０ ３ —
