Loading [MathJax]/jax/output/SVG/jax.js

桩锤作用下岸坡悬挂钢板桩的动力响应特性

任俊刚 高永胜 赵博文 蒋楠

奚小波, 缪宏, 赵荔, 朱河霖, 张瑞宏, 金亦富. 大耕深旋耕刀激光冲击强化残余应力研究[J]. 高压物理学报, 2015, 29(1): 52-58. doi: 10.11858/gywlxb.2015.01.009
引用本文: 任俊刚, 高永胜, 赵博文, 蒋楠. 桩锤作用下岸坡悬挂钢板桩的动力响应特性[J]. 高压物理学报, 2023, 37(6): 064103. doi: 10.11858/gywlxb.20230704
XI Xiao-Bo, MIAO Hong, ZHAO Li, ZHU He-Lin, ZHANG Rui-Hong, JIN Yi-Fu. Study on Residual Stress of Big Tilling Depth Rotary Blade for Laser Shock Peening[J]. Chinese Journal of High Pressure Physics, 2015, 29(1): 52-58. doi: 10.11858/gywlxb.2015.01.009
Citation: REN Jungang, GAO Yongsheng, ZHAO Bowen, JIANG Nan. Dynamic Response Characteristics of Double-Row Suspended Steel Sheet Piles on Nearshore Slope under the Impact of Pile Hammer[J]. Chinese Journal of High Pressure Physics, 2023, 37(6): 064103. doi: 10.11858/gywlxb.20230704

桩锤作用下岸坡悬挂钢板桩的动力响应特性

doi: 10.11858/gywlxb.20230704
详细信息
    作者简介:

    任俊刚(1979-),男,高级工程师,主要从事能源类工程土建结构研究. E-mail:zbwmail@icloud.com

    通讯作者:

    高永胜(1991-),男,工程师,主要从事电气工程及自动化研究. E-mail:1202120576@cug.edu.cn

  • 中图分类号: O347

