• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

基于荷载时程分析法的钢板混凝土结构墙的抗冲击性能敏感性分析

朱秀云 林皋 潘蓉 路雨

引用本文:
Citation:

基于荷载时程分析法的钢板混凝土结构墙的抗冲击性能敏感性分析

    作者简介: 朱秀云(1985—),女,博士研究生, 高级工程师, lyzhuxiuyun@163.com;
  • 基金项目: 国家科技重大专项基金项目 2011ZX06002-10

  • 中图分类号: O347.1;TL371

Sensitivity analysis for impact resistance of steel plate concrete walls based on force vs. time-history analysis

  • CLC number: O347.1;TL371

  • 摘要: 选取由表面钢板、拉结筋、剪力钉及混凝土组成的钢板混凝土结构墙为研究对象,运用经典显式非线性动力分析软件ANSYS/LS-DYNA,基于荷载时程分析法进行了一系列影响钢板混凝土结构墙抗冲击性能的参数敏感性分析。这些参数包括:墙体厚度、钢板厚度、拉结筋直径与间距等。分析结果表明,以上参数均会影响墙体的抗冲击性能,尤其是墙体与钢板的厚度以及拉结筋的间距。本文的研究工作对于核电厂核岛厂房钢板混凝土结构外墙的设计具有一定的指导与参考意义。
  • 图 1  波音707-320的撞击荷载时程曲线

    Figure 1.  Curve of impact force vs. time-history

    图 2  波音707-320的撞击面积时程曲线

    Figure 2.  Curve of impact area vs. time-history

    图 3  机身和机翼的撞击荷载时程曲线

    Figure 3.  Curves of impact force vs. time-history for fuselage and wings

    图 4  机身和机翼的撞击面积时程曲线

    Figure 4.  Curves of impact area vs. time-history for fuselage and wings

    图 5  波音707-320机身和机翼的均布荷载时程曲线

    Figure 5.  Curves of evenly distributed impact pressure vs. time-history for fuselage and wings of Boeing 707-320

    图 6  波音707-320整个飞机的均布荷载时程曲线

    Figure 6.  Curve of evenly distributed impact pressure vs. time-history for the whole aircraft of Boeing 707-320

    图 7  典型的钢板混凝土结构墙示意图

    Figure 7.  Schematic illustration of the typical steel plate concrete (SC) wall

    图 8  钢板混凝土结构墙的有限元模型

    Figure 8.  FEM model of the steel plate concrete wall

    图 9  不同加载模式的墙体最大位移时程曲线

    Figure 9.  Curves of max displacement vs. time-history of SC wall for different loads

    图 10  不同加载模式的混凝土最大塑性应变时程曲线

    Figure 10.  Curves of max plastic strain vs. time-history of concrete for different loads

    图 11  不同加载模式下混凝土的最大塑性应变云图

    Figure 11.  Contour plot of max plastic strain of concrete for different loads

    图 12  不同厚度SC墙体的最大位移时程曲线

    Figure 12.  Curves of max displacement vs. time-history for SC wall with different thickness

    图 13  不同厚度SC墙体的混凝土最大塑性应变时程曲线

    Figure 13.  Curves of max plastic strain vs. time-history for concrete of SC wall with different thickness

    图 14  不同钢板厚度的墙体最大位移时程曲线

    Figure 14.  Curves of max displacement vs. time-history for SC wall of steel plates with different thickness

    图 15  不同钢板厚度的混凝土最大塑性应变时程曲线

    Figure 15.  Curves of max plastic strain vs. time-history for concrete for steel plates with different thickness

    图 16  不同钢板厚度的墙体最大位移时程曲线

    Figure 16.  Curves of max displacement vs. time-history for SC wall of steel plates with different thickness

    图 17  不同钢板厚度的混凝土最大塑性应变时程曲线

    Figure 17.  Curves of max plastic strain vs. time-history for concrete of steel plates with different thickness

    图 18  SC墙体14 mm厚背部钢板的等效应力云图

    Figure 18.  Contour plot of the equivalent stress of back face of 14mm steel plate for SC wall

    图 19  SC墙体18 mm厚背部钢板的等效应力云图

    Figure 19.  Contour plot of the equivalent stress of back face of 18mm steel plate for SC wall

    图 20  不同拉结筋直径的墙体最大位移时程曲线

    Figure 20.  Curves of max displacement vs. time-history for SC wall with different diameter of tie-bar

    图 21  混凝土的最大位移云图

    Figure 21.  Contour plot of max displacement of concrete

    图 22  SC墙体背部的残余位移云图

    Figure 22.  Contour plot of residual displacement of SC wall back surface

    图 23  SC墙体钢板的最大塑性应变云图

    Figure 23.  Contour plot of max plastic strain of SC wall steel plate

    图 24  不同拉结筋间距的墙体最大位移时程曲线

    Figure 24.  Curves of max displacement vs. time-history for differently spaced tie-bar

    表 1  钢板混凝土墙的材料参数

    Table 1.  Material parameters of steel plate concrete wall

    材料 ρ/(kg·m-3) E/GPa υ fc/MPa ft/MPa
    混凝土 2 400 35.2 0.2 46.0 4.2
    钢板 7 800 206.0 0.3 400.0 400.0
    拉结筋 7 800 206.0 0.3 330.0 330.0
    剪力钉 7 800 206.0 0.3 330.0 330.0
    下载: 导出CSV
  • [1] US Nuclear Regulatory Commission.10CFR50.150 Aircraft impact assessment[S]. Washington, DC: US Nuclear Regulatory Commission, 2009.
    [2] ERIN Engineering and Research, Inc. NEI 07-13, Revision 8P, methodology for performing aircraft impact assessments for new plant designs[S]. Palo Alto, 2011.
    [3] Tsubota H, Koshika N, Mizuno J, et al. Scale model tests of multiple barriers against aircraft impact: Part 1. Experimental program and test results[C]//Transactions of the 15th International Conference on Structural Mechanics in Reactor Technology (SMiRT-15). Seoul, Korea, 1999: Ⅶ137-Ⅶ144.
    [4] Mizuno J, Koshika N, Sawamoto Y, et al. Investigations on impact resistance of steel plate reinforced concrete barriers against aircraft impact part 1: Test program and results[C]//Transactions of the 18th International Conference on Structural Mechanics in Reactor Technology, 2005: 2566-2579.
    [5] Muhammad S, Zhu X Y, Pan R. Simulation analysis of impact tests of steel plate reinforced concrete and reinforced concrete slabs against aircraft impact and its validation with experimental results[J]. Nuclear Engineering and Design, 2014, 273:653-667. doi: 10.1016/j.nucengdes.2014.03.031
    [6] 朱秀云, 潘蓉, 林皋, 等.基于荷载时程分析法的钢筋混凝土和钢板混凝土墙的冲击响应对比分析[J].振动与冲击, 2014, 33(22):172-177.
    Zhu Xiuyun, Pan Rong, Lin Gao, et al. Comparative analysis of impact response with reinforced concrete and steel plate concrete walls based on force time-history analysis method[J]. Journal of Vibration and Shock, 2014, 33(22):172-177.
    [7] 潘蓉, 吴婧姝, 张心斌.钢板混凝土结构在核电工程中应用的发展状况[J].工业建筑.2014, 44(12):1-7.
    Pan Rong, Wu Jingshu, Zhang Xinbin. Application and development of steel plate reinforced concrete structure in nuclear power engineering[J]. Industrial Construction, 2014, 44(12):1-7.
    [8] Hallquist J Q. LS-DYNA keyword user's manual, revision 971[M]. California: Livermore Software Technology Corportion, 2007.
    [9] 中国国家核安全局.HAD101/04.核电厂厂址选择的外部人为事件[S].北京: 中国法制出版社, 1989.
    [10] Mullapudi T R S, Summers P, Hwan M. Impact analysis of steel plated concrete wall[C]//Structures Congress 2012.ASCE, 2012: 1881-1893.
    [11] Arros J, Doumbalski N. Analysis of aircraft impact to concrete structures[J]. Nuclear Engineering and Design, 2007, 237(12/13):1241-1249.
    [12] Abu-Odeh A. Modeling and simulation of bogie impacts on concrete bridge rails using LS-DYNA[C]//10th International LS-DYNA Users Conference. 2008.
    [13] Comite Euro-International du Beton. CEB-FIP model code 1990[M]. Trowbridge, Wiltshire, U K: Redwood Books, 1993.
    [14] Wu Y C, Crawford J E, Magallanes J M. Performance of LS-DYNA concrete constitutive models[C]//12th International LS-DYNA Users Conference, 2012.
    [15] 朱秀云, 潘蓉, 林皋, 等.基于ANSYS/LS-DYNA的钢板混凝土墙冲击实验的有限元分析[J].爆炸与冲击, 2015, 35(2):222-228.
    Zhu Xiuyun, Pan Rong, Lin Gao, et al. FEM analysis of impact experiments with steel plate concrete walls based on ANSYS/LS-DYNA[J]. Explosion and Shock Waves, 2015, 35(2):222-228.
    [16] Mizuno J, Koshika N, Morikawa H, et al. Investigations on impact resistance of steel plate reinforced concrete barriers against aircraft impact. Part 2: Simulation analyses of scale model impact tests[C]//Transactions of the 18th International Conference on Structural Mechanics in Reactor Technology, 2005: 2580-2590.
  • [1] 张振华钱海峰王媛欣牟金磊梅志远牛闯 . 球头落锤冲击下金字塔点阵夹芯板结构的动态响应实验. 爆炸与冲击, 2015, 35(6): 888-894. doi: 10.11883/1001-1455(2015)06-0888-07
    [2] 李砚召王肖钧郭晓辉曹海赵凯 . 部分预应力混凝土梁抗冲击性能试验研究. 爆炸与冲击, 2006, 26(3): 256-261. doi: 10.11883/1001-1455(2006)03-0256-06
    [3] 吴松林杜扬李国庆张培理 . C1~C4烷烃混合气体热爆燃的简化机理与分析. 爆炸与冲击, 2015, 35(5): 641-650. doi: 10.11883/1001-1455(2015)05-0641-10
    [4] 陈鹏程程欣弓磊路国运 . 不同宽厚比H形截面钢梁抗冲击性能. 爆炸与冲击, 2018, 38(6): 1378-1385. doi: 10.11883/bzycj-2017-0125
    [5] 赖建中孙伟戎志丹 . 活性粉末混凝土在多次冲击荷载下的力学行为. 爆炸与冲击, 2008, 28(6): 532-538. doi: 10.11883/1001-1455(2008)06-0532-07
    [6] 姜世平于海龙芮筱亭洪俊黎超 . 散体系统冲击破碎的动力学分析. 爆炸与冲击, 2014, 34(2): 247-251. doi: 10.11883/1001-1455(2014)02-0247-05
    [7] 宋延泽李志强赵隆茂 . 基于十四面体模型的闭孔泡沫材料动态力学性能的有限元分析. 爆炸与冲击, 2009, 29(1): 49-55. doi: 10.11883/1001-1455(2009)01-0049-07
    [8] 孙伟王建国姜绍飞 . 减小高速列车荷载对地面结构损害的数值分析. 爆炸与冲击, 2009, 29(2): 155-161. doi: 10.11883/1001-1455(2009)02-0155-07
    [9] 朱秀云潘蓉林皋李亮 . 基于ANSYS/LS-DYNA的钢板混凝土墙冲击实验的有限元分析. 爆炸与冲击, 2015, 35(2): 222-228. doi: 10.11883/1001-1455(2015)02-0222-07
    [10] 王志华马宏伟赵隆茂杨桂通 . 反分析法在泡沫金属材料动态性能实验中的应用. 爆炸与冲击, 2006, 26(5): 410-415. doi: 10.11883/1001-1455(2006)05-0410-06
    [11] 姚熊亮戴绍仕周其新张阿漫闻雪友 , . 船体与设备一体化抗冲击分析. 爆炸与冲击, 2009, 29(4): 367-374. doi: 10.11883/1001-1455(2009)04-0367-08
    [12] 刘逸平黄小清汤立群罗立峰 . 新型混凝土桥面铺装材料的冲击力学性能. 爆炸与冲击, 2007, 27(3): 217-222. doi: 10.11883/1001-1455(2007)03-0217-06
    [13] 范飞林许金余李为民杨进勇翟毅 . 矿渣-粉煤灰基地质聚合物混凝土的冲击力学性能. 爆炸与冲击, 2009, 29(5): 516-522. doi: 10.11883/1001-1455(2009)05-0516-07
    [14] 陶俊林秦李波李奎刘丹贾彬陈小伟陈刚 . 混凝土高温动态压缩力学性能实验. 爆炸与冲击, 2011, 31(1): 101-106. doi: 10.11883/1001-1455(2011)01-0101-06
    [15] 任兴涛周听清钟方平胡永乐王万鹏 . 钢纤维活性粉末混凝土的动态力学性能. 爆炸与冲击, 2011, 31(5): 540-547. doi: 10.11883/1001-1455(2011)05-0540-08
    [16] 王斌李夕兵 . 单轴荷载下饱水岩石静态和动态抗压强度的细观力学分析. 爆炸与冲击, 2012, 32(4): 423-431. doi: 10.11883/1001-1455(2012)04-0423-09
    [17] 王吉王肖钧卞梁 . 光滑粒子法与有限元的耦合算法及其在冲击动力学中的应用. 爆炸与冲击, 2007, 27(6): 522-528. doi: 10.11883/1001-1455(2007)06-0522-07
    [18] 李本平 . 冲击荷载作用下粘土孔隙水压的数值分析. 爆炸与冲击, 2005, 25(3): 281-284. doi: 10.11883/1001-1455(2005)03-0281-04
    [19] 张延昌杨代玉王自力 . 舱内液体对VLCC舷侧结构碰撞性能的影响. 爆炸与冲击, 2010, 30(5): 479-486. doi: 10.11883/1001-1455(2010)05-0479-08
    [20] 宋顺成王庭辉才鸿年王富耻 . 钢制弹丸冲击混凝土时绝热剪切局部化的数值计算. 爆炸与冲击, 2007, 27(3): 246-250. doi: 10.11883/1001-1455(2007)03-0246-05
  • 加载中
