考虑动态流体网格的颗粒-流体耦合算法

doi:10.16183/j.cnki.jsjtu.2020.317

上海交通大学学报 ›› 2021, Vol. 55 ›› Issue (6): 645-651.doi: 10.16183/j.cnki.jsjtu.2020.317

所属专题：《上海交通大学学报》2021年“土木建筑工程”专题；《上海交通大学学报》2021年12期专题汇总专辑

考虑动态流体网格的颗粒-流体耦合算法

何金辉, 李明广(), 陈锦剑, 夏小和

上海交通大学土木工程系, 上海 200240

收稿日期:2020-09-30 出版日期:2021-06-28 发布日期:2021-06-30
通讯作者: 李明广 E-mail:lmg20066028@sjtu.edu.cn
作者简介:何金辉(1995-),男,河南省洛阳市人,硕士生,从事岩土流固耦合方面研究
基金资助:
国家自然科学基金面上项目(41977216);国家自然科学基金面上项目(51978399);上海市科委项目(BI0100054)

Particle-Fluid Coupling Algorithm Considering Dynamic Fluid Mesh

HE Jinhui, LI Mingguang(), CHEN Jinjian, XIA Xiaohe

Department of Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

Received:2020-09-30 Online:2021-06-28 Published:2021-06-30
Contact: LI Mingguang E-mail:lmg20066028@sjtu.edu.cn

摘要/Abstract

摘要：

传统颗粒-流体耦合算法难以考虑流体动态边界问题,在模拟大变形案例时因流固边界不匹配易造成计算误差,影响结果的准确性.针对该问题,引入流体网格动态更新方法,推导动网格下的达西渗流方程和颗粒-流体相互作用方程,在离散元商业软件PFC^2D基础上,开发了考虑流体动网格的颗粒-流体耦合算法.将该算法用于模拟饱和土不排水剪切双轴试验,通过与常体积法的结果对比,验证了该方法的有效性.最后,采用该算法模拟了不同围压下的不排水双轴试验,计算结果的规律与室内试验具有较好的一致性.该算法考虑了流体动网格的问题,在模拟大变形下的三轴压缩试验、一维固结试验等案例时可获得相适应的流固边界,有助于提高模拟精度,因此可为类似研究提供参考.

关键词: 流固耦合, 动网格, 常体积法, 双轴试验

Abstract:

It is generally difficult to consider the fluid dynamic boundary problem in the traditional particle-fluid coupling algorithm, causing calculation errors owing to the mismatching of the fluid-solid boundary and affecting the accuracy of the results in the modeling of large deformation issues. In view of this problem, the dynamic updating method of fluid mesh is introduced and Darcy’s seepage equation and the particle-fluid interaction equation in dynamic mesh are derived. Based on the discrete element commercial software PFC^2D, the particle-fluid coupling algorithm considering dynamic fluid mesh is developed. The proposed algorithm is applied to simulate the undrained shear biaxial test of saturated soil. The comparison result with the constant volume method verifies the effectiveness of the developed algorithm. Finally, the algorithm is used to simulate the undrained biaxial tests at different confining pressures. The law of the calculated results agrees well with that of the laboratory tests. By considering the problem of fluid dynamic boundary, the developed algorithm can obtain the fluid-solid boundary matching in the simulation of triaxial compression test and one-dimensional consolidation test or other cases in large deformations, which can help to improve the simulation accuracy and offer a theoretical reference for similar studies.

Key words: fluid-solid coupling, dynamic mesh, constant volume method, biaxial test

中图分类号:

TU411

何金辉, 李明广, 陈锦剑, 夏小和. 考虑动态流体网格的颗粒-流体耦合算法[J]. 上海交通大学学报, 2021, 55(6): 645-651.

HE Jinhui, LI Mingguang, CHEN Jinjian, XIA Xiaohe. Particle-Fluid Coupling Algorithm Considering Dynamic Fluid Mesh[J]. Journal of Shanghai Jiao Tong University, 2021, 55(6): 645-651.