Dynamic Response Characteristics of Double-Row Suspended Steel Sheet Piles on Nearshore Slope under the Impact of Pile Hammer

  • 摘要: 岸坡锤击管桩施工过程中,为保证岸坡的安全稳定性,准确评估桩锤的动力荷载作用对岸坡支护钢板桩结构的影响非常重要。针对华容煤炭某码头一期工程,采用动力有限元数值模拟方法建立双排悬挂钢板桩结构计算模型;结合现场振动测试验证,分析桩锤振动作用下双排钢板桩的动力响应特性,实现桩锤振动影响下双排钢板桩的安全性评价。结果表明:钢板桩的位移变形与锤击位点相对位置的关系不大,上、下排钢板桩的最大位移分别为2.51和3.14 cm;最大主应力与锤击位点相对位置的关系不大,上、下排钢板桩的最大主应力均小于20 MPa;Mises应力最大值出现在钢板桩与锤击点垂直连线处,上、下排钢板桩的最大Mises应力分别为20.85和25.40 MPa。双排钢板桩的位移变形与应力处于安全控制范围内。

     

  • 旋耕刀是旋耕机作业最主要的受力部件, 也是最易受损的部件, 其性能直接影响旋耕机的作业效率和耕作质量。随着社会生产力的急剧发展, 大马力旋耕机的使用成为农业机械化发展的必然趋势, 而大耕深旋耕刀的加工制造成为关键技术难题。目前, 大耕深旋耕刀的制造工艺为同比例放大法, 即按比例将原有旋耕刀回转半径扩大至350 mm, 这就带来刀面是否加厚的问题。事实证明[1], 当刀面厚度不增加时, 机具耕作100 h内的刀具折断率高达30%以上; 当刀面厚度同比例增加时, 机具耕作100 h内的刀具折断率仍在20%左右, 且两种情况下后者的能耗增加了15%。为节省能耗、减少刀具材料使用, 不宜增加刀面厚度, 通过合适的表面处理提高刀具性能是一种有效方法[2-5]

    激光冲击强化(Laser Shock Peening, LSP)技术是一种新型的表面技术, 目前已广泛应用于材料表面改性的研究。激光冲击强化采用几十纳秒的短脉冲高峰值功率密度(> 109 W/cm2)激光辐射向金属表层, 使金属表面涂覆的保护层吸收激光能量并发生爆炸性气化蒸发, 产生大于1 GPa高压等离子体冲击波, 利用冲击波的力效应使表层材料微观组织发生变化引入残余压应力, 从而提高金属材料抗疲劳、耐磨损和防应力腐蚀等性能[6-8]

    旋耕刀切土时受到的土壤阻力大多集中在刀具的正切削刃上, 会导致旋耕刀刀柄处产生较大的集中应力, 刀柄在循环交变应力的作用下极易形成疲劳裂纹并失效[9-12]。为解决这一问题, 可在刀柄应力集中部位引入残余压应力以抵消旋耕刀作业时的外界应力, 抑制或减缓疲劳裂纹的产生, 提高刀具使用寿命。本研究借助ANSYS分析出旋耕刀应力集中区域, 并模拟激光冲击试验得出残余应力引入情况, 同时采用激光冲击强化技术处理应力集中部位, 利用X射线衍射法对残余应力进行测试, 讨论激光冲击强化处理后刀面的残余应力的影响。

    旋耕机作业时, 旋转刀对土壤进行切削、破碎及抛掷, 土壤对旋耕刀的反作用力构成了土壤阻力[13]。刀具在切削土壤时, 其耕作深度及切削土壤面均先由小到大, 后由大到小, 其在土壤中的位置也不断变化, 所以刀具受到土壤阻力的大小、方向和作用点在其切削土壤时是不断变化的。在切土过程中, 刀具基本不受轴向力的作用, 可将阻力F沿xz两个坐标轴分解成FxFz两分力, 则有

    {F=Fx+FzF=F2x+F2z (1)

    旋耕刀的回转半径为350 mm, 耕作深度为21 cm, 刀面切土宽度为6 cm, 土壤阻力计算公式为

    Rx=0.7abK (2)

    式中:Rx为土壤阻力, N;a为刀面切土宽度, cm;b为耕作深度, cm;K为土壤比阻, N·cm-2。一般土壤的K为4~5 N·cm-2, 粘土的K为6~8 N·cm-2, 本研究K取7 N·cm-2, 根据公式(2)可求出旋耕刀所受土壤阻力Rx约为620 N。

    在PRO/E 5.0中建立回转半径为350 mm、刀面厚度为7 mm的旋耕刀三维模型, 并将模型导入ANSYS中进行有限元分析。65Mn旋耕刀的弹性模量E为208 GPa, 泊松比n为0.288。采用Solid 168单元将模型自动划分网格, 共计86 021个节点、53 902个单元。在旋耕刀柄两个侧面及顶面与连接孔处施加约束, 并在受力处施加620 N的集中载荷, 最后运算求解。

    图 1(a)为旋耕刀的位移分布等值线图, 工作时旋耕刀的最大位移可以达到0.508 mm, 发生在旋耕刀切削土壤时的正切削刃和侧切削刃处, 而且位移变形是向两边等距递减的, 由于旋耕刀柄固定, 切土时正切削刃和侧切削刃受力, 因而此处变形较大。图 1(b)为旋耕刀的应力分布等值线图, 旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 旋耕刀工作时该处极易折断。

    图  1  (a) 位移分布等值线图(b)应力分布等值线图
    Figure  1.  (a) Isoline of displacement distribution (b) Isoline of stress distribution

    旋耕刀采用65Mn为材料, 其主要化学成分为(质量分数, %):C 0.65, Si 0.26, Mn 1.10, S 0.002, P 0.019, Cr 0.02, Ni 0.01, Cu 0.01, Fe余量。材料在激光冲击过程中产生动态变形, 因此在选择材料应力-应变本构方程时, 静态方程已经不能反映材料的真实响应, 必须用动态本构方程替代。一般材料的动态本构方程, 目前较多采用Johnson-Cook方程[14], 其本构方程为

    σ=(A+Bεn)[1+Cln(˙ε˙ε0)][1(θθrθmθr)m] (3)

    式中:εθ分别为应变和温度为应变率分别为参考应变率和参考温度; θm为材料熔点; ABnCm为待定系数, ABn表征材料应变强化项系数, C表征材料应变率强化项系数, m表征材料热软化系数。

    段志勇等[15]在实验基础上, 提出了一个与实际比较相符的半经验型激光冲击波估计模型公式

    pmax=0.8ρI0 (4)

    式中:pmax为激光峰值压力, GPa; ρ为折合密度, g/cm3; I0为激光功率密度, GW/cm2。折合密度公式为

    2ρ=1ρ1+1ρ2 (5)

    式中:ρ1为约束层密度, g/cm3; ρ2为靶材密度, g/cm3。激光功率密度公式为

    I0=Eπr2τ (6)

    式中:E为激光能量, J;τ为激光脉宽, ns;r为光斑半径, mm。

    已知水密度ρ1=1 g/cm3, 靶材密度ρ2=7.85 g/cm3, 算得折合密度ρ=1.774 g/cm3; 激光能量E=10 J, 脉冲宽度τ=30 ns, 光斑半径r=3 mm, 算得激光功率密度I0=1.179 GW/cm2, 激光峰值压力pmax=1.16 GPa。激光冲击波作用时间为激光脉宽的3倍左右甚至更高[16], 实验采用的激光脉宽为30 ns, 取激光冲击波作用时间为100 ns, 其冲击波作用压力与时间的关系如图 2所示, 在50 ns时的压力最大, 最大压力值为1.16 GPa。

    图  2  激光冲击波加载曲线
    Figure  2.  Shock wave of laser loading curve

    采用Solid 70单元划分网格并定义材料属性, 设定环境温度为22 ℃, 设定对流及辐射边界条件, 并在应力集中区域的旋耕刀另一侧面设定绝热边界条件。在应力集中部位施加高斯表面热源, 由于区域面积较小, 高斯热源固定施加在待冲击区域中心, 设定载荷步选项和分析选项运算求解, 得到温度场模型。

    进行应力分析时, 通过前处理模块定义温度场, 按图 2中激光冲击波加载曲线定义载荷函数, 采用Solid 45单元重新划分网格并定义材料属性, 在应力集中区域定义加载单元组元, 将定义的冲击压力载荷施加在定义的单元组元。同时在单元组元对应的模型各个表面施加透射条件将应力波透射, 防止过大的应力波使模型发生弯扭变形影响仿真结果。最后施加约束并设置求解选项进行运算, 输出表面残余应力结果。

    激光冲击强化前, 将激光冲击强化区域进行抛光, 利用超声波清洗器将试样在乙醇中清洗15 min, 除去表面油污等杂质, 最后烘干试样[17]。采用铝箔作为激光冲击的金属涂敷层, 铝箔涂层长为25 mm, 宽为10 mm, 厚为0.1 mm。将铝箔粘贴在旋耕刀应力最集中的刀柄外弯角处, 如图 3(a)所示。使用专用夹具将贴有铝箔的旋耕刀夹于激光冲击强化处理工作台。用水作约束层, 水流厚度为0.5 mm。试验采用江苏大学Gaia-R系列高能量脉冲灯抽运YAG激光器, 激光器调整冲击参数为:激光能量10 J, 波长1.06 μm, 脉冲宽度30 ns, 频率0.1 Hz, 光斑直径6 mm。

    图  3  (a) 旋耕刀冲击强化前准备(b)残余应力测点位置
    Figure  3.  (a) The preparation before LSP (b) Residual stress measuring point location

    利用超声波清洗器将试样在丙酮溶液中清洗15 min, 并用酒精反复冲洗, 最后吹干[18-19]。采用MSF-3M型X射线应力测定仪进行残余应力测量, 应力测试取点如图 3(b)所示。X射线分析测试条件为V靶材, 所用衍射晶面为α-Fe(211)衍射晶面, 交相关方法定峰[20]

    从ANSYS后处理输出各测试点的表面残余应力仿真值, 对比试验实测数据, 如图 4所示。各测点的仿真值与实测值曲线呈上下交替分布, 曲线拟合较好, 误差在±20 MPa内。仿真值呈下凹状态, 主要是因为ANSYS中施加的高斯热源为固定热源, 且集中施加在冲击区域的中心, 造成中心温度向边缘扩散, 结果表现为整个区域中心的残余压应力明显高于周边。同理, 实测值呈现下凹状态也是因为激光冲击时, 光斑对准区域中心, 中心受热过大, 从而残余应力的引入明显多于周边。

    图  4  激光冲击仿真与实测残余应力
    Figure  4.  Residual stress datum of simulation and experiment after LSP

    利用ANSYS模拟激光冲击强化, 其在材料表层引入的残余应力与试验实测的结果吻合得较好, 说明ANSYS可完成激光冲击强化试验, 且误差在允许范围内, 通过此类仿真分析研究, 可对激光冲击强化的工艺参数进行优化, 进一步提高材料性能。

    图 5(a)、图 5(b)分别是激光冲击前后旋耕刀表面材料残余应力沿xy轴的分布情况, 未经强化的平均残余压应力仅为146.90 MPa, 强化后的平均残余压应力达到390.70 MPa, 提高了166%。经激光冲击强化后, 最大残余压应力达到412.25 MPa, 在冲击区域的中心位置; 最小残余压应力为375.82 MPa, 在冲击区域的最边缘位置。同时, 冲击后沿x轴及y轴方向的残余应力误差大小在±20 MPa内, 说明整个冲击强化区域的残余应力分布均匀, 无应力集中现象, 说明激光冲击强化能有效提高材料表层残余应力。

    图  5  冲击前、后的残余应力
    Figure  5.  Residual stress before and after LSP

    由于旋耕刀工作环境恶劣, 且工作时间较长, 易存在应力集中, 引发疲劳裂纹, 激光冲击强化处理引入的残余压应力可消除应力集中的影响, 减慢或抑制裂纹的扩展, 能起到提高板料疲劳寿命的作用, 进而提高旋耕刀使用寿命。

    (1) 利用ANSYS对旋耕刀进行应力分析, 得出旋耕刀柄外弯角处应力最为集中, 最大拉应力达到303 MPa, 与日常刀具易受损折断处的位置基本吻合。

    (2) ANSYS可模拟激光冲击强化工艺, 其算得的残余应力引入值与试验实测的结果吻合很好, 误差在±20 MPa内, 通过仿真分析可优化相关工艺参数, 进一步提高材料性能。

    (3) 经激光冲击强化处理的旋耕刀表层材料残余压应力明显增大, 最大残余压应力达412.25 MPa, 可消除旋耕刀在恶劣工况下应力集中的影响, 减慢或抑制裂纹的扩展, 提高材料疲劳寿命。

  • 图  拟定边坡开挖及支护方案

    Figure  1.  Proposed slope excavation and support plan

    图  现场测试点布置

    Figure  2.  Field test point layout

    图  数值计算模型

    Figure  3.  Numerical calculation model

    图  锤击力荷载曲线

    Figure  4.  Hammer loading force curve

    图  监测点P1在x方向的振动速度对比

    Figure  5.  x-axis velocity comparison at monitoring point P1

    图  监测点P1的合振动速度对比

    Figure  6.  Resultant velocity comparison at monitoring point P1

    图  位移和应力测点的布置

    Figure  7.  Layout of displacement and stress measurement points

    图  钢板桩结构的峰值位移分布

    Figure  8.  Peak displacement distribution of steel sheet pile structure

    图  典型测点的最大主应力

    Figure  9.  Maximum principal stress at typical measuring points

    图  10  Mises应力沿水平方向的分布

    Figure  10.  Mises stress distribution along the horizontal direction

    表  1  现场振动测试数据

    Table  1.   Field vibration test data

    No. Test No. PPV/(cm·s−1) N Distance/m
    1 2.000 7 13
    2 0.632 3 21
    3 0.158 6 38
    1 2.045 7 13
    2 0.644 2 21
    3 0.160 6 38
    1 2.025 7 13
    2 0.639 3 21
    3 0.160 3 38
    下载: 导出CSV

    表  2  数值模拟参数

    Table  2.   Numerical simulation parameters

    Material ρ/(kg·m−3) E0/GPa μ A/MPa B/MPa n C σ0/GPa
    Steel 7830 205 0.30 792 510 0.26 0.014
    Clay 1900 11 0.20 0.19, 0.15
    Mudstone 2230 15 0.22 0.30
    下载: 导出CSV

    表  3  质点峰值振动速度对比

    Table  3.   Comparison of PPVs

    No.Burst distance/mTest PPV/(cm·s−1)Simulated PPV/(cm·s−1)Relative error/%
    1132.0231.8677.7
    2210.6380.5769.7
    3380.1590.20528.9
    下载: 导出CSV
  • [1] 吕鹏, 董胜, 焉振, 等. U型钢板桩横向承载性能模型试验研究[J]. 海洋湖沼通报, 2023, 45(1): 1−9.

    LYU P, DONG S, YAN Z, et al. Experimental study on lateral bearing performance model of U-shaped steel sheet piles [J]. Marine and Lake Sciences Bulletin, 2023, 45(1): 1–9.
    [2] 焉振, 王元战. 考虑软基不排水强度循环弱化的格型钢板桩防波堤动力有限元分析 [J]. 岩土力学, 2017, 38(5): 1454–1462.

    YAN Z, WANG Y Z. Dynamic finite element analysis of geogrid steel sheet pile breakwater considering undrained strength degradation of soft foundation [J]. Rock and Soil Mechanics, 2017, 38(5): 1454–1462.
    [3] 郑国兵, 黄朝煊, 袁文喜, 等. 深厚淤泥中双排钢板桩结构安全稳定性研究 [J]. 水利水电科技进展, 2020, 40(4): 71–76.

    ZHENG G B, HUANG C X, YUAN W X, et al. Research on the safety and stability of double-row steel sheet pile structure in deep and thick silt [J]. Advances in Science and Technology of Water Resources and Hydropower Engineering, 2020, 40(4): 71–76.
    [4] 杨熠. 双排钢板桩围堰结构受力变形及稳定性分析 [D]. 南昌: 南昌大学, 2020.

    YANG Y. Analysis of stress, deformation, and stability of double-row steel sheet pile caisson structure [D]. Nanchang: Nanchang University, 2020.
    [5] 曹胜敏. 高桩码头桩竖向荷载下静动力学行为研究 [D]. 成都: 西南交通大学, 2008.

    CAO S M. Study on the static and dynamic behavior of high pile wharf piles under vertical loads [D]. Chengdu: Southwest Jiaotong University, 2008.
    [6] 王娜娜. 打桩下沉过程中的动力响应研究 [D]. 天津: 天津大学, 2012.

    WANG N N. Study on the dynamic response during pile driving process [D]. Tianjin: Tianjin University, 2012.
    [7] XIA Y, JIANG N, ZHOU C, et al. Dynamic behaviors of buried reinforced concrete pipelines with gasketed bell-and-spigot joints subjected to tunnel blasting vibration [J]. Tunnelling and Underground Space Technology, 2021, 118: 104172. doi: 10.1016/j.tust.2021.104172
    [8] JIANG N, ZHU B, ZHOU C, et al. Blasting vibration effect on the buried pipeline: a brief overview [J]. Engineering Failure Analysis, 2021, 129: 105709. doi: 10.1016/j.engfailanal.2021.105709
    [9] HALLQUIST J O. LS-DYNA theory manual [M]. Livermore, USA: Livermore Software Technology Corporation, 2006.
    [10] HALLQUIST J O. LS-DYNA keyword user’s manual version 971 [M]. Livermore, USA: Livermore Software Technology Corporation, 2007.
    [11] ZHOU Y, ZHANG S. Rupture and perforation responses of pressurized tubular members subjected to medium-velocity transverse impact loading [J]. Engineering Failure Analysis, 2021, 127: 105387. doi: 10.1016/j.engfailanal.2021.105387
    [12] YAO S, ZHANG D, LU Z, et al. Experimental and numerical investigation on the dynamic response of steel chamber under internal blast [J]. Engineering Structures, 2018, 168: 877–888. doi: 10.1016/j.engstruct.2018.03.067
    [13] 邓灵敏. 基于LS-DYNA的动力沉桩全过程分析 [D]. 武汉: 武汉理工大学, 2013.

    DENG L M. Dynamic analysis of the entire process of pile driving based on LS-DYNA [D]. Wuhan: Wuhan University of Technology, 2013.
    [14] JIANG N, ZHU B, HE X, et al. Safety assessment of buried pressurized gas pipelines subject to blasting vibrations induced by metro foundation pit excavation [J]. Tunnelling and Underground Space Technology, 2020, 102: 103448. doi: 10.1016/j.tust.2020.103448
    [15] 杜闯, 丁红岩, 张浦阳, 等. 钢板桩围堰有限元分析 [J]. 岩土工程学报, 2014, 36(Suppl 2): 159–164.

    DU C, DING H Y, ZHANG P Y, et al. Finite element analysis of steel sheet pile caisson [J]. Chinese Journal of Geotechnical Engineering, 2014, 36(Suppl 2): 159–164.
  • 加载中
图(10) / 表(3)
计量
  • 文章访问数:  144
  • HTML全文浏览量:  68
  • PDF下载量:  17
出版历程
  • 收稿日期:  2023-08-08
  • 修回日期:  2023-09-07
  • 网络出版日期:  2023-12-07
  • 刊出日期:  2023-12-15

目录

/

返回文章
返回