基于近场动力学理论的混凝土重力坝地震响应数值模拟

赵凯丽, 邱流潮, 李敬军

PDF(3007 KB)
节水灌溉
中国农村水利水电 ›› 2021 ›› (7) : 139-145,150.
水电建设

基于近场动力学理论的混凝土重力坝地震响应数值模拟

作者信息 +

Numerical Simulation of Seismic Response in Concrete Gravity Dam Based on Peridynamics

Author information +
稿件信息 +

摘要

近场动力学基于非局部作用思想直接求解积分形式的运动方程，避免了传统的局部微分方程求解不连续问题时的奇异性，可以有效地求解连续和非连续变形问题，在众多领域已得到成功应用。本文建立了混凝土重力坝地震动力响应的键基近场动力学模型，采用FORTRAN90语言编程实现，并通过在三层虚拟节点施加地震加速度实现了地震波的输入。首先通过悬臂梁自由端受力和底部周期运动时的动力响应两个算例验证了模型和程序的精度和有效性，最后基于本文模型对混凝土重力坝的地震动力响应进行了数值模拟。结果表明：本文建立的键基近场动力学模型可以较好的模拟混凝土重力坝的地震动力响应，为混凝土坝地震动力响应数值模拟方法提供了一种新的思路和方法。

Abstract

Peridynamics is based on the idea of non-local action and directly solves the motion equation in the integral form， thus avoiding the singularity of traditional local differential equation in solving discontinuous problems. It can effectively solve continuous and discontinuous deformation problems and has been successfully applied in many fields. A bond-based peridynamics model is established to simulate the seismic dynamic response of a concrete gravity dam. It is implemented in FORTRAN90 programming language. In the proposed model， the input of seismic wave is realized by applying seismic acceleration to the three layer virtual node. The paper first verifies the accuracy and effectiveness of the model and the program through two examples of a force acts on the free end of the cantilever beam and the dynamic response under the input periodic motion at the bottom. Finally， the numerical simulation of the seismic dynamic response of the concrete gravity dam is performed based on the proposed model. The results show that it is feasible to simulate the seismic dynamic response of the concrete gravity dam with the bond-based peridynamics model and provides a new idea for the numerical simulation of seismic response of the concrete dam.

关键词

混凝土重力坝 / 地震输入 / 近场动力学 / 地震响应 / 数值模拟

Key words

concrete gravity dam / earthquake input / peridynamics / seismic response / numerical simulation

基金

国家自然科学基金项目(11772351)

引用本文

EndNote

Ris (Procite)

Bibtex

导出引用
赵凯丽 , 邱流潮 , 李敬军. 基于近场动力学理论的混凝土重力坝地震响应数值模拟[J].中国农村水利水电, 2021(7): 139-145,150
Kai-li ZHAO , Liu-chao QIU , Jing-jun LI. Numerical Simulation of Seismic Response in Concrete Gravity Dam Based on Peridynamics[J].China Rural Water and Hydropower, 2021(7): 139-145,150

0 引 言

