留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

P波斜入射下非基岩场地中核岛结构地震响应规律研究

赵密 王鑫 钟紫蓝 杜修力

赵密, 王鑫, 钟紫蓝, 杜修力. P波斜入射下非基岩场地中核岛结构地震响应规律研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.12.0744
引用本文: 赵密, 王鑫, 钟紫蓝, 杜修力. P波斜入射下非基岩场地中核岛结构地震响应规律研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.12.0744
Mi ZHAO, Xin WANG, Zi-lan ZHONG, Xiu-li DU. STUDY ON SEISMIC RESPONSES OF NUCLEAR ISLAND STRUCTURE IN NON-BEDROCK SITE UNDER OBLIQUELY-INCIDENCE OF P WAVES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.12.0744
Citation: Mi ZHAO, Xin WANG, Zi-lan ZHONG, Xiu-li DU. STUDY ON SEISMIC RESPONSES OF NUCLEAR ISLAND STRUCTURE IN NON-BEDROCK SITE UNDER OBLIQUELY-INCIDENCE OF P WAVES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.12.0744

P波斜入射下非基岩场地中核岛结构地震响应规律研究

doi: 10.6052/j.issn.1000-4750.2019.12.0744
基金项目: 国家自然科学基金项目(51738001,51421005);教育部创新团队发展计划项目(IRT_17R03)
详细信息
    作者简介:

    赵 密(1980−),男,吉林人,教授,博士,主要从事重大工程抗震研究(E-mail: zhaomi@bjut.edu.cn)

    王 鑫(1995−),男,河南人,硕士生,主要从事重大工程抗震研究(E-mail: wx5048@emails.bjut.edu.cn)

    杜修力(1963−),男,四川人,教授,博士,主要从事地震工程研究工作(E-mail: duxiuli@bjut.edu.cn)

    通讯作者: 钟紫蓝(1986−),男,江西人,教授,博士,主要从事土结构动力相互作用研究(E-mail: zilanzhong@bjut.edu.cn)
  • 中图分类号: P315.9

