EI、Scopus 收录

## 留言板 引用本文: 姚学昊, 陈丁, 武立伟, 黄丹. 流固耦合破坏分析的多分辨率PD-SPH方法. 力学学报, 2022, 54(12): 3333-3343 Yao Xuehao, Chen Ding, Wu Liwei, Huang Dan. A multi-resolution PD-SPH coupling approach for structural failure under fluid-structure interaction. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(12): 3333-3343 doi: 10.6052/0459-1879-22-268
 Citation: Yao Xuehao, Chen Ding, Wu Liwei, Huang Dan. A multi-resolution PD-SPH coupling approach for structural failure under fluid-structure interaction. Chinese Journal of Theoretical and Applied Mechanics, 2022, 54(12): 3333-3343 ## 流固耦合破坏分析的多分辨率PD-SPH方法

##### doi: 10.6052/0459-1879-22-268

###### 作者简介:黄丹, 教授, 主要研究方向: 计算力学与工程仿真研究. E-mail: danhuang@hhu.edu.cn
• 中图分类号: O353.4

## A MULTI-RESOLUTION PD-SPH COUPLING APPROACH FOR STRUCTURAL FAILURE UNDER FLUID-STRUCTURE INTERACTION

• 摘要: 流固耦合破坏是一类涉及结构变形与破坏以及复杂自由表面现象的强非线性力学问题. 结合近场动力学(peridynamics, PD)与光滑粒子流体动力学(smoothed particle hydrodynamics, SPH)各自的优势并考虑其计算效率问题, 提出一种适用于分析流−固耦合破坏问题的多分辨率PD-SPH混合方法. 分别采用SPH和PD方法以不同的空间和时间分辨率对流体和结构进行离散与求解, 利用具有与流体粒子相同光滑长度的虚粒子处理流−固界面, 以高精度满足界面边界条件. 通过两个经典算例: 液柱静压力下弹性板的变形和溃坝流体冲击弹性闸门的变形问题, 表明提出的多分辨率PD-SPH方法兼具较高的计算精度和计算效率; 对含裂缝的Koyna重力坝水力劈裂问题进行模拟计算, 所得裂缝扩展路径与文献结果吻合, 说明该方法适用于涉及结构破坏的流固耦合问题仿真. 最后尝试采用该方法进行流体冲击作用下含裂纹混凝土板崩塌过程数值仿真, 准确描述混凝土板的断裂破坏和全过程中的流体运动. 多分辨率PD-SPH混合方法或可为流−固耦合作用下的结构损伤破坏仿真提供一种新选择.

• 图  1  多分辨率PD−SPH耦合方案

Figure  1.  Multi-resolution PD−SPH coupling scheme

图  2  $\ell = 6$ 时多时间步长方案

Figure  2.  Multi-timestep scheme ($\ell = 6$ )

图  3  多分辨率PD-SPH耦合计算流程

Figure  3.  Flow chart for multi-resolution PD-SPH coupling calculation

图  4  液柱静压力下铝板变形的初始条件(单位: m)

Figure  4.  Initial condition for the aluminum plate under hydrostatic pressure (unit: m)

图  5  铝板跨中挠度时间历程

Figure  5.  Time histories for mid-span deflection of aluminum plate

图  6  t = 5 s时工况III的流体压力场与板的挠度场

Figure  6.  Fluid pressure and plate deflection field in case-III at t = 5 s

图  7  FSI系统归一化总能量时间历程

Figure  7.  Time histories for normalized total energy of FSI system

图  8  不同工况的计算时间

Figure  8.  Computational times for different cases

图  9  溃坝流体冲破弹性闸门的初始条件示意图(单位: mm)

Figure  9.  Schematic sketch for initial conditions in the dam-break flow through an elastic gate (unit: mm)

图  10  不同时刻弹性闸门变形和自由液面演变

Figure  10.  Elastic gate deformation and change of fluid free surface

图  11  弹性闸门自由端水平位移时程

Figure  11.  Time histories of horizontal displacements at the free end of elastic gate

图  12  3种工况的计算总时长

Figure  12.  Total calculation time in three cases

图  13  Koyna重力坝模型(单位: m)

Figure  13.  Koyna gravity dam model (unit: m)

图  14  Koyna大坝裂缝扩展过程

Figure  14.  Crack propagation process in Koyna dam