随着国际能源发展战略方向的转移，大力发展水电等可再生能源成为了许多国家战略发展的新方向。我国如三峡水电站、糯扎渡水电站等大、中型水电站多数位于西部或西南部高烈度地震区［1］，地震作用下的安全问题是影响大坝可持续发展的重要因素，强烈地震作用下水工结构可能会产生结构破坏甚至垮塌的工程事故，从而造成人员伤亡、资金损失等一系列的灾难性后果。此外，我国大多数大坝都是在地震现场数据有限的情况下设计的［2］，大坝的除险加固研究也显得尤为重要。
随着计算机技术的高速发展，相较物理模型试验而言，数值模拟方法具有成本低、效率高、简单易行、可研究材料微观特性等一系列优点［3］，为大坝的地震动力响应研究提供了一种可行可靠的研究手段。在各种数值计算方法中，有限元方法的理论框架经过几十年的发展已相对较为成熟，在大坝的地震动力响应研究中被广泛应用［4-7］。然而，有限元方法是一种基于连续介质力学的方法，在求解非连续问题时仍存在局限性［5］，求解断裂破坏等非连续问题时会产生奇异性。针对上述不足，相关学者提出了一些改进方法，例如，Zhang等［8］基于扩展有限元方法研究重力坝在地震作用下的裂纹扩展情况，Fang等［9］应用扩展有限元方法模拟Koyna大坝地震裂缝，该方法相对于有限元法有所改进，可以解决弱不连续问题，但是难以模拟多裂纹扩展问题。Das等［10］基于无网格法预测地震作用下大坝溃坝前破坏机制，该方法在处理网格畸变问题时具有很大优势，克服了对网格的依赖性，然而其计算量较大且在模拟多裂纹扩展时受到限制［11］。Pekau等［12］基于非连续介质力学的离散元法对大坝进行裂缝分析，成功模拟了开裂问题，但是需要明确指出裂纹萌生点，在复杂的地震荷载作用下难以实现，并且计算量较大，因此需要一种新的方法来克服上述问题。
本文采用近场动力学（peridynamics，PD）方法进行混凝土重力坝地震响应数值模拟。该方法由美国Sandia实验室的 Silling［13］于2000年提出，黄丹等［14］2010年首次在国内应用。近场动力学法［15，16］是一种联系经典连续介质力学和分子动力学的数学力学理论，是分子动力学的连续版本，同时兼有分子动力学和无网格法的优点，在解决断裂等非连续大变形问题方面具有明显优势。近场动力学理论包括键基（bond-based）近场动力学理论和态基（state-based）近场动力学理论［17］，该理论已被应用到了许多研究领域，包括边坡稳定性分析［18］、非均质材料的水力压裂破坏模拟［19］、冰块撞击破坏［20］、岩土类材料［21］的大变形研究以及混凝土结构破坏过程模拟［22-24］等。Gu等［25］运用近场动力学键基理论研究了混凝大坝冲击破坏；张钰彬等［26］模拟了高水头作用下混凝土重力坝的水力劈裂过程，为研究大坝破坏机理提供了新方法，具有重要的工程意义。
本文应用键基近场动力学理论对混凝土重力坝的地震动力响应进行了数值研究。论文的余下部分安排如下：在第一节中简要介绍了键基近场动力学理论及数值计算方法；在第二节中首先通过悬臂梁自由端施加集中力和底部输入周期运动的正弦波时的动力响应两个算例来验证模型和程序的精度和有效性，随后基于验证过的模型对混凝土重力坝的地震动力响应进行数值模拟；在第三节中列出了本文研究得出的一些结论和进一步研究的方向。

1 键基近场动力学方法

1.1 运动方程

近场动力学理论不再基于连续性假设以及微分方程求解力学问题［27］，而是采用非局部作用的积分形式的运动方程，避免了非连续奇异性问题的发生。非局部思想反映在物质的相互作用方面，如图1所示，在某一时刻，空间域R内物体的任一物质点xi 与以其为中心，近场范围δ（horizon）为半径的圆形平面或球体范围内的其他物质点xj 之间产生相互作用，即当|xi -xj |≤δ时，两物质点之间产生相互作用力，而在近场范围内的物质点xj 构成xi 的家族（familiy），用H表示。
图1 物质点xi 和xj 变形过程中的相互作用

Fig.1 Interaction between material points xi and xj during deformation

Full size|PPT slide

根据牛顿第二定律可以得到物质点xi 的基于键基近场动力学理论的运动方程［28］：
ρxiu¨xi,t=∫HT̲xi,txj-xi-T̲xj,txi-xjdV+bxi,t
式中：ρ表示密度； u¨表示加速度；dV是物质点xj 处的体积微元； b 表示单位体积物质上施加的外部荷载，即体力密度； T̲xi,t为力的密度矢量，表示在t时刻物质点之间的相互作用。
对于键基近场动力学理论，物质点之间的相互作用力大小相等，方向相反，如图1所示。力密度矢量形式［29］如下：
T̲xi,txj-xi=12Cxj'-xi'xj'-xi'=12fuj-ui,xj-xi,t
T̲xj,txi-xj=-12Cxj'-xi'xj'-xi'=-12fuj-ui,xj-xi,t
式中：C是一个与物质点 xi和 xj之间相对伸长量以及近场动力学范围有关的辅助参数； xi'和 xj'表示物质点 xi和 xj变形后的位置。则有物质点xi 的运动方程为：
ρxiu¨xi,t=∫Hfuj-ui,xj-xi,tdV+bxi,t
式中：f表示xi 和xj 之间的相互作用，用“键”来表示，由其构成的函数称作本构力函数。
两物质点之间的相对位置用 ξ=xj-xi来表示，相对位移表示为： η=uj-ui， 式（4）可进一步写为：
ρxiu¨xi,t=∫Hfξ,ηdV+bxi,t
对于各向同性材料，可以假设力密度矢量f与物质点之间的拉伸成线性相关［29，30］：
fξ,η=c1sξ,ηxj'-xi'xj'-xi'
其中伸长率计算可以看做经典连续介质力学中的应变值计算如下［29］：
sξ,η=xj'-xi'-xj-xixj-xi
常数c 1的计算公式如下［29］：
c1=c=18 kπ δ4
式中：c为键常数微观模量；体积模量 k=E31-2 υ，E是弹性模量，υ是泊松比。

