基于指定应力方法的混凝土开裂的三维数值模拟
3D Numerical Simulation of Concrete Cracking Based on Specified Stress Method
张晓庆 1王家林 1易志坚 1张拓 1王敏1
作者信息
- 1. 重庆交通大学土木工程学院,重庆 400074
- 折叠
摘要
基于指定应力方法,根据线弹性理论推导出了一种新的混凝土开裂的空间有限元列式,根据该有限元公式编制了相应的C++计算程序.通过3个算例,将该文算法与理论结果以及ABAQUS中XFEM计算结果进行对比,验证了该文开裂算法的正确性.与常规开裂算法相比,该文开裂算法开裂积分点的应力一旦指定为零(开裂状态),在后续计算中将不会出现非零,不需要迭代过程去调整为零,大大减少了迭代次数和每次迭代过程中需要处理的数据量;相较于ABAQUS中的XFEM算法只能使用一阶单元,本文开裂算法可以使用二阶单元进行开裂计算,在其他计算环境相同的情况下对开裂区域和开裂状态的判断更准确,这为商业有限元软件采用二阶单元进行更精细的开裂计算提供了一种新的途径和算法.
Abstract
Based on the specified stress method,a new spatial finite element formulation for concrete cracking was derived according to the linear elasticity theory.This formulation was used to create a calculation program with the C++language.The accuracy of the proposed cracking algorithm was validated through 3 numerical ex-amples,where the theoretical results were compared with the ABAQUS XFEM calculations.Beyond convention-al cracking algorithms,the proposed cracking algorithm has the advantage that once the stress at the cracking integration point is specified as zero(in the cracking state),it will remain zero in subsequent calculations.There is no need for an iterative process to adjust it to zero,which means significant reduction of the number of iterations and the amount of data processing required in each iteration.In comparison to the ABAQUS XFEM al-gorithm,which is limited to the 1st-order elements,the proposed cracking algorithm can utilize the 2nd-order elements for crack calculation,and allows for a more accurate determination of the cracked regions and states under the same computational conditions.This work provides a new approach and algorithm for commercial fi-nite element software to conduct more refined crack calculations with the 2nd-order elements.
关键词
指定应力方法/有限元/混凝土开裂/三维数值模拟Key words
specified stress methods/finite element/concrete cracking/3D numerical simulation引用本文复制引用
出版年
2024
国家科技期刊平台
NSTL
万方数据