图  15  重力坝裂缝扩展路径

Figure  15.  Crack paths in the gravity dam

图  16  水流冲击混凝土板模型(单位: m)

Figure  16.  Model of water impacting on a concrete slab (unit: m)

图  17  水流冲击作用下混凝土板崩塌过程

Figure  17.  Collapse process of concrete slab subjected to impact of water flow

图  18  A点竖向位移时程曲线

Figure  18.  Time history of vertical displacement at point A

表  1  液柱静压力下铝板变形问题工况布置

Table  1.   The setup of different cases for aluminum plate deformation with hydrostatic pressure

 Cases $\dfrac{{d{x_{\text{F}}}}}{{d{x_{\text{S}}}}}$ $d{x_{\text{S}}}{\text{/m}}$ $\dfrac{{\Delta {t_{\text{F}}}}}{{\Delta {t_{\text{S}}}}}$ $\Delta {t_{\text{S}}}{\text{/s}}$ I 1.0 6.25 × 10−3 5.0 1.0 × 10−6 II 2.0 6.25 × 10−3 10.0 1.0 × 10−6 III 4.0 6.25 × 10−3 20.0 1.0 × 10−6

表  2  溃坝流体冲破弹性闸门问题工况布置

Table  2.   The setup of different cases for dam-break flow through an elastic gate

 Cases $\dfrac{{d{x_{\text{F}}}}}{{d{x_{\text{S}}}}}$ $d{x_{\text{S}}}{\text{/m}}$ $\dfrac{{\Delta {t_{\text{F}}}}}{{\Delta {t_{\text{S}}}}}$ $\Delta {t_{\text{S}}}{\text{/s}}$ I 1.0 1.0 × 10−3 1.0 5.0 × 10−6 II 2.0 1.0 × 10−3 1.0 5.0 × 10−6 III 2.0 1.0 × 10−3 4.0 5.0 × 10−6
• 点击查看大图
##### 计量
• 文章访问数:  518
• HTML全文浏览量:  187
• PDF下载量:  145
• 被引次数: 0
##### 出版历程
• 收稿日期:  2022-06-14
• 录用日期:  2022-08-10
• 网络出版日期:  2022-08-11
• 刊出日期:  2022-12-15

### 目录 / 下载:  全尺寸图片 幻灯片
• 分享
• 用微信扫码二维码

分享至好友和朋友圈
{"ft_lang_label":"__label__en","ft_lang_prob":0.6342759,"math_prob":0.9750764,"size":11858,"snap":"2023-14-2023-23","text_gpt3_token_len":4295,"char_repetition_ratio":0.14332715,"word_repetition_ratio":0.0397878,"special_character_ratio":0.3304942,"punctuation_ratio":0.22167689,"nsfw_num_words":0,"has_unicode_error":false,"math_prob_llama3":0.967559,"pos_list":[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18],"im_url_duplicate_count":[null,null,null,null,null,null,null,3,null,3,null,null,null,3,null,3,null,null,null],"WARC_HEADER":"{\"WARC-Type\":\"response\",\"WARC-Date\":\"2023-06-07T06:09:50Z\",\"WARC-Record-ID\":\"<urn:uuid:f93b2d60-c868-4d50-a993-c44ba84c7013>\",\"Content-Length\":\"183732\",\"Content-Type\":\"application/http; msgtype=response\",\"WARC-Warcinfo-ID\":\"<urn:uuid:13da0bbf-c34d-4c56-be02-f83789ab14b4>\",\"WARC-Concurrent-To\":\"<urn:uuid:e54a7778-70d9-4605-8368-3d9277a93e73>\",\"WARC-IP-Address\":\"123.56.203.87\",\"WARC-Target-URI\":\"https://lxxb.cstam.org.cn/cn/article/doi/10.6052/0459-1879-22-268\",\"WARC-Payload-Digest\":\"sha1:AZFONDFWA75Z5A64RCVDIOI66P45E2RL\",\"WARC-Block-Digest\":\"sha1:4FWVTK73KEEFWTGPJUGCD32YWQMBGM24\",\"WARC-Identified-Payload-Type\":\"text/html\",\"warc_filename\":\"/cc_download/warc_2023/CC-MAIN-2023-23/CC-MAIN-2023-23_segments_1685224653608.76_warc_CC-MAIN-20230607042751-20230607072751-00105.warc.gz\"}"}