CN118153468B 一种用于桥梁冲刷sph多相流模型的单层边界仿真方法 （东南大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN118153468B(45)授权公告日2025.07.08

(21)申请号202311848820.1

(22)申请日2023.12.29

(65)同一申请的已公布的文献号申请公布号CN118153468A

(43)申请公布日2024.06.07

(73)专利权人东南大学

地址211189江苏省南京市江宁区东南大

学路2号

(72)发明人熊文张嵘钊

(74)专利代理机构南京苏高专利商标事务所(普通合伙)32204

专利代理师王安琪

(51)Int.CI.

GO6F30/28(2020.01)

GO6F30/25(2020.01)

GO6F30/13(2020.01)

GO6F113/08(2020.01)

G06F119/14(2020.01)

GO6F111/10(2020.01)

(56)对比文件

CN115859432A,2023.03.28

朱晓临；汪欢欢；殷竞存.SPH流固耦合模拟自适应采样固壁虚粒子边界处理方法.图学学报.2018,(03),全文.

审查员洪必文

权利要求书3页说明书6页附图2页

(54)发明名称

一种用于桥梁冲刷SPH多相流模型的单层边界仿真方法

(57)摘要

CN118153468B本发明公开了一种用于桥梁冲刷SPH多相流模型的单层边界仿真方法，包括如下步骤：建立SPH多相流模型，建立冲刷墩单层粒子边界模型；根据上一步骤所得冲刷墩单层粒子边界，构建幽灵粒子映射边界；基于上一步骤所得映射边界，利用核函数完成所有幽灵粒子周围流体粒子检索；基于上一步骤所得幽灵粒子及对应检索粒子，逆向得出冲刷墩边界粒子密度；根据上一步骤所得更新密度后的冲刷墩边界粒子，对于任一冲刷墩边界粒子b检索半径内的流体粒子i,根据相对距离施加人工斥力；根据上一步骤所得人工斥力，计算流体粒子所受冲刷墩边界合力。本发明能够实现桥梁冲刷SPH

CN118153468B

CN118153468B权利要求书1/3页

2

1.一种用于桥梁冲刷SPH多相流模型的单层边界仿真方法，其特征在于，包括如下步

骤：

步骤1、建立SPH多相流模型，确定流体粒子参数，设定粒子间距，确立核函数形式，建立目标边界单层粒子模型；具体包括如下步骤：

步骤11、设定粒子初始间距d,确定SPH算法采用的核函数形式，采用Wendland方程：

其中，W为粒子j相对于粒子i的核函数值，α为维度系数，取21/16πh3,h为核函数搜索

半径，取值2d,,r;j=r;-r;,r;为粒子i的坐标；

步骤12、确定低水粒子的密度p,质量m,确定泥沙粒子的密度p,质量m,确定粒子体积V;流体粒子质量及体积计算如下：

步骤13、根据粒子初始间距d划分冲刷墩单层粒子边界；

步骤2、根据步骤1所得冲刷墩单层粒子边界，构建幽灵粒子映射边界；

步骤3、基于步骤2所得映射边界，循环构建所有冲刷墩边界粒子对应的映射幽灵粒子，利用核函数完成所有幽灵粒子周围流体粒子检索；

步骤4、基于步骤3所得幽灵粒子及对应检索粒子，引入等效密度法更新幽灵粒子密度，并逆向得出冲刷墩边界粒子密度；

步骤5、根据步骤4所得更新密度后的冲刷墩边界粒子，对于任一冲刷墩边界粒子b检索半径内的流体粒子i,根据相对距离r,施加人工斥力F,以避免流体粒子出现穿透现象；

步骤6、根据步骤5所得人工斥力F,结合粒子梯度压力及冲刷墩边界粒子粘性系数，计

算流体粒子i所受冲刷墩边界合力Fb。

2.如权利要求1所述的用于桥梁冲刷SPH多相流模型的单层边界仿真方法，其特征在

于，步骤2中，幽灵粒子映射边界应构建在边界平面法向距离d/2处，边界法向向量应指向计算域中流体方向。

3.如权利要求1所述的用于桥梁冲刷SPH多相流模型的单层边界仿真方法，其特征在于，步骤3中，幽灵粒子g周围检索所得流体粒子i应满足：

rgil≤h。

4.如权利要求1所述的用于桥梁冲刷SPH多相流模型的单层边界仿真方法，其特征在于，步骤4中，基于步骤3所得幽灵粒子及对应检索粒子，引入等效密度法更新幽灵粒子密度，并逆向得出冲刷墩边界粒子密度具体包括如下步骤：

步骤41、引入等效密度法，对每一个幽灵粒子检索到的泥沙粒子进行质量折算，得到等效质量m;

步骤42、根