大坝结构静动力分析精细化模型.docVIP

大坝结构静动力分析精细化模型 　　摘要：阐述近期大连理工大学抗震研究所在大坝抗震分析方面所取得的研究成果，为了提高大坝结构静动力分析的计算精度与效率，提出多边形单元的计算技术与奇异应力场计算方法，该方法使应力计算成果较有限元法有很大的改进，且计算工作量比较节省。 　　关键词：静动力分析；多边形单元；奇异应力场 　　中图分类号：TU45 文献标识码：A 文章编号：1000-0666（2016）01-0001-09 　　0 前言 　　大坝结构的安全对城市抗震防灾有重要影响。2008年5月12日发生的汶川8.0级特大地震对位于成都市西北60km岷江上游的紫坪铺大坝造成了一定的地震损伤，威胁到成都市的安全。吉林市的安全与丰满大坝的安全息息相关，密云水库、官厅水库的安全也对北京市的安全产生一定影响。研究城市防灾需要关注周边大坝的安全设防。 　　新中国成立以来，特别是改革开放以来，我国水利水电建设迅速发展。进入21世纪，世界大坝建设中心已经转向中国。一大批世界级的高坝已经或正在中国大地上兴建，坝高屡破世界纪录，其成就举世瞩目。我国地处世界两大活跃地震带――环太平洋地震带和欧亚地震带的交汇地段，地震活动频繁。为适应地震危险区大坝建设的需要，我国的大坝抗震技术也蓬勃发展，并进入世界先进行列。 　　自20世纪50年代以来，大连理工大学工程抗震研究所一直根据我国大坝抗震的需要，致力于研究如何提高大坝抗震分析的水平。大坝抗震的特点在于大坝所处环境复杂（图1），参照世界坝工专家Chopra（1992，2008）的意见，大坝抗震的关键技术包括坝与库水的动力相互作用分析，坝与无限地基的动力相互作用分析，筑坝材料的动态特性及大坝地震损伤发展与风险分析。在2013年中国工程院举办的“中国工程科技论坛――水安全与水利水电可持续发展”会议上，林皋（2014）提出了大坝抗震模型的改进意见。关于坝与库水的动力相互作用分析需要考虑水的可压缩性以及动水压力波在水库边界的吸收作用，Chopra（2008）采用有限元法（FEM）进行求解，动水压力按坝的振动模态展开，计算复杂、难以实际应用。Lin等（2007a，2012a）采用比例边界有限元法（SBFEM）进行求解，按库水的模态进行展开，计算简便、精度高、便于实际应用。关于坝与无限地基的动力相互作用，由于计算模型比较复杂，文献中一般都将无限地基简化为均质地基进行考虑，难以反映实际情况。我们将无限地基分为层状分布、分块不均质和近场杂乱不规则等几种类型，分别提出了不同的计算模型（Linet al，2015，2007b；Yin et al，2013），计算结果都具有较好的精度与效率。为了提高坝-库水-地基系统时域分析的计算精度与效率，Lin等（2012b）提出以等几何分析方法（isogeometric a-nalysis，简称IGA）代替FEM来进行坝的应力分析，并将IGA与SBFEM相结合来进行坝与库水和坝与无限地基的动力相互作用分析。关于坝材料的动力特性方面，混凝土是速率敏感性材料。自从第一颗原子弹1945年在日本广岛爆炸以后，为了防卫核爆炸造成的灾难，许多国家都加强了混凝土材料抗击核冲击的研究。但这种研究偏重单调瞬时冲击，而且抗压性能居多，抗拉性能偏少。地震荷载的特点是不规则的循环往复，而且加载速率也偏低。为此，我们开展了2000多试件混凝土在地震作用下的动态特性的研究，取得了一批有价值的成果（Yan，Lin，2006；Lin et al，2007c）。 　　笔者认为除了计算模型以外，数值计算方法也对大坝抗震分析的计算精度与效率产生影响。大坝地震响应分析中，有限元方法的应用比较普遍。有限元法可以适应各种不规则的几何边界条件和分布不规律的单元材料特性，可以应用于各种类型物理问题的求解，计算能力强，不受问题求解自由度的限制。但有限元法也有其不足之处，其主要采用低阶插值函数，连续性差，应力计算精度低。但在工程设计中应力是对大坝安全性评价的重要依据，所以，有必要寻求更为有效的其它数值计算方法。根据我们研究，提出以下计算模型和方法供参考，欢迎批评指正。 　　1 多边形单元技术 　　多边形单元（PE-Polygon Element）不是一个新概念，但经过澳大利亚新南威尔士大学（UN-SW）宋崇民教授团队的发展，建立在SBFEM基础上的多边形单元已经成为与FEM一样具有广泛适应性的网格离散技术，对数值计算具有强大的功能。PE可称为广义的FE。 　　大连理工大学发展了PE技术在结构静动力分析中的应用，主要工作为：（1）单元细化与连续性升阶可以只在PE的边界上进行，不改变网格的拓朴结构，使应力计算精度显著提高，但计算工作量则很节约；（2）将PE与曲线边界相结合，既可构造纯直线边的PE，也可构造曲线边的PE，或直线边与曲