田荣 作品数:12 被引量:32 H指数:3 供职机构: 北京应用物理与计算数学研究所 更多>> 发文基金: 国家自然科学基金 国家自然科学基金委员会-广东省人民政府联合基金 更多>> 相关领域: 理学 一般工业技术 自动化与计算机技术 电子电信 更多>>
改进型XFEM进展 被引量:9 2016年 扩展有限元法(XFEM)在诞生后的十几年时间里,引起学术界和工业界的广泛关注,并已经成为目前裂纹分析的主流数值方法。然而,在实际应用中该法一直受到两方面的困扰,(1)总体方程高度病态;刚度阵条件数随网格尺寸呈h-6变化(普通有限元为h-2)。(2)裂尖强化插值由于能量一致性问题无法直接推广应用于动力学计算。前者表现在XFEM稳态问题的迭代求解收敛慢或难以收敛,后者长期以来导致XFEM裂纹扩展动力学计算实施困难。本文认为XFEM目前遇到的种种困难,均与单位分解引入的额外自由度相关。为此,提出了无额外自由度的单位分解插值格式,基于此格式,进一步构造出改进型扩展有限元方法。改进型XFEM具有如下特点,(1)可以消除原有XFEM的线性依赖性和总体方程病态的问题。(2)避免动力学问题中额外自由度引起的质量集中、零临界时间步长问题以及裂纹扩展过程中的能量一致性问题。本文结合静动力学测试问题综述上述改进。 田荣 文龙飞关键词:XFEM 砌体结构长脉宽爆炸荷载损伤等级评估方法 被引量:1 2021年 随着百千吨级当量爆炸工业事故的频繁发生,建筑结构的损伤评估和抗爆安全性更受到关注。目前,构件级的评估方法相对成熟,而大当量冲击波作用下的建筑结构整体毁伤评估依旧是个开放性问题。本文中,面向结构级的毁伤评估,提出了一种新的评估方法−构件损伤加权。该方法以承重构件损伤程度为基础,通过基于应变能的构件权系数加权,进而评估结构级损伤破坏程度。为了验证评估方法的有效性,以典型砌体结构为例,利用自主研发的冲击波结构毁伤模拟有限元程序,开展了百毫秒脉宽爆炸冲击波荷载下结构动力学响应数值模拟。根据数值模拟结果,结合构件损伤加权的评估方法,获取砌体结构损伤等级与冲击波超压的关系。预测的超压值的相对误差为−16.9%~26.2%,验证了评估方法的有效性。该评估方法为获取砌体结构的超压-冲量曲线提供了可行的途径,可为结构的抗爆安全设计提供参考。 曾繁 肖桂仲 冯晓伟 黄超 田荣关键词:砌体结构 超压 钢筋混凝土结构改进型分离式数值模型 2022年 提出一种考虑粘结滑移效应的钢筋混凝土改进型分离式数值模型。在混合物理论基础上,该模型兼顾混凝土基体和钢筋的力学行为,且基于钢筋混凝土界面粘结滑移模型,获得了钢筋等效模型。改进型分离式数值模型由于对钢筋及其界面无显式离散要求,使得钢筋的运用完全独立于其几何形状,同时对混凝土网格没有约束,并且不增加计算成本,因此该模型可适用于钢筋混凝土宏观结构层面分析。通过钢筋混凝土构件-结构的爆炸实验,对改进型分离式数值模型进行层次化验证。对比结果表明,考虑粘结滑移效应的有限元模型能够更好地预测钢筋混凝土结构的宏观力学行为。 曾繁 冯晓伟 黄超 徐权 徐权 田荣关键词:混凝土 钢筋 粘结-滑移 爆炸荷载 非连续问题中单元分割的模板方法 被引量:2 2021年 扩展有限元法(extended finite element method, XFEM)因具有裂纹几何独立于模拟网格、裂纹扩展时无需网格重分重映、计算精度高等优点,成为裂纹分析的主流数值方法之一.但该方法在工程实践中存在单元被裂纹分割的几何困难——现有精确几何分割方法实现复杂、计算量大、鲁棒性差.为克服这一困难,本文提出一种基于单元水平集的模板分割方法,用于非连续单元子剖分和数值积分.首先,遍历单元水平集值所有形态并建立标准单元分割模板库;然后,根据单元水平集值,对非标准单元进行形态查询和模板插值;最后,套用标准单元分割模板实现单元高效分割和子剖分.将该方法与常规XFEM、改进型XFEM进行结合,从而应用于孔洞、夹杂、裂纹等非连续问题分析中.算例分析表明,本文提出的模板分割方法具有较高计算精度.由于不引入复杂几何操作,该模板分割方法同时具有较高计算效率和鲁棒性,故可为XFEM类方法在实际工程应用中提供有效支撑. 王理想 文龙飞 肖桂仲 田荣关键词:水平集 隐格式近似不可压四面体4节点大变形单元 2017年 由于在处理体积自锁方面的优势,近似不可压问题的大变形求解多采用六面体单元/网格,但对于复杂工程问题,由于网格剖分上的限制,往往更需要一种可以很好解决体积自锁的四面体单元。Bonet和Burton的平均节点压力4节点四面体单元是为数不多能够较好处理体积自锁问题的四面体单元之一,但是该单元目前主要用于显式计算。利用单元平均压力对位移增量的精确方向导数,得到了严格的一致切线阵,保证了Newton-Raphson迭代的二阶收敛,从而使得该单元可以用于隐式计算。该单元的压力平均计算会耦合相邻单元的节点自由度,从而增加切线刚度阵的非零带宽,但不增加自由度总数。分别采用线性六面体选择缩减积分单元、标准线性四面体单元和本文的单元计算了3个近似不可压的典型算例。算例表明,本文推导的单元可以有效克服体积自锁,达到与常用六面体单元相近的效果,使得四面体网格可以方便地用于不可压问题的大变形隐式求解。 肖桂仲 陈成军 肖世富 田荣关键词:大应变 一种数理方程解算构件及并行软件研发方法与系统 本文涉及高性能数值模拟技术领域,尤其涉及一种数理方程解算构件及并行软件研发方法与系统。其中数理方程解算构件基于并行计算构件研发,并通过数理方程解算构件接口将物理建模插件单元,工程模板插件单元,数值离散解算单元三类单元隔离... 莫则尧 张爱清 王景焘 程汤培 刘娜 李刚 陈晓洁 李光荣 田荣 刘青凯 许开龙文献传递 一种数理方程解算构件及并行软件研发方法与系统 本文涉及高性能数值模拟技术领域,尤其涉及一种数理方程解算构件及并行软件研发方法与系统。其中数理方程解算构件基于并行计算构件研发,并通过数理方程解算构件接口将物理建模插件单元,工程模板插件单元,数值离散解算单元三类单元隔离... 莫则尧 张爱清 王景焘 程汤培 刘娜 李刚 陈晓洁 李光荣 田荣 刘青凯 许开龙文献传递 基于改进型XFEM的裂纹分析并行软件实现 被引量:5 2018年 扩展有限元法(XFEM)目前已成为裂纹分析的主流数值方法.然而,在实际应用中该方法一直受到两方面的困扰:总体刚度矩阵高度病态以及动力学计算时额外自由度上的能量无法正确传递.前者导致在以迭代法求解为主的大规模计算中难以收敛,后者导致动力学求解精度低、实施困难.为解决这两大困扰,基于无额外自由度单位分解插值格式,提出改进型扩展有限元法(IXFEM).基于该方法,在高性能数值模拟软件开发框架JAUMIN上开发裂纹分析并行软件—PANDA_Fracture.实际算例表明, PANDA_Fracture软件具有如下特色:1)允许引入高精度裂尖加强,支持裂纹任意长度扩展; 2)同时支持动力学裂纹问题隐式求解与显式求解; 3)支持任意多裂纹问题高效求解; 4)支持上千CPU核高性能可扩展计算. 王理想 文龙飞 王景焘 田荣关键词:PANDA 动载下裂纹应力强度因子计算的改进型扩展有限元法 被引量:19 2018年 相较于常规扩展有限元法(extended finite element method,XFEM),改进型扩展有限元法(improved XFEM)解决了现有方法线性相关与总体刚度矩阵高度病态问题,在数量级上提升了总体方程的求解效率,克服了现有方法在动力学问题中的能量正确传递、动态应力强度因子数值震荡、精度低下问题.本文基于改进型XFEM,采用Newmark隐式时间积分算法,重点研究了动载荷作用下扩展裂纹尖端应力强度因子的求解方法,与静力学方法相比,增加了裂纹扩展速度项与惯性项的贡献.通过数值算例研究了网格单元尺寸、质量矩阵、时间步长、裂尖加强区域、惯性项、扩展速度项及相互作用积分区域J-domain的网格与单元尺寸对动态应力强度因子求解精度的影响,验证了改进型XFEM计算动态裂纹应力强度因子方法的有效性.针对文献中具有挑战性的"I型半无限长裂纹先稳定后扩展"问题,改进型XFEM给出目前为止精度最好的动态应力强度因子数值解. 文龙飞 王理想 田荣关键词:动态应力强度因子 线弹性断裂力学 改进型XFEM 扩展有限元法(XFEM)在诞生后的十几年时间里,引起学术界和工业界的广泛关注、成为目前裂纹分析的主流数值方法。然而,在实际应用中该法也同时一直受到两方面算法上的困扰:(1)总体方程高度病态;刚度阵条件数随网格尺寸呈h变化... 田荣关键词:XFEM 文献传递