Earthquake Response Analysis of Soil-reactor Plant-nuclear power Auxiliary Plant Interaction System
-
摘要: 土-结构相互作用分析是核电结构抗震设计的重要环节,考虑到附属厂房可能导致反应堆厂房处于最不利工况状态,对地基土-反应堆厂房-核电辅助厂房结构相互作用体系地震响应进行研究。基于PASSI算法,提出显-隐式单元层计算方法,实现显、隐式交替计算,保持高效性的同时,提高计算稳定性。通过场地分析算例和土-结构相互作用分析算例,与ABAQUS软件计算结果进行对比,验证计算方法的可行性。以某核电站为对象,分析同一基础上相邻厂房对反应堆厂房地震响应的影响。研究结果表明,在基岩场地上,安全厂房通过基础和场地对反应堆厂房的作用较小,对反应堆厂房地震响应的影响较小;地震动输入下,燃料厂房和电器厂房加大了反应堆厂房位移峰值,减小了反应堆厂房加速度反应谱峰值,并使反应堆厂房顶部点加速度反应谱峰值向高频移动;与反应堆厂房共用同一基础的辅助厂房,应与反应堆厂房作为整体进行地震响应分析,至少应将与反应堆厂房相连的辅助厂房作为整体进行分析。Abstract: Soil-structure interaction analysis is an important step in seismic design, Considering that the existence of the auxiliary plant may cause the reactor plant in the most unfavorable condition, In the paper the coupled dynamic response of soil- reactor plant-auxiliary plant structure system to earthquake excitation is studied. Based on PASSI, a computation scheme of explicit and implicit element layer is proposed, which can realize explicit and implicit computation alternately, and improve the computation stability with high efficiency. The site analysis and soil-structure interaction analysis examples were compared with the results of ABAQUS commercial software, and the feasibility of the calculation scheme was verified. On this basis, the influence of adjacent plants on the seismic response of a nuclear power plant on the same basis was analyzed. The results show that on the bedrock site, the safety plant has little effect on the reactor plant through the foundation and site, It also has little effect on the seismic response of the reactor plant. Under the seismic wave input, the fuel and electrical plant increased the displacement peak of the reactor plant, reduced the peak of the acceleration response spectrum, and made the peak of the acceleration response spectrum of the top point of the reactor plant move to high frequency. The auxiliary plant sharing the same foundation with the reactor plant should be analyzed holistically, and at least the auxiliary plant connected with the reactor plant should be analyzed holistically.
-
Key words:
- Nuclear power plant /
- SSI /
- Explicit-implicit element method /
- PASSI /
- Adjacent structure
-
引言
汶川地震和雅安地震的严峻灾情和惨重损失充分反映出我国城镇及农村地区在抗御地震灾害能力方面的欠缺(张培震等,2013)。地震中75%的伤亡是由建筑结构的倒塌破坏造成(Coburn等,2002)。《国家中长期科学和技术发展规划纲要(2006—2020)》中提出:“发展完善地震参数速报与灾情信息快速收集发布技术,为决策层和社会公众提供快速准确的地震及相关灾情信息服务”。承灾体数据的准确性和时效性直接决定了震前预防监控、震时救援和震后重建工作的质量与效率(Sahar等,2010),在防灾与应急处置一体化体系中处于最基础和最关键的地位,是连接地震防御和辅助决策等模块的纽带和核心动力。
在承灾体数据采集领域有3个问题长期困扰地震工作人员:
(1)建筑抗震数据标准的设立。《震害防御信息服务数据采集规范》给出了国家参考标准(国家地震社会服务工程震害防御规范编写组,2011),但随着不同区域承灾体及其脆弱性呈多元化趋势,建筑抗震数据标准需要因地制宜得到进一步地改进(安基文等,2015)。定制化的数据采集标准可以实现更加精确的灾害预测(陈洪富,2012)。因此,本文根据西南地区县城实际环境预设采集项目,定制承灾体表单,使得承灾体数据采集的过程更加科学、规范,使其更适用于我国西南乡镇地区的实际环境,如本系统的示范区四川省丹棱县。
(2)数据采集手段效率的提升。现有的以调查人员现场调查填写表格为主的承灾体信息采集手段,对调查人员专业水平要求较高,且往往难以获得标准化的数据,后期数据汇总耗时耗力(FEMA等,2002)。基于卫星图像的采集方法成本较高,间隔较长,不够灵活(Satio等,2011;Wieland等,2012)。随着便捷智能终端的普及和互联网通讯技术的发展,手机的作用不再限于通信工具,而是更进一步成为采集人员的工作助手。在终端选择上,有定制PDA终端(张翼等,2014)、基于SETP的GNSS终端(黄经国等,2017)、基于iOS系统的iPad终端(徐柳华等,2012)、基于Windows FMC的ROVER终端(Annunziato等,2010)、基于Android的IDCT终端(龙立等,2016)等。本文的创新之处在于,结合西南县城地区实际情况,针对Android与iOS平台开发APP。使用者无需对硬件和操作系统进行定制化修改,只需下载采集APP即可使用,且APP兼顾Android与iOS两个平台,保证了终端使用的兼容性和便捷性。
(3)数据时效性与有效性的保证。承灾体数据难以有效实时更新,导致震害预测和地震风险评估工作的时效性和准确性难以保障(王铮,2015)。承灾体数据缺乏有效的认证确认机制,难以发现外业人员疏忽下的数据录入错误。因此本文设计新的承灾体数据更新机制,以现有的城镇规划数据为基础,结合数据采集终端软件传回的数据进行迭代和更新,最终将该数据更新机制集成在后台数据更新与数据管理系统中。该数据更新系统有效解决了城镇规划数据体量大、更新慢的问题(张云龙,2006),城镇规划数据和采集终端数据二者互为补充,辅以自动化交叉验证,保障了承灾体数据准确及时地入库。手机端数据采集终端软件和后台数据更新与数据管理系统相辅相成,构成了一套西南地区承灾体数据采集与动态更新系统,使得承灾体数据从采集到入库更加智能和易于管理,为地震风险评估预测、灾后应急提供强力支持。
1. 系统架构和系统功能
“西南地区承灾体数据采集与动态更新系统”是集数据采集、校验、动态更新为一体的综合系统,该系统采用C/S架构,分为数据采集终端软件(Android及iOS版本)和后台数据更新与数据管理系统两部分,系统架构图见图 1。
数据采集终端软件部署在外业人员手持的Android或iOS系统移动终端上,对接图层数据库,提供基于地图的标准化采集功能,以确保采集数据的准确性和一致性。
后台数据更新与数据管理系统在以数据采集终端软件为主要数据来源的同时,对接城镇规划数据库,拓展第二数据来源渠道,二者交叉验证,相互补充,提高数据的可靠性和准确性。所有数据将于后台系统进行格式化、标准化处理和自动校验,剔除采集重复或不合格数据,将异常数据标注,通过人工审核进行删除或修改,减少采集误差,以保证数据的质量。最终承灾体数据库对通过自动校验和人工审核的数据进行更新,对接图层数据库,使图层数据得到同步更新。同时,数据库将为灾情综合分析处理系统提供接口,为震害预测和风险评估系统提供及时标准的承灾体信息。
承灾体数据采集更新的业务流程为:系统首先尝试对接城镇规划数据,若该地区拥有城镇规划数据,则用其进行数据初始化。之后以智能终端采集和城镇规划数据更新记录为信息渠道,进行数据采集。系统在接受到新采集的数据后,对数据进行标准化、数据校验、人工审核等处理,最终将实时数据录入数据库,并生成更新记录,与图层数据库动态更新。
2. 关键技术点
2.1 采集标准化
西南乡镇地区承灾体属性与其他地区相比有差异性,并且经历多次地震灾后重建,承灾体种类及其脆弱性呈多元化趋势,传统的通用数据标准已经无法适用于西南地区(张斌等,2010)。为了对承灾体数据采集标准进行因地制宜地改进,同时最大程度采集承灾体的各方面特性,本系统借助承灾体信息采集终端对采集流程和采集项目进行了规范和定制优化,解决了以往采集工作与当地情况结合不紧密的问题。
2.1.1 采集项目标准化
承灾体数据包括:位置信息、建筑信息、其他3大类。
位置信息包括:坐标(经度和纬度)、所属区(县)、所属乡镇(街道办、社区)、行政村、门牌号。
建筑信息包括:名称、层数、结构类型、平立面规整性、建设年代、设防标准、结构现状、用途、是否为重要建筑。
其他信息包括:系统编码、承灾体编码、备注、承灾体图片(正面照、侧面照、反面照以及认为有必要补充的照片)、采集时间。
作为承灾体数据的重点部分,建筑信息的具体采集字段和取值如表 1所示。
表 1 建筑物信息采集内容Table 1. Information collected for buildings序号 采集字段 字段取值 1 名称 手动输入 2 层数 手动输入 3 结构类型 砖混结构、框架结构、框剪(框筒)结构、砖木结构、土木结构、钢结构、混凝土小砌块结构、木结构、穿斗木结构、生土结构、石结构、其它 4 平立面 规整,不规整 5 建设年代 1979年以前、1980年至1989年、1990年至2000年、2001年至2008年、2009至2013年、2013年以后 6 设防标准 未设防、Ⅵ度设防、Ⅶ度设防、Ⅷ度设防、Ⅸ度设防 7 结构现状 完好、开裂、腐蚀、不均匀沉陷 8 用途 住宅、办公、商住、商业、其它、已拆 9 重要建筑 政府部门、学校、医院、公安(消防、武警)、大型公共场所 在制定承灾体数据采集标准时,本文在国家地震社会服务工程项目震害防御信息服务系统数据采集规范的基础上,进行城镇地震防灾与应急处置一体化服务系统数据需求分析,做出多项改进:
(1)针对西南地区承灾体结构特征,在结构类型中添加了我国西南山区常见的穿斗木结构。
(2)将震害防御信息服务系统数据采集规范中年代划分的“2001年以后”细分为“2001年至2008年”、“2009年至2013年”和“2013年以后”,该做法主要是考虑2008年汶川地震和2013年芦山地震后西南地区承灾体分别有过大规模的更新。
(3)在重要建筑中,系统表单保留了“政府部门”、“学校”、“医院”、“公安(消防、武警)”和“大型公共场所”5类,参考地震灾害预测重要建筑中的“体育馆”、“大型影剧院”、“大型百货商场”、“大型酒店”等都统一为“大型公共场所”,而震害防御信息服务系统数据采集规范中的“金融、通信、交通、能源、水利”等项,则考虑到县域地区的实际情况,删除或合并在“政府部门”该项中,这样的处理不仅进一步简化了数据采集过程中的分析判断过程,也减少了判断失误的可能性。
2.1.2 采集流程标准化
首先,外业采集人员通过专有账号和密码登录数据采集终端软件,查看基于地图的整体承灾体采集状况。账号密码由后台系统添加,保障访问权限安全。地图以蓝色标识已采集承灾体,红色标识未采集承灾体,橘色标识本次作业采集还未上传的承灾体。通过地图,外业采集人员可以迅速了解采集区域情况,制定采集策略和路线。采集人员可以选择已采集承灾体对信息进行更新,也可选择未采集承灾体添加新的数据。
其次,外业人员通过在APP地图上点击相应承灾体进入标准化采集表单,对信息进行录入。数据采集后缓存在本地,当网络通畅时,数据将通过GSM/GPRS/3G/4G移动网络或Internet上传至后台数据更新与数据管理系统,向管理人员实时反馈采集员的在线情况、位置以及数据的采集情况。
最后,后台数据更新与数据管理系统对进入承灾体数据库的数据初始化,随后以来自采集终端的承灾体数据和城镇规划数据的更新记录为凭证,对新采集的承灾体数据进行数据校验和人工审核,最终数据入库,完成数据采集。
具体采集流程如图 3所示。
2.2 校对自动化
管理人员可以根据当地承灾体分布情况,因地制宜地设定校验条件。当后台数据更新与数据管理系统接收到来自信息采集终端或城镇规划中心的承灾体数据时,会按照设定的关联规则自动对全体数据进行校验处理(史培军等,2014),校验通过的数据将自动添加或更新覆盖进数据库,校验失败的数据将转人工处理,由专业工作人员进行检查和更新。
相比以纸质表格为载体的信息校对方式,自动校对模块降低了人工纸质校对的时间成本和工作量,并且提供更加快捷的后续删改功能,使得系统能够根据建筑物数据的变更动态更新。
2.3 离线缓存数据
2.3.1 地图数据离线缓存
在登录成功后,APP中需要设置数据传输的服务器IP和端口号以连接后台数据更新与数据管理系统。考虑到该套系统示范应用的地区属于县城,当手机无法连接网络时,还可以通过离线模式登录,节省流量、避免并行压力,提高系统效率。
该地图的最底层是天地图,叠加有丹棱县承灾体基础数据图层。地图文件会在网络通畅的情况下自动进行更新和缓存,离线模式下外业采集人员使用缓存的地图数据进行采集工作。
2.3.2 采集数据离线缓存
若数据采集终端软件处于离线模式,则采集完毕的承灾体信息会被保存在手机中,当网络质量好或者采集完成后,进行统一上传。同时,离线采集完毕的承灾体会在地图上以橘色表示。
3. 实例演示
3.1 数据采集终端软件
(1)登录与退出系统
打开软件,输入用户名密码,设置数据传输的服务器IP和端口号与管理端连接(图 4)。
连接成功,用户进入地图界面,软件自动进行地图数据和用户数据的离线缓存,并将用户的实时位置、操作记录等信息传输至管理端。
在采集过程中,采集员数次模仿网络不畅通情景,关闭Wi-Fi和流量,进入离线模式。离线模式下,采集员使用的是基于最新缓存数据的地图界面。
(2)现场采集
地图上,已采集的承灾体以蓝色标识,而未采集的承灾体以红色标识,采集员通过观察地图制定采集策略与路线。
采集员到达龙鹄村,在地图上选中承灾体,点击右上角采集按钮,进入采集页面。
采集员在采集界面填写表单录入承灾体信息,点击“正面”,“反面”,“侧面”自动开启系统相机界面对承灾体进行拍照,最终点击保存完成采集(图 5)。
(3)数据的统计与缓存
采集若干建筑物后,采集员进入采集统计页面查看承灾体的采集概况。采集统计界面将承灾体数据分为全部、已采集、更新、新增4类,操作员选择不同的选项,软件同步显示不同的分类列表以及数量统计。
离线采集的承灾体以本地缓存数据的方式,呈黄色显示在离线缓存地图上,避免重复采集。
在网络通畅的情况下,承灾体数据会自动在后台被上传到服务器。
在网络堵塞的情况下,采集员可在缓存提交界面提交或清除缓存的采集数据(图 6)。
3.2 后台数据更新与数据管理系统
(1)数据初始化
操作员控制系统连接城镇中心规划数据库,以城镇规划数据为基础对承灾体数据库进行初始化,并且根据城镇规划数据更新记录对承灾体数据库进行更新。
(2)实时监控
当采集员使用连线模式时,后台操作人员借助后台数据更新与数据管理系统对采集员工作状况进行实时监控(图 7)。
操作员于页面左侧显示采集员列表与基本信息,于地图查看采集员当前位置,于地图下方两个显示框分别查看承灾体数据及用户状态更新提示。于页面右下角查看当前在线人数。
(3)采集数据的管理
该模块实现对采集数据的显示、查询、删除、编辑等管理功能。其中,数据编辑功能可以对采集数据的值进行修改,自定义条件查询功能可以通过设置字段值范围查询所需的数据。操作员点击“导出当前列表”,将数据以CSV格式导出。
(4)审核管理
审核管理模块具有自动校验及人工审核两大功能。自动校验通过预先设定的规则,对所有采集的数据或已校验的数据进行自动校验,如有不合格的数据,点击校验结果中的“不合格”,弹出错误报告,即可得知错误的信息。例如输入校验条件“1979年以前的建筑物不能高于3层”,则所有不满足条件的建筑会被标为不合格,那么采集员误输入的数据就可被检查出来。
后台操作员可在人工审核模块对通过自动校验的数据进行二次确认,可修改不合格的数据,点击审核通过按钮,提交所有已审核完成的数据,最终将数据录入数据库(图 8)。
(5)用户管理
后台操作人员在系统中操作生成采集员第一次的用户信息,之后用户则可通过用户名密码登录编辑个人信息(图 9)。
操作员于后台系统查看并修改用户基本信息,从而保证前端APP的用户可控和系统安全,防止非授权用户非法入侵和授权用户的越权访问。
(6)数据维护
操作员于后台系统查看自动生成的数据库备份记录、业务统计报表以及数据库日志,以了解业务进展概况。
3.3 示范应用总结
在丹棱县的示范应用中,数据采集终端软件单次单栋建筑物的平均采集时间控制在3分钟以内,单次作业(2小时)平均采集建筑物25座,连线连通率100%,极大地提高了外出采集人员的工作效率和上报数据的时效性、规范性。后台数据分析管理中心与地震部门日常工作紧密结合,通过用户管理和用户监控,显著地改善和提升了部门的管理服务水平。在数据校验环节,承灾体数据自动化校验通过率为98%,数据一致性100%,千组数据平均耗时2秒,显著地减轻了审核人员的工作负担,使得数据汇总工作更加方便、高效、智能。
4. 结论
本文系统地介绍了西南地区承灾体数据采集与动态更新系统的设计方法与实例演示,得出了以下主要结论:
(1)准确、实时的承灾体数据是震害预测、防御工作的基石。强震区县市地震工作部门传统的手工表单填写采集及人工档案整理的方式耗时耗力,数据一致性和实效性差,致使后续灾情预测工作缺乏可靠的数据支持,服务能力弱。
(2)本系统是“通用终端+专用数据采集软件+标准化采集流程”与城镇规划数据联网更新有机结合而开发的专用数据采集平台,旨在借助标准化采集项目增强数据质量,保障数据一致性;通过数据采集终端软件和优化采集流程加快数据采集速度,提高采集效率;通过对接城镇规划数据库提高数据体量;通过定制化的交叉检验审核提高数据准确性。
(3)示例应用数据表明,本系统显著提高了承灾体数据的采集、更新速度及数据质量,降低了人力和时间的消耗,实现了在最短时间用最小消耗的情况下采集最实时准确的承灾体数据的目标。切实增强了县市防震减灾工作科技水平和支撑服务能力,为县市震害防御、地震应急、政府规划工作提供实时可靠的基础数据和辅助决策信息支撑。
(4)系统设计时,充分考虑了产品与其他系统的集成及调用规划数据,为未来继续开发的功能提供了清晰的二次开发环境与接口。同时,系统的多个可定制模块决定了其优秀的可拓展性,能够因地制宜适应当地承灾体环境。
(5)示例应用的成功,在完善更新机制、规范采集标准、提高数据质量、保障数据时效性、提高系统针对性等方面为后续系统在西部地震重点监视防御区县市的推广提供了宝贵经验,加快了系统在西部地震重点监视防御区县市的推广进程。在不久的将来,有望在平台上集成整个西南地区承灾体数据,以打通不同地区信息壁垒,保障重点灾区的重点监控,切实提高西部强震区县市防震减灾能力。
-
表 1 土体参数
Table 1. Soil parameters
材料 厚度/m 弹性模量/GPa 泊松比 密度/kg·m−3 剪切波速/m·s−1 压缩波速/m·s−1 软土 50 0.108 0.35 1 000 200 416 表 2 显-隐式单元层计算方法和ABAQUS软件计算效率
Table 2. PASSI and ABAQUS calculation efficiency
算例 单元数/个 节点数/个 自由度数/个 显-隐式单元层计算方法
计算时间/minABAQUS软件
计算时间/min场地 80 000 85 731 257 193 4 106 土-结构相互作用 80 010 85 771 257 313 33 106 注:自由度数=节点数*3(每个节点自由度个数) 表 3 核电站结构材料参数
Table 3. Material parameters of nuclear power plant
编号 材料 结构 弹性模量/GPa 泊松比 密度/kg·m−3 1 C30 厂房/基础 32.5 0.2 2 400 2 C40 内/外壳 36.0 0.2 2 450 表 4 核电站厂址参数
Table 4. Soil parameters of nuclear power plant site
材料 厚度/m 弹性模量/GPa 泊松比 密度/kg·m−3 剪切波速/m·s−1 压缩波速/m·s−1 岩石 60 46.9 0.26 2 650 2 673 4 639 表 5 模态分析
Table 5. Modality analysis
振型 Case1 Case2 Case3 Case4 频率/Hz 振型参与系数 频率/Hz 振型参与系数 频率/Hz 振型参与系数 频率/Hz 振型参与系数 1 3.006 69 0.186 26 3.006 69 0.186 19 4.016 96 1.000 00 3.726 48 0.426 64 2 3.007 49 1.000 00 3.007 49 1.000 00 4.068 74 0.277 89 3.861 54 0.413 38 3 4.078 57 0.540 83 3.725 20 0.677 04 4.114 50 0.094 62 3.885 43 0.002 48 4 4.096 08 0.381 57 3.858 00 0.683 88 4.210 24 0.482 59 4.017 18 1.000 00 5 5.488 44 0.002 04 3.885 28 0.009 75 4.739 49 0.029 30 4.068 75 0.275 37 6 5.513 01 0.001 43 4.078 57 0.540 69 5.472 46 0.169 67 4.114 53 0.097 12 7 5.540 46 0.003 31 4.096 08 0.381 50 5.566 56 0.002 17 4.210 38 0.480 55 8 5.572 03 0.004 84 4.559 21 0.068 57 5.610 47 0.012 11 4.560 35 0.045 69 9 5.629 19 0.003 16 5.439 43 0.064 35 5.762 07 0.080 10 4.739 62 0.029 71 10 5.658 22 0.010 71 5.488 44 0.002 05 6.163 02 0.037 35 5.441 09 0.042 47 11 6.426 34 0.004 54 5.513 01 0.001 43 6.470 85 0.006 88 5.472 66 0.167 66 12 6.471 34 0.015 75 5.540 46 0.003 31 6.496 26 0.389 72 5.566 56 0.002 12 13 6.496 20 0.598 36 5.572 03 0.004 82 6.615 10 0.003 04 5.610 47 0.011 96 14 6.614 21 0.001 58 5.629 19 0.003 17 6.875 02 0.134 96 5.702 48 0.044 00 15 6.875 13 0.207 35 5.658 22 0.010 72 7.121 85 0.000 71 5.762 18 0.079 50 16 7.117 11 0.000 87 5.698 29 0.070 07 7.183 56 0.005 66 6.163 16 0.037 35 17 7.171 63 0.005 92 6.426 34 0.004 55 7.186 37 0.071 54 6.470 85 0.006 79 18 7.186 59 0.070 52 6.471 34 0.015 75 7.228 46 0.359 01 6.496 26 0.386 15 19 7.218 46 0.014 67 6.496 20 0.598 31 7.349 15 0.546 56 6.615 10 0.003 03 20 7.235 54 0.000 05 6.614 21 0.001 59 7.778 06 0.042 94 6.875 02 0.133 73 21 7.778 05 0.071 15 6.875 13 0.207 34 7.984 62 0.026 00 7.121 85 0.000 71 22 7.984 60 0.042 40 7.117 11 0.000 87 8.112 55 0.020 30 7.183 56 0.005 58 23 8.112 54 0.033 21 7.171 63 0.005 93 8.241 69 0.004 93 7.186 37 0.070 85 24 8.241 67 0.008 19 7.186 59 0.070 51 8.355 48 0.012 78 7.228 55 0.356 06 25 8.355 47 0.020 32 7.218 46 0.014 67 8.603 98 0.003 84 7.349 13 0.541 36 26 8.457 00 0.005 00 7.235 54 0.000 04 8.879 18 0.000 60 7.778 06 0.042 54 27 8.603 97 0.006 27 7.778 05 0.071 15 8.980 48 0.011 40 7.984 62 0.025 76 28 8.977 70 0.014 41 7.984 60 0.042 39 9.092 81 0.031 82 8.112 55 0.020 11 29 9.092 80 0.051 76 8.112 54 0.033 21 9.161 85 0.066 65 8.241 69 0.004 88 30 9.229 79 0.118 33 8.241 67 0.008 19 9.229 87 0.071 67 8.355 48 0.012 67 31 9.414 48 0.015 02 8.355 47 0.020 32 9.414 48 0.009 21 8.603 98 0.003 80 32 9.518 82 0.077 79 8.457 00 0.005 04 9.450 98 0.026 21 8.879 18 0.000 58 33 9.698 92 0.086 25 8.603 97 0.006 27 9.518 87 0.045 52 8.980 48 0.011 29 34 9.776 22 0.456 36 8.977 70 0.014 42 9.787 15 0.014 03 9.092 81 0.031 53 35 9.841 21 0.024 00 9.092 80 0.051 76 9.841 17 0.020 96 9.161 81 0.065 98 -
曹金凤, 王旭春, 孔亮, 2011. Python语言在Abaqus中的应用. 北京: 机械工业出版社. 韩冰, 陈少林, 梁建文, 2019. 结构-土-结构动力相互作用对结构系统频率的影响. 地震工程学报, 41(6): 1574—1580Han B. , Chen S. L. , Liang J. W. , 2019. Effects of structure-soil-structure dynamic interaction on the frequency of structural systems. China Earthquake Engineering Journal, 41(6): 1574—1580. (in Chinese) 何涛, 姜南, 2020. 土-相邻结构相互作用子结构振动台试验研究. 振动与冲击, 39(4): 207—214He T. , Jiang N. , 2020. A substructure shaking table test of soil-adjacent structure interaction. Journal of Vibration and Shock, 39(4): 207—214. (in Chinese) 李培振, 严克非, 徐鹏, 2014. 地震下考虑群体效应的高层建筑土-结构相互作用研究. 土木工程学报, 47(S1): 1—5Li P. Z. , Yan K. F. , Xu P. , 2014. Study on dynamic interaction between soil and group of high-rise buildings under seismic excitation. China Civil Engineering Journal, 47(S1): 1—5. (in Chinese) 廖振鹏, 2002. 工程波动理论导论. 2版. 北京: 科学出版社.Liao Z. P. , 2002. Introduction to wave motion theories in engineering. 2 nd ed. Beijing: Science Press. (in Chinese) 刘晶波, 王振宇, 杜修力等, 2005. 波动问题中的三维时域粘弹性人工边界. 工程力学, 22(6): 46—51 doi: 10.3969/j.issn.1000-4750.2005.06.008Liu J. B. , Wang Z. Y. , Du X. L. , et al. , 2005. Three-dimensional visco-elastic artificial boundaries in time domain for wave motion problems. Engineering Mechanics, 22(6): 46—51. (in Chinese) doi: 10.3969/j.issn.1000-4750.2005.06.008 刘晶波, 杜义欣, 闫秋实, 2007. 粘弹性人工边界及地震动输入在通用有限元软件中的实现. 见: 第三届全国防震减灾工程学术研讨会论文集. 南京: 中国土木工程学会, 43—48. 刘晶波, 宝鑫, 李述涛等, 2022. 采用黏弹性人工边界时显式算法稳定性条件. 爆炸与冲击, 42(3): 034201Liu J. B. , Bao X. , Li S. T. , et al. , 2022. Stability conditions of explicit algorithms when using viscoelastic artificial boundaries. Explosion and Shock Waves, 42(3): 034201. (in Chinese) 柳玉印, 尹训强, 2018. 考虑结构—土—结构动力相互作用的重力坝地震响应分析. 水利水电技术, 49(10): 52—58Liu Y. Y. , Yin X. Q. , 2018. Structure-soil-structure dynamic interaction-considered seismic response analysis of gravity dam. Water Resources and Hydropower Engineering, 49(10): 52—58. (in Chinese) 苏景鹤, 江丙云, 2016. ABAQUS Python二次开发攻略. 北京: 人民邮电出版社. 田彼得, 俞载道, 1987. 结构-土-结构相互作用体系的动力分析. 同济大学学报, 15(2): 157—168Tian B. D. , Yu Z. D. , 1987. Dynamic analysis of structure-soil-structure interaction. Journal of Tongji University, 15(2): 157—168. (in Chinese) 王飞, 宋志强, 刘昱杰等, 2018. 基于ABAQUS无限元的静-动力统一人工边界研究. 水资源与水工程学报, 29(6): 170—177Wang F. , Song Z. Q. , Liu Y. J. , et al. , 2018. Research on static-dynamic unified artificial boundary based on ABAQUS infinite element. Journal of Water Resources and Water Engineering, 29(6): 170—177. (in Chinese) 王国波, 袁明智, 苗雨, 2018. 结构-土-结构相互作用体系地震响应研究综述. 岩土工程学报, 40(5): 837—847Wang G. B. , Yuan M. Z. , Miao Y. , 2018. Review of seismic response of structure-soil-structure interaction system. Chinese Journal of Geotechnical Engineering, 40(5): 837—847. (in Chinese) 应稼年, 王荣昌, 俞载道, 1995. 地基土-桩基础-核电站辅助厂房结构相互作用体系的地震响应分析. 地震工程与工程振动, 15(1): 44—52 doi: 10.13197/j.eeev.1995.01.006Ying J. N. , Wang R. C. , Yu Z. D. , 1995. Earthquake response analysis of soil-pile-nuclear power station auxiliary workshop structure interaction system. Earthquake Engineering and Engineering Vibration, 15(1): 44—52. (in Chinese) doi: 10.13197/j.eeev.1995.01.006 中华人民共和国住房和城乡建设部, 国家市场监督管理总局, 2019. GB 50267—2019 核电厂抗震设计标准. 北京: 中国计划出版社.Ministry of Housing and Urban Rural Development of the People's Republic of China, State Administration of Market Supervision and Administration, 2019. GB 50267—2019 Standard for seismic design of nuclear power plants. Beijing: China Planning Press. (in Chinese) Aldaikh H. , Alexander N. A. , Ibraim E. , et al. , 2015. Two dimensional numerical and experimental models for the study of structure–soil–structure interaction involving three buildings. Computers & Structures, 150: 79—91. Aldaikh H. , Alexander N. A. , Ibraim E. , et al. , 2016. Shake table testing of the dynamic interaction between two and three adjacent buildings (SSSI). Soil Dynamics and Earthquake Engineering, 89: 219—232. doi: 10.1016/j.soildyn.2016.08.012 ASCE, Seismic analysis of safety-related nuclear structures, ASCE/SEI 4-16. 2016, American Society of Civil Engineers: New York Belytschko T. , Mullen R. , 1978. Stability of explicit-implicit mesh partitions in time integration. International Journal for Numerical Methods in Engineering, 12(10): 1575—1586. doi: 10.1002/nme.1620121008 Belytschko T., Yen H. J., Mullen R., 1979. Mixed methods for time integration. Computer Methods in Applied Mechanics and Engineering, 17—18: 259—275. Belytschko T. , Smolinski P. , Liu W. K. , 1985. Stability of multi-time step partitioned integrators for first-order finite element systems. Computer Methods in Applied Mechanics and Engineering, 49(3): 281—297. doi: 10.1016/0045-7825(85)90126-4 Chen S. L. , Lv H. , Zhou G. L. , 2022. Partitioned analysis of soil-structure interaction for Nuclear Island Buildings. Earthquake Engineering & Structural Dynamics, 51(10): 2220—2247. Ghandil M. , Behnamfar F. , Vafaeian M. , 2016. Dynamic responses of structure–soil–structure systems with an extension of the equivalent linear soil modeling. Soil Dynamics and Earthquake Engineering, 80: 149—162. doi: 10.1016/j.soildyn.2015.10.014 Kitada Y., Hirotani T., Iguchi M., 1999. Models test on dynamic structure–structure interaction of nuclear power plant buildings. Nuclear Engineering and Design, 192(2—3): 205—216. Lee T. H. , Wesley D. A. , 1973 a. Soil-structure interaction of nuclear reactor structures considering through-soil coupling between adjacent structures. Nuclear Engineering and Design, 24(3): 374—387. doi: 10.1016/0029-5493(73)90007-1 Lee T. H., Wesley D. A., 1973b. Influence of through-soil coupling between adjacent structures on seismic response of nuclear reactors. In: Proceedings of the Second International Conference on Structural Mechanics in Reactor Technology. Berlin: IASMiRT, 2—9. Lehmann L. , Antes H. , 2001. Dynamic structure-soil-structure interaction applying the Symmetric Galerkin Boundary Element Method (SGBEM). Mechanics Research Communications, 28(3): 297—304. doi: 10.1016/S0093-6413(01)00177-X Lou M. L. , Wang H. F. , Chen X. , et al. , 2011. Structure-soil-structure interaction: literature review. Soil Dynamics and Earthquake Engineering, 31(12): 1724—1731. doi: 10.1016/j.soildyn.2011.07.008 Luco J. E. , Contesse L. , 1973. Dynamic structure-soil-structure interaction. Bulletin of the Seismological Society of America, 63(4): 1289—1303. doi: 10.1785/BSSA0630041289 Matthees W. , Magiera G. , 1982. A sensitivity study of seismic structure-soil-structure interaction problems for nuclear power plants. Nuclear Engineering and Design, 73(3): 343—363. doi: 10.1016/0029-5493(82)90011-5 Murakami H. , Luco J. E. , 1977. Seismic response of a periodic array of structures. Journal of the Engineering Mechanics Division, 103(5): 965—977. doi: 10.1061/JMCEA3.0002286 Wang J. S. , Guo T. , Du Z. Y. , 2022. Experimental and numerical study on the influence of dynamic structure-soil-structure interaction on the responses of two adjacent idealized structural systems. Journal of Building Engineering, 52: 104454. doi: 10.1016/j.jobe.2022.104454 Warburton G. B. , Richardson J. D. , Webster J. J. , 1971. Forced vibrations of two masses on an elastic half space. Journal of Applied Mechanics, 38(1): 148—156. doi: 10.1115/1.3408735 Whitman R. V. , 1969. The current status of soil dynamics. Applied Mechanics Reviews, 22(1): 1—8. Wong H. L. , Trifunac M. D. , 1975. Two-dimensional, antiplane, building-soil-building interaction for two or more buildings and for incident planet SH waves. Bulletin of the Seismological Society of America, 65(6): 1863—1885. -