图(24)表(1)
计量
  • 文章访问数:  625
  • HTML全文浏览量:  183
  • PDF下载量:  816
  • 被引次数: 0
出版历程
  • 收稿日期:  2015-02-04
  • 录用日期:  2015-04-21
  • 刊出日期:  2016-09-25

基于荷载时程分析法的钢板混凝土结构墙的抗冲击性能敏感性分析

    作者简介:朱秀云(1985—),女,博士研究生, 高级工程师, lyzhuxiuyun@163.com
  • 1. 大连理工大学海岸与近海国家重点实验室抗震分室,辽宁 大连 116024
  • 2. 环境保护部核与辐射安全中心厂址与土建部, 北京 100082
基金项目:  国家科技重大专项基金项目 2011ZX06002-10

摘要: 选取由表面钢板、拉结筋、剪力钉及混凝土组成的钢板混凝土结构墙为研究对象,运用经典显式非线性动力分析软件ANSYS/LS-DYNA,基于荷载时程分析法进行了一系列影响钢板混凝土结构墙抗冲击性能的参数敏感性分析。这些参数包括:墙体厚度、钢板厚度、拉结筋直径与间距等。分析结果表明,以上参数均会影响墙体的抗冲击性能,尤其是墙体与钢板的厚度以及拉结筋的间距。本文的研究工作对于核电厂核岛厂房钢板混凝土结构外墙的设计具有一定的指导与参考意义。

