相关文章:廉艳平, 张雄, 周旭, 马志涛. A FE-MP method and its application in modeling dynamic response of reinforced concrete subjected to impact loading. Computer Methods in Applied Mechanics and Engineering, 200(17-20): 1659-1670, 2011. (DOI: 10.1016/j.cma. 2011.01.019).
钢筋混凝土中钢筋的直径远远小于混凝土的尺寸,这种尺寸上的差异给离散带来了很大的困难。若将钢筋和混凝土采用同样的质点间距离散,则需要大量的物质点,计算量将不可接受;不考虑钢筋存在又无法准确把握钢筋混凝土的特性。为此,我们提出了一种新的杂交物质点有限元法。我们将有限元的索单元引入到物质点法中,构造了一种杂交杆单元。在每个时间步中,单元的节点如质点一样,其动量方程在背景网格上积分,但其本构方程在单元的高斯点上进行计算。基于该方法,我们用物质点离散素混凝土,杂交杆单元离散混凝土中的钢筋,如图1所示。通过背景网格,杆单元代表的钢筋和物质点代表的混凝土实现了相互作用,从而较好的模拟了钢筋混凝土材料。
图 1 钢筋混凝土离散示意图。其中混凝土离散成物质点(用圆圈表示),钢筋离散成有限杆单元(实心点及其连线表示)。
我们用该方法模拟了钢筋混凝土侵彻实验,实验中卵形弹以不同的初始速度正侵彻钢筋混凝土靶。我们分别考虑了弹体未打中钢筋和弹体打中钢筋两种情况,所得相应的剩余弹速如图2所示。在未打中钢筋的情况下,钢筋对剩余弹速的影响较小,与实验结果吻合较好;在打中钢筋的情况下,钢筋对剩余弹速影响显著,降低了剩余弹速。图3给出了弹体打中钢筋的情况下,靶体和钢筋的变形损伤图。
基于该思想,可以构造其它类型的有限单元如薄膜单元,以解决由于几何尺寸的巨大差异带来的离散困难的问题。
图2 不同初始速度下卵形弹的剩余速度实验值与模拟结果比较
图3 弹体打中钢筋工况下的靶体损伤图和钢筋的变形图
创建: Apr 09, 2018 | 16:20
相关文章:廉艳平,张雄,刘岩. Coupling of finite element method with material point method by local multi-mesh contact method. Computer Methods in Applied Mechanics and Engineering 200: 3482–3494, 2011. (DOI: 10.1016/j.cma.2011.07.014)
物质点法在模拟材料的大变形甚至破碎现象方面有着天然的优势,但是在模拟变形较小的材料时,其精度和效率不如显式有限元方法。如果将有限元和物质点法耦合起来模拟冲击侵彻和流固耦合类等问题,无疑将提高求解结果的精度和效率。为解决显式有限元和物质点法的耦合问题,我们提出了一种基于局部多重背景网格的物质点有限元耦合算法(CFEMP)。在该算法中,变形较小的物体采用有限元方法模拟,变形较大的物体采用物质点法模拟。两者之间的相互作用通过接触算法完成。如图1所示,物体r采用有限元模拟,物体s采用物质点模拟。在接触面上,将有限元材料区域的a、b、c结点的物理量映射到背景网格上。采用基于局部多重背景网格的接触方法,计算出接触力并将力作用两个物体上。
图1 物质点有限元耦合示意图
通过两个物体的对撞、斜坡上球体的滚动问题解析解与模拟结果的对比验证了该算法的效率和精度均高于物质点法。然后,我们利用该方法模拟了卵形钢弹斜侵彻铝靶的问题,其中弹体用有限元模拟,靶体用物质点法模拟。弹体的剩余弹速模拟结果与实验结果吻合较好,如表1所示,弹体的变形过程同实验结果相吻合,如图所示2。
表1:不同初速下弹体剩余弹速实验值同模拟结果的对比
(a)弹体初速为 400m/s的工况 (b) 弹体初始为575m/s的工况
图3 不同时刻弹体变形图实验图片同模拟结果的对比
我们还用该方法模拟了一流固耦合问题:水柱在重力作用下的垮塌并冲击弹性障碍物,其中水柱用物质点法模拟,障碍物用有限元模拟。模拟的结果同国际上其他的学者采用的其它方法模拟的结果一致,如图3所示。
(a)水柱垮塌问题尺寸示意图 (b) 障碍物左上角点挠度时间历程曲线
图 3 水柱垮塌并冲击弹性障碍物问题模型尺寸及模拟结果
创建: Apr 09, 2018 | 16:24