时间相关的连续非连续数法非全局接触识别方法及系统技术方案

技术编号:41527336 阅读:27 留言:0更新日期:2024-06-03 23:02
本发明专利技术公开了一种时间相关的连续非连续数法非全局接触识别方法及系统,所述方法包括以下步骤:创建几何体数值模型,施加边界条件,计算单元力,对处于边界的实体单元和离散断裂网络单元列表中可能产生接触的单元进行识别,存储并更新潜在接触式实体单元列表,然后在潜在接触式实体单元列表中进行接触检测,更新各单元状态并重复附图中S300‑S600的步骤,直至模拟完成。本发明专利技术提供的一种时间相关的连续非连续数法非全局接触识别方法,可实时反映潜在接触对随裂纹扩展的逐步增加,能有效继承失效节理单元的单元力,维持数值模型的数值稳定,具有充分的可并行能力,对提高连续非连续数值模型运算效率和提高模拟精度都具有重要意义。

【技术实现步骤摘要】

本专利技术属于岩石力学和岩石工程领域,更具体地,涉及一种时间相关的连续非连续数法非全局接触识别方法及系统


技术介绍

1、在实际工程应用中,通常需要根据具体地质条件和构造环境设计合理的施工方案,并通过现场试验对施工方案进行调整优化。然而,受试验条件限制和实验成本的影响,很多情况很难直接开展大尺度或等尺度的物理模型实验,如采矿、基坑开挖及隧道开挖等。数值模拟可以通过数值计算的方法对现场试验过程进行仿真,减少成本,还能模拟物理试验无法实现的边界条件,有着重要的参考价值。结合了连续法和非连续法的优势,连续非连续数值方法能够真实模拟岩体结构从连续向不连续的演变,描述岩体的真实断裂、破碎过程,近年来一直是岩土工程计算方法的研究热点和重点。

2、当前大部分连续非连续数法(fdem)软件的接触检测部分是通过空间哈希接触检测方法完成的,例如cundall的fcp方法、munjiza的nbs方法以及刘宏元的潜在接触对方法,空间哈希接触检测方法会对离散体周围一定区域(圆形或方形)进行均匀的网格划分,将该区域内包括的每个离散体都作为潜在的接触对,并根据网格计算和比较本文档来自技高网...

【技术保护点】

1.一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,包括以下步骤:

2.根据权利要求1中所述的一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,步骤S300中所述数值模拟计算结果包括接触力、节点力、节点位移、节点速度、实体单元应力、实体单元应变、节理单元作用力、节理单元张开度、节理单元滑移位移以及节理单元状态。

3.根据权利要求1或2中所述的一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,步骤S400中根据失效的节理单元将共享该节理单元的实体单元对存储于潜在接触式实体单元列表中,包括:

4.根据权利要求3所述的一种时间...

【技术特征摘要】

1.一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,包括以下步骤:

2.根据权利要求1中所述的一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,步骤s300中所述数值模拟计算结果包括接触力、节点力、节点位移、节点速度、实体单元应力、实体单元应变、节理单元作用力、节理单元张开度、节理单元滑移位移以及节理单元状态。

3.根据权利要求1或2中所述的一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,步骤s400中根据失效的节理单元将共享该节理单元的实体单元对存储于潜在接触式实体单元列表中,包括:

4.根据权利要求3所述的一种时间相关的连续非连续数法非全局接触识别方法,其特征在于,所述s402步骤中,实体单元继承失效的节理单元的单元内力的计算方法为:

5.根据...

【专利技术属性】
技术研发人员:李晓锋李潜李海波傅帅旸
申请(专利权)人:中国科学院武汉岩土力学研究所
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1