Dynamic response of flowing ice colliding with a sluice pier under hydrodynamic action
-
摘要: 高寒地区河冰撞击河道的闸墩结构会产生极端冰载荷和冰激振动，水的动力效应使得碰撞过程更加复杂。采用任意拉格朗日-欧拉流固耦合方法，考虑作用在流冰和闸墩表面的流体力，建立了水-冰-闸墩耦合模型，探究了偶然极端条件下冰-闸墩碰撞的力学特性，设计了冰-砼碰撞实验。结果表明：冰-砼碰撞实验中，撞击力的模拟结果与实验结果吻合良好；对流固耦合的水动力效应分析发现，水-冰-闸墩耦合模型能够体现水的流体特性，在流冰撞击闸墩近场逼近过程中，初始时刻水的动力效应能够增加流冰的动能，撞击楔入闸墩过程中，水介质形成一个瞬态高压力场，产生水垫效应吸收冰体部分动能，从而抑制流冰运动；在不同流冰体积和压缩强度工况下，闸墩结构所承受的冰力随着流冰体积的增大而增大，流冰压缩强度对冰力的影响较小，流冰损伤与闸墩结构响应主要集中在碰撞接触区，流冰撞击闸墩结构引起冰激振动，流冰体积对闸墩振动加速度的影响较大，相同体积的流冰随着压缩强度的增大，振动幅值差异不明显，表明流冰体积是影响冰-闸墩碰撞的关键参数。Abstract: In frigid regions, the construction of sluice pier structures within river systems is confronted with considerable challenges arising from the presence of severe ice loads and ice-induced vibrations. The collision process between ice and sluice piers is further complicated due to the intricate hydrodynamic effects exerted by water. The arbitrary Lagrangian-Eulerian (ALE) fluid-structure interaction (FSI) method is employed in this research to meticulously account for the fluid forces acting upon both the ice and sluice pier surfaces. A comprehensive coupled model encompassing the interactions among water, ice, and sluice piers is established to thoroughly investigate the mechanical characteristics associated with ice-sluice pier collisions under highly unpredictable conditions. Corresponding ice-concrete collision tests are meticulously designed and conducted, revealing an exemplary concurrence between the simulated impact forces and the values obtained from experimental observations. Upon analyzing the fluid-structure interaction and hydrodynamic effects, the present study demonstrates that the water-ice-sluice pier coupled model adeptly captures the fluid characteristics inherent to water. During the approach of an ice mass towards a sluice pier, the initial hydrodynamic effects initiated by the water medium effectively augment the kinetic energy possessed by the ice. As the ice forcefully interacts with the sluice pier, the water medium swiftly generates a transient high-pressure field, thereby establishing a phenomenon colloquially referred to as the water cushion effect. This effect is manifested by absorbing a portion of the ice’s kinetic energy, effectively dampening its movement. Distinctive scenarios characterized by varying ice volumes and compression strengths elucidate that the ice forces exerted upon the sluice pier structure directly correlate with the magnitude of the ice volume, while the influence of ice compression strength on said forces is relatively negligible. The consequential damages inflicted upon the ice and the response exhibited by the sluice pier structure primarily manifest within the contact area at the moment of collision. Consequently, the collisions between ice and the sluice pier structure induce vibrations that are uniquely attributed to ice-related factors. The volume of ice significantly influences the acceleration of sluice pier vibrations. Furthermore, under the condition of maintaining a consistent ice volumes, an increase in compression strength yields only marginal discrepancies in vibration amplitude. This finding convincingly substantiates the critical role played by ice volume as the paramount parameter governing ice-sluice pier collisions.
-
Key words:
- fluid-structure interaction /
- hydrodynamic /
- sluice pier /
- collision forces /
- ice-induced vibration
-
多孔材料由胞孔组成, 不同冲击速度下胞孔的变形特性不同, 导致多孔材料宏观动态力学性能也不同。S.Lee等[1]和P.J.Tan等[2-3]指出, 冲击速度高时多孔材料以冲击波模式变形。Zheng Zhi-jun等[4]和刘耀东等[5]通过数值模拟总结出随着速度提高, 多孔材料分别呈现出准静态模式、过渡模式和冲击模式等3种变形模式。随着速度的提高, 试件的变形均匀性和两端应力均匀性都越来越差, 如何测量变形和应力的不均匀性以及讨论两者的关系成为了解多孔材料动态力学性能的难题。王鹏飞等[6]利用改进的Hopkinson杆实验装置测得了不同冲击速度下泡沫铝试件两端的应力。在变形不均匀性方面, Z.Zou等[7]和S.Pal等[8]根据数值模拟结果直接给出了冲击模式下试件中一维和二维的应变分布, 但用此法得出的结果误差较大。Liao Shen-fei等[9]将数值模拟结果结合最佳局部变形梯度原理推导出了更精确的应变场。如何在实验中准确测量泡沫铝试件压缩过程中的变形特性成了难题, 虽然通过数值模拟能够从规律上进行解释, 但是与实验仍有一定差距。采用常规的电测方法仅能测试试件的平均应变, 无法测量其中的不均匀性。高速摄影只能提供一种直观判断, 是一种定性的测量手段, 而且对于微小的变形无法判断。W.H.Peters等[10]提出的数字图像相关法(digital image correlation method, DICM)是光测力学中的一种简单方便精度高的测试手段, 经过多年的发展, 已经在生物力学、岩石力学、纳米力学等众多领域得到广泛应用。近年来, 这一方法也用于多孔材料研究, 但主要集中在准静态实验方面[11-12], 而动态性能方面研究较少, 主要有:Wang Li-li等[13]利用DICM测试了泡沫铝的泰勒杆实验中试件的速度随撞击时间的变化; S.Lee等[1]和I.Elnasri等[14]分别研究了开孔和闭孔泡沫铝高速变形过程中的全场应变; H.Luo等[15]在试件表面粘上带有格子的标签, 测量了泡沫非晶金属材料。然而上述研究并未对3种模式下的变形场作比较, 也没有给出全场应变发展的过程, 更没有进一步讨论变形均匀性与应力均匀性的关系。基于此, 本文中, 将结合Hopkinson压杆实验技术和高速摄影技术, 研究不同冲击速度下泡沫铝的全场应变及其变化趋势, 并依此研究泡沫铝试件中的应力不均匀性。
1. 冲击压缩实验
试件为直径为32 mm、高度为32 mm的闭孔泡沫铝, 密度为0.320~0.340 g/cm3。实验在直径为37 mm的Hopkinson压杆装置上进行, 实验系统示意图如图 1所示, 杆材为铝。共进行了3组实验, 子弹速度分别为12、50和110 m/s, 对应3种不同的变形模式:即准静态模式、过渡模式和冲击模式。为了研究准静态模式下试样的应力均匀性, 利用常规Hopkinson杆装置进行准静态模式下的实验, 见图 1(a)。由于子弹速度过高会导致打击杆屈服, 采用改进的Hopkinson杆装置(图 1(b))进行后2种模式的实验, 即将试件放在支撑杆与子弹之间, 让子弹直接撞击试件。图中的Wheatstone电桥接法中含有2套单独的应变片系统, 每个系统中含有2片应变片且使用的电桥接法为半桥接法。
实验所用高速摄影相机为Phantom V12.1, 全画幅分辨率为1 280 Pixel×800 Pixel, 全画幅拍摄速率可达6 242 s-1, 最高拍摄速率可达到1 000 000 s-1, 3组实验的拍摄频率皆为74 000 s-1, 此时帧间隔时间为13.5 μs, 3种模式对应的实际图像分辨率见表 1。实验中使用了2个脉冲氙灯来提高曝光量。
表 1 实验参数Table 1. Experimental parameters模式 v/(m·s-1) ρ/(g·cm-3) 实验方法 图像分辨率 拍摄速率/s-1 准静态模式 12 0.325 常规SHPB 256 Pixel×200 Pixel 74 000 过渡模式 50 0.340 直接撞击 304 Pixel×168 Pixel 74 000 冲击模式 110 0.338 直接撞击 288 Pixel×232 Pixel 74 000 2. 变形均匀性
根据高速摄影图像可以分别得到准静态模式、冲击模式和过渡模式下试件的变形过程, 如图 2~4所示。其中高速摄影的第1张图像都对应着试件刚开始发生变形的前一时刻, 并将该时刻定为0 μs。
图 2 准静态模式下试件的变形过程Figure 2. Deformation progress of the specimenin the quasi-static mode
图 4 过渡模式下试件的变形过程Figure 4. Deformation progress of the specimen in the transition mode采用自编的程序处理高速摄影图像, 即可得到变形过程中每一个时刻相对于上一时刻试件的增量应变场, 如处理0和13.5 μs的2张图像即可得到0 μs时刻试件的增量应变场。本文中高速摄影的图像间隔时间均为13.5 μs, 所以本文中采用DICM得到的计算结果的间隔时间也是13.5 μs。对于冲击模式下的图像, 在冲击波波阵面附近的区域由于应变过大, 同时伴有翻转和弯折, 区域内大部分点的灰度值在相邻2张图像中已不存在联系。这就违背了数字图像相关方法的假定, 使该区域无法与变形前图像进行相关计算。因此在计算过程中只计算冲击波波阵面前方的区域, 如图 4(g)所示。图 5~7分别表示3种模式下试件沿着加载方向全场应变的变化过程。本方法的标定可参见文献[16]。
图 5 准静态模式下试件变形过程中的应变场Figure 5. Strain field of deformation progress under quasi-static mode
图 6 冲击模式下试件变形过程中的应变场Figure 6. Strain field of deformation progress under dynamic mode
图 7 过渡模式下试件变形过程中的应变场Figure 7. Strain field of deformation progress under transition mode从高速摄影图像上来看, 准静态模式下(图 2)试件变形较均匀, 在变形的前半段试件的变形以均匀变形为主, 没有产生局部变形。在变形的后半段, 试件开始产生局部变形并发展为如图 2(h)虚线所示的局部变形带。冲击模式下(图 3), 试件从一开始就变形不均匀。由于冲击端速度很高, 泡沫铝试件中与冲击端紧挨着的胞孔很快被压溃, 随后压溃区前方的胞孔也被压溃并发生逐层垮塌, 以一种压实波[17]的形式向支撑端发展。表现在图像中为试件从靠近打击杆开始变形, 并快速形成了密实带, 然后密实带(图 3中箭头所指)快速向前传播, 直至整个试件压实。过渡模式下(图 4), 试件的变形模式介于准静态模式和冲击模式之间, 先在靠近冲击端的区域以密实带的形式发展, 随后密实带前方的区域的变形以局部变形带为主。
采用DICM得到的计算结果反映了与高速摄影图像完全一样的变形趋势, 然而计算结果是量化的, 并且揭示了更多变形过程的细节。高速摄影图像中相邻时刻的图像之间变形较小, 无法判断其应变变化。采用DICM得到的计算结果则可以弥补这个缺陷, 能够展示整个变形过程中应变场的变化。
准静态模式下试件的变形过程如图 5所示, 按变形均匀性可以将整个过程分为2部分:均匀变形阶段A和不均匀变形阶段B。在阶段A中整个试件的变形较均匀, 塑形变形遍布整个试件而且较均匀。随着变形的发展, 试件的变形逐渐局部化, 塑形应变开始集中。在162 μs时, 即图中A、B区域的临界点, 塑形应变的集中最终形成了局部变形带。随后试件的变形主要集中在局部变形带上, 形状和方向与图 2(h)中的局部变形带一致。在压缩过程中, 试件中的局部变形带不断发展, 在试件变形的最后阶段试件的变形逐渐变小直至不变。值得一提的是, 试件中的局部变形带并不是一直发展的, 其变形到一定密实程度后又会逐渐停止发展, 并在试件的其他薄弱区域产生新的局部变形带[16]。在整个变形过程中, 开始段和结束段的应变率都较小, 而中间段的应变率较高。这是由于采用了常规SHPB装置, 冲击端处的入射波是一个梯形波, 应力幅值在加载和卸载段均较低造成的。
冲击模式下应变场变化见图 6, 从一开始试件的变形就集中在靠近冲击端附近, 形成一个很小的局部变形区域, 而其他区域的应变很小。随着压缩的发展, 局部变形区快速向前发展, 直至整个试件压溃。
为了更好地观测冲击波波阵面前方的应变场, 对应变场作进一步处理。将图 7中增量应变场换算成全量应变场, 再沿着y方向求和并求平均, 即可得到整个试件沿着加载方向的平均应变分布如图 8所示。结果表面:波阵面前方仍存在变形区且应变较大, 应变沿加载方向的分布与冲击波波阵面附近的应变分布类似。这与文献[7-9]的结果相似但仍存在区别, 这是因为本文中的试件较短, 在变形的开始段和结束段应变分布在很大程度上受到边界条件的影响。在起始段(0~67.5 μs), 应力波还未作用到冲击波前端所在地区, 此时冲击波前端的应变较小; 中间段(67.5~121.5 μs), 冲击波前端的应变分布与文献[7-9]的结果一致, 即只在冲击波前端约一个胞孔距离的区域有应变, 而远离冲击端的地方无应变; 末尾段(135μs), 此时冲击波距离支撑端很近, 整个试件都产生了较大的应变。图 8中变形中间段存在较清晰的应变下降段, 可认为是冲击波波阵面。根据其数据可以算出压实波波阵面分别相差17.9和19.6像素, 本实验中图像与真实试件的比例为5.625像素/毫米, 即可得波阵面的传播速度为236和258 m/s, 此结果略大于理论值, 其原因还有待研究。
图 8 冲击模式下沿着冲击方向的应变分布Figure 8. Strain distribution along the impact direction under dynamic mode过渡模式下的应变场变化如图 7所示。试件的变形过程可以分为2部分:局部致密区变化阶段和局部变形带变化阶段。在试件变形的前段, 由于冲击端速度较高但又未达到冲击模式的临界速度, 靠近冲击端部分的胞孔很快被压溃并形成了致密区, 但没有形成冲击波。试件的变形集中在靠近冲击端的位置, 并逐渐形成了致密区, 此时远离冲击端部分未发生变形。在变形后半段, 由于子弹速度降低和试件中存在缺陷等原因, 试件进入了一个局部变形带变形期。此时的变形由试件中的缺陷主导, 在有缺陷处产生局部变形带直到整个试件压实。
泡沫铝试件由于其本身的结构特性在冲击载荷下会产生局部化变形, 不同冲击速度下对应着不同的变形模式。本文利用数字图像相关方法, 从实验中测试了3种不同模式下试件的变形特征。准静态模式下, 试件先经历一个较长的变形均匀区, 随着外加载荷超过试件临界失稳载荷, 试件在薄弱处发生失稳坍塌, 形成了局部变形带。冲击模式下, 整个变形过程都是不均匀的, 试件的变形主要以压实波的形式发展。过渡模式下试件的变形过程是冲击模式和准静态模式变形过程的综合, 前半部分与冲击模式相似, 是一个局部致密化的过程, 后半部分与准静态模式相似, 变形集中在局部变形带上。
综上所述, 可知:(1)低速冲击时, 试件的局部化变形源于试件中含初始缺陷区域的失稳。在此情况下, 试件两端的载荷可基本达到均衡, 可认为整个试件受力均匀。试件中含缺陷区域最薄弱, 一定载荷下, 该区域开始产生局部化变形, 此时局部化变形以局部变形带为主。该区域的失稳发展到一定程度后, 由于变形导致该区域强度增大而停止失稳。此时变形转移到试件中其他薄弱的区域并形成新的局部变形带, 这可由文献[16]中给出的多个局部变形带证实。(2)高速冲击时, 试件的局部化变形源于冲击端区域的压缩失稳。当冲击速度达到一定值时, 试件两端的载荷很难达到均衡。此情况下, 冲击波的强度即可使材料发生破坏, 局部化变形发生在靠近冲击端的位置。因此, 对应试件的局部化变形由惯性效应决定, 表现为剧烈体积压缩的致密区。由此可见, 低速和高速冲击下材料的破坏分别对应着泡沫铝材料的整体失稳和局部失稳两种不同的机制。在这里, 应力的均匀性有重要作用。
3. 应力均匀性
为了讨论试件的应力均匀性, 在准静态模式下做了另外一组试件密度接近, 子弹速度、试件厚度等其他参数都相同的实验。利用石英片测出冲击端和支撑端的应力-时间曲线, 如图 9(a)所示。可以看到试件变形过程中存在2个区域:应力均匀区A和应力不均匀区B。在A区中冲击端应力大于支撑端应力, 而在B区中两端应力平衡。
图 9 泡沫铝试件两端的应力均匀性Figure 9. Stress uniformity of the two ends of the aluminum foam specimen在A区中两端应力不均匀的原因是2条曲线存在时间差, 即应力波在试件中传播的时间, t=L/c。现在将2条曲线对齐来判断应力波在试件中传播的耗散性。图 9(b)中曲线可以分为2部分:非耗散区C和耗散区D。非耗散区C中, 前后端应力相同, 说明应力在试件中的传播没有耗散。应力耗散区D内, 冲击端的应力明显大于支撑端的应力, 即应力发生了耗散。
同时可以看到, 应力开始发生耗散的时刻(165 μs)恰好是试件中产生局部变形带的时刻(162 μs), 即非耗散区对应着变形均匀区, 耗散区对应着变形非均匀区。当试件均匀变形时, 应力波相当于在一个均匀的物体中传播, 此时应力不发生耗散。而当胞孔发生了坍塌, 试件中形成了局部变形带。局部变形带所在区域的波阻抗大于其附近区域的波阻抗, 当应力波传播到这个界面上时会发生反射, 导致透射应力低于入射应力, 即应力发生耗散。这就说明变形的不均匀性影响到应力的不均匀性。
为了更好地说明上述问题, 建立了如下的简化模型, 即假设在试件中存在一个垂直加载方向的局部变形带。假设弹性波在局部变形带中传播速度、密度和弹性模量分别为c2、ρ2和E2, 在试件其他部分的参数对应为c1、ρ1和E1, 显然有ρ1c1 < ρ2c2。根据弹性波在不同介质界面上的传播规律, 有:
σr=Fσi,σt=Tσi (1) 式中:
和σt分别为入射波、反射波和透射波的应力。
则图 10中1、2、3区中弹性波存在关系:
。
由于n < 1, 即可得σ3/σ1 < 1。随着局部变形带的发展, 波阻抗不匹配度增大, n越小时, σ3/σ1越小, 即应力波的耗散越大。这就解释了随着局部变形带的发展, 两端应力不均匀性越来越大的原因。
随着速度的增大, 变形模式变成过渡模式和冲击模式, 过渡模式和冲击模式下两端的应力不均匀性显著增高[6]。从上文中所得到的应变场分析可以看到, 过渡模式和冲击模式下试件没有均匀变形的时期, 而是在靠近冲击端位置直接产生了密实区。局部变形带和密实区对于波传播的影响是一致的, 而且密实区相对密度更大, 与其余部分的波阻抗不匹配度更高。应力波经过界面后产生了更低的透射应力。此时冲击端和吸收端的应力差距更大, 这也从另一角度解释了随着速度提高, 试件两端的应力不均匀性更大的原因。
与图 10(a)类似, 高速冲击下试件产生致密区后的模型如图 10(b)所示, 参数定义与上文中相同。根据公式(1)可得:
σ1。
因为n=ρ1c1/(ρ2c2)≪1, 所以σ2 < σ1, 即密实区和试件其他部分的波阻抗不匹配度很高时, 支撑端应力远小于冲击端的应力。
综上所述:当试件变形均匀时, 其两端应力是均匀的; 当局部变形区域产生时, 试件变形开始不均匀, 此时试件两端应力也开始不均匀。这是由于局部变形不均匀的区域与其他部分波阻抗不匹配, 导致波在试件中的传播发生了反射, 即波在试件中传播发生了耗散, 使冲击端和吸收端的应力不均匀。
4. 结论
利用Hopkinson杆装置对泡沫铝试件进行了不同冲击速度的压缩实验, 利用高速摄影机对所有实验过程进行了跟踪拍摄。采用数字图像相关方法程序对高速摄影图像进行处理, 得到了不同变形模式下试件的应变场发展过程。
(1) 不同冲击速度下试件中的变形模式完全不同:准静态模式下, 试件开始段变形较均匀, 随后在试件中产生了局部变形带; 过渡模式下, 开始段变形集中在靠近冲击端的区域并形成了局部致密区, 结束段试件中的变形以局部变形带的方式发展; 冲击模式下, 试件变形主要以胞孔压实的模式向前发展。
(2) 不同速度下试件的变形机理也不同:低速下, 产生局部变形带的机理是试件中含缺陷胞孔的失稳; 高速下, 产生局部致密区主要是由惯性效应引起的。
(3) 材料不均匀性所导致的失稳(整体效应)和惯性效应(局部效应)是应力不均匀和变形不均匀的原因。在泡沫铝试件均匀变形时, 应力波在试件中的传播没有耗散。当试件中产生了局部变形带或致密区这样的局部变形区时, 由于局部变形带和致密区对波的反射作用, 使应力波产生了耗散。而且随着速度增高, 局部变形区与均匀变形区波阻抗不匹配程度增大, 反射波也越来越大, 这就导致了应力不均匀的加剧。
-
图 6 模拟和实验得到的撞击过程的应力时程曲线
Figure 6. Stress time history curves of impact processobtained by simulation and experiment
图 7 不同冰厚下模拟与实验最大应力的对比
Figure 7. Comparison of simulated and experimental maximum stresses for different ice thicknesses
图 8 不同时刻水流形态变化和y方向应变
Figure 8. Changes in water flow pattern and y-strain of water mediam at different times
图 10 不同工况下冰-闸墩撞击力曲线
Figure 10. Ice-sluice pier impact force curves under different working conditions
图 12 不同体积和压缩强度下的力-时间历程
Figure 12. Force-time histories under different volumes and compression strengths
图 13 不同体积和压缩强度与撞击力的关系
Figure 13. Relationships of the volume and compressive strength of ice with collision force
图 14 不同体积和压缩强度工况下冰激振动的加速度-时间历程
Figure 14. Acceleration-time histories of ice-excited vibration under different volume and compression strength conditions
图 15 流冰损伤和闸墩响应等效应力
Figure 15. Flowing ice damage and contours of pier response equivalent stress
混凝土材料参数 密度/（kg·m–3） 弹性模量/GPa 泊松比 初始抗拉极限/MPa 抗剪极限/MPa 断裂韧度/（N·m–1） 剪切保持力 2 500 30 0.2 4.02 21 0.14 0.03 混凝土材料参数 钢筋材料参数 体积黏度 压屈应力/MPa 弹性模量/GPa 屈服应力/MPa 硬化模量/GPa 失效应变 0.72 42 200 335 10 0.75 
下载: 导出CSV
表 2 冰材料模型参数
Table 2. Material parameters of ice
密度/（kg·m–3） 剪切模量/GPa 屈服应力/MPa 塑性硬化模量/GPa 体积模量/GPa 失效应变 截断应力/MPa 910 2.2 2.1 4.26 5.26 7.69×10–4 –4.0
下载: 导出CSV
表 3 水和空气介质材料参数
Table 3. Material parameters of water and air media
流体介质 密度/（kg·m–3） 截断压力/Pa 黏度系数/（N·s·m–2） C0 C1 C2 C3 C4 C5 E0/MPa V0 空气 1.184 5 –10 1.844×10−5 0 0 0 0 0.4 0.4 0.253 1.0 水 998.21 –1.0×10−5 1.790×10−3 1.0133×105 2.25×109 1.0
下载: 导出CSV
表 4 流冰-闸墩碰撞工况
Table 4. Flowing ice -pier collision conditions
工况 冰厚/m 冰温/℃ 冰速/（m·s–1） 冰体积/m3 冰压缩强度/MPa 1 0.3 –8 1.5 7.2 2.186 2 0.3 –8 1.5 14.4 2.186 3 0.3 –8 1.5 28.8 2.186 4 0.3 –8 1.5 64.8 2.186 5 0.3 –8 1.5 115.2 2.186 6 0.3 –2 1.5 64.8 1.123 7 0.3 –5 1.5 64.8 1.825 8 0.3 –8 1.5 64.8 2.186 9 0.3 –14 1.5 64.8 2.615 10 0.3 –20 1.5 64.8 2.889
下载: 导出CSV
-
[1] YANG X, PAVELSKY T M, ALLEN G H. The past and future of global river ice [J]. Nature, 2020, 577: 69–73. DOI: 10.1038/s41586-019-1848-1. [2] ROKAYA P, BUDHATHOKI S, LINDENSCHMIDT K E. Trends in the timing and magnitude of ice-jam floods in Canada [J]. Scientific Reports, 2018, 8(1): 5834. DOI: 10.1038/s41598-018-24057-z. [3] 郭新蕾, 王涛, 付辉, 等. 河渠冰水力学研究进展和趋势 [J]. 力学学报, 2021, 53(3): 655–671. DOI: 10.6052/0459-1879-20-407.GUO X L, WANG T, FU H, et al. Progress and trend in the study of river ice hydraulics [J]. Chinese Journal of Theoretical and Applied Mechanics, 2021, 53(3): 655–671. DOI: 10.6052/0459-1879-20-407. [4] 贡力, 杨腾腾, 靳春玲, 等. 水-空气耦合介质中冰载荷对闸墩的撞击影响研究 [J/OL]. 工程力学, (2023-04-03)[2022-11-28]. DOI: 10.6052/j.issn.1000-4750.2022.05.0491.GONG L, YANG T T, JIN C L, et al. Research on collision of ice loads on the sluice pier in water-air coupling medium [J/OL]. Engineering Mechanics, (2023-04-03)[2022-11-28]. DOI: 10.6052/j.issn.1000-4750.2022.05.0491. [5] WU C G, WEI Y M, JIN J L, et al. Comprehensive evaluation of ice disaster risk of the Ningxia-Inner Mongolia reach in the upper Yellow River [J]. Natural Hazards, 2015, 75(2): 179–197. DOI: 10.1007/s11069-014-1308-z. [6] 杨开林. 长距离输水水力控制的研究进展与前沿科学问题 [J]. 水利学报, 2016, 47(3): 424–435. DOI: 10.13243/j.cnki.slxb.20150824.YANG K L. Review and frontier scientific issues of hydraulic control for long distance water diversion [J]. Journal of Hydraulic Engineering, 2016, 47(3): 424–435. DOI: 10.13243/j.cnki.slxb.20150824. [7] 王娟, 黄樾, 邓宇, 等. 基于数字图像相关方法的黄河冰断裂性能研究 [J]. 水利学报, 2021, 52(9): 1036–1046. DOI: 10.13243/j.cnki.slxb.20201064.WANG J, HUANG Y, DENG Y, et al. Research on ice fracture of the Yellow River performance based on digital image correlation method [J]. Journal of Hydraulic Engineering, 2021, 52(9): 1036–1046. DOI: 10.13243/j.cnki.slxb.20201064. [8] WANG Z, SHI H B, LIU X M, et al. Analysis on the ice regime change characteristics in the Inner Mongolia reach of the Yellow River from 1950 to 2010 [J]. Journal of Coastal Research, 2020, 115(S1): 405. DOI: 10.2112/JCR-SI115-115.1. [9] 余同希, 朱凌, 许骏. 结构冲击动力学进展(2010–2020) [J]. 爆炸与冲击, 2021, 41(12): 121401. DOI: 10.11883/bzycj-2021-0113.YU T X, ZHU L, XU J. Progress in structural impact dynamics during 2010−2020 [J]. Explosion and Shock Waves, 2021, 41(12): 121401. DOI: 10.11883/bzycj-2021-0113. [10] JEON S, KIM Y. Numerical simulation of level ice-structure interaction using damage-based erosion model [J]. Ocean Engineering, 2021, 220: 108485. DOI: 10.1016/j.oceaneng.2020.108485. [11] CAI W, ZHU L, YU T X, et al. Numerical simulations for plates under ice impact based on a concrete constitutive ice model [J]. International Journal of Impact Engineering, 2020, 143: 103594. DOI: 10.1016/j.ijimpeng.2020.103594. [12] GONG L, DONG Z Q, JIN C L, et al. Flow-solid coupling analysis of ice-concrete collision nonlinear problems in the Yellow River Basin [J]. Water, 2023, 15(4): 643. DOI: 10.3390/w15040643. [13] ZHOU L, ZONG Z, LI J N. A numerical study of hydrodynamic influence on collision of brash ice with a structural plate [J]. Journal of Hydrodyn, 2022, 34: 43–51. DOI: 10.1007/s42241-022-0004-9. [14] YU Z L, AMDAHL J. A numerical solver for coupled dynamic simulation of glacial ice impacts considering hydrodynamic-ice-structure interaction [J]. Ocean Engineering, 2021, 226: 108827. DOI: 10.1016/j.oceaneng.2021.108827. [15] KIM J H, KIM Y. Numerical simulation on the ice-induced fatigue damage of ship structural members in broken ice fields [J]. Marine Structures, 2019, 66: 83–105. DOI: 10.1016/j.marstruc.2019.03.002. [16] 王帅霖, 刘社文, 季顺迎. 基于GPU并行的锥体导管架平台结构冰激振动DEM-FEM耦合分析 [J]. 工程力学, 2019, 36(10): 28–39. DOI: 10.6052/j.issn.1000-4750.2018.10.0560.WANG S L, LIU S W, JI S Y. Coupled discrete-finite element analysis for ice-induced vibration of conical jacket platform based on GPU-BASED parallel algorithm [J]. Engineering Mechanics, 2019, 36(10): 28–39. DOI: 10.6052/j.issn.1000-4750.2018.10.0560. [17] 蒋文灿, 程祥珍, 梁斌, 等. 一种组合药型罩聚能装药战斗部对含水复合结构毁伤的数值模拟及试验研究 [J]. 爆炸与冲击, 2022, 42(8): 083303. DOI: 10.11883/bzycj-2021-0389.JIANG W C, CHENG X Z, LIANG B, et al. Numerical simulation and experimental study on the damage of water partitioned structure by a shaped charge warhead with a combined charge liner [J]. Explosion and Shock Waves, 2022, 42(8): 083303. DOI: 10.11883/bzycj-2021-0389. [18] ERCEG S, ERCEG B, VON BOCK UND POLACH F, et al. A simulation approach for local ice loads on ship structures in level ice [J]. Marine Structures, 2022, 81: 103117. DOI: 10.1016/j.marstruc.2021.103117. [19] LSTC. LS-DYNA keyword user’s manual [M]. CA: Livermore Software Technology Corporation (LSTC), 2014. [20] CHIQUITO. M, CASTEDO. R, SANTOS. A. P, et al. Numerical modelling and experimental validation of the behaviour of brick masonry walls subjected to blast loading [J]. International Journal of Impact Engineering, 2021, 148: 103760. DOI: 10.1016/j. ijimpeng. DOI: 10.1016/j.ijimpeng. [21] 崔堃鹏. 汽车撞击荷载及其作用下高速列车与桥梁系统动力响应与列车运行安全研究 [D]. 北京: 北京交通大学, 2015.CUI K P. Research of motor collision loads and dynamic responses of highspeed train-bridge system and running safety evaluation of trains subjected to motor collision loads [D]. Beijing: Beijing Jiaotong University, 2015. [22] DENG Y, LI C J, LI Z J, et al. Dynamic and full-time acquisition technology and method of ice data of Yellow River [J]. Sensors, 2021, 22(1): 176. DOI: 10.3390/s22010176. [23] SONG M, MA J, HUANG Y. Fluid-structure interaction analysis of ship-ship collisions [J]. Marine Structures, 2017, 55: 121–136. DOI: 10.1016/j.marstruc.2017.05.006. [24] YE X D, FAN W, SHA Y Y, et al. Fluid-structure interaction analysis of oblique ship-bridge collisions [J]. Engineering Structures, 2023, 274: 115129. DOI: 10.1016/j.engstruct.2022.115129. [25] INCE S T, KUMAR A, PAIK J K. A new constitutive equation on ice materials [J]. Ships and Offshore Structures, 2017, 12(5): 610–623. DOI: 10.1080/17445302.2016.1190122. [26] 王庆凯, 张宝森, 邓宇, 等. 黄河冰单轴压缩强度的试验与影响因素探究 [J]. 水利水电技术, 2016, 47(9): 90–94. DOI: 10.13928/j.cnki.wrahe.2016.09.018.WANG Q K, ZHANG B S, DENG Y, et al. Study on test of uniaxial compressive strength of ice in Yellow River and its influencing factors [J]. Water Resources and Hydropower Engineering, 2016, 47(9): 90–94. DOI: 10.13928/j.cnki.wrahe.2016.09.018. [27] 张健, 王甫超, 刘海冬, 等. 水介质对船冰碰撞结构响应的影响 [J]. 船舶工程, 2019, 41(7): 12–15, 22. DOI: 10.13788/j.cnki.cbgc.2019.07.03.ZHANG J, WANG F C, LIU H D, et al. The influence of water medium on the structure response of ship ice collision [J]. Ship Engineering, 2019, 41(7): 12–15, 22. DOI: 10.13788/j.cnki.cbgc.2019.07.03. [28] 王鸿, 贡力, 王忠慧, 等. 基于不同碰撞模型的流冰-输水隧洞碰撞动态响应研究 [J]. 水资源与水工程学报, 2021, 32(1): 164–171. DOI: 10.11705/j.issn.1672-643X.2021.01.24.WANG H, GONG L, WANG Z H, et al. Dynamic response of drift ice-water tunnel collision based on different collision models [J]. Journal of Water Resources and Water Engineering, 2021, 32(1): 164–171. DOI: 10.11705/j.issn.1672-643X.2021.01.24. [29] SONG M, KIM E, AMDAHL J, et al. A comparative analysis of the fluid-structure interaction method and the constant added mass method for ice-structure collisions [J]. Marine Structures, 2016, 49: 58–75. DOI: 10.1016/j.marstruc.2016.05.005. [30] LIU Y Z, SHI W, WANG W H, et al. Dynamic analysis of monopile-type offshore wind turbine under sea ice coupling with fluid-structure interaction [J]. Frontiers in Marine Science, 2022, 9: 12. DOI: 10.3389/fmars.2022.839897. -
下载:
计量
- 文章访问数: 255
- HTML全文浏览量: 67
- PDF下载量: 107
- 被引次数: 0