English Abstract

  • 9.11事件后,核电厂在大型商用飞机撞击下的安全问题成为公众的焦点。在美国, 将大型商用飞机撞击作为核电厂一种超设计基准事件,以美国联邦法规的新条款10CFR50.150[1]形式颁布, 要求新设计的核动力堆需就抵御大型商用飞机恶意撞击进行专门的评价。其中,文件NEI07-13[2]提供了美国电力研究院发展的一套评价大型商用飞机撞击的方法,对于构筑物整体破坏的评估,推荐了荷载时程分析法和飞射物-靶体相互作用分析法。在我国,核安全法规中还没有针对核电厂抗大型商用飞机撞击评估提出相关要求,但从国际核行业安全评价的发展趋势看,掌握和发展核电厂构筑物抵御大型商用飞机撞击技术,建立和完善相关法规,具有重大意义。

    关于钢板混凝土结构与钢筋混凝土结构的抗冲击性能对比,日本Kobori综合研究所分别对1/7.5缩尺飞机模型垂直撞击不同厚度钢筋混凝土(RC)墙[3]及钢板混凝土(SC)墙[4]进行了实验研究;M.Sadiq等[5]和朱秀云等[6]基于非线性有限元软件ANSYS/LS-DYNA对以上实验分别进行了飞射物-靶体相互作用分析法和荷载时程分析法的有限元数值模拟,并对不同厚度的钢筋混凝土(RC)墙与钢板混凝土(SC)墙的抗冲击性能进行了一系列的对比分析,其结论是钢板混凝土结构墙的抗冲击性能优于钢筋混凝土结构墙,尤其是背面钢板能够有效地约束混凝土在撞击方向上的运动,并限制混凝土碎片的飞溅,因此,用于抗飞机撞击的钢板混凝土结构墙体的厚度可以较钢筋混凝土结构适当减薄。针对大型商用飞机恶意撞击事件的补充考虑,美国西屋公司将AP1000屏蔽厂房结构类型由美国核管会批准的DCD第15版钢筋混凝土结构变更为DCD第16版的钢板混凝土结构;国内新设计的大型先进压水堆CAP1400为了抵抗大型商用飞机的撞击,同样将屏蔽厂房设计为钢板混凝土结构。当然,潘蓉等[7]指出钢板混凝土结构由于钢板暴露在外,没有混凝土面层的保护,受火时比钢筋混凝土结构更易于受到影响而丧失强度,需要进一步研究其防火方法。所以在抗商用飞机撞击分析中,一般只允许钢板混凝土结构墙体发生少量变形,不允许被穿透以防止飞机燃油燃烧对钢板混凝土结构墙体造成破坏影响。

    为了探讨影响钢板混凝土结构墙体抗冲击性能的参数,本文中基于非线性有限元动力分析软件ANSYS/LS-DYNA[8],采用HAD101/04《核电厂厂址选择的外部人为事件》[9]提供的波音707-320撞击荷载函数,进行一系列影响钢板混凝土结构墙体抗冲击性能的参数敏感性研究,其参数包括:混凝土厚度、钢板厚度、拉结筋直径与间距等。

    • 对于荷载时程分析法,首先是确定飞机撞击荷载时程函数。撞击荷载选用导则HAD101/04《核电厂厂址选择的外部人为事件》[9]附录Ⅰ中列出的波音707-320在典型起落速度370 km/h的荷载时程曲线,如图 1所示;撞击面积时程曲线,如图 2所示。可考虑2种不同的均布荷载时间函数加载模式。第1种模式是根据商用飞机的特征,其撞击作用分为机身和机翼部分,其机身和机翼的撞击荷载及撞击面积时程曲线分别如图 3~4所示;由此相同时刻对应的撞击荷载除以撞击面积,可得到分别作用于墙体的机身与机翼区域的均布荷载时程曲线,如图 5所示。第2种模式是如图 1~2所示的撞击荷载和撞击面积时程,两者直接相除得到作用于墙体整个飞机区域的一致均布荷载时间函数曲线,如图 6所示。

      图  1  波音707-320的撞击荷载时程曲线

      Figure 1.  Curve of impact force vs. time-history

      图  2  波音707-320的撞击面积时程曲线

      Figure 2.  Curve of impact area vs. time-history

      图  3  机身和机翼的撞击荷载时程曲线

      Figure 3.  Curves of impact force vs. time-history for fuselage and wings

      图  4  机身和机翼的撞击面积时程曲线

      Figure 4.  Curves of impact area vs. time-history for fuselage and wings

      图  5  波音707-320机身和机翼的均布荷载时程曲线

      Figure 5.  Curves of evenly distributed impact pressure vs. time-history for fuselage and wings of Boeing 707-320

      图  6  波音707-320整个飞机的均布荷载时程曲线

      Figure 6.  Curve of evenly distributed impact pressure vs. time-history for the whole aircraft of Boeing 707-320

    • 钢板混凝土结构(SC)墙一般由表面钢板、拉结筋、剪力钉以及混凝土等部分组成,其典型的结构形式示意如图 7所示。

      图  7  典型的钢板混凝土结构墙示意图

      Figure 7.  Schematic illustration of the typical steel plate concrete (SC) wall

      研究对象为长36 m、高10 m、4个侧立面固定约束的钢板混凝土结构墙体,厚度为0.75~1.0 m,内外层钢板厚12~20 mm,剪力钉直径12 mm、长15 cm,剪力钉间距0.2~0.3 m,拉结筋直径12~20 mm,拉结筋长度与墙体厚度一致,拉结筋间距为剪力钉间距的2倍。通过改变上述参数,进行钢板混凝土墙抗撞击性能的敏感性分析。

      钢板混凝土结构墙体的钢板、拉结筋、剪力钉和混凝土分离建模。混凝土、钢板的单元类型分别为solid 164、shell 163,拉结筋和剪力钉的单元类型为beam 161,各个部件的有限元模型如图 8所示。混凝土沿墙体的厚度共划分为5层单元,墙体中心区域单元长度约0.2 m,墙体边缘单元长度约为0.4 m,从中间区域向外侧网格渐变;钢板的网格尺寸与剪力钉的间距保持一致,整个墙体模型的单元数为50 721,节点数为70 687。拉结筋、剪力钉及钢板与混凝土之间的连接通过ANSYS/LS-DYNA软件中自带的耦合约束*CONSTRAINED_LAGRANGE_IN_SOLID[8]。撞击荷载作用于墙体中心区域的钢板上,加载区域如图 8(b)所示。为有效传递荷载,钢板部件与混凝土部件之间定义为自动面-面接触。

      图  8  钢板混凝土结构墙的有限元模型

      Figure 8.  FEM model of the steel plate concrete wall

      在高速碰撞过程中,钢材和混凝土的强度会随着材料的高应变率有所提高,采用NEI07-13报告[2]中推荐的动力强化系数来考虑此特性。考虑动力强化系数后钢板混凝土结构墙的基本材料参数见表 1。表中,ρ为材料密度,E为弹性模量,υ为泊松比,fc为抗压强度, ft为抗拉强度。

      材料 ρ/(kg·m-3) E/GPa υ fc/MPa ft/MPa
      混凝土 2 400 35.2 0.2 46.0 4.2
      钢板 7 800 206.0 0.3 400.0 400.0
      拉结筋 7 800 206.0 0.3 330.0 330.0
      剪力钉 7 800 206.0 0.3 330.0 330.0

      表 1  钢板混凝土墙的材料参数

      Table 1.  Material parameters of steel plate concrete wall

    • 由于金属材料在碰撞时材料应变率效应十分明显,故选用考虑屈服、强化及应变率效应的Cowper-Symonds本构模型,该模型为碰撞分析中常用的金属本构。在ANSYS/LS-DYNA软件中Cowper-Symonds本构模型对应于*MAT_PLASTIC_KINEMATIC[8],其屈服函数为:

      $ \sigma_{\rm{y}}=\left[1+\left(\frac{\dot{\varepsilon}}{C}\right)^{\frac{1}{p}}\right]\left(\sigma_{0}+\beta E_{\rm{p}} \varepsilon_{\rm{p}, \rm{eff}}\right) $

      式中:σy为考虑应变率效应的屈服强度;σ0为初始屈服应力; $\dot{\varepsilon}$为应变率;Cp为应变率常数;εp, eff为等效塑性应变;β为硬化参数,β=0代表随动强化,β=1代表各向同性强化;Ep是塑性硬化模量。该模型通过等效塑性应变来定义单元的失效准则,删除失效单元。

    • 进行飞机撞击的瞬时响应分析,需要研究混凝土的材料非线性本构关系、应变率效应及失效准则。在前人研究[10-12]和运用的基础上,混凝土材料选用了模型参数相对较少且应用较为广泛的Winfrith本构模型,并考虑了应变率效应($\dot{\varepsilon}$=0)[8],其中断裂能的取值参考文献[13]。此模型[14]的屈服函数为:

      $ Y\left(I_{1}, J_{2}, J_{3}\right)=a J_{2}+\lambda \sqrt{J_{2}}+b I_{1}-1 $

      式中:I1为应力张量第一不变量;J2J3分别为应力偏张量第二不变量、应力偏张量第三不变量:

      $ \lambda=\left\{\begin{array}{ll}{k_{1} \cos \left[\frac{1}{3} \cos ^{-1}\left(k_{2} \cos (3 \theta)\right)\right]} & {\cos (3 \theta) \geqslant 0} \\ {k_{1} \cos \left[\frac{\pi}{3}-\frac{1}{3} \cos ^{-1}\left(-k_{2} \cos (3 \theta)\right)\right]} & {\cos (3 \theta)<0}\end{array}\right. $

      $ \cos (3 \theta)=\frac{3 \sqrt{3}}{2} \frac{J_{3}}{J_{2}^{3 / 2}} $

      式中:参数abk1k2均是混凝土抗拉强度与抗压强度比值(ft/fc)的函数。

      通过关键字*MAT_ADD_EROSION[8]控制材料失效,由于冲击荷载作用下,混凝土材料强度会随应变率变化,且在冲击作用下混凝土主要是受压破坏,因此选取主应变作为混凝土材料失效准则。

    • 钢板混凝土结构材料本构模型、失效准则以及荷载时程分析方法的合理性验证工作主要基于日本Kobori综合研究所对1/7.5缩尺飞机模型垂直撞击不同厚度钢板混凝土(SC)墙的实验研究[4]。文献[6]中将荷载时程分析法的计算结果与飞射物-靶体相互作用方法计算结果[5, 15-16]以及实验结果[4]进行了对比分析,结果表明,荷载时程分析法能够较好地模拟此冲击实验,且计算结果略偏大,验证了钢板混凝土结构墙材料本构模型选取的有效性,以及由于靶体完全刚性假设造成的此方法具有一定的保守性。

    • 基于荷载时程分析法进行一系列影响钢板混凝土墙体抗冲击性能的参数敏感性分析,其中参数包括混凝土厚度(0.75、0.80、0.85、0.90、1.0 m),钢板厚度(12、14、16、18、20 mm),拉结筋直径(12、16、18、20 mm),拉结筋间距(0.4、0.5、0.6 m)。

    • 为了对比如图 5~6所示的2种不同加载模式的保守性,选取钢板厚度为16 mm、拉结筋直径为12 mm、拉结筋间距为0.5 m、墙体厚度为0.90 m的钢板混凝土结构墙为研究对象,进行2种加载模式下的墙体冲击响应对比分析,其墙体背部的最大位移时程曲线与撞击区域混凝土的最大塑性应变时程曲线如图 9~10所示。其中,图例“Load case 1”、“Load case 2”分别表示如图 5~6所示的加载模式。此外,2种不同的加载模式下混凝土的最大塑性应变云图见图 11

      图  9  不同加载模式的墙体最大位移时程曲线

      Figure 9.  Curves of max displacement vs. time-history of SC wall for different loads

      图  10  不同加载模式的混凝土最大塑性应变时程曲线

      Figure 10.  Curves of max plastic strain vs. time-history of concrete for different loads

      图  11  不同加载模式下混凝土的最大塑性应变云图

      Figure 11.  Contour plot of max plastic strain of concrete for different loads

      图 11可见,在不同的加载模式下,其撞击区域混凝土的塑性分布是不同的。第1种加载方式即机身和机翼区域分别加载,由于作用于机翼区域的荷载值远大于机身区域,导致机翼区域的混凝土进入塑性区;而第2种加载方式,即整个飞机撞击区域一致均布加载,则混凝土的塑性区集中于机身及机翼的根部区域。由图 9~10中的曲线对比可见,第2种加载模式下墙体的位移响应值和撞击区域混凝土的塑性应变值均比第1种加载模式下的响应值偏大。可见,对于墙体的整体破坏效应评估,其作用于整个撞击区域的较长持时的一致均布加载比单独作用于机身和机翼区域的加载起控制作用。综上所述,鉴于第2种加载模式的保守性,在下文的墙体参数敏感性分析中均采用此整个撞击区域一致均布加载的模式。

    • 对钢板厚度16 mm,拉结筋直径12 mm,拉结筋间距0.5 m,墙体厚度分别选0.75、0.80、0.85、0.90、1.00 m的钢板混凝土结构墙体进行抗冲击性能的对比分析。在如图 6所示的冲击荷载作用下,以上不同厚度墙体背部的撞击区域处节点沿冲击方向的最大位移时程曲线和撞击区域混凝土的最大塑性应变时程曲线分别如图 12~13所示。

      图  12  不同厚度SC墙体的最大位移时程曲线

      Figure 12.  Curves of max displacement vs. time-history for SC wall with different thickness

      图  13  不同厚度SC墙体的混凝土最大塑性应变时程曲线

      Figure 13.  Curves of max plastic strain vs. time-history for concrete of SC wall with different thickness

      图 12~13可见,厚度为0.75、0.80、0.85、0.90、1.0 m的墙体,背部撞击区域的最大位移值分别为46.5、38.2、23.6、6.67、3.13 cm;撞击区域混凝土的最大塑性应变值分别为1.0×10-2、3.68×10-3、2.84×10-3、1.96×10-3、1.24×10-3。可见,墙体厚度对位移响应起到显著的控制作用。比如,墙体厚度从0.8 m增大到0.85 m时,最大位移响应由38.2 cm减小为23.6 cm,相当于减少了38.2%;当墙体厚度从0.8 m增大到0.9 m时,最大位移响应由38.2 cm减为6.67 cm,相当于减少了82.5%;当墙体厚度从0.9 m增大到1.0 m时,最大位移响应由6.67 cm减为3.13 cm。从以上位移对比可得,墙体厚度对于减少位移响应起到显著的控制作用,但当墙体厚度达到一定值(比如0.9 m),其结构的响应已经较小时,再进一步增大截面的厚度,对于减小位移响应影响不是很显著。同样,当墙体厚度为0.75 m时,撞击区域混凝土的塑性应变达到1.0%,当墙体厚度为0.8 m时,混凝土的塑性应变骤减至3.68×10-3,随着墙体厚度的增大,撞击区域混凝土的塑性应变随之减少。所以,在其他参数不变的情况下,增大墙体的厚度能够非常有效地减小冲击作用下结构的响应。

    • 对墙体厚度为0.85 m,拉结筋直径为12 mm、间距0.5 m,钢板厚度分别取14、16、18和20 mm的钢板混凝土结构墙,进行墙体的抗冲击性能对比分析。

      墙体背部撞击部位区域沿冲击方向的最大位移时程曲线与混凝土的最大塑性应变时程曲线分别见图 14~15。从图中可见,随着钢板厚度的增大,墙体背部的最大位移响应以及混凝土的最大塑性应变是减小的;对应钢板厚度为14、16、18和20 mm,其最大位移值分别为24.4、23.6、21.7和19.8 cm,但墙体背部钢板的最终残余位移值相当,可见增大钢板的厚度,可以减小冲击作用下结构的位移响应;撞击区域混凝土的最大塑性应变值分别为2.73×10-3、2.71×10-3、2.65×10-3和2.63×10-3,可见随着钢板厚度的增加,撞击区域混凝土的最大塑性应变值有所减小,但效果并不显著。

      图  14  不同钢板厚度的墙体最大位移时程曲线

      Figure 14.  Curves of max displacement vs. time-history for SC wall of steel plates with different thickness

      图  15  不同钢板厚度的混凝土最大塑性应变时程曲线

      Figure 15.  Curves of max plastic strain vs. time-history for concrete for steel plates with different thickness

    • 对墙体厚度为0.80 m,拉结筋直径为12 mm、间距为0.5 m,钢板厚度分别为12、14、18和22 mm的钢板混凝土结构墙进行冲击响应的对比分析。

      撞击区域墙体背部沿冲击方向的最大位移时程曲线以及混凝土的最大塑性应变时程曲线分别见图 16~17。从图中可见,随着钢板厚度的增大,墙体背部的最大位移响应值、最终残余位移值以及混凝土的最大塑性应变均是减小的。钢板厚度为12、14、18和22 mm的墙体,其最大位移值分别为44.8、39.3、31.3和27.0 cm,背部钢板的残余位移值分别为34.3、27.6、15.8、10.0 cm;撞击区域混凝土的最大塑性应变值分别为3.05×10-3、2.84×10-3、2.60×10-3和2.51×10-3。可见,对于相对较薄的钢板混凝土结构墙,增大钢板的厚度,可以有效地减小冲击作用下墙体的位移响应值,以14 mm厚钢板为基准,18、22 mm厚度的钢板墙的最大位移分别减少20.4%与31.3%,残余位移分别减小42.7%与63.8%。同样,随着钢板厚度的增加,撞击区域混凝土的最大塑性应变值也有效地减小。在冲击荷载作用下,钢板厚度分别为14、18 mm的墙体,其最大位移时刻对应的背部钢板的Mises等效应力云图分别见图 18~19。背部钢板中心区域的红色云图表示其应力已达到屈服值,通过对比可见,14 mm厚钢板的中间撞击部位区域均已屈服,而18 mm厚钢板的中间撞击部位区域,只有小部分区域达到屈服,大部分区域的应力小于或接近屈服应力。综上可得,对于厚度为0.80 m的钢板混凝土结构墙体,增大钢板的厚度能够有效地减小结构的最大位移、残余位移、混凝土的塑性应变以及钢板的等效应力等。

      图  16  不同钢板厚度的墙体最大位移时程曲线

      Figure 16.  Curves of max displacement vs. time-history for SC wall of steel plates with different thickness

      图  17  不同钢板厚度的混凝土最大塑性应变时程曲线

      Figure 17.  Curves of max plastic strain vs. time-history for concrete of steel plates with different thickness

      图  18  SC墙体14 mm厚背部钢板的等效应力云图

      Figure 18.  Contour plot of the equivalent stress of back face of 14mm steel plate for SC wall

      图  19  SC墙体18 mm厚背部钢板的等效应力云图

      Figure 19.  Contour plot of the equivalent stress of back face of 18mm steel plate for SC wall

      通过以上厚度为0.85、0.80 m的钢板混凝土结构墙体在不同厚度钢板情况下的冲击响应对比分析可知,总体来讲,增大钢板的厚度能够约束混凝土在撞击方向上的运动,有效地减小冲击作用下结构的响应,对于防护飞机撞击起到良好的作用。特别对于相对较薄的墙体,增加钢板的厚度在减小墙体的残余变形、混凝土的塑性应变以及钢板的有效应力等方面,其效果更显著。

    • 以墙体厚度为0.85 m、钢板厚度为16 mm、拉结筋间距为0.5 m的墙体作为研究对象,其拉结筋直径分别取为12、16、18和20 mm时,进行墙体结构响应的对比分析。通过如图 20所示的撞击区域墙体背部沿冲击方向的最大位移时程曲线对比可得,拉结筋直径为12、16、18和20 mm的墙体,其最大位移值分别为23.6、22.7、21.4和20.4 cm。可见,随着拉结筋直径的增大,墙体的最大响应略有减小,其效果不是很显著。可见增大拉结筋直径可有效地抵抗平面外剪力,但在顺着钢筋长度方向施加冲击荷载,抵抗飞机撞击的效果不明显。

      图  20  不同拉结筋直径的墙体最大位移时程曲线

      Figure 20.  Curves of max displacement vs. time-history for SC wall with different diameter of tie-bar

    • 以墙体厚度为0.85 m、钢板厚度为16 mm、拉结筋直径12 mm的墙体作为研究对象,其拉结筋间距分别取为0.4、0.5、0.6 m,进行墙体结构响应的对比分析。当拉结筋的间距取为0.6 m时,墙体混凝土部分的最大位移云图见图 21、其墙体背部的残余位移云图见图 22、墙体前部和背部钢板的最大塑性应变云图见图 23。从图 21可见,墙体中心撞击区域的混凝土发生较大位移且部分被压碎,最大位移达到2.137 m,且此区域的剪力钉和拉结筋发生部分脱落。从图 22可见,墙体背部中心区域的钢板被撕裂,最后混凝土与钢板的残余变形达到1.824 m;由于具有良好延展性的钢板的存在,约束了混凝土在撞击方向上的运动,虽然位移响应较大,仍限制了混凝土碎片的飞溅。从图 23可知,前部和背部钢板均发生了局部破坏,最大塑性应变达到9.95%。当拉结筋间距取为0.5、0.4 m时,墙体背部节点的最大位移时程对比见图 24,其最大位移分别为23.6、17.9 cm。墙体的前部和背部钢板仅局部发生塑性变形,没有单元发生失效。

      图  21  混凝土的最大位移云图

      Figure 21.  Contour plot of max displacement of concrete

      图  22  SC墙体背部的残余位移云图

      Figure 22.  Contour plot of residual displacement of SC wall back surface

      图  23  SC墙体钢板的最大塑性应变云图

      Figure 23.  Contour plot of max plastic strain of SC wall steel plate

      图  24  不同拉结筋间距的墙体最大位移时程曲线

      Figure 24.  Curves of max displacement vs. time-history for differently spaced tie-bar

      通过对比以上不同拉结筋间距的墙体响应可知,拉结筋间距对墙体的抗冲击性能影响显著。若拉结筋间距过大,墙体的抗冲击性能较差,会发生整体破坏。拉结筋的加密有效地增强了钢板与混凝土之间的连接,提高了墙体构件的整体刚度,能够有效地提高墙体的抗冲击性能,以减小冲击作用下的结构响应,对于防护飞机撞击起到良好效果。

    • 建立了三维钢板混凝土结构墙体的精细化有限元模型,基于荷载时程分析方法对影响钢板混凝土结构墙抗冲击性能的一系列参数进行了敏感性分析。这些参数包括混凝土厚度、钢板厚度和拉结筋直径与间距。得出结论如下:

      (1) 总体来讲,以上参数均会影响钢板混凝土结构墙体的抗冲击性能。增大混凝土与钢板厚度、增大拉结筋直径以及减小拉结筋与剪力钉间距,均有利于墙体抗冲击性能的提高。

      (2) 在其他参数不变的情况下,单方面的增大墙体混凝土厚度,即增大墙体的截面尺寸,对于提高钢板混凝土结构墙体的抗冲击性能效果显著,起到决定性作用。

      (3) 通过不同钢板厚度的结构响应敏感性分析,可以得出,增大钢板厚度能够有效地减小冲击作用下结构的响应,对于防护飞机撞击起到良好的作用。特别对于相对较薄墙体,增加钢板厚度能够有效的减小墙体的残余变形、混凝土的塑性应变以及钢板的有效应力等,其效果更显著。

      (4) 适当减小拉结筋间距,相对于增大拉结筋的直径,在减小结构的冲击响应、提高墙体的抗冲击性能方面,其效果更显著。

参考文献 (16)

目录

    /

    返回文章
    返回