贮液结构是一种极为重要的结构,与人民群众日常生产和生活息息相关。随着贮液不同,贮液结构广泛应用于生命线工程和能源与化工领域,如水塔、水箱、储油罐、储酸罐、消化池、降温池、冷凝器等,其结构形式也多种多样,包括半地下式、立式、卧式、支承式等。由于存在液固耦合作用,这种结构动力反应极为复杂,尤其是强地震动作用下,引起的震害现象千差万别。历次震害表明,高烈度区大量贮液结构受到破坏,如浮放支座移位、支承结构移位、基底开裂、结构壁断裂等等,这些地震灾害直接导致贮液结构损坏,其原有功能丧失,生产和生活被迫中断,造成极大的经济损失。此外,贮液结构内部的贮液有时也是极大危险源,贮液结构受到破坏后,容易引起贮液泄露,进而引发火灾、爆炸、环境污染等次生灾害。2011年3月11日日本东海岸发生9.0级地震,地震和随后的海啸造成福岛核电站中的贮液结构系统受损,致使降温系统失去功能,引发了重大核泄漏事故。从上世纪三十年代开始,大批学者针对贮液结构地震反应进行了相关研究,其中包括建立理论分析模型、液固耦合作用效应、简化方法工程应用、地震反应规律分析、抗震设计方法和数值分析技术等等。但是从工程应用角度,贮液结构地震反应仍然有很多课题需要去研究,包括有限元建模方法、非线性算法、地震反应规律分析、工程抗震分析应用方法、优化设计等等。本文基于此背景,针对贮液结构地震反应及相关数值分析方法进行了研究,主要工作和取得的成果如下：1.总结了贮液结构解析解难点是求解给定边界条件和初始条件下的Laplace方程,弹性壁解析解和Housner简化计算均是建立在脉冲压力和对流压力理论假设基础上的；Haroun-Housner模型和有限元数值模型分别是建立在Housner刚性壁理论和速度势理论基础上的,本文对四种方法进行了对比,说明了有限元数值模型的可靠性。2.基于数值迭代理论,进行了贮液结构非线性数值迭代计算方法的研究。利用经典的N-R数值迭代理论,提出改进N-R算法,考虑切线刚度和法线刚度共同协作的方式,增加单步迭代步长,提高非线性数值迭代求解效率；基于改进N-R算法,提出超N-R算法,建立数值迭代和材料本构之间的关联性,进一步提高非线性数值迭代求解效率,并通过大量试算给出对于钢材,比例系数77较为稳定的取值范围在1.46左右；并针对贮液结构非线性数值分析方法,提出采用外挂开发模式和生死单元技术相结合的方法进行数值模拟,通过举例,验证该方法的可行性和有效性。这两种算法极大地提高了贮液结构非线性数值反应分析的计算效率。3.基于APDL语言,进行了有限元建模数值二次开发研究。提出外挂模式的有限元二次开发方法,成功实现了在ANSYS平台上进行用户自定义子程序编制和调用；基于该方法,有效地简化了利用ANSYS进行复杂贮液结构地震反应建模、计算、数据处理等工作,便于工程技术人员应用。本文的有限元二次开发工作包括：提出针对大型复杂的有限元模型,进行三维梯度网格优化划分方法,以提高网格划分质量；提出编制针对有限元结构信息存储、读取、批处理等程序,利用波前法技术,实现通过刚度矩阵修改等多种方法进行数据存储量和计算量的压缩,以提高有限元分析的高效性和实用性。4.基于模态影响线控制网格划分技术、ATOS转换程序的开发、圣维南原理子模型技术等多方面说明了贮液结构有限元分析模型的有效性；基于此,对贮液结构基本地震反应规律进行了研究,获取了贮液结构液固耦合体系特性,以及地震动输入、壁厚、贮液量、径高比等因素对贮液结构地震反应的影响；并研究了贮液结构弹性变形对液动压力、结构有效应力、基底弯矩、基底剪力等的影响,为贮液结构抗震分析和设计提供了理论依据。5.基于有限元数值分析和插值理论,对常规立式贮液结构“象足”效应问题进行了研究,通过编制了样条插值分析和拟合程序,提出对“象足”效应易发位置进行快速确定的函数形式,并通过大量计算,给出了拟合函数中各系数取值参照表。利用该方法,实现了快速确定“象足”效应位置,据此,可以一定程度上针对“象足”效应进行高效的抗震设计应用。6.基于ESO形状拓扑理论,进行了贮液结构轮廓形状优化设计方法的研究,介绍了基于Mises应力、刚度约束、频率控制的ESO优化算法；提出将单元非连续性分类方式引入到ESO理论中,可以有效地解决在贮液结构轮廓拓扑计算过程中,为提高效率而增加删除率BDR,时引起的结构刚度矩阵奇异、拓扑过程不连续问题,提高了贮液结构拓扑分析的稳定性。7.基于智能优化算法,进行了贮液结构参数优化设计方法的研究,对穷举法、ACO算法、PSO算法和DE算法进行了介绍和分析；提出将ACO方法和子结构方法进行互相嵌入使用,通过该方法可以有效减小贮液结构参数优化计算的计算量；提出了PSO-DE联合优化方法对,利用飞行、变异、交叉理论,扩大了全局搜索范围,避免了优化计算过早陷入局优,提高了优化算法效率,并加速了迭代收敛的速度。本文的研究工作为贮液结构地震反应研究提供了有效的数值分析方法,为贮液结构地震反应研究结果在工程中的应用提供了极大地便利性,为贮液结构工程抗震设计提供了重要的应用方法,为贮液结构防灾减灾工作提供了重要的理论和应用基础。