1.2 数值离散与时间积分

近场动力学方法求解积分形式的运动方程，通常将所研究物体离散为一系列物质点，于是得到控制方程（4）在空间域的离散表达式为：
ρ u¨xi,t=∑jfuj-ui,xj-xi,tVj+bxi,t
式中： Vj表示物质点 xj代表的体积。
本文对（9）式采用显式的向前差分时间积分求解，即按下列步骤计算：
u¨xi,tn+1=1ρ∑jfuj-ui,xj-xi,tnVj+bxi,tn
u˙xi,tn+1=u˙xi,tn+u¨xi,tn+1Δt
uxi,tn+1=uxi,tn+u˙xi,tn+1Δt
式中：n表示时间步数；Δt是时间积分步长，tn +1=tn +Δt，对于显式时间积分法，时间步长的选取需要满足稳定条件［30］。

2 数值算例与结果

本文通过使用FORTRAN 90语言编程实现了近场动力学的数值计算。首先通过悬臂梁自由端受集中力响应以及在底部输入周期运动时动力响应两个算例验证来模型和程序的精度和有效性，接着对混凝土重力坝的地震动力响应进行了数值模拟。

2.1 集中力作用下的悬臂梁动力响应

这里首先对悬臂梁自由端在集中力作用下的动力响应进行模拟。悬臂梁的长度L=0.02 m，宽度W=0.002 m和高度H=0.002 m，如图2所示。悬臂梁自由端施加动荷载P随时间变化的关系如图3所示。
图2 悬臂梁尺寸及示意图（单位：mm）

Fig.2 Dimensions and schematic diagram of cantilever beam

Full size|PPT slide

图3 荷载与时间之间的关系图

Fig.3 The relation diagram between load and time

Full size|PPT slide

数值计算中，悬臂梁一共离散为337 50个物质点，如图4所示。时间积分步长大小设置为Δt=1×10-8 s，仿真总时间为0.02 s。本算例中使用的材料参数如下：密度为ρ=1 000 kg/m3，弹性模量E=1 GPa，泊松比为 υ=0.25。本算例将集中力作为体力密度施加在悬臂梁最外层单元上，并且暂时不考虑阻尼 影响。
图4 悬臂梁初始时刻的离散效果图

Fig.4 Discrete rendering of the cantilever at the initial moment

Full size|PPT slide

图5给出了t=5×10-4 s、t=2×10-3 s、t=5.5×10-3 s和t=2×10-2 s 四个时刻悬臂梁的竖向位移云图。图6给出了悬臂梁自由端的竖向位移时程图，由于不考虑阻尼影响，1 ×10-3 s之后集中力大小不变时，自由端部数值模拟结果在1.5 mm与2.5 mm之间波动，即在解析解2 mm上下波动，与实际相符，并将本文键基近场动力学模拟结果与有限元软件ABAQUS的计算结果进行了比较，图6表明二者结果非常吻合，验证了本文键基近场动力学模型和程序的正确性。
图5 悬臂梁在不同时刻z方向上的位移

Fig.5 The displacement of the cantilever beam in the z direction at differenct moments

Full size|PPT slide

图6 自由端位移曲线比较

Fig.6 The comparison of free end displacement curves

Full size|PPT slide

