基于振型叠加法电站进水塔基础应力三维有限元分析.docVIP

基于振型叠加法电站进水塔基础应力三维有限元分析 　　【摘要】电站进水塔为水电站最常见至水工永久建筑结构，在许多长引水电站中，常采用独立岸塔式进水塔。进水塔基础应力分析是进水塔设计中比较重要的一部分。而在分析进水塔基础应力时，库区正常蓄水位时发生地震的情况乃重点关注之工况。简单的拟静力解析法计算分析基础应力仿真度不高，近年来在工程中多采用线弹性动力三维有限元法进行计算。线弹性动力有限元通过弹性力学变分原理建立动力问题有限单元表达格式――结构运动方程，在地震分析的数值解法中，可采用振型叠加法、直接积分法和振型分解反应谱法求解。对于大型复杂的进水塔模型，采用直接积分法较费时间和资源，通常可采用振型叠加法或振型分解反应谱法。本文以云南保山某电站进水塔地震工况为例，采用基于振型叠加法的线弹性动力三维有限元分析方法，对进水塔基础应力进行分析研究，为类似工程的进水塔设计提供有益思路。 　　【关键词】电站进水塔；基础应力；振型分解反应谱法；线弹性动力三维有限元； 　　1. 电站进水塔布置概况 　　该电站进水口采用岸塔式，底板高程1837.00m；塔顶高程为1900.00m，塔高63m。塔宽最宽处14.0m。塔前正常蓄水位EL.1895.000m，设计洪水位EL.1895.429m，校核洪水位EL.1896.301m，死水位EL.1848.000m。 　　地震设防烈度为8度，a=0.19g。 　　进水塔混凝土强度等级为C25 　　2. 计算荷载及工况 　　（1）进水口结构及设备自重 　　进水口自重按结构的体积计算，在大气中混凝土容重 =24KN/m3，钢筋混凝土容重 =25KN/m3。当自重效应对结构不利时作用分项系数取1.05，对结构有利时作用分项系数取0.95。 　　（2）静水压力 　　按实际水深进行计算，水的容重γ=9.8KN/m3，作用分项系数1.0。 　　（3）扬压力 　　按塔体与地基接触面处的相应水头进行计算，作用分项系数1.0，可仅考虑浮托力。 　　（4）浪压力 　　按《水工建筑物荷载设计规范》（DL 5077-1997）计算，作用分项系数1.2，本计算中未考虑。 　　（5）风荷载 　　按《水工建筑物荷载设计规范》（DL 5077-1997）计算，作用分项系数1.3，本计算中未考虑。 　　（6）地震作用 　　设防烈度为8度，a=0.19g，荷载分项系数取1.0。 　　（7）计算工况 　　正常蓄水位1895.000m运行+地震工况，荷载组合为：塔体结构自重+静水压力+扬压力+地震力 　　3. 计算理论及方法 　　线弹性动力有限元通过弹性力学变分原理建立动力问题有限单元表达格式――运动方程。动力学分析与静力学分析的主要区别是基本方程中的各参数都是时空的函数，因而方程中将有加速度、速度、初始条件，且引入质量矩阵和阻尼矩阵等内容。动力有限元法不同于静力有限元法的思想是：将结构只对空间域离散为一系列连续分布的单元，仍以单元结点位移为基本变量，由平衡方程及力的边界条件的等效积分形式的Galerkin法得到单元的运动方程，集成得结构运动方程，方程为常微分方程组，在地震分析的数值解法中，可采用振型叠加法、直接积分法和振型分解反应谱法求解。其它计算过程同静力有限元思想相同。 　　3.1 运动方程 　　形成单元阻尼矩阵和质量阻尼矩阵，以及结构运动方程，是不同于静力分析的内容，也是动力分析中的最重要部分，简述如下： 　　由平衡方程及力的边界条件的等效积分形式的Galerkin法得到的运动方程为： 　　（3.1-1） 　　式中：[M]―结构质量矩阵；[C]―结构阻尼矩阵；[K]―结构刚度矩阵； ―结点位移向量； ―结点速度向量； ―结点加速度向量。 　　结构质量矩阵[M]、阻尼矩阵[C]、刚度矩阵[K]和结点荷载向量分别由各自的单元矩阵和向量集成。 　　（3.1-2） 　　（3.1-3） 　　式（3.1-1）右项取地震载荷时，经过对上式变形整理可得类似方程： 　　（3.1-4） 　　式中： ―结点相对地面位移向量； ―结点相对地面速度向量； ―结点相对地面加速度向量； ―地面加速度向量。 　　3.2 计算方法 　　直接积分法是指在积分运动方程之前不进行方程形式的变换，而直接进行逐步数值积分。对于大型复杂的模型，采用直接积分法较费时间和资源，通常可采用振型叠加法或振型分解反应谱法。振型叠加法是在积分运动方程以前，对自由振动方程进行模态分析，获得结构的固有振型和固有频率，利用系统自由振动的固有振型将方程组转换为n个不相互耦合的方程，对这种方程可以解析或数值地进行积分。对时间历程较长，只需少数较低阶振型的结果时，振型叠加法是十分有效的。 　　地震分析中，由于其偶然性和地域性特点，荷载随