图/表 11

图1

图2

图3

图4

表1

图5

图6

图7

图8

图9

图10

参考文献 15

[1]	蒋明镜, 张望城. 一种考虑流体状态方程的土体CFD-DEM耦合数值方法[J]. 岩土工程学报, 2014, 36(5):793-801.
	JIANG Mingjing, ZHANG Wangcheng. Coupled CFD-DEM method for soils incorporating equation of state for liquid[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(5):793-801.
[2]	刘广, 荣冠, 侯迪, 等. 颗粒-流体耦合算法与饱和土不排水剪切特性研究[J]. 岩土力学, 2016, 37(1):210-218.
	LIU Guang, RONG Guan, HOU Di, et al. Fluid-particle coupled model and a numerical investigation on undrained shear behavior of saturated soil[J]. Rock and Soil Mechanics, 2016, 37(1):210-218.
[3]	ZHANG Z L, SHAO Z S, FANG X B, et al. Research on the fracture grouting mechanism and PFC numerical simulation in loess[J]. Advances in Materials Science and Engineering, 2018, 2018:1-7.
[4]	ZENG Y W, XIA L. Numerical simulation of hydraulic fracturing in transversely isotropic rock masses based on PFC-2D[J]. Journal of Vibroengineering, 2019, 21(4):833-847. doi: 10.21595/jve.2018.19962 URL
[5]	WANG Y, DONG Q, CHEN Y. Seepage simulation using pipe network flow model in a discrete element system[J]. Computers and Geotechnics, 2017, 92:201-209. doi: 10.1016/j.compgeo.2017.08.010 URL
[6]	王学红. 高频液压振动锤沉桩机理的颗粒离散元模拟[D]. 福州: 福州大学, 2011.
	WANG Xuehong. Discrete element modelling of mechanisms of high frequency hydraulic vibratory pile driving[D]. Fuzhou: Fuzhou University, 2011.
[7]	ZEGHAL M, EL SHAMY U. A continuum-discrete hydromechanical analysis of granular deposit liquefaction[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2004, 28(14):1361-1383. doi: 10.1002/(ISSN)1096-9853 URL
[8]	TAO H, TAO J. Quantitative analysis of piping erosion micro-mechanisms with coupled CFD and DEM method[J]. Acta Geotechnica, 2017, 12(3):573-592. doi: 10.1007/s11440-016-0516-y URL
[9]	ZHANG A, JIANG M, THORNTON C. A coupled CFD-DEM method with moving mesh for simulating undrained triaxial tests on granular soils[J]. Granular Matter, 2020, 22(1):1-13. doi: 10.1007/s10035-019-0969-4 URL
[10]	ZHANG F, WANG T, LIU F, et al. Modeling of fluid-particle interaction by coupling the discrete element method with a dynamic fluid mesh: Implications to suffusion in gap-graded soils[J]. Computers and Geotechnics, 2020, 124:103617. doi: 10.1016/j.compgeo.2020.103617 URL
[11]	SHAFIPOUR R, SOROUSH A. Fluid coupled-DEM modelling of undrained behavior of granular media[J]. Computers and Geotechnics, 2008, 35(5):673-685. doi: 10.1016/j.compgeo.2007.12.003 URL
[12]	LIU G, RONG G, PENG J, et al. Numerical simulation on undrained triaxial behavior of saturated soil by a fluid coupled-DEM model[J]. Engineering Geology, 2015, 193:256-266. doi: 10.1016/j.enggeo.2015.04.019 URL
[13]	KHALILI Y, MAHBOUBI A. Discrete simulation and micromechanical analysis of two-dimensional sa-turated granular media[J]. Particuology, 2014, 15:138-150. doi: 10.1016/j.partic.2013.07.005 URL
[14]	OKADA Y, OCHIAI H. Coupling pore-water pre-ssure with distinct element method and steady state strengths in numerical triaxial compression tests under undrained conditions[J]. Landslides, 2007, 4(4):357-369. doi: 10.1007/s10346-007-0092-1 URL
[15]	KEISHING J, HANLEY K J. Improving constant-volume simulations of undrained behaviour in DEM[J]. Acta Geotechnica, 2020, 15(9):2545-2558. doi: 10.1007/s11440-020-00949-1 URL

颗粒参数	取值	颗粒参数	取值
颗粒密度/(kg·m^-3)	2500	切向刚度/(MN·m^-1)	100
最小颗粒直径/mm	3.6	摩擦因数	0.5
最大颗粒直径/mm	5.58	渗透性系数/(μm·s^-1)	10
法向刚度/(MN·m^-1)	100	试样初始孔隙率	0.16

考虑动态流体网格的颗粒-流体耦合算法

Particle-Fluid Coupling Algorithm Considering Dynamic Fluid Mesh

RichHTML

PDF (PC)

可视化

摘要/Abstract

引用本文

使用本文

图/表 11

参考文献 15

相关文章 15

编辑推荐

Metrics

本文评价

[1]	吴利红, 封锡盛, 叶作霖, 李一平. 自主水下机器人强制自航下潜的类物理模拟[J]. 上海交通大学学报, 2021, 55(3): 290-296.
[2]	张宇, 王晓亮. 基于显式动力学的软式飞艇流固耦合计算框架[J]. 上海交通大学学报, 2021, 55(3): 311-319.
[3]	田新亮. “软尾减阻”述评[J]. 上海交通大学学报, 2021, 55(2): 213-214.
[4]	郑启宇，夏小和，李明广，张扬清. 深基坑降承压水对墙体变形和地表沉降的影响[J]. 上海交通大学学报, 2020, 54(10): 1094-1100.
[5]	于鹏垚，赵勇，王天霖，甄春博，苏绍娟. 基于半解析砰击模型的弹性楔形体入水冲击分析[J]. 上海交通大学学报（自然版）, 2019, 53(2): 179-187.
[6]	刘广, 尚宇晴, 许泉, 俞刘建, 樊浩. 单侧提拉发射系统流固耦合动力学仿真[J]. 空天防御, 2019, 2(1): 14-21.
[7]	秦浩，唐文勇，薛鸿祥. 非线性畸形波所致的平台底部砰击载荷及结构响应数值模拟[J]. 上海交通大学学报（自然版）, 2018, 52(9): 1009-1016.
[8]	郭泽宇，陈作钢. 基于复杂柔性运动流场数值求解的动网格方案及应用[J]. 上海交通大学学报（自然版）, 2018, 52(4): 403-409.
[9]	单铁兵，沈志平. 深水半潜式平台的涡激运动I.关键特性研究[J]. 上海交通大学学报（自然版）, 2017, 51(12): 1493-1503.
[10]	王小庆a,b，金先龙a,b，曹源a. 大规模输水隧道水锤效应三维数值模拟[J]. 上海交通大学学报（自然版）, 2016, 50(01): 98-102.
[11]	刘晗a，马宁a,b*，邵闯a，顾解忡a,b. 限宽水域中船舶平面运动机构试验及水动力导数数值模拟[J]. 上海交通大学学报（自然版）, 2016, 50(01): 115-122.
[12]	程演, 张璐璐, 张磊, 王建华. 基于随机场的非饱和土固结分析[J]. 上海交通大学学报, 2014, 48(11): 1528-1535.
[13]	王宏志a，邹早建a,b. 航行船与停泊船水动力相互作用的数值分析[J]. 上海交通大学学报（自然版）, 2014, 48(04): 445-449.
[14]	刘鹏，苏玉民，刘焕兴，廖煜雷. 串列异步拍动翼推进性能分析[J]. 上海交通大学学报（自然版）, 2014, 48(04): 457-463.
[15]	胡世良a，鲁传敬a,b，何友声a,b. 平板流固耦合振动的数值分析[J]. 上海交通大学学报（自然版）, 2013, 47(10): 1487-1493.