一种基于面-面局部搜索的接触算法

接触-碰撞算法是影响动态接触问题数值分析精度的重要因素,发展健壮、精确的局部搜索算法对提高数值模拟精度具有重要意义。为解决点-面接触搜索算法存在的盲区问题,提出了一种基于面-面的局部搜索算法,并发展了相应的接触力计算方法。新的接触算法以面心坐标与特征长度表征接触片,进行预搜索快速排除不会发生接触的潜在接触对;利用主从接触片的投影与侵入深度关系确定接触片间的接触状态;利用等参逆变换解析给出接触点,避免了方程组的迭代求解。典型算例与工程实例计算表明,本文算法消除了接触搜索的盲区,具有很好的健壮性与计算精度。     

非连续问题中单元分割的模板方法

扩展有限元法 (extended finite element method, XFEM) 因具有裂纹几何独立于模拟网格、裂纹扩展时无需网格重分重映、计算精度高等优点,成为裂纹分析的主流数值方法之一. 但该方法在工程实践中存在单元被裂纹分割的几何困难 —— 现有精确几何分割方法实现复杂、计算量大、鲁棒性差. 为克服这一困难, 本文提出一种基于单元水平集的模板分割方法, 用于非连续单元子剖分和数值积分. 首先, 遍历单元水平集值所有形态并建立标准单元分割模板库; 然后, 根据单元水平集值, 对非标准单元进行形态查询和模板插值; 最后, 套用标准单元分割模板实现单元高效分割和子剖分. 将该方法与常规XFEM、改进型XFEM进行结合,从而应用于孔洞、夹杂、裂纹等非连续问题分析中. 算例分析表明, 本文提出的模板分割方法具有较高计算精度. 由于不引入复杂几何操作, 该模板分割方法同时具有较高计算效率和鲁棒性, 故可为XFEM类方法在实际工程应用中提供有效支撑.      

钢筋混凝土结构改进型分离式数值模型

提出一种考虑粘结滑移效应的钢筋混凝土改进型分离式数值模型。在混合物理论基础上,该模型兼顾混凝土基体和钢筋的力学行为,且基于钢筋混凝土界面粘结滑移模型,获得了钢筋等效模型。改进型分离式数值模型由于对钢筋及其界面无显式离散要求,使得钢筋的运用完全独立于其几何形状,同时对混凝土网格没有约束,并且不增加计算成本,因此该模型可适用于钢筋混凝土宏观结构层面分析。通过钢筋混凝土构件-结构的爆炸实验,对改进型分离式数值模型进行层次化验证。对比结果表明,考虑粘结滑移效应的有限元模型能够更好地预测钢筋混凝土结构的宏观力学行为。     

隐格式近似不可压四面体4节点大变形单元

由于在处理体积自锁方面的优势,近似不可压问题的大变形求解多采用六面体单元/网格,但对于复杂工程问题,由于网格剖分上的限制,往往更需要一种可以很好解决体积自锁的四面体单元。Bonet和Burton的平均节点压力4节点四面体单元是为数不多能够较好处理体积自锁问题的四面体单元之一,但是该单元目前主要用于显式计算。利用单元平均压力对位移增量的精确方向导数,得到了严格的一致切线阵,保证了Newton-Raphson迭代的二阶收敛,从而使得该单元可以用于隐式计算。该单元的压力平均计算会耦合相邻单元的节点自由度,从而增加切线刚度阵的非零带宽,但不增加自由度总数。分别采用线性六面体选择缩减积分单元、标准线性四面体单元和本文的单元计算了3个近似不可压的典型算例。算例表明,本文推导的单元可以有效克服体积自锁,达到与常用六面体单元相近的效果,使得四面体网格可以方便地用于不可压问题的大变形隐式求解。     

砌体结构长脉宽爆炸荷载损伤等级评估方法

随着百千吨级当量爆炸工业事故的频繁发生,建筑结构的损伤评估和抗爆安全性更受到关注。目前,构件级的评估方法相对成熟,而大当量冲击波作用下的建筑结构整体毁伤评估依旧是个开放性问题。本文中,面向结构级的毁伤评估,提出了一种新的评估方法?构件损伤加权。该方法以承重构件损伤程度为基础,通过基于应变能的构件权系数加权,进而评估结构级损伤破坏程度。为了验证评估方法的有效性,以典型砌体结构为例,利用自主研发的冲击波结构毁伤模拟有限元程序,开展了百毫秒脉宽爆炸冲击波荷载下结构动力学响应数值模拟。根据数值模拟结果,结合构件损伤加权的评估方法,获取砌体结构损伤等级与冲击波超压的关系。预测的超压值的相对误差为?16.9%～26.2%,验证了评估方法的有效性。该评估方法为获取砌体结构的超压-冲量曲线提供了可行的途径,可为结构的抗爆安全设计提供参考。     

一种基于面-面局部搜索的接触算法

接触-碰撞算法是影响动态接触问题数值分析精度的重要因素,发展健壮、精确的局部搜索算法对提高数值模拟精度具有重要意义。为解决点-面接触搜索算法存在的盲区问题,提出了一种基于面-面的局部搜索算法,并发展了相应的接触力计算方法。新的接触算法以面心坐标与特征长度表征接触片,进行预搜索快速排除不会发生接触的潜在接触对;利用主从接触片的投影与侵入深度关系确定接触片间的接触状态;利用等参逆变换解析给出接触点,避免了方程组的迭代求解。典型算例与工程实例计算表明,本文算法消除了接触搜索的盲区,具有很好的健壮性与计算精度。