STUDY ON SEISMIC RESPONSES OF NUCLEAR ISLAND STRUCTURE IN NON-BEDROCK SITE UNDER OBLIQUELY-INCIDENCE OF P WAVES

  • 摘要: 随着我国核电产业飞速发展,核电站厂址的选择将不可避免地遇到非基岩场地。针对地震动斜入射下非基岩场地中核电结构抗震问题,建立了包括CAP1400核岛结构和岩土场地的土-结构系统三维有限元模型,基于黏弹性人工边界和斜入射地震动输入方法,研究了P波斜入射角度和非基岩分层场地对核岛结构地震响应的影响规律。研究表明:非基岩场地将增大核岛结构加速度响应和楼层反应谱,多数条件下也将增大楼层相对位移;P波入射角度增大,核岛结构竖向响应减小、水平响应增大,楼层反应谱的频谱特性发生改变。
  • 图  1  核岛结构有限元模型

    Figure  1.  Finite element model of nuclear island structure

    图  2  场地有限元模型 /m

    Figure  2.  Finite element model of site

    图  3  地震动时程及反应谱

    Figure  3.  Time history and response spectrum of seismic motion

    图  4  土-结构相互作用系统x-y平面示意图 /m

    Figure  4.  x-y plane view of soil-structure interaction system

    图  5  黏弹性边界示意图

    Figure  5.  Diagram of viscous-elastic artificial boundary

    图  6  参考点位置示意图

    Figure  6.  Locations of selected reference points

    图  7  参考点N1的加速度时程

    Figure  7.  Acceleration time histories of reference point N1

    图  8  参考点的峰值加速度

    Figure  8.  Peak accelerations of reference points

    图  9  楼层加速度反应谱

    Figure  9.  Floor acceleration response spectra

    图  10  楼层相对位移

    Figure  10.  Relative floor displacements

    图  11  关键楼层相对位移

    Figure  11.  Relative displacements of key floors

    表  1  场地材料参数

    Table  1.   Material properties of site

    类型厚度/m密度/(kg/m3)剪切波速/(m/s)压缩波速/(m/s)泊松比
    分层场地20195032020900.488
    20160055414920.420
    302800160038260.394
    硬岩场地702800160038260.394
    下载: 导出CSV

    表  2  场地和结构的阻尼系数

    Table  2.   Damping coefficients of the site and structure

    类型第一阶频率/Hz质量系数α刚度系数β
    分层场地2.691.110430.00203
    硬岩场地5.711.701970.00146
    核岛结构4.092.004470.00241
    下载: 导出CSV
  • [1] 王涛, 王飞, 侯钢领, 等. 核电厂隔震结构的振动台试验研究[J]. 工程力学, 2014, 31(10): 62 − 68, 84.

    Wang Tao, Wang Fei, Hou Gangling, et al. Shaking table tests on base-isolated nuclear power plant [J]. Engineering Mechanics, 2014, 31(10): 62 − 68, 84. (in Chinese)
    [2] 黄欢, 丁文杰, 郭海兵. 影响我国内陆核电发展的关键性问题分析[J]. 南华大学学报(社会科学版), 2019, 20(3): 9 − 15.

    Huang Huan, Ding Wenjie, Guo Haibing. Analysis of the key factor affecting the development of Chinese inland nuclear power [J]. Journal of University of South China (Social Science Edition), 2019, 20(3): 9 − 15. (in Chinese)
    [3] 胡进军, 李琼林, 邬迪, 等. 基于PSHA的核电厂近断层抗震设计谱构建方法[J]. 工程力学, 2019, 36(5): 192 − 199.

    Hu Jinjun, Li Qionglin, Wu Di, et al. Construction method of the near fault seismic design spectrum for nuclear power plants based on probabilistic seismic hazard analysis [J]. Engineering Mechanics, 2019, 36(5): 192 − 199. (in Chinese)
    [4] Leonardo T S, Richard S O, Sener T, et al. Finite element modeling of the AP1000 nuclear island for seismic analyses at generic soil and rock sites [J]. Nuclear Engineering and Design, 2007, 237(12): 1474 − 1485.
    [5] Elkhoraibi T, Hashemi A, Ostadan F. Probabilistic and deterministic soil structure interaction analysis including ground motion incoherency effects [J]. Nuclear Engineering and Design, 2014, 269: 250 − 255. doi:  10.1016/j.nucengdes.2013.08.036
    [6] Politopoulos I, Sergis I, Wang F. Floor response spectra of a partially embedded seismically isolated nuclear plant [J]. Soil Dynamics and Earthquake Engineering, 2015, 78: 213 − 217. doi:  10.1016/j.soildyn.2015.06.017
    [7] Farahani R V, Dessalegn T M, Vaidya N R, et al. Seismic soil–structure interaction analysis of a nuclear power plant building founded on soil and in degraded concrete stiffness condition [J]. Nuclear Engineering and Design, 2016, 297: 320 − 326. doi:  10.1016/j.nucengdes.2015.12.002
    [8] 李小军, 王晓辉, 王贵英. 非基岩场地核电厂CAP1400结构地震反应分析[J]. 应用基础与工程科学学报, 2017, 25(3): 584 − 594.

    Li Xiaojun, Wang Xiaohui, Wang Guiying. Seismic response analysis of nuclear power plant CAP1400 structure in non-rock site [J]. Journal of Basic Science and Engineering, 2017, 25(3): 584 − 594. (in Chinese)
    [9] 王贵英. 非基岩场地核电厂CAP1400结构—地基相互作用效应分析 [D]. 北京: 北京工业大学, 2016.

    Wang Guiying. Interaction effect analysis of nuclear power plant CAP1400 structure-foundation in non-rock site [D]. Beijing: Beijing University of technology, 2016. (in Chinese)
    [10] 李小军, 王晓辉, 贺秋梅, 等. 非基岩核电厂结构地震响应振动台试验研究[J]. 核动力工程, 2017, 38(04): 31 − 35.

    Li Xiaojun, Wang Xiaohui, He Qiumei, et al. Shaking table test for evaluation of seismic behavior of nuclear power plants on non-rock site [J]. Nuclear Power Engineering, 2017, 38(04): 31 − 35. (in Chinese)
    [11] 王晓辉. 非基岩场地核电厂地震反应试验与数值模拟分析 [D]. 北京: 北京工业大学, 2017.

    Wang Xiaohui. Shaking table tests and numerical simulation analysis of nuclear power plants in non-rock site [D]. Beijing University of Technology, 2017. (in Chinese)
    [12] 高永武, 王涛, 戴君武, 等. 不同场地条件下某新型核电厂房的地震响应试验研究[J]. 振动与冲击, 2017, 36(18): 214 − 222.

    Gao Yongwu, Wang Tao, Dai Junwu, et al. Experimental research on seismic responses of a new type of nuclear power plant under different site conditions [J]. Journal of Vibration and Shock, 2017, 36(18): 214 − 222. (in Chinese)
    [13] 陈少林, 常梦利. SH波斜入射时三维结构-土-结构相互作用分析[J]. 地震工程与工程振动, 2017, 37(1): 81 − 90.

    Chen Shaolin, Chang Mengli. Analysis of three-dimensional structure-soil-structure interaction subjected to inclined SH plane wave [J]. Earthquake Engineering and Engineering Dynamics, 2017, 37(1): 81 − 90. (in Chinese)
    [14] 杜修力, 陈维, 李亮, 等. 斜入射条件下地下结构时域地震反应分析初探[J]. 震灾防御技术, 2007(3): 290 − 296. doi:  10.3969/j.issn.1673-5722.2007.03.009

    Du Xiuli, Chen Wei, Li Liang, et al. Preliminary study on time-domain seismic response analysis of underground structures under oblique incident [J]. Technology for Earthquake Disaster Prevention, 2007(3): 290 − 296. (in Chinese) doi:  10.3969/j.issn.1673-5722.2007.03.009
    [15] 徐海滨, 杜修力, 赵密, 等. 地震波斜入射对高拱坝地震反应的影响[J]. 水力发电学报, 2011, 30(6): 159 − 165.

    Xu Haibin, Du Xiuli, Zhao Mi, et al. Effect of oblique incidence of seismic waves on seismic responses of high arch dam [J]. Journal of Hydroelectric Engineering, 2011, 30(6): 159 − 165. (in Chinese)
    [16] 孙建超. 地震波斜入射条件下核安全壳的地震响应分析 [D]. 哈尔滨: 哈尔滨工业大学, 2015.

    Sun Jianchao. Seismic response analysis of nuclear containment under the condition of oblique incidence of seismic wave [D]. Harbin: Harbin Institute of Technology, 2015. (in Chinese)
    [17] 宗娟. 地震波空间斜入射时核电结构土—结动力相互作用分析 [D]. 南京: 南京航空航天大学, 2018.

    Zong Juan. Soil-structure interaction analysis of nuclear reactor structure subjected to obliquely incident seismic wave [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2018. (in Chinese)
    [18] ABAQUS. Standard user’s manual (Version 2017) [M]. Providence: Dassault Systemes Simulia Corporation, 2016.
    [19] GB 50267−97, 核电厂抗震设计规范 [S]. 北京: 中国计划出版社, 1998.

    GB 50267−97, Code for seismic design of nuclear power plant [S]. Beijing: China Planning Press, 1998. (in Chinese)
    [20] 李忠献, 李忠诚, 梁万顺. 考虑土-结构相互作用和岩土参数不确定性的核电厂结构地震响应分析[J]. 地震工程与工程振动, 2006(2): 143 − 148. doi:  10.3969/j.issn.1000-1301.2006.02.023

    Li Zhongxian, Li Zhongcheng, Liang Wanshun. Seismic response analysis of nuclear power plant structures considering soil-structure interaction and uncertainty of soil properties [J]. Earthquake Engineering and Engineering Vibration, 2006(2): 143 − 148. (in Chinese) doi:  10.3969/j.issn.1000-1301.2006.02.023
    [21] Lysmer J, Kuhlemeyer R L. Finite dynamic model for infinite media [J]. Journal of the Engineering Mechanics Division, 1969, 95(4): 859 − 878.
    [22] Brinkgreve R B J, Vermeer P A. Finite element code for soil and rock analysis [M]. Rotterdam: A. A. Balkema Publishers, 1998.
    [23] 廖振鹏. 工程波动理论导论 [M]. 第二版. 北京: 科学出版社, 2002.

    Liao Zhenpeng. Introduction to wave motion theory in engineering [M]. 2nd ed. Beijing: Science Press, 2002. (in Chinese)
    [24] Du X L, Zhao M. A local time-domain transmitting boundary for simulating cylindrical elastic wave propagation in infinite media [J]. Soil Dynamics and Earthquake Engineering, 2010, 30(10): 937 − 946. doi:  10.1016/j.soildyn.2010.04.004
    [25] 刘晶波, 王振宇, 杜修力, 等. 波动问题中的三维时域粘弹性人工边界[J]. 工程力学, 2005, 22(6): 46 − 51. doi:  10.3969/j.issn.1000-4750.2005.06.008

    Liu Jingbo, Wang Zhenyu, Du Xiuli, at al. Three-dimensional visco-elastic artificial boundaries in time domain for wave motion problems [J]. Engineering Mechanics, 2005, 22(6): 46 − 51. (in Chinese) doi:  10.3969/j.issn.1000-4750.2005.06.008
    [26] 杜修力, 赵密, 王进廷. 近场波动模拟的人工应力边界条件[J]. 力学学报, 2006, 38(1): 49 − 56. doi:  10.3321/j.issn:0459-1879.2006.01.007

    Du Xiuli, Zhao Mi, Wang Jinting. A stress artificial boundary in for near-field wave problem [J]. Chinese Journal of Theoretical and Applied Mechanics, 2006, 38(1): 49 − 56. (in Chinese) doi:  10.3321/j.issn:0459-1879.2006.01.007
    [27] 赵密, 孙文达, 高志懂, 等. 阶梯地形成层场地的斜入射地震动输入方法[J]. 工程力学, 2019, 36(12): 62 − 68.

    Zhao Mi, Sun Wenda, Gao Zhidong, at al. Input method of obliquely incident earthquake for step-shaped layered site [J]. Engineering Mechanics, 2019, 36(12): 62 − 68. (in Chinese)
    [28] Zhao M, Gao Z D, Wang L T, et al. Obliquely incident earthquake input for soil-structure interaction in layered half space [J]. Earthquakes and Structures, 2017, 13(6): 573 − 588.
    [29] 尹侯权. 地震波斜入射时成层半空间场地反应分析方法及其应用 [D]. 北京: 北京工业大学, 2015.

    Yin Houquan. A method for layered half- space site response under obliquely incident earthquake and its application [D]. Beijing: Beijing University of Technology, 2015. (in Chinese)
  • [1] 李志远, 钟红, 胡志强, 林皋.  局部褶皱对层状地基马蹄形孔洞散射的影响分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2019.09.0512
    [2] 韩润波, 许成顺, 许紫刚, 蒋家卫, 杜修力.  对称地下结构抗震分析的边界强制反应位移法 . 工程力学, doi: 10.6052/j.issn.1000-4750.2020.02.0075
    [3] 巴振宁, 张家玮, 梁建文, 吴孟桃.  地震波斜入射下层状TI饱和场地地震反应分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2019.07.0363
    [4] 王丕光, 赵密, 李会芳, 杜修力.  一种高精度圆柱形人工边界条件:水-柱体相互作用问题 . 工程力学, doi: 10.6052/j.issn.1000-4750.2017.10.0787
    [5] 赵密, 高志懂, 杜修力, 王君杰.  地震作用下深厚土层-结构相互作用的高效分析方法 . 工程力学, doi: 10.6052/j.issn.1000-4750.2018.04.0245
    [6] 张奎, 夏佩林, 李伟华, 赵成刚.  斜入射下水下地基场地地震输入的一维化时域方法 . 工程力学, doi: 10.6052/j.issn.1000-4750.2018.09.0523
    [7] 张立红, 胡晓, 曾迪, 周德才, 毛宇, 吕玮.  基于抗震性能的高烈度区高端阀厅选型研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2017.05.S061
    [8] 王笃国, 赵成刚.  地震波斜入射下考虑场地非线性、地形效应和土结动力相互作用的大跨连续刚构桥地震响应分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.03.0180
    [9] 侯钢领, 李祯, 贾晓飞, 宋天舒.  储液结构隔震优化模型及其影响参数分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.04.S052
    [10] 张季, 梁建文, 巴振宁.  水平层状饱和场地地震响应分析的等效线性化方法 . 工程力学, doi: 10.6052/j.issn.1000-4750.2015.03.0215
    [11] 黄福云, 李建中, 陈宝春, 吴庆雄.  钢管混凝土单圆管拱结构振动台阵试验研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.11.0857
    [12] 刘昊苏, 董军, 杨昀.  不同设计参数下矮塔斜拉桥的抗震动力学分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.05.S017
    [13] 刘晶波, 王文晖, 赵冬冬.  地下结构横截面地震反应拟静力计算方法对比研究 . 工程力学, doi: 10.6052/j.issn.1000-4750.2011.05.0316
    [14] 张海顺, 姜忻良, 张亚楠.  高架桥-地铁站-桩-土复杂结构体系地震反应分析 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.06.S041
    [15] 闫维明, 李晰, 陈彦江, 李勇, 顾大鹏.  钢管混凝土拱桥台阵试验研究:场地条件的影响 . 工程力学, doi: 10.6052/j.issn.1000-4750.2012.03.0178
    [16] 赵密, 杜修力, 刘晶波.  一种高阶精度人工边界条件:出平面外域波动问题 . 工程力学,
    [17] 郜新军, 赵成刚, 刘秦.  地震波斜入射下考虑局部地形影响和土结动力相互作用的多跨桥动力响应分析 . 工程力学,
    [18] 赵 密, 杜修力.  时间卷积的局部高阶弹簧-阻尼-质量模型 . 工程力学,
    [19] 刘晶波, 刘祥庆, 薛颖亮.  地下结构抗震分析与设计的Pushover方法适用性研究 . 工程力学,
    [20] 苏 亮, 董石麟.  竖向多点输入下两种典型空间结构的抗震分析 . 工程力学,
  • 加载中
计量
  • 文章访问数:  49
  • HTML全文浏览量:  22
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-12-15
  • 修回日期:  2020-04-13
  • 网络出版日期:  2020-06-24

P波斜入射下非基岩场地中核岛结构地震响应规律研究

doi: 10.6052/j.issn.1000-4750.2019.12.0744
    基金项目:  国家自然科学基金项目(51738001,51421005);教育部创新团队发展计划项目(IRT_17R03)
    作者简介:

    赵 密(1980−),男,吉林人,教授,博士,主要从事重大工程抗震研究(E-mail: zhaomi@bjut.edu.cn)

    王 鑫(1995−),男,河南人,硕士生,主要从事重大工程抗震研究(E-mail: wx5048@emails.bjut.edu.cn)

    杜修力(1963−),男,四川人,教授,博士,主要从事地震工程研究工作(E-mail: duxiuli@bjut.edu.cn)

    通讯作者: 钟紫蓝(1986−),男,江西人,教授,博士,主要从事土结构动力相互作用研究(E-mail: zilanzhong@bjut.edu.cn)
  • 中图分类号: P315.9

摘要: 随着我国核电产业飞速发展,核电站厂址的选择将不可避免地遇到非基岩场地。针对地震动斜入射下非基岩场地中核电结构抗震问题,建立了包括CAP1400核岛结构和岩土场地的土-结构系统三维有限元模型,基于黏弹性人工边界和斜入射地震动输入方法,研究了P波斜入射角度和非基岩分层场地对核岛结构地震响应的影响规律。研究表明:非基岩场地将增大核岛结构加速度响应和楼层反应谱,多数条件下也将增大楼层相对位移;P波入射角度增大,核岛结构竖向响应减小、水平响应增大,楼层反应谱的频谱特性发生改变。

English Abstract

赵密, 王鑫, 钟紫蓝, 杜修力. P波斜入射下非基岩场地中核岛结构地震响应规律研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.12.0744
引用本文: 赵密, 王鑫, 钟紫蓝, 杜修力. P波斜入射下非基岩场地中核岛结构地震响应规律研究[J]. 工程力学. doi: 10.6052/j.issn.1000-4750.2019.12.0744
Mi ZHAO, Xin WANG, Zi-lan ZHONG, Xiu-li DU. STUDY ON SEISMIC RESPONSES OF NUCLEAR ISLAND STRUCTURE IN NON-BEDROCK SITE UNDER OBLIQUELY-INCIDENCE OF P WAVES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.12.0744
Citation: Mi ZHAO, Xin WANG, Zi-lan ZHONG, Xiu-li DU. STUDY ON SEISMIC RESPONSES OF NUCLEAR ISLAND STRUCTURE IN NON-BEDROCK SITE UNDER OBLIQUELY-INCIDENCE OF P WAVES[J]. Engineering Mechanics. doi: 10.6052/j.issn.1000-4750.2019.12.0744
  • 近年来频繁发生的核电事故,如柏崎·刈羽核电站核泄漏事故、福岛核电站核泄漏事故、北安娜核电站停机事件等,表明地震已经成为威胁核电安全的重要因素,地震及其次生灾害会对核电厂的正常运行和安全产生重大的影响[1]。在已有的核电工程建设中,厂址选择非常苛刻,通常为整块硬岩场地。随着我国核电发展,沿海硬岩场地逐渐稀缺,在内陆非基岩场地上建设核电工程不可避免[2-3]

    屏蔽厂房的地震反应分析往往需要考虑土-结构相互作用的影响[4-7]。李小军等[8-11]开展了CAP1400核岛结构的振动台试验和非基岩场地核岛结构体系地震响应数值模拟研究,发现非基岩场地对核岛结构地震响应影响显著,建议在进行核电结构地震反应分析时应考虑场地条件。高永武等[12]针对某新型屏蔽厂房及安全壳内部结构进行了振动台试验,建议《核电厂抗震设计规范》(GB 50267−1997)应进一步提高不考虑土-结构相互作用的场地条件。以上研究表明,非基岩场地条件对结构响应可能产生不利影响,开展非基岩场地下核岛结构地震响应研究十分必要。

    地震动的倾斜入射对结构响应会产生不利影响[13],特别是对于较大尺寸的工程结构,例如大坝[14]、大型地铁车站[15]、核电[16-17]等。杜修力等[14-15]研究发现地震动倾斜输入下大型结构响应与垂直入射下的响应具有明显差异,地震动倾斜入射可能显著增加结构关键点的频谱幅值放大系数。孙建超[16]基于黏弹性人工边界研究了CPR1000核安全壳在地震作用下的响应,认为按照地震波垂直入射下的计算分析结果对核安全壳结构进行设计是偏于不安全的。宗娟[17]针对某一核电模型,使用传递矩阵法计算场地反应,采用并行技术对空间斜入射时土-结构动力相互作用进行分析,讨论了地震波空间入射角度对土体、基础、结构响应的影响。综上,关于地震动斜入射下非基岩场地对结构响应的影响规律,仍然缺乏相关研究。

    本文针对我国具有自主知识产权的核岛结构“国和一号”(CAP1400),利用有限元软件ABAQUS建立土-核岛结构系统的有限元模型,进行斜入射下不同场地上结构地震响应数值模拟,探讨地震动斜入射下非基岩场地对核电结构地震响应的影响规律。

    • “国和一号”(CAP1400)是以浙江三门和山东海阳核电站为依托,在引进消化吸收AP1000三代核电技术基础上,通过再创新形成的具有自主知识产权、功率更大的大型先进非能动压水堆核电型号。与同属第三代结构的西屋公司AP1000核岛结构相比,CAP1400核岛结构整体尺寸有所加大。CAP1400核岛结构主要分为4部分,即屏蔽厂房、安全壳、内部结构和辅助厂房。结构总高度约87.750 m(含地下室,标高−12.200 m~75.550 m),平面最大尺寸为57.055 m×91.400 m。屏蔽厂房为圆形筒状斜屋顶结构,筒体底部外径为49.970 m、顶部标高为75.550 m、筒体壁厚1.100 m。屋顶为型钢混凝土结构,设有钢筋混凝土结构的环形冷却系统储水箱。辅助厂房为钢筋混凝土剪力墙结构,三面围绕屏蔽厂房布置,其屋顶标高呈阶梯状,为18.050 m、20.350 m和35.300 m,辅助厂房还有2层地下室,地下室底层标高为−12.200 m[8]

    • 参照CAP1400核岛结构的相关图纸,使用有限元软件ABAQUS[18]建立土-核岛结构系统的有限元模型。土-核岛结构系统包括场地和核岛结构,采用NEWMARK隐式积分进行计算。

      核岛结构有限元模型如图1所示,包括辅助厂房、屏蔽厂房、基础底板三部分,x轴正方向为东,z轴正方向为南,y轴为竖直方向。由于实际结构非常复杂,为了简化建模和计算过程,内部结构以增加楼板质量代替,屏蔽厂房上的进气孔、设备孔、屏蔽厂房内部的钢安全壳以及组合屋架处含有的水箱都未考虑。屏蔽厂房和辅助厂房采用壳单元,基础底板采用实体单元。核岛结构有限元模型共87057个单元,92689个节点。

      图  1  核岛结构有限元模型

      Figure 1.  Finite element model of nuclear island structure

      场地有限元模型如图2所示,均为实体单元。尺寸为500 m×500 m×70 m,共155992个单元,184214个节点。核岛屏蔽厂房埋深12.2 m,基础埋深18.6 m。

      图  2  场地有限元模型 /m

      Figure 2.  Finite element model of site

    • 选取基岩地表处的神户大学地震记录作为输入地震动,其加速度时程及对应的加速度反应谱如图3所示,地震动总持时20 s,卓越周期为1.22 s。加速度峰值按照极限安全地震动[19]的标准调幅到0.15 g,考虑地表放大效应折半后在场地底部进行输入。

      图  3  地震动时程及反应谱

      Figure 3.  Time history and response spectrum of seismic motion

      为研究入射角度的影响,入射角度分别取竖向剖面x-y平面内的0°、±10°、±20°、±30°,共7种工况,如图4所示。

    • 结构模型参数均采用CAP1400结构设计参数,弹性模量取34.5 GPa,密度取2350 kg/m3,泊松比取0.15。出于安全裕度考虑,核电厂结构往往非常厚重,刚度很大,且核电厂经常选址在刚度较高、地基条件较好的地基上,所以结构和地基均采用弹性本构[20]

      为了对比不同场地的影响,分别选取简化后的三层水平成层场地[17]和均质硬岩场地。场地材料参数见表1

      图  4  土-结构相互作用系统x-y平面示意图 /m

      Figure 4.  x-y plane view of soil-structure interaction system

      表 1  场地材料参数

      Table 1.  Material properties of site

      类型厚度/m密度/(kg/m3)剪切波速/(m/s)压缩波速/(m/s)泊松比
      分层场地20195032020900.488
      20160055414920.420
      302800160038260.394
      硬岩场地702800160038260.394

      计算分析中,选用Rayleigh阻尼模拟材料引起的能量消散,计算公式如下:

      $$C = \alpha { M} + \beta { K}$$ (1)

      式中:MK分别为质量矩阵和刚度矩阵;αβ分别为质量系数和刚度系数。

      αβ分别由下式计算确定:

      $$\alpha = \frac{{2{\omega _i}{\omega _j}({\omega _j}{\zeta _i} - {\omega _i}{\zeta _j})}}{{\omega _j^2 - \omega _i^2}}$$ (2)
      $$\beta = \frac{{2({\omega _i}{\zeta _i} - {\omega _j}{\zeta _j})}}{{\omega _i^2 - \omega _j^2}}$$ (3)

      式中,ζiζjωiωj分别为第ij阶模态的阻尼比和圆频率。

      本文计算中,场地的阻尼比均取为0.05,核岛结构的取0.07[19]。使用ABAQUS的模态分析功能可得到不同场地及核岛结构的固有频率,分析得到的第一阶固有频率和阻尼系数见表2

      表 2  场地和结构的阻尼系数

      Table 2.  Damping coefficients of the site and structure

      类型第一阶频率/Hz质量系数α刚度系数β
      分层场地2.691.110430.00203
      硬岩场地5.711.701970.00146
      核岛结构4.092.004470.00241
    • 在土-结构相互作用体系分析模型中,为了模拟无限地基对体系的辐射阻尼效应,需要对截取的土-结构体系模型施加人工边界条件,如黏性边界[21-22]、透射边界[23]、黏弹性边界[24-27]等。其中,黏弹性边界在人工边界处设置远端固定的并联弹簧-阻尼器系统,可以更好的模拟半无限介质的弹性恢复能力,稳定性良好,不存在低频飘移失稳、高频振荡失稳问题,其概念清晰,实现简单,广泛应用于各种结构计算中。

      图5为黏弹性边界示意图,其中弹簧-阻尼元件参数为[26]

      $$ {\text{法向}}:\;{K_{{\rm{BN}}}}{\rm{ = }}\frac{{\rm{1}}}{{{\rm{1 + }}A}} \cdot \frac{{\lambda + 2G}}{r},\;{C_{{\rm{BN}}}}{\rm{ = }}B\rho {c_{\rm p}} $$ (4)
      $$ {\text{切向}}:\;{K_{{\rm{BT}}}}{\rm{ = }}\frac{{\rm{1}}}{{{\rm{1 + }}A}} \cdot \frac{G}{r},\;{C_{{\rm{BT}}}}{\rm{ = }}B\rho {c_{\rm s}}\qquad\;\; $$ (5)

      式中:KBNKBT分别表示法向和切向弹簧系数;CBNCBT分别为法向和切向阻尼系数;ρ为介质密度;AB为无量纲参数,建议取值为0.8、1.1;r取为近场结构几何中心到该人工边界面的距离;cpcs分别是地震动的P波波速和S波波速。

      图  5  黏弹性边界示意图

      Figure 5.  Diagram of viscous-elastic artificial boundary

      通过将地震场地反应和黏弹性边界相结合,实现将场地地震反应转化为截断边界面上等效荷载的地震动输入[27-29]

    • 为研究不同场地和不同地震动入射角度对结构响应的影响,本文选取核岛结构屏蔽厂房、辅助厂房上共8个参考点进行分析。参考点位置如图6所示,N1~N3位于屏蔽厂房南侧,N4~N7位于屏蔽厂房南侧的辅助厂房上,N8位于基础底板上。

      图  6  参考点位置示意图

      Figure 6.  Locations of selected reference points

    • 通过计算,分别得到不同地震动入射角度和不同场地下屏蔽厂房上参考点N1~N3的加速度时程曲线。图7为参考点N1在分层场地下不同入射角度时加速度时程曲线,图8为分层场地和硬岩场地下不同地震动入射角度时参考点N1、N2和N3的峰值加速度结果。

      图  7  参考点N1的加速度时程

      Figure 7.  Acceleration time histories of reference point N1

      图  8  参考点的峰值加速度

      Figure 8.  Peak accelerations of reference points

      图8可知,分层场地条件下各参考点加速度峰值均大于均质硬岩场地,即分层场地的场地放大效应显著,符合场地放大效应规律,且3个方向分层场地条件和均质硬岩场地条件下加速度峰值随角度的变化规律基本一致。

      y向的峰值加速度最大,x向次之,z向最小。x方向上,各参考点的峰值加速度随着地震动入射的角度增大而增大;y向则相反,即各参考点峰值加速度随着角度增大而减小;z方向峰值加速度随入射角度的变化则没有明显规律。根据P波斜入射时的传播特性,随着入射角度的增大,振动方向在y方向的分量逐步减小,在x方向的分量逐步增大,所以上述参考点峰值加速度的规律符合P波的传播特性。

    • 取辅助厂房上的参考点N4和N7计算5%阻尼比下的楼层加速度反应谱,如图9所示。

      图  9  楼层加速度反应谱

      Figure 9.  Floor acceleration response spectra

      图9可知,y向楼层反应谱加速度峰值最大,x向次之,z向最小,与峰值加速度响应规律一致。除z向外,分层场地条件下楼层反应谱峰值均大于均质硬岩场地。分层场地除了会增大反应谱峰值外,也改变了反应谱的形状,体现在峰值出现的频率不同。z向上,分层场地的反应谱相较硬岩场地还会在中低频区域上多出一个较明显的峰值,出现在周期为0.32 s时。

      x向上,反应谱随着地震动入射角度增大而增大,不同入射角度之间相差较大,且随着角度的变化反应谱形状也发生了变化;y向上,整体上反应谱随着地震动入射角度增大而减小,但不同入射角度之间相差较小,不同入射角度下的反应谱形状也基本一致;z向上,不同地震动入射角度下楼层反应谱无明显规律,且小于xy向反应。

      综上,大多数情况下分层场地会增大反应谱,反应谱的峰值随着入射角度的变化规律仍然符合P波的传播特性,入射角度和场地类型的变化会改变反应谱的频谱特性。

    • 取辅助厂房上的参考点N4~N7,基础底板上的参考点N8,计算参考点N4~N7相对N8的峰值位移差,结果如图10所示,相对高度取各参考点相对N8的高度。

      图  10  楼层相对位移

      Figure 10.  Relative floor displacements

      图10可知,x向上分层场地上楼层相对位移均大于硬岩场地上的结果,分层场地的放大效应明显。y向上,地震动入射角为正时,分层场地上楼层相对位移大于硬岩场地的结果;地震动入射角为负时,分层场地上楼层相对位移小于硬岩场地的结果。z向上分层场地对楼层相对位移影响没有明显规律。硬岩场地上考虑地震动斜入射带来的影响小于分层场地,表现为不同入射角度之间差距相对分层场地的情况小。

      x向上,随着入射角度的增大,楼层的相对位移也随之增大,且正负入射角度的结果基本一致;y向上,正角度入射时楼层的相对位移随入射角度增大而增大,负角度入射时楼层的相对位移随入射角度增大而减小;z向上,随着入射角度的增大,楼层的相对位移也随之增大。

      为了便于分析,取出参考点N4和参考点N6所在楼层的相对位移结果,其随入射角度的变化如图11所示。

      图  11  关键楼层相对位移

      Figure 11.  Relative displacements of key floors

      图11所示,参考点N4和参考点N6所在楼层的相对位移随入射角度变化的规律基本一致,且参考点N4所在楼层的相对位移变化较参考点N6幅度更大。x向上,都表现为随着入射角度的增大,相对位移也增大。y向上,整体上表现为随着入射角度从−30°向30°变化,相对位移逐步增大,但硬岩场地上不明显。z向上,整体上相对位移随着入射角度的增大而增大。

    • 本文基于斜入射地震动作用下土-结构系统的整体有限元方法,分析了非基岩分层场地和均质硬岩场地条件下P波斜入射对“国和一号”核岛结构地震反应的影响。研究表明:

      (1) 不同场地对核岛结构的影响明显。非基岩场地将增大核岛结构加速度响应和楼层反应谱,多数条件下也将增大楼层相对位移。两种场地下的结构加速度响应、楼层加速度反应谱和楼层相对位移随入射角度变化的规律基本一致。

      (2) 整体来看,随着入射角度的增大,核岛结构水平响应增大,竖向响应减小,楼层反应谱的频谱特性也发生改变。

目录

    /

    返回文章
    返回