地下洞室三维弹塑性动力显式有限元锚杆模拟

摘要

锚杆的数值模拟方法是分析地震灾变中锚杆的作用机理及锚杆支护对地下洞室围岩抗震的加固效果的关键。基于地震动荷载作用下锚杆围岩联合受力机理,在对FLAC3D锚杆模拟算法进行一定改进的基础上,提出了一种适用于抗震分析的全长粘结式锚杆模拟方法,并将该方法应用到自主开发的地下洞室三维显式动力有限元程序中,取得较好效果。该方法继承了FLAC3D中对锚杆、锚杆与砂浆接触面、锚固砂浆、砂浆与围岩接触面及围岩本身共5种结构材料的受力破坏情况的详细模拟,但不再进行相邻锚杆计算段之间的受力平衡迭代,而是将各锚杆计算段的受力转移到单元结点上,将锚杆求解与单元求解相统一,从而降低动力求解中锚杆的计算耗时。同时,详细推导了锚杆各结构材料的屈服极限表达式,缩短了锚杆结构材料屈服判断及应力修正耗时。该锚杆模型与FLAC3D中的模型假定不同,认为锚杆变形由围岩变形引起,各锚杆计算段间受力相对独立,受力并不连续,与地震中围岩锚杆的受力情况较吻合。采用该方法对某地下厂房洞室群进行地震动力时程分析,表明该模拟方法提高了求解速度,充分反映了锚杆对围岩的支护效果。