2.2 正弦加速度激励下垂直方柱的动力响应

在结构的地震动力响应分析中，需要输入地震波，本文计算中假设地基是刚性的，因而在结构基底输入地震加速度。本文近场动力学模拟中将已知加速度施加在三层虚拟边界层的物质点上实现地震波输入，虚拟边界层的定义如图7所示。
图7 虚拟边界层的定义

Fig.7 The definition of virtual boundary layer

Full size|PPT slide

为了验证虚拟边界层物质点施加加速度模拟地震波的适用性和正确性，对垂直方柱在其底部正弦波激励作用下的动力响应进行数值模拟计算。方柱高度L=1.0 m，宽度W=0.02 m，厚度H=0.02 mm，如图8所示。在本算例中，在方柱的底部虚拟边界三层粒子的y方向上施加总时长为4 s的正弦加速度如图9所示。
图8 垂直悬臂柱体尺寸图（单位：m）

Fig.8 The dimension of vertical cantilever column

Full size|PPT slide

图9 正弦加速度

Fig.9 Sinusoidal acceleration

Full size|PPT slide

数值计算中，方柱一共离散为10 800个物质点，如图10所示。时间积分步长大小设置为 Δt=1×10-5 s，仿真总时间为4 s。本算例中使用的材料参数为：ρ=7 780 kg/m3，弹性模量E=1 GPa，泊松比为 υ=0.25。计算中不考虑阻尼影响。
图10 垂直悬臂柱体离散图

Fig.10 The discrete diagram of vertical cantilever column

Full size|PPT slide

图11给出了t=1 s、t=2 s、t=3 s和t=4 s共4个时刻方柱在y向位移云图，显示了柱体在施加周期波后在y方向上的摆动情况，且从图像可以看出悬臂柱体在以上4个时刻顶部位移要大于底部位移，与图12同一时刻相对应。图12给出了方柱顶部和底部y向的相对位移时程图，更直观的显示出了柱体的动力响应，并将本文键基近场动力学模拟结果与有限元软件ABAQUS的计算结果进行了比较，二者结果非常吻合，验证了本文基于键基近场动力学模型的虚拟边界层物质点施加加速度的方法实现地震波输入是可行的。
图11 悬臂梁在不同时刻y方向上的位移

Fig.11 The displacement of the cantilever in the y direction at different moments

Full size|PPT slide

图12 近场动力学模型与ABAQUS软件计算得到的柱顶部和柱底部相对位移

Fig.12 The relative displacement of column top and column bottom was calculated by the peridynamics model and ABAQUS software

Full size|PPT slide

2.3 混凝土重力坝地震动力响应

为了验证本文键基近场动力学模型模拟混凝土重力坝地震动力响应的可行性，本小节对Koyna混凝土重力坝在水平及竖向地震波共同作用下的动力响应进行模拟。研究坝断面高度103 m，底宽70 m，顶宽14.8 m，原坝型形状及尺寸如图13所示。
图13 大坝基本尺寸示意图（单位：m）

Fig.13 The schematic diagram of basic dam dimensions

Full size|PPT slide

数值计算中按照原坝型进行建模，为简化计算，将厚度方向取0.5 m，为一个单元高度大小，混凝土重力坝一共离散为 28 840个物质点，如图14所示。时间积分步长大小设置为 Δt=1×10-5 s，仿真总时间为10 s。本算例中使用的材料参数为：ρ=2 643 kg/m3，弹性模量E=31.02 GPa，泊松比为 υ=0.25。在本算例中，在坝的底部虚拟边界三层粒子的x及y方向上施加图15所示的地震加速度。计算不考虑库水及阻尼的影响，同时假设坝基为刚体。
图14 坝体离散图

Fig.14 The discrete diagram of the dam

Full size|PPT slide

图15 地震加速度

Fig.15 seismic acceleration

Full size|PPT slide

