扩展有限元法及其在钢筋混凝土结构中的应用综述.pdf

水利水电科技进展,2015,35(3)摇 Tel:[phone]摇 E鄄mail:jz@ hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 第 35 卷第 3 期 Vol. 35 No. 3 水 利 水 电 科 技 进 展 Advances in Science and Technology of Water Resources 2015 年 5 月 May 2015 基金项目:国家大坝安全工程技术研究中心开放基金(NDSKFJJ1101) 作者简介:李俊(1990—),男,江苏沛县人,硕士研究生,主要从事混凝土结构研究。 E鄄mail:lijunhhucn@163. com DOI:10. 3880/ j. issn. 10067647. 2015. 03. 020 扩展有限元法及其在钢筋混凝土结构中的应用综述 李摇 俊,汪基伟,冷摇 飞 (河海大学土木与交通学院,江苏 南京摇 210098) 摘要:简述扩展有限元法的基本思想和应用价值,以钢筋混凝土的裂纹问题(强不连续问题)为例 对扩展有限元法中描述裂纹的水平集法、不连续位移场富集策略、积分方案进行了介绍,并给出相 关领域近年来的研究成果。 概述扩展有限元法在断裂力学领域中的应用动态。 在此基础上,分别 从混凝土开裂和钢筋与混凝土两种材料的黏结滑移两个方面介绍扩展有限元法在钢筋混凝土结构 开裂分析中的研究与应用进展,指出为将该方法应用于水工钢筋混凝土结构等规模巨大、配筋和受 力复杂的实际工程中,需对单元模型、裂缝处理方式以及裂缝处钢筋应力等问题进行深入研究。 关键词:扩展有限元法;钢筋混凝土结构;断裂力学;综述 中图分类号:TV313摇 摇 摇 文献标志码:A摇 摇 摇 文章编号:10067647(2015)03010608 Review of extended finite element method and its application in reinforced concrete structures/ / LI Jun, WANG Jiwei, LENG Fei (College of Civil and Transportation Engineering, Hohai University, Nanjing 210098, China) Abstract: The basic theory and application value of extended finite element method(XFEM) are introduced briefly. The level set method for describing crack, the enrichment strategy of discontinuous displacement field, and the integration scheme of XFEM are presented with an example of crack problem(strong discontinuity) of reinforced concrete(RC), and recent achievements in these fields are exposed. Thereby, a summary of the situation of applications of XFEM in fracture machanics is given. On the basis of the foregoing, the research and application development of XFEM for the analysis of reinforced concrete structure crack are presented from two aspects: concrete cracking, and bond鄄slip between reinforcement and concrete. It is pointed out that element models, crack treatment, steel stress at crack section and other issues have to be further investigated in order to meet the requirements of characteristics of practical projects(e. g. , hydraulic reinforced concrete structures), such as large scale, complex reinforcement, and multi鄄axial force state. Key words: extended finite element method; reinforced concrete structures; fracture machanics; review 摇 摇 20 世纪末期,Mo觕s 等[1]在标准有限元法的基 础上,结合单位分解法、水平集法和断裂力学研究成 果,提出了扩展有限元法(extended finite element method, XFEM),为不连续问题提供了新的解决方 法。 该方法通过在标准位移形函数中增加具有不连 续性质的富集函数,在单元内部描述位移间断(强 不连续)、应变间断(弱不连续),使得有限元网格独 立于裂纹几何路径,在追踪裂纹扩展时亦无需重剖 网格,提高了裂纹问题的求解效率。 在水利水电工程中,钢筋混凝土结构应用广泛, 其限裂分析对结构的正常使用性能及耐久性有重要 意义,重要的大体积非杆系混凝土结构的限裂分析 通常需要借助非线性有限元进行数值模拟[2]。 传 统钢筋混凝土有限元中,弥散裂缝模型应用广 泛[3],该模型不需随开裂区扩展而更新网格,适用 于大型工程,但裂缝分布及裂缝宽度有一定网格依 赖性,且无法准确模拟裂缝处的钢筋应力。 事实上, 钢筋混凝土结构开裂分析时,不仅要考虑混凝土基 体的位移不连续,还要考虑钢筋的加强作用及其与 混凝土之间的黏结滑移,较为复杂。 将扩展有限元 法应用于钢筋混凝土结构开裂分析,能调和传统有 限元计算效率和计算精度之间的矛盾,有较高的理 论价值和广阔的应用前景。 1摇 扩展有限元法基本原理 1. 1摇 水平集法 水平集法(level set method, LSM)可以从几何 上描述静态或动态裂纹,以及孔洞、夹杂等不连续界 601 水利水电科技进展,2015,35(3)摇 Tel:[phone]摇 E鄄mail:jz@ hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 面,它将界面的变化表示为比界面高一维(时间 t) 的水平集曲线[4],设不连续界面为 祝(t),任意考察 点 x沂祝(t)满足 准(x,t) = 0(1) 式中:准(x,t)为考察点水平集函数。 对于裂纹问题,扩展有限元中常用符号距离函 数作为其水平集函数,设 x祝为不连续界面上距 x 最 近的点,有 准(x,t) = 依 min x祝沂祝(t)椰x - x祝椰 (2) 摇 摇 Sukumar 等[5]介绍了如何在扩展有限元中采用 水平集法描述孔洞、夹杂等弱不连续界面。 Ventura 等[6]采用水平集法描述裂纹扩展过程,提出一种改 进的向量水平集公式,使得在描述裂纹扩展的同时 不改变已有裂纹的水平集值,并且裂纹沿原裂纹面 法向扩展时,新裂纹水平集值也不会改变。 水平集 法可以精确刻画不连续面的几何性质及其动态变 化,易于从平面问题推广到三维问题[7]。 1. 2摇 裂纹问题的扩展有限元富集策略 Melenk 等[8]研究了有限元形函数的单位分解 法,单位分解法的基本思想是任意函数 F(x)都可 以在计算域内表示为如下形式: F(x) =移 I [NI(x)f(x)](3) 其中 NI(x)满足单位分解,即 移 I NI(x) = 1(4) 式中:NI(x)为标准有限元形函数;I 为求解域节点 集;f(x)为 F(x)在计算域内的局部近似函数。 单位分解法为改进标准有限元形函数提供了理 论依据。 基于单位分解法,可以根据不连续问题的 特点,在标准有限元位移近似场的基础上增加富集 项,用来模拟单元内部存在的强、弱不连续性。 扩展有限元中,改进后的位移场 uh(x)由标准 有限元部分和富集项叠加得到[9]: uh(x) =usta(x) + uenr(x) =移 i沂I [Ni(x)ui,sta] + 移 i沂I* [Ni *(x)追(x)a i] (5) 式中:usta(x)为标准有限元位移;uenr(x)为富集位 移; ui,sta为标准有限元位移未知量;Ni(x)为标准有 限元形函数;Ni *(x)为单位分解函数;追(x)为富集 函数;ai为附加自由度;I*为富集节点集。 Ni *(x) 可以和 Ni(x)相同,也可不同,为了编制程序方便, 可取为相同。 不连续问题包括强不连续问题(如裂纹问题中 的位移不连续)、弱不连续问题(如夹杂问题和复合 材料问题中的应变不连续)等,不同类型的问题需 要不同的富集位移构建方法。 以下以二维裂纹问题 为例,对扩展有限元的位移富集思路进行说明。 图 1 为包含任意形状裂纹的二维规则化网格,裂纹所 在单元的相关节点为富集节点,裂纹影响单元分为 裂纹贯穿单元(图 1 中 2 ~ 6 号单元)和裂尖单元 (图 1 中 1、7 号单元)两部分,相应富集节点亦分为 两部分,即裂尖单元富集节点和裂纹贯穿单元富集 节点,若某节点同时属于裂纹贯穿单元和裂尖单元, 为保证裂尖位移场计算精度,令其优先属于裂尖富 集节点[1]。 富集单元和标准单元之间由混合单元 (图 1 中阴影区单元)连接,混合单元中同时包含富 集节点和标准节点(不参与富集的节点)。 图 1摇 包含任意形状裂纹的二维网格 对裂纹贯穿单元,采用阶跃函数构建富集函数, 以反映裂纹两侧的位移不连续;对裂尖单元,采用断 裂力学中裂尖位移场解析解的各项作为基函数构建 富集函数[10]。 参照式(5),二维裂纹问题的位移模 式可统一写为 uh(x) =移 i沂I [Ni(x)ui,sta] +移 j沂N祝[N j(x)H(x)aj] + 移 k沂N撰 Nk(x)移 4 w = 1 [Fw(x)bw,k{}] (6) 其中 H(x) = sign[准(x,t)] = 1摇 准(x,t) > 0 - 1摇 准(x,t) < { 0 (7) Fw(x) = rsin 兹 2 摇 w = 1 rcos 兹 2 摇 w = 2 rsin 兹 2 sin兹摇 w = 3 rcos 兹 2 sin兹 摇 w = 4 (8) 式中:N祝为裂纹贯穿单元富集节点集;N撰为裂尖单 元富集节点集;aj和 bw,k分别为上述两类富集节点 的附加自由度;Nj(x)、Nk(x) 为单位分解函数; H(x)为阶跃函数;Fw(x)为裂尖函数;r、兹 为裂尖极 701 水利水电科技进展,2015,35(3)摇 Tel:[phone]摇 E鄄mail:jz@ hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 坐标系所定义的位置参数。 Belytschko 等[11]对裂纹贯穿单元富集函数做出 修正,在断裂问题中应用较多: 追i(x) = H[准(x)] - H[准(xi)](9) 摇 摇 将 追i(x)代替式(5)中 追(x),并令 x=xj(xj为 混合单元节点坐标),有 uh(xj)= uj,sta,式中 uj,sta为 j 节点的标准有限元位移未知量。 可以看到,改进后 得到的节点位移为标准有限元形函数的系数,这一 点和标准有限元一致,便于位移后处理,其原因在于 对富集函数进行修正后,混合单元内位移函数富集 项取值为零,即混合单元中富集节点的附加自由度 不对混合单元产生影响,只对裂纹贯穿单元内部产 生影响。 对于裂尖置于单元内部的情况,可在富集 函数中引入斜坡函数 R(x) [12],以提高混合单元收 敛速度,此时用式(10) 中的 追i(x) 代替式(5) 中 的 追(x): 追i(x) = Ni(x)追(x)R(x)(10) 其中 R(x)取值为移 i沂I* Ni *(x) ,可以反映附加自由 度在混合单元内的渐变影响。 1. 3摇 扩展有限元的积分方案 扩展有限元中形函数具有单位分解属性,使得 刚度矩阵保持了对称、稀疏的特性,这一点和传统有 限元一致。 但对于不连续单元(如裂尖单元和裂纹 贯穿单元),在获得刚度矩阵时,其积分不能依照传 统连续单元进行。 不连续单元的积分方案中,分区 域积分较为常用[13]。 对于裂纹贯穿单元,由于裂纹两侧子区域的位 移场是连续的,所以可以在子区域分别积分,然后将 子区域的积分叠加得到单元刚度矩阵。 对于裂尖单 元,当以裂纹尖端为顶点分成若干小三角形进行分 区域积分时,裂尖处会出现奇异性,Laborde 等[14]介 绍了通过坐标变换消除奇异性的方法。 由于分区域积分方案计算量较大,也有学者提 出采用连续积分方案作为替代。 Ventura[15]针对简 单富集函数提出一种可将分区域积分转化为连续高 斯积分的方法,但当单元形态不是平行四边形时有 明显误差。 Song 等[16]对扩展有限元中位移基函数和 自由度进行重置,采用叠合单元和虚拟节点表示裂 缝,在此基础上提出简易积分方案,该方案每个单元 设置一个积分点,为保证矩阵非奇异,需采取控制措 施避免零能模式。 喻葭临等[17]结合文献[16]的思路 和传统双线性四边形单元的四节点积分方案,提出同 时适用于含裂缝单元与标准单元(不含裂缝)的简易 积分方案,避免了分区域积分的复杂计算,亦无需对 零能模式采取控制措施,有较高的实用价值。 2摇 扩展有限元法应用概况 扩展有限元允许在单元内部描述裂缝,裂缝扩 展时无需重剖网格,同时保留了传统有限元求解非 线性问题的优势[18],被广泛应用于断裂力学领域。 Sukumar 等[19鄄20]针对准静态裂纹扩展问题对扩 展有限元数值实现步骤进行说明。 董玉文等[21]采 用扩展有限元法计算应力强度因子,指出其精度和 采用 J 积分计算得到的应力强度因子精度相当。 除 静态和准静态裂纹问题外,扩展有限元法亦被应用 于动态裂纹问题中[22鄄23]。 Belytschko 等[22]采用渐进 的动态裂纹尖端解作为基函数,把 I 型和 II 型应力 强度因子作为附加自由度,使用显示积分方案求解 弹性动态裂纹扩展问题。 Menouillard 等[23]对显示 积分系统的稳定性进行研究,并给出质量矩阵的对 角化方法。 扩展有限元的优势之一在于不预设开裂 路径,尽可能地还原裂缝的真实形态,这一点在分支 裂纹、交叉裂纹[24鄄25]的模拟中得到展现。 真实裂纹之间可能存在摩擦接触,在水工结构 中,裂纹中还可能存在渗透压力作用,是否考虑裂缝 内的作用对裂缝预测有较大影响。 Dolbow 等[26]采 用扩展有限元法研究裂纹面间的摩擦接触,采用非 线性本构关系来描述裂纹面的接触,并采用 LATIN 迭代法求解;余天堂[27]采用线性互补法求解裂纹面 的非线性接触,能够避免复杂的迭代求解过程;李建 波等[28鄄29]采用虚功原理推导了考虑单元内部裂纹 面上分布荷载及缝内粘连的扩展有限元基本公式, 提出和标准有限元协调一致的扩展有限元刚度矩阵 形成模式。 三维断裂问题一直是研究的热点和难点, Sukumar 等[30]首次利用扩展有限元法研究了三维 裂纹扩展问题,其裂尖富集函数建立在与裂纹尖端 垂直的平面内,和平面问题类似,其函数形式仍采用 极坐标表示。 三维裂纹扩展问题的难点在于裂纹面 的连续性与光滑性较难模拟,Duan 等[31]利用单元 水平集描述三维裂缝,采用最小二乘法改进了原裂 缝面方向和新预测开裂方向的一致性,从而改善了 三维裂缝面的光滑性。 除断裂力学领域之外,扩展有限元在复合材料 界面失效(同时包含强、弱不连续) [32]以及变形局部 化(如剪切带、损伤过程区) [16,33] 等领域也有应用。 值得注意的是,扩展有限元研究集中于线弹性断裂 力学领域,当材料属性发生变化时,其断裂性能也会 发生变化,相关参数如断裂能和富集函数构造形式 也会发生变化,应根据所求问题具体分析。 801 水利水电科技进展,2015,35(3)摇 Tel:[phone]摇 E鄄mail:jz@ hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 3摇 扩展有限元法在钢筋混凝土结构开裂分 析中的应用 摇 摇 钢筋混凝土结构开裂分析时存在两种不连续问 题:一是混凝土开裂导致的位移不连续;二是钢筋与 混凝土两种材料的黏结滑移。 扩展有限元法在上述 两种不连续模拟中均有应用,现简述如下。 3. 1摇 在混凝土开裂分析中的应用 混凝土断裂力学认为在裂缝尖端区域存在断裂 过程区,断裂过程区是微细裂缝发生、发展并逐步转 化为宏观裂纹的区域。 黏聚裂纹模型(又称凝聚力 模型) [34]采用宏观裂纹尖端的黏聚区表示断裂过程 区,黏聚区的传力机理通过应力-张开位移曲线进 行描述,该曲线与坐标轴所围面积代表相应材料的 断裂能。 黏聚裂纹模型能够较真实地反映裂缝尖端受力 状态,在包括混凝土在内的准脆性材料断裂分析中 得到较多应用。 Mo觕s 等[35]把扩展有限元法引入黏 聚裂纹分析中,认为稳定裂纹的裂尖 I 型应力强度 因子保持为零,并以此作为裂纹是否扩展的评判标 准,针对黏聚裂纹尖端位移场和线弹性断裂问题的 位移场有所不同,提出采用 rsin 兹 ( 2 或 r3/2sin 兹 2 或 r2sin 兹 ) 2 对线弹性扩展有限元的裂尖场富集项进行 修正。 Mariani 等[36]利用扩展有限元位移富集技术 构建三次裂尖位移场,从而可以描述黏聚裂纹模型 的尖角状裂尖。 实际混凝土开裂分析中,需要确定裂缝开展准 则,包括开展方向和裂缝扩展步长。 Unger 等[37]把 扩展有限元法和自适应网格加密算法相结合以追踪 混凝土结构黏聚裂纹的扩展,将采用不同裂纹扩展 准则的数值算例计算结果与试验结果进行对比研 究,结果表明基于线弹性断裂力学的最大能量释放 率准则在预测开裂方向时有较好的稳定性。 方修君 等[38]在研究混凝土梁的复合型开裂时,推荐采用简 化的最大切向应力准则确定裂缝扩展方向。 裂缝开 展步长尚缺乏明确的理论研究,多根据具体问题选 择一个较小值。 黏聚裂纹模型可以较为真实地模拟裂尖受力状 态,但需预先设定起裂位置,杜效鹊等[39]结合扩展 有限元和黏聚裂纹模型,研究了预制缝重力坝模型 的断裂特性,得出与试验结果一致的荷载响应曲线 和裂缝扩展路径;张晓东等[40]利用黏聚裂纹模型研 究了带初始边缘裂纹的混凝土板在单向拉伸作用下 的裂纹扩展过程,效果良好。 然而实际工程中并不 是总能预先获得开裂位置,而且一般有多条裂纹发 生,此时过分关注裂尖稳定是不必要的。 陈胜宏 等[41]采用扩展有限元法对小湾拱坝坝踵开裂进行 分析,采用最大拉应力准则判断裂缝是否扩展,并且 认为裂缝发生后即贯穿整个单元,该处理方法便于 处理大型复杂结构,但计算结果表明裂缝扩展范围 与网格尺寸有一定关系。 在细观层次上,混凝土被认为是由骨料、硬化砂 浆和二者之间的过渡区组成的复合材料,其中骨料 和硬化砂浆交界处易产生微细裂纹,微细裂纹的发 展和贯通导致宏观裂缝的出现。 采用扩展有限元法 模拟混凝土细观断裂过程需要解决以下问题:淤不 同材料间的弱不连续模拟,即复合材料界面(过渡 区)的模拟;于微细裂纹产生、发展导致的强不连续 模拟;盂复合材料界面附近断裂参数的确定。 江守 燕等[42]采用扩展有限元法对包含圆形随机骨料的 多夹杂问题进行弱不连续分析,用水平集法描述材 料界面,采用标准网格进行离散,减轻了网格剖分的 负担,但该研究假定材料为线弹性,不涉及微细裂纹 的产生和扩展。 杜修力等[43]采用扩展有限元法模 拟了混凝土微细裂纹的产生和发展过程,将骨料周 围 2mm 厚度范围内设置为过渡区,给出骨料、过渡 区、砂浆的力学参数(开裂前假定三者都是线弹性 材料,开裂后采用线性的应力-裂纹宽度关系表征 混凝土软化特性)。 由于假定裂纹产生后随即贯穿 整个单元,为保证收敛性和求解精度,该研究网格剖 分较为精细,对过渡区需单独考虑,计算成本较大。 于红军[44]研究了含复杂界面的非均匀材料的断裂 参数确定方法,采用相互作用积分计算位于界面附 近或界面上的裂尖应力强度因子,验证了相互作用 积分的稳定性和区域无关性。 上述研究将扩展有限元法应用到素混凝土的开 裂模拟中,可为扩展有限元法应用于钢筋混凝土结 构提供参考,但由于没有考虑钢筋在限制混凝土开 裂中起到的作用,因而无法直接应用于钢筋混凝土 结构开裂分析。 3. 2摇 在钢筋与混凝土界面模拟中的应用 不同于素混凝土结构,钢筋混凝土结构需要考 虑钢筋的加强作用以及钢筋与混凝土两种材料的相 互作用,其中钢筋与混凝土的黏结滑移是关键问题, 是否考虑黏结滑移对裂缝分布及裂缝宽度有明显影 响。 茹忠亮等[45]用扩展有限元法对预设初始裂缝 的钢筋混凝土梁复合断裂过程进行模拟,假定钢筋 与混凝土之间没有滑移,得到的极限承载力较试验 值偏大,钢筋布置区的裂缝分布和试验情况也有较 大差别。 传统钢筋混凝土有限元模型依据钢筋模拟 901 水利水电科技进展,2015,35(3)摇 Tel:[phone]摇 E鄄mail:jz@ hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 方法的不同分为 3 种:分离式、组合式、整体式,其中 分离式和组合式模型可以模拟钢筋与混凝土的黏结 滑移,在裂缝分析中应用较多[46],也可在扩展有限 元分析中借鉴采用。 Simone 等[47]把扩展有限元的思想引入钢筋与 混凝土界面的模拟,利用有限元形函数的单位分解 性质,对钢筋与混凝土界面上的节点预设附加自由 度以考虑二者的相对滑移。 Deb 等[48]在研究灌浆 锚杆时,定义含锚杆的实体单元为富集单元,富集单 元有附加自由度,用于确定锚杆的位移、应力,富集 单元的刚度矩阵由岩体、 砂浆、 锚杆共同组成。 Simone 等[47]的研究基于分离式模型,Deb 等[48]的 研究基于组合式模型,但二者在分析中没有考虑混 凝土 开 裂 或 岩 体 开 裂 导 致 的 不 连 续 位 移。 Ibrahimbegovic 等[49]基于平面组合式模型,在包含 钢筋的混凝土单元位移场中加入黏结滑移附加项, 基体开裂则采用内嵌不连续模型进行模拟,然后把 位移分解为不计滑移时钢筋混凝土的位移与黏结滑 移引起的位移之和,该研究优点在于同时考虑了基 体开裂和黏结滑移两种不连续。 何学[50]建立了分 离式钢筋混凝土扩展有限元模型,采用四节点扩展 有限元等参单元模拟混凝土基体,采用杆单元模拟 钢筋,并对传统四边形黏结单元进行改进,与钢筋单 元相连的部分保持位移连续,与混凝土相连的部分 添加由裂缝产生的附加位移,使之可以考虑裂缝引 起的混凝土位移间断。 文献[49-50]保留了传统钢 筋混凝土有限元的架构,对混凝土、钢筋及二者的黏 结滑移分别考虑,适用于不同配筋及复杂受力状况, 为扩展有限元法应用于实际钢筋混凝土工程作出了 有益探索。 非杆系钢筋混凝土结构常通过控制钢筋应力间 接控制裂缝宽度[2],因而准确得到裂缝处的钢筋应 力有重要意义。 Contrafatto 等[51]研究了规则配筋轴 心受拉构件的钢筋应力分布,用变分原理推导了考 虑基体开裂、黏结滑移的基本方程,探讨了混凝土形 函数和钢筋形函数的改进方法,从而能够描述由裂 缝、滑移引起的局部大梯度应力(应变)场,该研究 对准确求得复杂配筋及复杂受力状态下裂缝处钢筋 应力有参考意义。 另外,通用有限元软件 ABAQUS 在原有的断裂 力学分析功能基础上增加了扩展有限元的分析模 块,许多学者基于此进行了素混凝土构件和钢筋混 凝土构件的断裂分析。 可以充分利用通用有限元程 序的非线性求解功能[52],也可通过加入弹簧单元的 方法考虑钢筋与混凝土的黏结滑移[53],目前该软件 的扩展有限元分析功能仍处在不断发展和丰富中。 3. 3摇 需要继续研究的问题 由于钢筋混凝土结构开裂机理的复杂性,多数 研究侧重于两种不连续问题中的一个方面,少数研 究对两种不连续问题进行统一考虑[49鄄51,53]。 扩展有 限元的思想可以和传统钢筋混凝土有限元结合(文 献[50]称之为钢筋混凝土扩展有限元),为同时解 决钢筋混凝土结构中的两类不连续问题提供思路。 随着结构设计中正常使用性能及耐久性的重要性逐 渐增大,需要更高效的裂缝分布计算方法,既能保证 一定的计算精度,其计算成本又不是太高,扩展有限 元在该领域有广阔的应用前景。 但是,目前钢筋混 凝土扩展有限元还处于研究阶段,尚有许多问题需 要深入研究: a. 构建合适的单元模型,能够模拟钢筋与混凝 土的黏结滑移,同时满足实际工程中复杂配筋的需 要。 传统钢筋混凝土有限元中,组合式模型可以考 虑任意方向、任意数量的配筋,通过界面单元模拟钢 筋与混凝土的黏结滑移,可借鉴文献[50]的方法对 传统界面单元进行改进,在界面单元与混凝土相连 处增加附加位移项,使之能够考虑混凝土开裂带来 的影响。 b. 构建合适的裂缝模型,既能满足裂缝分布计 算精度要求,又不过分复杂,符合工程需求。 断裂力 学裂缝模型需要较多参数,在分析、设计时存在较多 不确定性,而且钢筋混凝土结构中由于钢筋参与工 作,混凝土裂缝的开展受到限制,裂缝不会因裂尖处 应力集中导致不稳定扩展,过分关注裂尖应力状态 是不必要的,因而传统的强度理论(如最大拉应力 准则)在裂缝分析、限裂设计中仍有应用价值。 此 外,通常大体积混凝土结构裂缝开展深度较大,在较 小的网格尺寸下,仍可假定裂缝产生后随即贯穿整 个单元。 c. 多裂纹问题及复合型开裂问题。 当采用断 裂力学裂缝模型时,单一裂纹的 I 型开裂已有较多 研究,但多裂纹、交叉裂纹以及复合型开裂问题研究 较少。 多裂纹以及复合型开裂问题需要建立相应的 位移模式、裂纹扩展准则和裂纹汇合准则,随着裂纹 数目增多,裂纹间相互作用更趋复杂。 钢筋混凝土 结构一般会产生较多裂缝,因而要反映其真实开裂 情况仍有一定难度。 如采用传统强度理论作为开裂 准则,可简化分析,但仍需对位移模式以及裂缝间作 用进行研究。 d. 基体开裂情况下钢筋与混凝土界面滑移的 高精度模拟以及裂缝处钢筋应力的高精度算法。 裂 缝处钢筋应力和裂缝宽度具有正相关性,同时钢筋 应力也是结构(构件)承载力的重要指标。 钢筋一 011 水利水电科技进展,2015,35(3)摇 Tel:[phone]摇 E鄄mail:jz@ hhu. edu. cn摇 http:/ / www. hehaiqikan. cn 般可采用杆单元进行模拟,杆单元的平均应力即代 表该单元内钢筋应力,而事实上由于黏结切应力的 存在,单元内钢筋应力分布是不均匀的,当钢筋横穿 裂缝时,裂缝处钢筋应力有突跃现象。 可考虑利用 扩展有限元位移富集思想构建能够描述裂缝处局部 变形集中的钢筋形函数,非均匀有理 B 样条曲线有 良好的局部修改能力,且有良好的可导性,可用来准 确模拟裂缝处钢筋应力突跃。 4摇 结摇 论 a. 扩展有限元法充分利用有限元形函数的单 位分解属性,使得不连续界面的描述独立于网格剖 分,在追踪不连续界面扩展时亦无需更新网格,同时 由于其保留了标准有限元求解非线性问题的优势,因 而在多种材料的强、弱不连续问题研究中应用广泛。 b. 将扩展有限元法应用于钢筋混凝土结构开 裂分析,有理论价值和应用前景。 同时考虑混凝土 开裂和黏结滑移两种不连续问题的研究较少,钢筋 混凝土扩展有限元可为此提供新的思路,但该方法 尚处于研究阶段,其理论体系有待于进一步完善。 参考文献: [ 1 ] MO魮S N,DOLBOW J,BELYTSCHKO T. A finite element methodforcrackgrowthwithoutremeshing [ J ]. InternationalJournalforNumericalMethodsin Engineering,1999,46(1): 131鄄150. [ 2 ] SL91—2008摇 水工混凝土结构设计规范[S]. [ 3 ] 江见鲸,陆新征. 混凝土结构有限元分析[M]. 2 版. 北 京: 清华大学出版社,2013. [ 4 ] OSHERS, SETHIANJA. Frontspropagatingwith curvature鄄dependentspeed:algorithmsbasedon Hamilton鄄Jacobiformulations[ J ].Journalof Computational Fhysics,1988,79(1): 12鄄49. [ 5 ] SUKUMAR N, CHOPP D L, MO魮S N, et al. Modeling holes and inclusions by level sets in the extended finite鄄 element method [ J ]. Computer MethodsinApplied Mechanics and Engineering,2001,190(46): 6183鄄6200. [ 6 ] VENTURA G,BUDYN E,BELYTSCHKO T. Vector level sets for description of propagating cracks in finite elements [ J]. International Journal for Numerical Methods in Engineering,2003,58: 1571鄄1592. [ 7 ] STOLARSKA M,CHOPP D L,MO魮S N,et al. Modelling crack growth by level sets in the extended finite element method[J]. International Journal for Numerical Methods in Engineering,2001,51(8): 943鄄960. [ 8 ] MELENK J M,BABU譒KA I. The partition of unity finite element method: basic theory and applications [ J ]. ComputerMethodsinAppliedMechanicsand Engineering,1996,139(1/2/3/4): 289鄄314. [ 9 ] CHESSAJ, WANGH, BELYTSCHKOT.Onthe construction of blending elements for local partition of unity enriched finite elements[J]. International Journal for Numerical Methods in Engineering,2003,57(7): 1015鄄 1038. [10] BELYTSCHKO T,BLACK T. Elastic crack growth in finite elements with minimal remeshing[J]. International Journal for Numerical Methods in Engineering,1999,45(5): 601鄄 620. [11] BELYTSCHKO T,CHEN H,XU J,et al. Dynamic crack propagation based on loss of hyperbolicity and a new discontinuous enrichment [ J]. International Journal for Numerical Methods in Engineering,2003,58(12): 1873鄄 1905. [12] FRIES T P. A corrected XFEM approximation without problems in blending elements[J]. International Journal for Numerical Methods in Engineering,2008,75(5): 503鄄 532. [13] 余天堂. 扩展有限元法的数值方面[J]. 岩土力学, 2007,28(增刊 1): 305鄄310. (YU Tiantang. Numerical aspects of the extended finite element method[J]. Rock and Solid Mechanics,2007,28 ( Sup1): 305鄄310. ( in Chinese)) [14] LABORDE P,POMMIER J,RENARD Y,et al. High鄄order extended finite element method for cracked domains[J]. InternationalJournalforNumericalMethodsin Engineering,2005,64(3): 354鄄381. [15] VENTURA G. On the elimination of quadrature subcells for description of propagating cracks in finite elements [ J]. International Journal for Numerical Methods in Engineering,2006,66: 761鄄795. [16] SONG J H,AREIAS P,BELYTSCHKO T. A method for dynamic crack and shear band propagation with phantom nodes[J]. International Journal for Numerical Methods in Engineering,2006,67(6): 868鄄893. [17] 喻葭临,张其光,于玉贞,等. 扩展有限元中非连续区域 的一种积分方案[J]. 清华大学学报:自然科学版,2009 49 (3): 351鄄354. ( YU Jialin, ZHANG Qiguang, YU Yuzhen,et al. An integration scheme for discontinuities in the extended finite element method [ J ]. Journalof Tsinghua University: Science and Technology,2009,49 (3): 351鄄354. (in Chinese)) [18] KARIHALOO B L,XIAO Q Z. Modelling of stationary and growing cracks in FE framework without remeshing: a state鄄of鄄the鄄art review[J