混凝土拱坝地震动应力的精细化分析与抗震策略研究

混凝土拱坝地震动应力的精细化分析与抗震策略研究一、引言1.1研究背景与意义我国地处环太平洋地震带与欧亚地震带之间，是世界上地震灾害最为严重的国家之一。在过去的几十年里，我国发生了多次强烈地震，如1976年的唐山大地震、2008年的汶川大地震以及2010年的玉树地震等。这些地震不仅造成了大量的人员伤亡和财产损失，也对各类基础设施，尤其是水利工程设施，构成了严重威胁。随着我国经济的快速发展和能源需求的不断增长，水利水电工程建设在我国能源战略中占据着越来越重要的地位。高拱坝作为一种高效的挡水建筑物，因其具有良好的经济性和稳定性，在我国西南、西北等地区得到了广泛应用。这些地区往往处于强震频发的高地震烈度区，如金沙江、雅砻江、大渡河等流域，在建和拟建的高拱坝众多，像小湾拱坝（坝高294.5米）、溪洛渡拱坝（坝高285.5米）、白鹤滩拱坝（坝高289米）等。由于高拱坝一旦在地震中失事，将会引发下游地区的洪水泛滥，对人民生命财产安全造成毁灭性打击，因此，其抗震安全问题至关重要，已成为大坝设计中的关键控制因素。混凝土拱坝在地震作用下的响应极其复杂，地震动应力是衡量拱坝抗震性能的关键指标。当拱坝遭受地震时，坝体将承受随时间快速变化的地震惯性力、动水压力和动土压力等荷载，这些荷载会在坝体内产生复杂的应力分布，包括拉应力、压应力和剪应力。过高的拉应力可能导致坝体混凝土开裂，削弱坝体的强度和整体性；过大的压应力则可能引发混凝土的局部破坏和失稳；剪应力过大还可能导致坝体出现剪切破坏。准确分析混凝土拱坝在地震作用下的动应力分布规律和变化特征，对于评估拱坝的抗震安全性、优化拱坝的抗震设计以及采取有效的抗震加固措施具有重要意义。从工程应用角度来看，地震动应力分析为拱坝的抗震设计提供了关键的技术支撑。通过对地震动应力的准确计算，可以合理确定拱坝的体型、尺寸和材料强度，优化坝体的结构布置，从而提高拱坝的抗震能力。例如，在设计过程中，根据地震动应力分析结果，对坝体应力集中区域进行加强处理，如增加混凝土厚度、配置钢筋或采用预应力技术等，以增强坝体的抗裂性能和承载能力。此外，地震动应力分析结果还可以为拱坝的运行管理和维护提供科学依据，通过实时监测坝体的应力状态，及时发现潜在的安全隐患，采取相应的措施进行处理，确保拱坝的长期安全运行。从学术研究角度来讲，混凝土拱坝地震动应力分析是水工抗震学科的重要研究内容之一，其涉及到固体力学、结构动力学、地震工程学等多个学科领域，研究过程中需要综合考虑坝体-库水-地基的动力相互作用、材料非线性、几何非线性以及地震波的传播特性等诸多复杂因素。对这些复杂问题的深入研究，不仅有助于推动水工抗震学科的发展，完善拱坝抗震理论体系，也能够为其他大型工程结构的抗震分析提供有益的借鉴和参考。在我国大力推进清洁能源开发和基础设施建设的背景下，开展混凝土拱坝地震动应力分析与研究具有重要的现实意义和理论价值，它是保障高拱坝安全运行、促进水利水电事业可持续发展的关键环节。1.2国内外研究现状1.2.1混凝土拱坝地震动应力分析方法研究早期，由于计算能力的限制，混凝土拱坝地震动应力分析主要采用简化方法。如纯拱法，将拱坝视为一系列独立的水平拱圈，忽略了拱圈之间的相互作用以及梁向的作用，虽计算简便，但无法准确反映拱坝的实际受力状态，仅适用于初步估算和小型拱坝分析。随着对拱坝受力机理认识的深入，拱梁分载法应运而生，该方法将拱坝的荷载通过拱和梁两种作用进行分配，考虑了拱坝的空间受力特性，比纯拱法更为合理，在很长一段时间内成为拱坝应力分析的主要方法，如在早期的一些中低高度拱坝设计中得到广泛应用。随着计算机技术的飞速发展，有限元法逐渐成为混凝土拱坝地震动应力分析的主流方法。有限元法能够将复杂的拱坝结构离散为众多小单元，通过求解单元的平衡方程，进而得到整个结构的应力分布。它可以方便地考虑坝体-库水-地基的动力相互作用、材料非线性和几何非线性等复杂因素。例如，在对小湾拱坝的地震动应力分析中，利用有限元软件建立了详细的坝体-库水-地基耦合模型，考虑了库水的可压缩性和地基的辐射阻尼，准确地模拟了地震作用下拱坝的应力响应。除了常规有限元法，一些改进的有限元方法也不断涌现，如无单元法、流形元法等，这些方法在处理复杂边界条件和大变形问题时具有独特优势，为拱坝地震动应力分析提供了新的思路和手段。在时域分析方面，逐步积分法是常用的方法之一，如Newmark法、Wilson-θ法等，它们通过将地震过程划分为一系列微小的时间步，逐步求解结构在每个时间步的动力响应，能够准确地反映结构在地震过程中的动态变化。在频域分析中，振型分解反应谱法应用广泛，该方法通过将结构的地震响应分解为各个振型的贡献，利用反应谱来确定每个振型的最大响应，进而得到结构的总响应。例如，在对溪洛渡拱坝的抗震分析中，采用振型分解反应谱法计算了坝体在不同地震工况下的动应力，为拱坝的抗震设计提供了重要依据。1.2.2混凝土拱坝地震动应力影响因素研究坝体-库水-地基的动力相互作用是影响混凝土拱坝地震动应力的重要因素之一。库水的存在会增加坝体的附加质量，改变坝体的自振特性，同时库水与坝体之间的相互作用还会产生动水压力，对坝体的应力分布产生影响。地基的刚度和阻尼特性也会对拱坝的地震动响应产生显著影响，地基的柔性会使坝体的地震反应增大，而地基的辐射阻尼则会消耗地震能量，降低坝体的地震响应。例如，研究表明，考虑库水可压缩性和地基辐射阻尼后，拱坝的动拉应力会显著降低。材料非线性对混凝土拱坝地震动应力的影响也不容忽视。混凝土在地震作用下会表现出非线性特性，如弹性模量的降低、开裂和损伤等，这些非线性行为会导致坝体的刚度退化和应力重分布。例如，采用混凝土塑性损伤模型研究发现，在强震作用下，坝体混凝土会出现开裂和损伤，导致坝体的刚度降低，动应力重新分布，损伤区域的应力集中现象明显。几何非线性也是影响拱坝地震动应力的重要因素，尤其是对于高拱坝，在地震作用下坝体的大变形可能会导致几何形状的改变，从而引起应力的变化。如在对超高拱坝的研究中发现，考虑几何非线性后，坝体的应力分布更加复杂，局部区域的应力集中现象加剧。1.2.3混凝土拱坝抗震措施研究为了提高混凝土拱坝的抗震性能，国内外学者提出了多种抗震措施。在坝体结构设计方面，优化拱坝的体型是一种有效的方法，通过合理设计拱坝的曲率、厚度和坝肩形状等参数，可以改善坝体的受力状态，降低地震动应力。例如，采用优化算法对拱坝体型进行优化设计，使坝体在地震作用下的应力分布更加均匀，减小了应力集中区域。在坝体配筋方面，通过在坝体内部配置钢筋，可以提高坝体的抗拉强度和抗裂性能，减少地震作用下坝体的开裂和损伤。如印度Koyna坝在强震开裂后，采用垂直预应力锚索加固，有效地限制了裂缝的扩展；我国新丰江坝在地震后，对坝体裂缝进行插筋加固并加大坝体剖面，提高了坝体的抗震能力。设置抗震构造措施也是提高拱坝抗震性能的重要手段，如在拱坝横缝处设置插筋、阻尼器等，能够限制横缝的张开度，增强坝体的整体性。例如，前苏联英古里拱坝采用横缝插筋和坝面钢筋网，智利拱坝采用水平连续式钢筋加固带，都取得了较好的抗震效果；国内学者研究了横缝配阻尼器对横缝张开度、坝体加速度及应力与位移的影响，发现横缝配筋和阻尼器对横缝张开度限制能够起到一定的作用。1.2.4研究不足与展望尽管国内外在混凝土拱坝地震动应力分析方面取得了丰硕的成果，但仍存在一些不足之处。在分析方法方面，虽然有限元法能够考虑多种复杂因素，但对于一些特殊的工程问题，如坝体-库水-地基系统的流固耦合问题、混凝土材料的多尺度非线性问题等，现有的分析方法还存在一定的局限性，需要进一步改进和完善。在影响因素研究方面，虽然对坝体-库水-地基动力相互作用、材料非线性和几何非线性等因素进行了大量研究，但这些因素之间的耦合作用机制尚未完全明确，需要深入研究它们之间的相互关系，以更准确地评估其对拱坝地震动应力的影响。在抗震措施研究方面，目前的抗震措施主要是基于经验和工程类比提出的，缺乏系统的理论支撑和优化设计方法。对于一些新型抗震措施的效果评估和作用机理研究还不够深入，需要进一步开展相关研究，以开发出更加有效的抗震措施。未来的研究可以朝着以下几个方向展开：一是进一步发展和完善混凝土拱坝地震动应力分析方法，结合多物理场耦合理论、人工智能技术等，建立更加精确、高效的分析模型；二是深入研究各种影响因素之间的耦合作用机制，揭示拱坝在地震作用下的破坏机理，为抗震设计提供更坚实的理论基础；三是开展抗震措施的优化设计研究，通过数值模拟和试验研究相结合的方法，对现有抗震措施进行优化，并探索新型抗震措施，提高拱坝的抗震性能；四是加强对实际工程的监测和反馈分析，通过对实际拱坝在地震作用下的响应监测，验证和改进理论分析结果，为拱坝的抗震安全提供更可靠的保障。1.3研究内容与方法1.3.1研究内容混凝土拱坝地震动应力分析方法研究：对现有混凝土拱坝地震动应力分析方法进行系统梳理和对比，包括振型分解反应谱法、时程分析法等。深入研究有限元法在拱坝地震动应力分析中的应用，建立考虑坝体-库水-地基动力相互作用、材料非线性和几何非线性的精细化有限元模型。结合实际工程案例，利用有限元软件对不同分析方法的计算结果进行验证和分析，评估各方法的优缺点和适用范围，为后续研究选择合适的分析方法提供依据。混凝土拱坝地震动应力影响因素分析：全面研究坝体-库水-地基动力相互作用对拱坝地震动应力的影响机制。通过理论推导和数值模拟，分析库水可压缩性、库水附加质量、地基刚度和阻尼等因素对坝体自振特性和地震动应力分布的影响规律。深入探讨材料非线性和几何非线性对拱坝地震动应力的影响。采用混凝土塑性损伤模型，研究混凝土在地震作用下的开裂、损伤和刚度退化等非线性行为对坝体应力分布和重分布的影响；考虑拱坝在地震作用下的大变形效应，分析几何非线性对坝体应力状态的改变。通过参数化分析，确定各影响因素的敏感程度，明确在拱坝抗震设计中需要重点关注的因素。混凝土拱坝抗震策略研究：基于地震动应力分析结果，研究优化拱坝体型对降低地震动应力的作用。采用优化算法，对拱坝的曲率、厚度、坝肩形状等关键参数进行优化设计，通过数值模拟分析优化后拱坝在地震作用下的应力分布情况，验证优化效果，总结优化规律。研究坝体配筋和设置抗震构造措施等抗震策略对提高拱坝抗震性能的作用。分析不同配筋方案对坝体抗拉强度和抗裂性能的提升效果，以及设置插筋、阻尼器等抗震构造措施对限制横缝张开度、增强坝体整体性和减小地震动应力的影响，通过数值模拟和对比分析，提出合理的配筋和抗震构造措施方案。案例分析与工程应用：选取实际的混凝土拱坝工程作为案例，运用前面研究的地震动应力分析方法和抗震策略，对该拱坝进行详细的地震动应力分析和抗震性能评估。结合工程实际情况，如坝址地质条件、地震设防标准等，对分析结果进行深入讨论，评估拱坝在设计地震作用下的抗震安全性。根据分析评估结果，针对拱坝存在的抗震薄弱环节，提出具体的抗震加固建议和措施，为实际工程的抗震设计和加固提供参考依据，验证研究成果的工程实用性和有效性。1.3.2研究方法数值模拟方法：利用大型通用有限元软件，如ANSYS、ABAQUS等，建立混凝土拱坝的三维有限元模型。在模型中，合理模拟坝体、库水和地基的材料特性和力学行为，考虑坝体-库水-地基的动力相互作用、材料非线性和几何非线性等复杂因素。通过施加不同的地震波输入，进行地震动应力分析和抗震性能模拟，得到坝体在地震作用下的应力、位移、加速度等响应结果。利用数值模拟方法，可以灵活地改变模型参数，进行参数化分析，深入研究各种因素对拱坝地震动应力的影响规律。理论分析方法：运用结构动力学、固体力学、地震工程学等相关理论，对混凝土拱坝在地震作用下的动力响应进行理论推导和分析。建立拱坝的简化力学模型，推导其动力平衡方程，求解结构的自振频率、振型和地震响应。通过理论分析，深入理解拱坝地震动应力的产生机制和分布规律，为数值模拟提供理论基础和验证依据。同时，理论分析方法还可以用于对一些复杂问题进行定性分析，指导数值模拟的方案设计和结果解释。案例研究方法：收集国内外实际的混凝土拱坝工程案例，对其设计资料、地震监测数据、运行情况等进行详细调研和分析。通过对实际工程案例的研究，了解不同拱坝在地震作用下的实际响应情况，验证理论分析和数值模拟结果的准确性和可靠性。同时，从实际工程案例中总结经验教训，发现存在的问题和不足，为进一步改进和完善混凝土拱坝的抗震设计和分析方法提供参考。二、混凝土拱坝地震动应力分析方法2.1拟静力法2.1.1原理与计算步骤拟静力法是一种用静力学方法近似解决动力学问题的简易方法，其基本思想是在静力计算的基础上，将地震作用简化为一个惯性力系附加在研究对象上。在混凝土拱坝地震动应力分析中，该方法的核心在于设计地震加速度的确定。从原理角度来看，拟静力法基于达朗贝尔原理，将结构在地震作用下所受的惯性力视为等效的静力荷载。假设拱坝的质量为m，在地震作用下产生的加速度为a，则等效的惯性力F=ma。通过将此惯性力施加在拱坝结构上，按照静力分析的方法来求解坝体的应力和变形。在计算步骤方面，首先需要确定设计地震加速度。这通常根据工程场地的地震危险性分析结果，结合相关的抗震设计规范来取值。例如，根据我国《水工建筑物抗震设计规范》，设计地震加速度一般根据场地的地震基本烈度、设计地震分组以及场地类别等因素来确定。在确定设计地震加速度后，计算拱坝各部分的惯性力。对于拱坝这样的复杂结构，可将其离散为若干个质点，每个质点的惯性力等于该质点的质量与设计地震加速度的乘积。然后，将这些惯性力作为静力荷载施加在拱坝的相应位置上，采用常规的结构力学方法，如材料力学公式、有限元法（在简化有限元模型情况下也可用于拟静力分析）等，计算拱坝在这些荷载作用下的应力分布。以材料力学方法为例，对于拱坝的拱圈，可将其视为受弯构件，根据弯矩和轴力计算公式，结合惯性力产生的荷载效应，计算拱圈截面的正应力和剪应力；对于拱坝的梁向结构，也采用类似的方法计算其应力。2.1.2应用范围与局限性拟静力法适用于结构简单、地震响应较规则的拱坝。当拱坝的体型较为规则，坝体-库水-地基系统的动力相互作用相对较弱，且地震作用强度较低时，拟静力法能够给出较为合理的结果。例如，在一些中低高度的小型拱坝设计中，由于其结构相对简单，地震响应相对规则，拟静力法能够满足工程设计的基本要求。然而，拟静力法存在明显的局限性。在考虑动力特性方面，它只能在有限程度上反映荷载的动力特性，不能反映各种材料自身的动力特性以及结构物之间的动力响应。例如，混凝土材料在地震作用下的动态力学性能与静态力学性能存在差异，拟静力法无法准确考虑这种差异。对于结构物之间的动力耦合关系，如坝体与库水之间的动水压力相互作用、坝体与地基之间的动力传递等，拟静力法也不能准确反映。这是因为拟静力法将地震作用简化为等效静力荷载，忽略了地震动的时间历程和结构在地震过程中的动态响应。当拱坝结构复杂，如具有不规则的体型、复杂的坝肩地形或地基条件时，拟静力法的计算结果往往与实际情况偏差较大。在强震作用下，结构的非线性行为（如混凝土的开裂、塑性变形等）显著，拟静力法无法考虑这些非线性因素对结构响应的影响。例如，在强震作用下，拱坝坝体混凝土可能出现开裂，导致结构刚度降低，应力重分布，而拟静力法基于弹性假设，无法准确模拟这种情况。因此，拟静力法不能用于地震时土体刚度有明显降低或者产生液化的场合，也只适用于设计加速度较小、动力相互作用不甚突出的结构抗震设计。在实际工程中，对于高拱坝以及处于强震区的拱坝，拟静力法通常作为一种初步的分析方法，为后续更精确的分析提供参考，而不能作为主要的设计依据。2.2反应谱法2.2.1振型分解反应谱法原理振型分解反应谱法基于结构动力学理论，其核心是将复杂的多自由度结构体系在地震作用下的响应，分解为多个独立的单自由度体系的响应。在地震作用下，多自由度结构体系的运动方程为：[M]\{\ddot{x}\}+[C]\{\dot{x}\}+[K]\{x\}=-[M]\{1\}\ddot{x}_{g}其中，[M]为质量矩阵，[C]为阻尼矩阵，[K]为刚度矩阵，\{\ddot{x}\}、\{\dot{x}\}、\{x\}分别为加速度、速度和位移向量，\ddot{x}_{g}为地面加速度，\{1\}为元素全为1的向量。利用振型正交性原理，引入主振型向量\{\Phi\}，将位移向量\{x\}展开为振型的线性组合，即\{x\}=\sum_{j=1}^{n}\{\Phi_{j}\}q_{j}(t)，其中q_{j}(t)为第j阶振型的广义坐标。将其代入运动方程，并利用振型正交性条件\{\Phi_{i}\}^{T}[M]\{\Phi_{j}\}=0(i\neqj)和\{\Phi_{i}\}^{T}[K]\{\Phi_{j}\}=0(i\neqj)，可将多自由度体系的运动方程解耦为n个独立的单自由度体系的运动方程：\ddot{q}_{j}+2\xi_{j}\omega_{j}\dot{q}_{j}+\omega_{j}^{2}q_{j}=-\gamma_{j}\ddot{x}_{g}其中，\omega_{j}为第j阶自振频率，\xi_{j}为第j阶阻尼比，\gamma_{j}为第j阶振型参与系数。对于单自由度体系，在给定的地震作用下，其最大反应（加速度、速度或位移）与体系的自振周期T（T=2\pi/\omega）之间存在一定的关系，这种关系通过反应谱来描述。反应谱是根据大量的地震记录，对不同自振周期的单自由度体系进行动力分析，得到其最大反应后绘制而成的曲线。例如，加速度反应谱S_{a}表示单自由度体系在地震作用下的最大加速度反应与自振周期的关系。根据反应谱理论，可得到第j阶振型的最大地震作用F_{j}为：F_{j}=\alpha_{j}\gamma_{j}G_{j}其中，\alpha_{j}为第j阶振型的地震影响系数，可由反应谱查得；G_{j}为第j阶振型的等效重力荷载。最后，通过一定的组合规则（如平方和开方SRSS法），将各阶振型的地震作用效应组合起来，得到结构的总地震作用效应。即结构的总地震作用效应S为：S=\sqrt{\sum_{j=1}^{n}S_{j}^{2}}其中，S_{j}为第j阶振型的地震作用效应。通过这种方式，振型分解反应谱法将地震荷载分解到各个分量并施加到结构上，再通过叠加得到最终的地震响应值。2.2.2计算流程与关键参数确定反应谱法的计算流程如下：结构离散与模型建立：利用有限元软件或其他数值方法，将混凝土拱坝离散为多个单元，建立结构的有限元模型。在模型中，准确定义坝体的几何形状、材料属性、边界条件等参数。例如，对于坝体材料，根据混凝土的类型确定其弹性模量、泊松比等参数；对于边界条件，考虑坝体与地基的连接方式，通常采用固定约束或弹簧约束来模拟地基对坝体的作用。自振特性计算：求解结构的特征方程[K-\omega^{2}M]\{\Phi\}=0，得到结构的自振频率\omega_{j}和振型\{\Phi_{j}\}（j=1,2,\cdots,n）。自振频率反映了结构振动的快慢，振型则描述了结构在振动时的变形形态。在实际计算中，通常采用子空间迭代法、Lanczos法等数值方法来求解特征方程。例如，子空间迭代法通过不断迭代，逐步逼近结构的真实自振特性，具有较高的计算精度和收敛速度。振型参与系数计算：根据公式\gamma_{j}=\frac{\{\Phi_{j}\}^{T}[M]\{1\}}{\{\Phi_{j}\}^{T}[M]\{\Phi_{j}\}}计算各阶振型的参与系数\gamma_{j}。振型参与系数表示第j阶振型在地震响应中的贡献程度，其值越大，说明该振型对结构地震响应的影响越大。地震影响系数确定：根据场地类别、设计地震分组和结构的自振周期，从规范规定的反应谱曲线中查取各阶振型的地震影响系数\alpha_{j}。我国《水工建筑物抗震设计规范》给出了不同场地条件下的反应谱曲线，其形状与地震动的频谱特性、场地土的性质等因素有关。例如，对于坚硬场地土，反应谱曲线在短周期段幅值较高，随着周期的增大而迅速衰减；对于软弱场地土，反应谱曲线在长周期段幅值相对较大。地震作用计算：利用公式F_{ji}=\alpha_{j}\gamma_{j}X_{ji}G_{i}计算第j阶振型第i质点的水平地震作用F_{ji}，其中X_{ji}为第j阶振型第i质点的相对水平位移，G_{i}为第i质点的重力荷载代表值。通过对各质点的地震作用进行叠加，可得到各阶振型的地震作用效应。效应组合：采用合适的组合规则，如SRSS法，将各阶振型的地震作用效应组合起来，得到结构的总地震作用效应。对于一些特殊结构或地震作用较为复杂的情况，也可采用其他组合规则，如完全二次型组合CQC法，以更准确地考虑各振型之间的相关性。在上述计算流程中，自振频率、振型和振型参与系数等关键参数的准确确定至关重要。自振频率和振型决定了结构的振动特性，直接影响到地震作用的计算结果。振型参与系数则反映了各阶振型在地震响应中的相对重要性，对于合理选择参与组合的振型数量具有指导意义。在实际工程计算中，通常需要根据结构的复杂程度和精度要求，合理确定参与计算的振型数量。一般来说，对于大多数混凝土拱坝，前几阶振型对结构的地震响应起主要作用，选取前5-10阶振型进行计算即可满足工程精度要求。但对于一些体型复杂或高地震烈度区的拱坝，可能需要选取更多的振型参与计算。2.3时程分析法2.3.1基本原理与数值积分方法时程分析法是一种直接求解动力平衡方程以计算结构地震反应的方法，在数学上属于步步积分法，抗震设计中也被称为“动态设计”。该方法通过输入与结构所在场地相应的地震波作为地震作用，从初始状态开始，一步一步地逐步积分，直至地震作用终了，从而求得整个时间历程的地震反应。在混凝土拱坝地震动应力分析中，时程分析法的基本原理基于结构动力学的基本方程。对于一个多自由度的拱坝结构体系，其动力平衡方程可表示为：[M]\{\ddot{x}\}+[C]\{\dot{x}\}+[K]\{x\}=\{F(t)\}其中，[M]为质量矩阵，[C]为阻尼矩阵，[K]为刚度矩阵，\{\ddot{x}\}、\{\dot{x}\}、\{x\}分别为加速度、速度和位移向量，\{F(t)\}为随时间变化的地震作用荷载向量。由于地震作用是随时间快速变化的动力荷载，该方程是一个非线性的常微分方程，一般难以获得解析解，因此需要采用数值积分方法进行求解。常用的数值积分方法有Newmark法、Wilson-θ法、中心差分法等。以Newmark法为例，其基本思想是假设在时间间隔\Deltat内加速度呈线性变化，通过对动力平衡方程进行离散化处理，将其转化为一系列的代数方程进行求解。在t+\Deltat时刻，位移、速度和加速度的近似计算公式如下：\{x\}_{t+\Deltat}=\{x\}_{t}+\Deltat\{\dot{x}\}_{t}+(\frac{1}{2}-\beta)\Deltat^{2}\{\ddot{x}\}_{t}+\beta\Deltat^{2}\{\ddot{x}\}_{t+\Deltat}\{\dot{x}\}_{t+\Deltat}=\{\dot{x}\}_{t}+(1-\gamma)\Deltat\{\ddot{x}\}_{t}+\gamma\Deltat\{\ddot{x}\}_{t+\Deltat}其中，\beta和\gamma是与积分精度和稳定性相关的参数，通常取\beta=1/4，\gamma=1/2时，Newmark法具有无条件稳定性。通过迭代求解上述方程，可以逐步得到结构在各个时刻的位移、速度和加速度响应，进而计算出结构的内力和应力。Wilson-θ法与Newmark法类似，也是一种逐步积分法，它通过引入一个大于1的系数\theta（通常取\theta=1.4），对加速度进行线性外推，从而提高积分的稳定性。中心差分法是一种显式积分方法，它直接利用相邻时刻的位移差来近似计算速度和加速度，计算过程相对简单，但稳定性条件较为苛刻，只适用于一些小型问题或对计算精度要求不高的情况。2.3.2地震波选取与输入地震波的选取是时程分析法中的关键环节，其对结构的地震反应影响非常大。合理选取地震波需要综合考虑拱坝场地条件和地震特性等多方面因素。从场地条件来看，主要考虑场地的地质条件和地震动参数。场地的地质条件决定了地震波在传播过程中的衰减和频谱特性。例如，对于坚硬场地土，地震波的高频成分相对较多，而对于软弱场地土，地震波的低频成分会得到放大。因此，应根据场地的类别（如I类、II类、III类、IV类场地）选择与之相匹配的地震波。我国《建筑抗震设计规范》给出了不同场地类别的特征周期，在选取地震波时，应使所选地震波的频谱特性与场地的特征周期相接近。地震特性方面，包括地震的震级、震中距、地震波的类型等。震级和震中距决定了地震波的强度和频谱特性，一般来说，震级越大、震中距越小，地震波的能量越强，对结构的破坏作用也越大。地震波的类型主要有天然地震波和人工合成地震波。天然地震波是实际地震记录得到的，能够真实反映地震的特性，但由于其数量有限，且每个地震记录都具有特殊性，不一定能完全满足工程分析的要求。人工合成地震波则是根据地震动的统计特性和场地条件，利用数学模型合成的地震波，其优点是可以根据需要调整地震波的频谱特性、峰值加速度等参数，以满足不同工程的需求。在实际工程中，为了更全面地考虑地震的不确定性，通常选取多条地震波进行计算，然后取其平均值作为计算结果。例如，我国《水工建筑物抗震设计规范》规定，进行时程分析时，应至少选取3条实际强震记录和1条人工合成地震波。在选择实际强震记录时，应尽量选择与工程场地条件相近、震级和震中距与工程所在地可能遭遇的地震情况相似的地震波。国内外进行结构时程分析时所经常采用的几条实际强震记录主要有适用于I类场地的滦河波、适用于II、III类场地的El-Centrol波(1940，N-S)和Taft波（1952，E-w）、适用于IV类场地的宁河波等。地震波的输入方式也会影响拱坝的地震反应。常见的输入方式有一致激励输入和多点激励输入。一致激励输入假定地震波在整个地基表面同时到达且具有相同的幅值和相位，这种输入方式计算简单，在早期的工程分析中应用广泛。然而，实际地震中，由于地基的不均匀性和地震波的传播特性，地震波在地基表面各点的到达时间和幅值存在差异，因此一致激励输入与实际情况存在一定偏差。多点激励输入则考虑了地震波的空间变化特性，通过在地基表面不同位置输入不同的地震波时程，更真实地模拟了地震对拱坝的作用。对于大型拱坝，尤其是坝基范围较大、地基条件复杂的拱坝，多点激励输入能够更准确地反映坝体-库水-地基系统的动力相互作用，得到更合理的地震反应结果。但多点激励输入的计算过程较为复杂，需要更多的计算资源和更精确的地基参数。三、影响混凝土拱坝地震动应力的因素3.1坝体材料特性3.1.1混凝土动态力学性能混凝土作为混凝土拱坝的主要建筑材料，其动态力学性能对拱坝在地震作用下的应力响应有着显著影响。混凝土是一种典型的率敏感材料，在不同应变率下，其强度、刚度和延性等力学性能会发生明显变化。在强度方面，大量试验研究表明，随着应变率的增加，混凝土的强度呈现出增大的趋势。如邹笃建等人通过混凝土柱的轴心动态抗压试验，在10^{-4}/s～10^{-2}/s应变率范围内研究发现，当应变率为10^{-4}/s、10^{-3}/s、10^{-2}/s时，混凝土的抗压强度相对准静态抗压强度（应变率为10^{-6}/s）分别增加了7.45\%、19.51\%和29.23\%。这是因为在高应变率下，混凝土内部的微裂缝发展受到抑制，材料的变形来不及充分发展，从而使混凝土能够承受更大的荷载。然而，当应变率过高时，混凝土的强度反而会下降。这是由于应变速度过快，混凝土无法快速适应外载荷的变化，导致内部结构迅速破坏，强度降低。在地震作用下，应变率的变化范围较大，坝体不同部位的混凝土可能处于不同的应变率状态，这就使得坝体的强度分布变得复杂。对于可能承受高应变率作用的坝体关键部位，如坝肩、坝踵等，混凝土强度的提高有助于增强坝体的承载能力，但过高的应变率也可能导致这些部位的混凝土强度下降，增加坝体的安全隐患。混凝土的刚度在不同应变率下也有所变化。一般来说，随着应变率的增加，混凝土的弹性模量有一定程度的提高，但变化幅度相对较小。在一些研究中发现，在10^{-4}/s～10^{-2}/s应变率范围内，混凝土的弹性模量相对于准静态弹性模量变化不大。混凝土刚度的变化会影响拱坝的自振特性和地震响应。刚度的提高会使拱坝的自振频率增加，根据结构动力学原理，自振频率的改变会影响地震作用下结构的响应大小和分布。在地震作用下，若拱坝的自振频率与地震波的卓越频率接近，可能会发生共振现象，导致坝体的地震响应急剧增大。因此，混凝土刚度的变化通过影响拱坝的自振频率，间接影响着坝体在地震作用下的应力分布和大小。延性方面，混凝土的延性是指其在破坏前能够承受较大变形而不丧失承载能力的性能。随着应变率的增加，混凝土的延性通常会降低，表现出更明显的脆性。在高应变率下，混凝土内部的裂缝迅速发展，材料在较小的变形下就可能发生破坏。混凝土延性的降低对拱坝的抗震性能有不利影响。在地震作用下，坝体需要通过一定的变形来耗散地震能量，延性好的混凝土能够在变形过程中吸收更多的能量，从而保护坝体结构。而延性降低后，坝体在地震作用下更容易发生脆性破坏，导致结构的整体性丧失，增加坝体失事的风险。例如，在一些地震灾害中，由于混凝土延性不足，拱坝坝体出现突然的脆性断裂，造成了严重的后果。3.1.2材料非线性对动应力的影响混凝土材料在地震作用下会进入非线性阶段，其应力-应变关系呈现出复杂的非线性特征。当混凝土受到的荷载较小时，处于弹性阶段，应力与应变成正比，满足胡克定律。随着荷载的增加，当应力达到一定程度后，混凝土内部开始产生微裂缝，进入塑性阶段，应力-应变关系不再是线性的。在塑性阶段，混凝土的应力-应变曲线表现出明显的非线性。以混凝土棱柱体一次单调加载试验为例，当应力达到比例极限后，曲线开始偏离直线，混凝土呈现塑性性质，此时内部微裂缝逐渐发展。当应力继续增加达到峰值应力后，试件的承载力随应变增长逐渐减小，出现“应变软化”现象，裂缝迅速发展，混凝土内部结构的整体性受到越来越严重的破坏。不同强度等级的混凝土，其应力-应变曲线也存在差异。高强混凝土加载时的线性段范围增大，峰值应变略有增大，但过峰值后曲线骤然下跌，表现出更明显的脆性。材料非线性对拱坝地震动应力分布和大小有着重要影响。在地震作用下，坝体不同部位的混凝土可能处于不同的受力状态，部分区域会进入塑性阶段。当混凝土进入塑性阶段后，其刚度降低，导致坝体的整体刚度下降。根据结构力学原理，刚度的变化会引起应力的重分布。原本应力分布较为均匀的坝体，在混凝土进入塑性阶段后，会出现应力集中现象。在坝体的薄弱部位，如坝肩、坝踵以及孔洞周边等，由于应力集中，这些部位的混凝土更容易进入塑性阶段，形成塑性铰区域。塑性铰的出现改变了坝体的传力路径，使得坝体的受力状态更加复杂。混凝土的开裂和损伤也是材料非线性的重要表现。在地震作用下，当坝体混凝土的拉应力超过其抗拉强度时，混凝土会发生开裂。裂缝的出现不仅削弱了混凝土的承载能力，还会导致坝体的刚度进一步降低，从而引起应力的再次重分布。随着地震作用的持续，裂缝会不断扩展和贯通，使坝体的整体性受到严重破坏。混凝土的损伤还会导致其弹性模量降低，进一步影响坝体的力学性能。通过混凝土塑性损伤模型的模拟分析可以发现，在强震作用下，坝体混凝土会出现明显的开裂和损伤区域，这些区域的应力集中现象显著，动应力大幅增加。因此，在混凝土拱坝地震动应力分析中，必须充分考虑材料非线性的影响，以准确评估坝体在地震作用下的受力状态和抗震安全性。3.2坝体结构特征3.2.1拱坝体型参数拱坝的体型参数，如高宽比、厚高比等，对其自振特性和地震动应力分布有着重要影响。高宽比是指拱坝的坝高与坝顶宽度的比值，它反映了拱坝的纵向形态。一般来说，高宽比越大，拱坝的纵向刚度相对越小，自振频率越低。以某实际拱坝为例，当高宽比从3.0增加到3.5时，其自振频率下降了约10%。这是因为随着高宽比的增大，拱坝的坝体相对更加高耸，结构的柔性增加，在地震作用下更容易发生变形，从而导致自振频率降低。厚高比是拱坝坝体厚度与坝高的比值，它体现了拱坝的横向刚度。厚高比越大，拱坝的横向刚度越大。当厚高比增大时，拱坝的自振频率会相应提高。通过对不同厚高比的拱坝模型进行数值模拟分析发现，当厚高比从0.15增加到0.20时，拱坝的自振频率提高了约15%。这是因为坝体厚度的增加使得拱坝的抗变形能力增强，结构的整体刚度提高，从而自振频率增大。拱坝的曲率半径也是一个重要的体型参数。曲率半径的大小决定了拱坝拱圈的弯曲程度，进而影响拱坝的受力性能。较小的曲率半径会使拱坝的拱作用增强，在承受荷载时，拱圈能够更有效地将荷载传递到坝肩，减小坝体内部的应力。然而，曲率半径过小也会导致坝体应力集中现象加剧，尤其是在拱坝的拱端部位。在一些高拱坝设计中，为了优化坝体的受力状态，会采用变曲率的拱圈设计，即从坝顶到坝底，拱圈的曲率半径逐渐变化，以适应不同高程处的荷载分布特点。拱坝的坝肩形状对地震动应力分布也有显著影响。合理的坝肩形状能够改善坝肩的受力条件，减少坝肩处的应力集中。例如，采用扩宽坝肩的设计方式，可以增加坝肩的承载面积，降低坝肩处的应力水平。一些拱坝在坝肩部位采用了特殊的形状设计，如梯形坝肩、弧形坝肩等，通过数值模拟和工程实践验证，这些特殊形状的坝肩能够有效地改善坝体的应力分布，提高拱坝的抗震性能。3.2.2坝体结构布置坝体孔洞、廊道等结构布置会对地震动应力集中和传播产生重要影响。坝体孔洞，如泄洪孔、导流孔等，改变了坝体的连续性，容易在孔洞周边产生应力集中现象。当拱坝遭受地震时，地震波在传播过程中遇到孔洞，会发生反射和绕射，导致孔洞周边的应力分布变得复杂。通过有限元模拟分析发现，在孔洞边缘，动应力会显著增大，尤其是在孔洞的拐角处，应力集中系数可达到2-3倍。这是因为在这些部位，地震波的能量聚集，使得坝体材料承受的应力远高于平均应力水平。为了降低孔洞周边的应力集中，工程中常采用在孔洞周边设置加强环、加厚混凝土等措施。加强环可以增强孔洞周边的结构刚度，分散应力；加厚混凝土则可以增加材料的承载能力，减小应力集中的影响。坝体廊道，如灌浆廊道、排水廊道等，也会影响地震动应力的传播。廊道的存在会改变坝体内部的应力流路径，使得地震动应力在廊道附近发生重分布。在廊道与坝体的连接处，由于结构的不连续性，容易出现应力集中。研究表明，廊道的走向和位置对地震动应力的影响较大。当廊道与地震波传播方向平行时，应力集中现象相对较轻；而当廊道与地震波传播方向垂直时，应力集中现象较为明显。为了优化坝体结构布置，减小廊道对地震动应力的影响，可以合理调整廊道的位置和走向。在设计阶段，通过数值模拟分析不同廊道布置方案下的地震动应力分布情况，选择应力集中最小的方案。此外，在廊道与坝体的连接处采用渐变过渡的设计方式，也可以有效缓解应力集中现象。在坝体结构布置优化方面，除了上述针对孔洞和廊道的措施外，还可以综合考虑坝体的整体结构形式。采用合理的分缝方式，如诱导缝、收缩缝等，能够在一定程度上控制坝体的应力分布，减少地震作用下坝体的开裂风险。诱导缝是一种人为设置的薄弱部位，在地震作用下，诱导缝可以优先开裂，释放部分能量，从而保护坝体其他部位的完整性。收缩缝则主要用于适应混凝土在施工和运行过程中的收缩变形，避免因收缩应力导致坝体产生裂缝。在坝体的关键部位，如坝肩、坝踵等，合理配置钢筋或采用预应力技术，也可以提高坝体的抗震性能。钢筋能够增强坝体的抗拉能力，预应力技术则可以在坝体内部预先施加压应力，抵消部分地震作用产生的拉应力，从而提高坝体的抗裂性能和承载能力。3.3地基条件3.3.1地基土动力特性地基土的动力特性对混凝土拱坝在地震作用下的动应力有着显著影响。地基土的刚度是其重要的动力特性之一，它直接影响着拱坝的地震响应。地基土刚度的大小决定了其对拱坝的约束能力。当地基土刚度较大时，对拱坝的约束作用较强，能够限制拱坝的变形，使得拱坝在地震作用下的位移减小。在一些地基条件较好，地基土刚度较大的拱坝工程中，地震作用下坝体的位移相对较小。然而，这种较强的约束也会导致坝体与地基之间的相互作用力增大，从而使坝体的地震动应力增加。在有限元模拟分析中可以发现，当地基土刚度增大一倍时，坝体某些部位的动应力可能会增加20%-30%。这是因为地基土刚度的增加使得地震波在传播过程中更容易在坝体与地基的交界面处产生反射和折射，导致坝体所承受的地震荷载增大。相反，当地基土刚度较小时，对拱坝的约束作用较弱，拱坝在地震作用下的变形相对较大。此时，坝体与地基之间的相互作用力相对较小，坝体的地震动应力可能会有所降低。但是，地基土刚度过小也会带来问题，它可能导致拱坝的稳定性下降，增加坝体在地震作用下发生滑动、倾倒等破坏的风险。在一些地基土较为软弱的地区，如软黏土场地，由于地基土刚度小，拱坝在地震作用下的变形明显增大，虽然动应力有所降低，但坝体的整体稳定性受到威胁。地基土的阻尼特性也是影响拱坝地震动应力的重要因素。阻尼是指材料在振动过程中消耗能量的能力。地基土的阻尼可以分为材料阻尼和辐射阻尼。材料阻尼是由于地基土内部的摩擦和黏滞作用而消耗能量，辐射阻尼则是由于地震波向地基深部传播而导致的能量耗散。地基土的阻尼能够消耗地震能量，从而减小拱坝的地震响应。当地基土的阻尼较大时，地震波在地基中传播时的能量衰减较快，传递到坝体的地震能量相应减少，使得坝体的地震动应力降低。在实际工程中，一些地基土含有较多的黏性土或粉质土，这些土的阻尼较大，能够有效地降低拱坝的地震响应。通过数值模拟研究发现，当考虑地基土的阻尼后，拱坝的动应力峰值可能会降低10%-20%。这是因为阻尼的存在使得地震波的振幅减小，从而减小了坝体所承受的地震荷载。不同类型的地基土具有不同的动力特性，其对拱坝地震响应的影响也各不相同。岩石地基通常具有较高的刚度和较低的阻尼。由于岩石地基的刚度大，对拱坝的约束作用强，在地震作用下，拱坝的变形较小，但坝体的动应力相对较大。岩石地基的阻尼较小，地震能量在地基中的衰减较慢，更多的地震能量会传递到坝体，进一步增加了坝体的地震响应。在一些建在坚硬岩石地基上的拱坝，如花岗岩地基，地震动应力往往较大，对坝体的抗震性能提出了较高的要求。而土体地基的动力特性则较为复杂，其刚度和阻尼会受到土的类型、密度、含水量等多种因素的影响。砂土地基的刚度相对较小，但在密实状态下具有一定的抗剪强度。在地震作用下，砂土地基的变形较大，对拱坝的约束作用相对较弱，坝体的动应力可能相对较小。然而，砂土地基在饱和状态下，容易发生液化现象，导致地基的刚度和强度急剧降低，从而对拱坝的稳定性造成严重威胁。1964年日本新潟地震中，许多建筑物由于地基砂土液化而遭受严重破坏，这也为拱坝在砂土地基上的抗震设计敲响了警钟。黏性土地基的刚度和阻尼相对较大，其变形特性与砂土有很大不同。黏性土具有较强的黏聚力，在地震作用下，能够较好地保持地基的整体性，对拱坝的约束作用较为稳定。黏性土地基的阻尼较大，能够有效地消耗地震能量，减小拱坝的地震响应。在一些建在黏性土地基上的拱坝，如黏土含量较高的冲积层地基，地震动应力相对较小，坝体的抗震性能相对较好。3.3.2地基与坝体相互作用地基与坝体相互作用对地震波传播和坝体动应力分布有着重要影响。在地震作用下，地震波从地基传入坝体，由于地基和坝体的材料性质和力学特性不同，地震波在传播过程中会发生复杂的变化。地震波在地基与坝体的交界面处会发生反射和折射现象。当波从地基传入坝体时，由于坝体的刚度通常大于地基，波速会发生变化，导致波的传播方向改变。一部分地震波会被反射回地基，另一部分则会折射进入坝体。这种反射和折射现象会使地震波的能量分布发生改变，进而影响坝体的动应力分布。在坝体与地基的交界面附近，由于地震波的反射和折射，会出现应力集中现象。当反射波和折射波在交界面附近相互叠加时，会导致该区域的应力显著增大。通过数值模拟分析可以发现，在坝体与地基的交界面处，动应力可能会比坝体其他部位高出30%-50%。这种应力集中现象对坝体的抗震性能是不利的，容易导致坝体在该区域出现裂缝或破坏。地基与坝体相互作用还会改变坝体的自振特性。坝体与地基是一个相互耦合的系统，地基的存在会增加坝体的附加质量和约束条件，从而影响坝体的自振频率和振型。当考虑地基与坝体相互作用时，坝体的自振频率通常会降低。这是因为地基的附加质量增加了坝体的惯性，使得坝体振动时需要克服更大的阻力，从而导致自振频率下降。坝体的振型也会发生变化，原本较为规则的振型可能会变得更加复杂。坝体自振特性的改变会进一步影响其在地震作用下的响应。由于自振频率的降低，坝体在地震作用下更容易与地震波的卓越频率产生共振，从而增大坝体的地震响应。在地震动应力分析中，如果不考虑地基与坝体相互作用对自振特性的影响，计算结果可能会与实际情况存在较大偏差。为了准确分析地基与坝体相互作用对拱坝地震动应力的影响，需要采用合适的分析方法。有限元法是目前应用最为广泛的一种方法，它能够将地基和坝体离散为有限个单元，通过求解单元的平衡方程，得到整个系统的应力和位移分布。在有限元模型中，可以通过设置合适的边界条件和材料参数，来模拟地基与坝体的相互作用。对于地基与坝体的交界面，可以采用接触单元来模拟其接触状态，考虑界面的摩擦、滑移等非线性行为。边界元法也是一种常用的分析方法，它将问题的求解域边界离散化，通过求解边界积分方程来得到边界上的未知量，进而计算域内的物理量。边界元法在处理无限域问题时具有独特的优势，能够准确地模拟地基的无限远边界条件。在分析地基与坝体相互作用时，边界元法可以将地基视为无限域，通过边界元与有限元的耦合，实现对整个系统的分析。除了有限元法和边界元法，还有一些其他的方法，如子结构法、比例边界有限元法等。子结构法将整个结构划分为若干个子结构，分别对每个子结构进行分析，然后通过界面协调条件将子结构的结果进行组合。这种方法可以有效地减少计算量，提高计算效率。比例边界有限元法是一种半解析方法，它在边界上离散，通过引入比例坐标，将二维或三维问题转化为一维问题进行求解。该方法在处理地基与坝体相互作用问题时，能够准确地模拟地基的动力特性和辐射阻尼，具有较高的计算精度。在实际工程应用中，需要根据具体问题的特点和要求，选择合适的分析方法，以准确评估地基与坝体相互作用对拱坝地震动应力的影响。3.4地震动特性3.4.1地震波频谱特性地震波的频谱特性是影响混凝土拱坝地震响应的重要因素之一，其卓越周期、频谱组成等特性与拱坝自振频率密切相关，并对拱坝地震动应力产生显著影响。地震波的卓越周期是指在地震记录中，能量相对集中的某一周期。当拱坝的自振频率与地震波的卓越频率接近时，容易发生共振现象。共振会导致拱坝的地震响应急剧增大，坝体的地震动应力显著提高。以某实际拱坝为例，通过数值模拟分析发现，当拱坝的某一阶自振频率与输入地震波的卓越频率接近时，坝体在该阶振型下的地震动应力峰值相比非共振情况下提高了2-3倍。这是因为共振时，地震波的能量能够持续地输入到拱坝结构中，使得结构的振动不断加剧，从而导致应力大幅增加。地震波的频谱组成包含了不同频率成分的谐波。不同频率成分的地震波对拱坝的作用效果不同。高频成分的地震波能量相对集中在坝体的上部和表面，容易引起坝体上部的局部振动和应力集中。在地震作用下，坝体顶部的一些附属结构，如坝顶的启闭机室等，由于其自身的刚度和质量较小，对高频地震波较为敏感，容易在高频地震波的作用下产生较大的振动和应力。低频成分的地震波则更容易使坝体产生整体振动，对坝体的整体受力状态产生影响。低频地震波能够使坝体的主体结构产生较大的变形，从而导致坝体内部的应力分布发生改变。通过对不同频谱组成的地震波作用下拱坝地震响应的数值模拟研究发现，当输入的地震波中低频成分占比较大时，坝体的整体位移和应力分布相对较为均匀；而当高频成分占比较大时，坝体上部和表面的局部应力集中现象明显。为了避免共振现象的发生，在拱坝的设计阶段，需要合理调整拱坝的体型参数，以改变拱坝的自振频率，使其与可能遭遇的地震波卓越频率错开。通过优化拱坝的高宽比、厚高比等体型参数，可以有效地改变拱坝的自振频率。当高宽比减小、厚高比增大时，拱坝的自振频率会相应提高。在选择地震波进行地震动应力分析时，应充分考虑地震波的频谱特性，选择与工程场地条件和拱坝自振特性相匹配的地震波。对于场地条件较为复杂的拱坝，可能需要选取多条不同频谱特性的地震波进行分析，以更全面地评估拱坝在不同地震作用下的抗震性能。3.4.2地震波幅值与持时地震波幅值和持时对拱坝地震动应力大小和累积损伤有着重要影响，在抗震设计中需要采取相应的对策来保障拱坝的安全。地震波幅值直接决定了地震作用的强度，对拱坝地震动应力大小有着显著影响。一般来说，地震波幅值越大，拱坝所承受的地震惯性力、动水压力等荷载就越大，坝体的地震动应力也就越大。当地震波幅值增加一倍时，通过有限元模拟计算，拱坝坝体关键部位的动应力可能会增加1.5-2倍。在强震作用下，由于地震波幅值较大，坝体可能会出现超过材料强度极限的应力，导致混凝土开裂、结构局部破坏等情况。在一些地震灾害中，高幅值的地震波使得拱坝坝体出现严重的裂缝，甚至部分坝段发生坍塌。地震波持时是指地震波从开始到结束的持续时间。较长的地震波持时会使拱坝在较长时间内受到地震作用，从而导致坝体的累积损伤增加。随着地震波持时的延长，坝体混凝土在反复的地震作用下，内部微裂缝不断发展、扩展和贯通，材料的性能逐渐劣化。通过混凝土损伤力学模型的模拟分析发现，当地震波持时增加50%时，坝体混凝土的损伤程度会增加20%-30%。这种累积损伤会降低坝体的刚度和承载能力，进一步影响坝体的地震动应力分布。由于坝体刚度的降低，在后续的地震作用下，坝体的变形会增大，应力集中现象会更加明显。为了应对地震波幅值和持时对拱坝的影响，在抗震设计中可以采取一系列对策。在坝体材料选择方面，应选用高强度、高韧性的混凝土材料。高强度混凝土能够承受更大的应力，减少在高幅值地震波作用下的开裂风险；高韧性混凝土则具有较好的耗能能力，能够在地震作用下吸收更多的能量，降低累积损伤。在坝体结构设计上，优化拱坝的体型，合理调整拱坝的曲率、厚度等参数，使坝体的应力分布更加均匀，降低应力集中程度。采用合理的分缝和配筋措施，也可以提高坝体的抗裂性能和整体性。分缝可以释放坝体在温度变化和地震作用下产生的应力，配筋则能够增强坝体的抗拉能力。设置抗震构造措施，如在坝体关键部位设置阻尼器等，能够有效地耗散地震能量，减小地震波幅值和持时对坝体的影响。阻尼器可以在地震作用下产生阻尼力，消耗地震能量，降低坝体的振动响应。四、混凝土拱坝地震动应力分析案例研究4.1工程概况本案例研究选取的是位于我国西南地区某峡谷河段的混凝土拱坝。该地区处于欧亚板块与印度洋板块碰撞挤压的边缘地带，地质构造复杂，地震活动频繁，属于高地震烈度区，地震设防烈度为Ⅷ度。坝址所在河谷呈V形，谷坡陡峭，两岸地形基本对称。坝址区出露的基岩主要为花岗岩，岩石完整性较好，强度较高，但存在一些节理裂隙和小断层，对坝基和坝肩的稳定性有一定影响。该拱坝为双曲拱坝，坝高150m，坝顶高程为1250m，坝底高程为1100m。坝顶宽度为10m，坝底最大厚度为30m，厚高比为0.2。坝顶弧长为350m，拱冠梁上游面最大倒悬度为0.2。坝体混凝土采用C25混凝土，弹性模量为30GPa，泊松比为0.167。该拱坝是一座以发电为主，兼顾防洪、灌溉等综合利用的水利枢纽工程。水库正常蓄水位为1240m，相应库容为10.5亿m³；死水位为1200m，相应库容为4.5亿m³。电站装机容量为1000MW，年发电量为35亿kW・h。坝体设有5个表孔和4个中孔用于泄洪，表孔尺寸为12m×10m（宽×高），中孔尺寸为8m×8m（宽×高）。坝后式厂房布置在下游坝脚处，安装有4台单机容量为250MW的水轮发电机组。在工程建设过程中，对坝址区进行了详细的地质勘察工作，包括地质测绘、钻探、物探等。通过地质勘察，查明了坝址区的地层岩性、地质构造、水文地质条件等。在坝基处理方面，对坝基进行了固结灌浆和帷幕灌浆，以提高坝基的强度和防渗性能。在坝肩处理方面，对坝肩岩体进行了锚固和混凝土置换等措施，以增强坝肩的稳定性。4.2有限元模型建立4.2.1模型简化与离散化在建立混凝土拱坝有限元模型时，需要根据拱坝的实际结构和受力特点进行合理的模型简化与离散化。考虑到实际工程中，坝体与地基的连接部分在地震作用下的应力分布较为复杂，但在整体分析中，坝体的主要受力区域和变形特征对拱坝的抗震性能起关键作用。因此，对于坝体与地基连接部分的一些细微构造，如局部的小凸起、小凹槽等，在不影响整体力学性能的前提下，可以进行适当简化。坝体内部的一些小型孔洞，如用于观测的小孔洞，其对整体结构的力学性能影响较小，也可忽略不计。对于坝体和地基的离散化，采用有限元方法将其划分为有限个单元。在坝体部分，选用八节点六面体单元进行离散，这种单元具有较好的计算精度和稳定性，能够准确模拟坝体的复杂形状和受力状态。对于地基部分，由于其范围较大且远离坝体的区域对坝体地震响应的影响相对较小，采用四面体单元进行离散。四面体单元划分较为灵活，适用于复杂的几何形状，虽然其计算精度相对六面体单元略低，但在地基离散中能够在保证一定计算精度的前提下，有效减少计算量。在网格划分时，遵循一定的原则，在坝体应力集中区域，如坝肩、坝踵以及孔洞周边等部位，采用较小的单元尺寸进行加密，以提高计算精度，准确捕捉这些部位的应力变化。在坝体其他部位，根据结构的复杂程度和受力情况，合理调整单元尺寸，以平衡计算精度和计算效率。通过这种方式，既保证了模型能够准确反映坝体和地基的力学特性，又能在可接受的计算资源范围内完成分析。4.2.2材料参数与边界条件设定混凝土和地基等材料的力学参数取值直接影响有限元模型的计算结果，因此需要依据相关规范和试验数据进行准确确定。坝体混凝土采用C25混凝土，根据《混凝土结构设计规范》，其弹性模量取值为30GPa，泊松比为0.167。这些参数是通过大量的混凝土材料试验得到的，能够较好地反映C25混凝土的力学性能。对于地基材料，根据坝址区的地质勘察报告，地基主要为花岗岩，其弹性模量根据岩石的完整性和风化程度取值为60-80GPa，泊松比取0.25。不同风化程度的花岗岩力学性能存在差异，在取值时充分考虑了这些因素，以确保地基材料参数的准确性。模型边界条件的设定对计算结果也有重要影响。在地基底部，采用固定约束，限制其在三个方向的平动和转动自由度，模拟地基底部与基岩的紧密连接，使地基底部在地震作用下不会发生位移和转动。在地基侧面，为了模拟地基的无限域效应，采用粘性边界条件。粘性边界通过在边界节点上施加阻尼力，吸收地震波传播到边界时的能量，从而模拟地基的辐射阻尼，减少边界反射对计算结果的影响。在坝体与库水的交界面上，采用流固耦合边界条件，考虑库水与坝体之间的相互作用，包括动水压力的传递和库水对坝体的附加质量效应。通过准确设定材料参数和边界条件，能够建立起符合实际工程情况的有限元模型，为后续的地震动应力分析提供可靠的基础。4.3地震动输入与计算工况4.3.1地震波选择与处理根据工程场地地震危险性分析结果，该地区可能遭遇的地震动参数表明，地震波的频谱特性和幅值对拱坝的地震响应影响显著。为了准确分析拱坝在地震作用下的动应力，需要选择合适的地震波进行输入。在地震波选择过程中，考虑到该地区的地质条件和地震特性，从国内外地震数据库中选取了三条具有代表性的天然地震记录和一条人工合成地震波。三条天然地震记录分别为1940年El-Centro波（NS向）、1952年Taft波（EW向）以及1999年台湾集集地震波（TCU082站NS向）。这些地震波在不同频率范围内具有丰富的能量分布，能够较好地模拟该地区可能遭遇的地震情况。El-Centro波在高频段能量相对集中，Taft波在中频段能量较为突出，台湾集集地震波则具有较宽的频谱特性。人工合成地震波是根据该地区的地震危险性分析结果和场地特征，利用随机振动理论和反应谱拟合技术生成的。它具有与该地区设计反应谱相匹配的频谱特性和峰值加速度，能够更准确地反映该地区的地震动特性。对选取的地震波进行了一系列处理，以满足分析需求。采用Butterworth滤波器对地震波进行滤波处理，去除高频噪声和低频漂移，使地震波的频谱更加清晰。设置滤波器的截止频率，对于高频噪声，选择合适的高通截止频率，滤除高于该频率的成分；对于低频漂移，采用低通截止频率，去除低于该频率的成分。通过这种方式，确保地震波的主要频率成分能够真实反映地震动的特性。根据工程场地的设计地震动参数，对地震波的幅值进行调整。该拱坝的设计地震动峰值加速度为0.2g，将选取的地震波峰值加速度调整到0.2g。采用线性缩放的方法，根据设计峰值加速度与原始地震波峰值加速度的比值，对地震波的每个时间点的加速度值进行缩放。若原始地震波的峰值加速度为0.1g，而设计峰值加速度为0.2g，则将地震波的每个加速度值乘以2，以满足设计要求。为了考虑地震波的不确定性，对地震波进行了相位调整。通过随机改变地震波的相位，生成多条具有不同相位特征的地震波。在分析过程中，采用这些不同相位的地震波进行计算，然后对计算结果进行统计分析，以评估地震波相位对拱坝地震响应的影响。通过相位调整，可以更全面地考虑地震波的随机性，提高分析结果的可靠性。4.3.2计算工况设置为全面分析拱坝在不同地震工况下的地震响应，确定了多种计算工况。考虑不同地震波输入方向，设置了单向输入和双向输入两种情况。单向输入包括顺河向（X向）、横河向（Y向）和竖向（Z向）输入。双向输入考虑了顺河向与竖向（X+Z向）、横河向与竖向（Y+Z向）的组合。在单向输入工况下，分别将选取的地震波沿X向、Y向和Z向单独输入到有限元模型中，计算拱坝在单一方向地震作用下的地震响应。在X向输入工况下，分析拱坝在顺河向地震作用下的应力、位移和加速度分布情况。在双向输入工况下，同时将地震波沿两个方向输入，考虑两个方向地震作用的耦合效应。在X+Z向输入工况下，研究顺河向和竖向地震作用共同作用下拱坝的地震响应，分析两个方向地震作用的相互影响。针对不同的地震强度，设置了多遇地震、设防地震和罕遇地震三种工况。多遇地震的超越概率为63%，设防地震的超越概率为10%，罕遇地震的超越概率为2%。根据地震危险性分析结果，确定不同地震强度下的地震波峰值加速度。多遇地震下的峰值加速度为0.1g，设防地震下为0.2g，罕遇地震下为0.4g。在不同地震强度工况下，分别输入相应峰值加速度的地震波，分析拱坝在不同地震强度下的地震响应，评估拱坝在不同地震作用下的抗震性能。考虑地震波持时对拱坝地震响应的影响，设置了不同持时的计算工况。根据该地区的地震特性和相关研究成果，选取了10s、20s和30s三种持时。在每个持时工况下，分别输入不同地震波和不同地震强度的地震波，分析拱坝在不同持时地震作用下的累积损伤和应力变化情况。当持时为10s时，分析拱坝在较短时间地震作用下的响应；当持时为30s时，研究拱坝在长时间地震作用下的累积损伤和应力重分布情况。通过设置不同的计算工况，能够全面分析拱坝在各种地震条件下的地震响应，为评估拱坝的抗震安全性提供丰富的数据支持。4.4计算结果与分析4.4.1地震动应力分布规律通过有限元模型，对拱坝在不同地震工况下的地震动应力进行计算，得到了坝体的应力分布云图。在多遇地震工况下，当顺河向输入地震波时，坝体的应力分布相对较为均匀。从坝体的上游面到下游面，压应力逐渐减小，拉应力相对较小。坝体的最大压应力出现在坝底靠近坝肩的部位，约为5MPa，这是由于坝底承受了较大的水压力和地震惯性力，且坝肩对坝体的约束作用使得该部位的应力集中。在坝体的中部和上部，压应力在3-4MPa之间。坝体的拉应力主要分布在坝顶和坝体的边缘部位，最大拉应力约为0.5MPa，这是因为坝顶和边缘部位的约束相对较弱，在地震作用下容易产生拉应力。当横河向输入地震波时，坝体的应力分布与顺河向输入时有一定差异。最大压应力仍然出现在坝底靠近坝肩的部位，但数值略有增加，约为5.5MPa，这是因为横河向地震作用对坝肩的受力影响较大，使得坝肩部位的应力进一步集中。坝体的拉应力分布范围有所扩大，在坝体的两侧和坝顶部位，拉应力相对较大，最大拉应力达到0.8MPa，这是由于横河向地震作用导致坝体两侧和坝顶部位的受力状态发生改变，产生了较大的拉应力。在设防地震工况下，坝体的应力明显增大。顺河向输入地震波时，坝底靠近坝肩部位的最大压应力达到8MPa，坝体中部和上部的压应力在5-6MPa之间。坝体的拉应力也显著增加，坝顶和边缘部位的最大拉应力达到1.5MPa。在坝体的局部区域，如孔洞周边和坝体与地基的交界面附近，出现了明显的应力集中现象，最大应力超过10MPa。这是因为设防地震的强度较大，坝体的变形和受力更加复杂，孔洞周边和交界面附近的应力集中效应更加明显。横河向输入地震波时，坝底靠近坝肩部位的最大压应力达到9MPa，坝体两侧和坝顶部位的拉应力进一步增大，最大拉应力达到2MPa。在坝体的一些薄弱部位，如坝体的角点和裂缝附近，应力集中现象更为严重，最大应力超过12MPa。这些部位的应力集中可能导致坝体混凝土开裂和局部破坏，对坝体的抗震安全性构成威胁。在罕遇地震工况下，坝体的应力急剧增大。顺河向输入地震波时，坝底靠近坝肩部位的最大压应力超过12MPa，坝体中部和上部的压应力在8-10MPa之间。坝体的拉应力分布范围更广，坝顶和边缘部位的最大拉应力达到3MPa。在坝体的多个区域，出现了超过混凝土抗拉强度的拉应力，导致坝体混凝土开裂，裂缝从坝体表面向内部扩展。坝体与地基的交界面处也出现了较大的应力集中，可能导致坝体与地基的连接失效。横河向输入地震波时，坝底靠近坝肩部位的最大压应力超过15MPa，坝体两侧和坝顶部位的拉应力超过3.5MPa。坝体的薄弱部位，如角点和裂缝附近，应力集中现象极为严重，最大应力超过18MPa。坝体的混凝土开裂和破坏现象更加严重，部分区域出现了混凝土脱落和局部坍塌的情况。这些现象表明，在罕遇地震作用下，坝体的抗震性能受到了极大的挑战，可能会发生严重的破坏。4.4.2关键部位动应力时程响应选取拱坝的坝肩、坝踵和坝趾等关键部位，对其在地震过程中的动应力时程响应进行分析，以评估结构的抗震安全性。在坝肩部位，顺河向地震作用下，动应力时程曲线呈现出明显的波动特征。在地震波的峰值时刻，动应力迅速增大，达到最大值后又逐渐减小。在多遇地震工况下，坝肩顺河向的最大动压应力约为4MPa，最大动拉应力约为0.4MPa。随着地震强度的增加，在设防地震工况下，最大动压应力达到7MPa，最大动拉应力达到1.2MPa；在罕遇地震工况下，最大动压应力超过10MPa，最大动拉应力达到2.5MPa。横河向地震作用下，坝肩的动应力变化规律与顺河向类似，但数值略有不同。多遇地震工况下，最大动压应力约为4.5MPa，最大动拉应力约为0.6MPa；设防地震工况下，最大动压应力达到8MPa，最大动拉应力达到1.5MPa；罕遇地震工况下，最大动压应力超过12MPa，最大动拉应力达到3MPa。坝肩部位的动应力较大，且在地震作用下波动明显，这是因为坝肩是拱坝的主要支撑部位，承受着拱坝传来的大部分荷载，在地震作用下，坝肩与地基之间的相互作用复杂，导致动应力变化剧烈。如果坝肩的动应力超过其承载能力，可能会导致坝肩失稳，进而影响整个拱坝的安全。坝踵部位在地震作用下，主要承受拉应力。顺河向地震作用时，动拉应力时程曲线在地震波的作用下不断波动。多遇地震工况下，坝踵的最大动拉应力约为0.6MPa；设防地震工况下，最大动拉应力达到1.5MPa；罕遇地震工况下，最大动拉应力超过2.8MPa。横河向地震作用下，坝踵的动拉应力也有类似的变化趋势，多遇地震工况下最大动拉应力约为0.8MPa，设防地震工况下达到1.8MPa，罕遇地震工况下超过3MPa。坝踵部位的拉应力过大容易导致混凝土开裂，裂缝的发展会削弱坝体的强度和整体性。当坝踵的动拉应力超过混凝土的抗拉强度时，坝踵处会出现裂缝，随着地震作用的持续，裂缝可能会不断扩展，影响坝体的防渗性能和结构安全。坝趾部位在地震作用下主要承受压应力。顺河向地震作用时，动压应力时程曲线在地震波的作用下呈现出波动上升的趋势。多遇地震工况下，坝趾的最大动压应力约为4.5MPa；设防地震工况下，最大动压应力达到7.5MPa；罕遇地震工况下，最大动压应力超过11MPa。横河向地震作用下，坝趾的动压应力变化规律与顺河向相似，多遇地震工况下最大动压应力约为5MPa，设防地震工况下达到8.5MPa，罕遇地震工况下超过13MPa。坝趾部位承受的压应力过大可能会导致混凝土的局部破坏和失稳。当坝趾的动压应力超过混凝土的抗压强度时，混凝土会发生塑性变形，甚至出现局部破碎的情况，这会影响坝体的稳定性。通过对坝肩、坝踵和坝趾等关键部位动应力时程响应的分析可知，随着地震强度的增加，这些部位的动应力显著增大。在罕遇地震工况下，部分关键部位的动应力已经接近或超过混凝土的强度极限，坝体存在较大的安全隐患。在拱坝的抗震设计中，需要针对这些关键部位采取有效的加固措施，如增加混凝土厚度、配置钢筋或采用预应力技术等，以提高这些部位的承载能力和抗震性能。4.4.3与其他方法结果对比验证将有限元计算结果与拟静力法和反应谱法的计算结果进行对比，以验证有限元模型和计算方法的准确性。在多遇地震工况下，对于坝体的最大压应力，拟静力法计算结果为4.5MPa，反应谱法计算结果为4.8MPa，有限元法计算结果为5MPa。对于坝体的最大拉应力，拟静力法计算结果为0.4MPa，反应谱法计算结果为0.45MPa，有限元法计算结果为0.5MPa。从对比结果可以看出，三种方法计算得到的最大压应力和最大拉应力数值较为接近，但有限元法的计算结果相对较大。这是因为拟静力法将地震作用简化为等效静力荷载，忽略了地震动的时间历程和结构的动态响应，计算结果相对保守。反应谱法虽然考虑了结构的动力特性，但在计算过程中采用了一些简化假设，如振型正交性假设等，导致计算结果与实际情况存在一定偏差。有限元法能够更全面地考虑坝体-库水-地基的动力相互作用、材料非线性和几何非线性等复杂因素，因此计算结果更接近实际情况。在设防地震工况下，拟静力法计算得到的坝体最大压应力为7MPa，反应谱法计算结果为7.5MPa，有限元法计算结果为8MPa。坝体最大拉应力方面，拟静力法计算结果为1MPa，反应谱法计算结果为1.3MPa，有限元法计算结果为1.5MPa。随着地震强度的增加，有限元法与其他两种方法的计算结果差异逐渐增大。这是因为在设防地震工况下，地震作用的非线性效应更加明显，拟静力法和反应谱法的简化假