图16和17分别给出了t=2 s、t=4 s、t=8 s和t=10 s共4个时刻混凝土坝水平向和竖向位移云图。图18给出了混凝土重力坝顶水平向和竖向相对位移时程图的键基近场动力学模拟结果，数值计算中采用键基近场动力学模型在虚拟边界层物质点施加加速度的方法实现地震波输入，从图像可以看出，坝体的各个部位在不同方向上的动力响应情况随时间是变化的，由于本模拟厚度方向上取一个单元厚度，故暂不考虑z方向上的动力响应情况。从图18可以看出大约在3 s左右坝体横向和纵向位移开始有明显波动，与施加的地震波波动大小相对应，很好说明了本文键基近场动力学模型模拟混凝土坝地震动力响应是可行的。
图16 重力坝在不同时刻沿x方向的位移

Fig.16 The displacement of the gravity dam in the x direction at different moments

Full size|PPT slide

图17 重力坝在不同时刻沿y方向的位移

Fig.17 The displacement of the gravity dam in the y direction at different moments

Full size|PPT slide

图18 坝顶相对位移时程图

Fig.18 The relative displacement time history of dam crest

Full size|PPT slide

3 结 论

本文建立了模拟混凝土重力坝地震动力响应的键基近场动力学模型，并采用FORTRAN90语言编程实现。其中，地震波的输入是通过在三层虚拟节点施加地震加速度来实现的。本文通过对比悬臂梁自由端和底部受动力和周期运动正弦波的近场动力学数值模型与ABAQUS软件的模拟结果，证明该近场动力学模型适宜模拟动力响应问题。接着利用本文模型对混凝土重力坝的地震动力响应进行数值模拟，结果表明：本文键基近场动力学模型模拟混凝土重力坝的地震动力响应是可行的，为混凝土坝地震动力响应数值模拟方法提供了新的思路。
近场动力学需要计算近场范围内的所有物质点的相互作用，较有限元等传统的网格方法需要计算的物质点多，因此计算时间较长，今后有待于进一步开展并行算法研究以提高计算效率。本文选用较为简单、应用较广的近场动力学键基模型，存在泊松比的限制，一些情况下可能产生误差，后续可以采用态基近场动力学模型进行研究。另外，本文计算中没有考虑坝基的变形以及辐射阻尼的影响，将在后期研究中进一步考虑。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
1
王旭东，张立翔，段其品.基于ANSYS的碾压混凝土重力坝抗震稳定性分析［J］.中国水运（下半月），2018，18（3）：189-190.
本文引用 [1]
2
GAVABAR S G， ALEMBAGHERI M. A study on seismic performance and fragility of gravity dams with various monolith shapes using nonlinear IDA ［J］. Geotechnical and Geological Engineering， 2020，38（2）：1 133-1 150.
本文引用 [1]
3
张国新，金 峰，王光纶.用基于流形元的子域奇异边界元法模拟重力坝的地震破坏［J］.工程力学，2001（4）：18-27.
本文引用 [1]
4
朱伯芳.有限单元法原理与应用［M］.2版.北京：中国水利水电出版社，1998.
本文引用 [1]
5
房树财，何 佳，刘小军.强震作用下戈兰滩重力坝动力响应分析［J］. 福建建材， 2015（10）：1-4.
本文引用 [1]
6
WANG H， FENG M， YANG H. Seismic nonlinear analyses of a concrete gravity dam with 3D full dam model ［J］. Bulletin of Earthquake Engineering， 2012，10（6）：1 959-1 977.
7
邱流潮，刘 桦， 金 峰.混凝土重力坝地震响应数值分析方法［J］. 河海大学学报（自然科学版）， 2006（1）：65-68.
本文引用 [1]
8
ZHANG S， WANG G， YU X. Seismic cracking analysis of concrete gravity dams with initial cracks using the extended finite element method［J］. Engineering Structures， 2013，56：528-543.
本文引用 [1]
9
FANG X， JIN F， WANG J. Seismic fracture simulation of the Koyna gravity dam using an extended finite element method ［J］. Journal of Tsinghua University （Science and Technology）， 2008，48（12）： 2 065-2 069.
本文引用 [1]
10
DAS R， CLEARY P W. A mesh-free approach for fracture modelling of gravity dams under earthquake ［J］. International Journal of Fracture， 2013， 179（1-2）： 9-33.
本文引用 [1]
11
卢杰志.基于近场动力学理论的混凝土损伤断裂行为数值模拟［D］.武汉：华中科技大学，2019.
本文引用 [1]
12
PEKAU O A， YUZHU C. Failure analysis of fractured dams during earthquakes by DEM ［J］. Engineering Structures， 2004，26（10）： 1 483-1 502.
本文引用 [1]
13
SILLING S A. Reformulation of elasticity theory for discontinuities and long-range forces ［J］. Journal of the Mechanics and Physics of Solids， 2000，48（1）：175-209.
本文引用 [1]
14
黄 丹，章 青，乔丕忠，等.近场动力学方法及其应用［J］. 力学进展，2010，40（4）：448-459.
本文引用 [1]
15
SILLING S A， LEHOUCQ R B. Peridynamic Theory of Solid Mechanics ［M］. Advances in Applied Mechanics， Aref H， Giessen E V D， Elsevier， 2010， 44， 73-168.
本文引用 [1]
16
Silling S A， Demmie P N， Warren T L. Peridynamic Simulation of High-Rate Material Failure ［Z］. Presented at the Applied Mechanics and Materials Conference， 2007.
本文引用 [1]
17
SILLING S A， EPTON M， WECKNER O， et al. Peridynamic States and Constitutive Modeling ［J］. Journal of Elasticity， 2007，88（2）：151-184.
本文引用 [1]
18
廖 洋，刘立胜，刘齐文.近场动力学方法在边坡稳定性分析中的应用［J］.安全与环境工程，2018，25（1）：45-50.
本文引用 [1]
19
NADIMI S， MISCOVI I， MCLENNAN J. A 3D Peridynamic Simulation of Hydraulic Fracture Process in a Heterogeneous Medium［J］. Journal of Petroleum Science and Engineering， 2016，145：444-452.
本文引用 [1]
20
SONG Y， YAN J， LI S， et al. Peridynamic Modeling and Simulation of Ice Craters By Impact ［Z］. 2019，121（2）：465-492.
本文引用 [1]
21
SONG X， KHALILI N. A peridynamics model for strain localization analysis of geomaterials ［J］. International Journal for Numerical and Analytical Methods in Geomechanics， 2018，43（1）：77-96.
本文引用 [1]
22
DEMMIE P， SILLING S. Mechanics of materials and structures ［Z］. 2007： 2.
本文引用 [1]
23
沈 峰，章 青，黄 丹，等.冲击荷载作用下混凝土结构破坏过程的近场动力学模拟［J］.工程力学，2012，29（）：12-15.
摘要
增刊1
24
郁杨天，章 青，顾 鑫.含单边缺口混凝土梁冲击破坏的近场动力学模拟［J］.工程力学，2016，33（12）：80-85.
本文引用 [1]
25
GU X， ZHANG Q. A modified conjugated bond-based peridynamic analysis for impact failure of concrete gravity dam ［J］. Meccanica， 2020，55（3）：547-566.
本文引用 [1]
26
张钰彬，刘一鸣，黄 丹.水力劈裂问题的态型近场动力学建模［J］.计算力学学报，2019，36（1）：21-26.
本文引用 [1]
27
乔丕忠，张 勇，张 恒，等.近场动力学研究进展［J］.力学季刊，2017，38（1）： 1-13.
本文引用 [1]
28
WU L， HUANG D， XU Y， et al. A non-ordinary state-based peridynamic formulation for failure of concrete subjected to impacting loads［J］. CMES-Computer Modeling in Engineering & Sciences， 2019，118（3）： 561-581.
本文引用 [1]
29
ERDOGAN M， ERKAN O. Peridynamic theory and its applications springer link ［M］. Springer， New York， NY， 2014.
本文引用 [4]
30
SILLING S A， ASKARI E. A meshfree method based on the peridynamic model of solid mechanics ［J］. Computers & Structures， 2005，83（17-18）：1 526-1 535.
本文引用 [2]
PDF(3007 KB)

946

访问

0

引用

详细情况
段落导航
相关文章
鄂公网安备 42010602001274号　鄂ICP备14015464号-1
网站版权所有 © 《中国农村水利水电》编辑部
地址：武汉大学二区《中国农村水利水电》编辑部　邮编：430072
电话：(027)68776133　传真：xsdbjb@188.com　
本系统由北京玛格泰克科技发展有限公司设计开发　技术支持：support@magtech.com.cn

/

〈 〉