基于反应谱法的重力坝有限元等效动应力响应分析
基于状态空间理论的砼重力坝振动响应分析
基于状态空间理论的砼重力坝振动响应分析摘要：介绍了一种基于状态空间理论的多自由度结构系统动力响应分析算法，将结构动力学的平衡方程经过一定的变换得到状态空间方程，再将该方程嵌入卡尔曼滤波算法，最终得到多自由度系统的状态随着系统输入的实时更新。
将该算法运用于混凝土重力坝算例，分别进行不同观测噪声水平的动力响应计算，并且将得到的动力响应信息与有限元软件数值模拟的结果进行比较，均具有较高的相似度，说明该算法的有效性，并且抗噪能力强，具有一定的工程应用价值。
关键词：状态空间理论；重力坝；动力响应；振动一、引言结构的动力响应分析和求解方法主要有振型分解反应谱法[1-2]、底部剪力法[3-4]和时程分析法[5-6]。
振型分解反应谱法计算精度较高，但是其需要用到结构的前大多数阶振型，而实际结构的自由度相当大甚至无穷大，这必然导致计算量的剧增，实际应用时具有一定的不可操作性；而底部剪力法仅适用于以剪切变形为主的重量和刚度沿高程分布均比较均匀的结构（否者会出现鞭梢效应[7]），当建筑结构的高度在40m以上或者结构的基频较小时，地震作用下将不能忽略结构的高阶振型影响，这必然导致其计算精度的降低；理论上，时程分析法是最准确的结构动力响应分析算法，可以在全时域内对结构进行实时分析[8-9]，但是由于其求解的复杂性以及外界激励的随机性，在实际工程运用时具有一定的困难，但随着计算机技术的发展和控制论中状态空间理论的发展，为这一精确方法的实际运用提供了捷径。
二、结构动力学的状态空间方程由结构动力学的基本知识[18]写出N个自由度结构系统的平衡方程：d（1）式中：M，L和K分别为多自由度结构的质量矩阵、阻尼矩阵和刚度矩阵；g(t)为N维激励。
g定义状态矢量y，其包含有位移矢量和速度矢量：（2）故方程（1）可以转换成状态空间形式：（3）式中：系统矩阵和控制矩阵如下：（4）（5）式中：表示a×b的零矩阵，表示a×b的单位矩阵。
重力坝分析的改进有限元等效应力法
重力坝分析的改进有限元等效应力法杨会臣;贾金生;郑璀莹【摘要】本文提出了重力坝分析的改进有限元等效应力法,分析了改进前后应力计算结果的变化,研究表明,当网格尺寸较大时,改进的等效应力方法具有更高的精度.对比研究了采用有限元等效应力方法、材料力学方法计算应力在结果上的相同点和不同点.分析表明,有限元等效应力法在解决应力集中的同时,还能反映地基刚度对建基面竖向正应力分布的影响,从而解决了材料力学的一些限制.结合我国第一座胶结颗粒料永久工程—守口堡工程的设计,详细分析了地基刚度对坝基面应力分布的影响,取得了很好效果.【期刊名称】《中国水利水电科学研究院学报》【年(卷),期】2013(011)002【总页数】5页(P112-116)【关键词】重力坝;有限元等效应力方法;材料力学方法;守口堡工程【作者】杨会臣;贾金生;郑璀莹【作者单位】中国水利水电科学研究院流域水循环模拟与调控国家重点实验室,北京100038;中国水利水电科学研究院流域水循环模拟与调控国家重点实验室,北京100038;中国水利水电科学研究院流域水循环模拟与调控国家重点实验室,北京100038【正文语种】中文【中图分类】TV641.31 改进的有限元等效应力法有限元等效应力法由我国学者首先提出[1-4]，其基本思想是将某一截面上有限元计算的应力结果等效为作用在该截面上的一组内力，然后利用材料力学的方法计算截面上的应力分布，以避免角缘处有限元计算结果的应力集中。
有限元等效应力法分析拱坝的文献[5-7]很多，但用于研究重力坝应力分布的文章还比较少[8-9]。
图1 有限元应力计算网格对于图1所示的有限元计算网格，xi处的竖向正应力为，xi+1处的竖向正应力为，则该截面上的等效内力为[8]：式中：N为截面轴力；M为截面弯矩；x0为截面转动中心；li=xi+1-xi；B为截面宽度，取1。
根据材料力学偏心受压构件正应力计算公式，结点i处的有限元等效应力即为：其中：L为截面长度；I为截面绕x0的转动惯量。
第四节重力坝的应力分析
W 6 M
a
B
B2
b
12 M B2
坝体内部应力计算图
（三）内部应力的计算
2、坝体内剪应力τ。呈抛物线分布
a1 b1x c1x2
3、坝内水平正应力σx。呈三次分布
x a2 b2 x c2 x2 d2 x3 接 近 直
线，对中小型工程可近似假定σx呈直线分布
运用期坝踵垂直应力不应出现拉应力（计扬压力）， 可按下式计算：
WR M RTR 0
AR
JR
WR B
6M R B2
0
规范要求，坝踵和坝体上游面的垂直应力的核算应按作用的标准值分别 计算作用的短期组合和长期组合。
2)短期组合下游坝面的垂直拉应力核算
•施工期属短暂状况，坝体下游面的垂直拉应力应不大于
x a3 b3x
4、坝内主应力σ1和σ2。(见课本P56式(3-48))
(四)考虑扬压力时的计算方法。
1、边缘应力 1)计算W和M计入扬压力 2)计算u 、xu、d、xd时
pu用 pu-puu代入 pd 用pd-pdd代入 2、求解坝内应力 可先不计扬压力计算、x、y，然后再 叠加由扬压力引起的应力
1、水平截面上的正应力σyu、σyd。 2、剪应力τu和τd。 3、水平正应力σxu和σxd。 4、主应力σ1u，σ2u和σ1d，σ2d。
1、水平截面上的正应力σyu、σyd。
根据偏心受压公式，坝
体上下游边缘垂直正应
力：
yu
W B
6M B2
yd
W B
6M B2
2、剪应力
4、边缘主应力
由于两个主应力面互相正交，由微分体
有限元等效应力法在重力坝强度分析中的应用
有限元等效应力法在重力坝强度分析中的应用范书立;陈健云;郭建业【期刊名称】《水利学报》【年(卷),期】2007(038)006【摘要】采用有限元法进行重力坝坝体应力分析时,坝踵、坝址的应力计算结果一直无法作为设计坝体断面的依据.本文将有限单元法所求得的应力合成截面内力,用材料力学公式求出重力坝建基面上任意点对应的线性化等效应力,并分析了影响有限元等效应力的各种因素.研究结果表明,有限元等效应力法有效地避免了坝踵、坝址的应力集中现象,解决了坝踵应力随网格的变化而变化的问题,计算出的建基面应力具有数值稳定性,得出的应力结果可以按现行的重力坝设计规范规定的应力标准评价重力坝坝体安全度.因此,结合有限元结果和等效应力法计算的坝踵、坝址应力结果,可为重力坝坝体强度的综合评定提供相应的依据.【总页数】7页(P754-759,766)【作者】范书立;陈健云;郭建业【作者单位】大连理工大学,海岸与近海工程国家重点实验室,辽宁,大连,116023;大连理工大学,海岸与近海工程国家重点实验室,辽宁,大连,116023;大连理工大学,海岸与近海工程国家重点实验室,辽宁,大连,116023【正文语种】中文【中图分类】TV64.3【相关文献】1.有限元等效应力法在拱坝设计中的应用 [J], 肖伟荣;苏志敏;唐涛2.重力坝分析的改进有限元等效应力法 [J], 杨会臣;贾金生;郑璀莹3.等效结构应力法在高强钢疲劳寿命预测中的应用 [J], 邹艳妮;李耀;罗东4.等效结构应力法在塔机关键焊缝疲劳强度研究中的应用 [J], 付玲;李耀;罗东;吴达鑫;佘玲娟5.等效结构应力法原理及其在转向架焊接构架疲劳寿命分析中的应用 [J], 薛俊谦;李向伟因版权原因，仅展示原文概要，查看原文内容请购买。
第四节--重力坝的应力分析
Gravity Dam Located onBatholith
第四节重力坝的应力分析
目的：
1、为了检验大坝在施工期和运用 期是否满足强度要求；
2、为解决设计和施工中的某些问 题，如砼分区，某些部位的配 筋等提供依据。
应力分析的过程：
1、进行荷载计算及荷载组合 2、选择合适的方法进行应力计算 3、检验大坝各部位的应力是否满
位等）的局部应力、个别部位（如宽缝重力坝的头部、闸 墩、导墙等）的应力等步骤，必要时分析坝基的上、下游 局部应力及内部应力。
三、材料力学方法
（一）基本假定 1、坝体砼为均质，连续各向同性
的弹性材料。 2、取单宽坝体作为固结在地基上
的悬臂梁计算，且不受两侧坝体的影 响。
3、水平断面上的垂直正应力σy是 直线分布。 （二）**边缘应力的计算
2)短期组合下游坝面的垂直拉应力核算
•施工期属短暂状况，坝体下游面的垂直拉应力应不大于
100kPa，其计算式为
Wc M cTc 100 (kPa)
Ac
Jc
第四节 重力坝的应力分析
其他： 坝体内一般不容许出现主拉应力，但以下情况例外：①宽
缝重力坝离上游面较远的局部区域，可出现拉应力，但不 得超过混凝土的容许拉应力；②当溢流坝堰顶部位出现拉 应力时，可考虑配置钢筋；③廊道及其它孔洞周边的拉应 力区域，宜配置钢筋，以承受拉应力。
坝体主应力 分布示意图
作业2
图与荷载同作业1，砼强度等级C10，标准抗压 强度fkc=10MPa，坝基为较完整的微风化花岗片 麻岩，标准抗压强度fkR=80MPa 试核算基本组合的设计洪水位情况下 (1)计算坝基面A、B及折坡处水平面C、D点的应 力x、y、及主应力； (2)坝趾B抗压强度和坝踵A应力是否满足要求； (3)根据所计算的稳定安全系数及应力情况．讨论 此坝断面设计是否得当；
基于反应谱法的重力坝有限元等效动应力响应分析
Wj 6 M j B - B2
E E 式中: Rj, yu , Rj, yd 分别为上下游边缘正应力; Wj 为计算截面上与第 j 阶振型对应的轴向约束力; Mj 为计算
摘 要: 提出了由有限元法得到的振型求解振型内力及振型等效 应力, 并根据反应 谱法理论 求解重力 坝等效 动应力 响
应的 分析方法。本文对四个在建或拟建的重力坝, 用悬臂梁反应谱 法和有限元等效应力法分别进行了计算, 比较了坝 体
几个 典型部位等效动应力响应。结果表明, 在坝体体型有突变的地 方, 有限元等效应力法相较悬臂梁法所 得值有明显 的 增高 , 更敏锐的反映出了应力集中对坝体垂直动应力的影 响。本文 提出的方法简单易行, 比悬臂梁反应谱 法更能反映 坝
N2j ) C2X
第 5期
李同春等: 基于反应 谱法的重力坝有限元等效动应力响应分析
21
式中: SE 为地震作用效应; Si 、Sj 分别为第 i 阶、第 j 阶振型的地震作用效应; m 为计算采用的振型数; Qij 为第 i
阶、第 j 阶振型的振型相关系数; Ni , Nj 分别为第 i 阶、第 j 阶振型的阻尼比; CX 为圆频率比, CX = XjPXi ; Xi , Xj 分
由弯曲变形和剪切变形产生的水平振动, 忽略坝体的竖向振动, 也忽略转动惯量 的影响[ 3] 。如图 1 所示, 将重力坝简化为一固结于地基的悬臂梁, 取单位厚度重
力坝坝段为分析对象, 将整个坝段沿高度按任意间距划分成 N - 1 段, 结点号自
下至上为 1 至 N, 结点 1 固定于地基上。取 y 轴向上, 当只考虑坝体由弯曲变形
体体型的影响。
关键词: 水工结构; 等效动应力; 反应谱法; 悬臂梁法; 有限元法; 重力坝
基于重力坝应力计算及稳定分析的优化设计
基于重力坝应力计算及稳定分析的优化设计重力坝是一种常见的水利工程结构，其稳定性是设计中需要考虑的重要问题。
在设计重力坝时，需要对其应力进行计算和稳定性进行分析，并进行优化设计。
首先，重力坝的应力计算需要考虑以下几个方面。
首先是坝体自重的计算，包括坝体上升水压力和上升地下水压力。
其次是坝顶压力的计算，包括抗倾覆稳定和抗滑移稳定的力学分析。
还需要考虑水侧坝体的压力计算，包括水压头的作用和大坝的承压强度。
最后是岩质坝体的应力分析，考虑岩性、节理的影响及坝体的变形与稳定性。
为了保证重力坝的稳定性，需要进行稳定分析。
稳定分析主要包括抗倾覆稳定和抗滑移稳定两个方面。
抗倾覆稳定分析是为了防止重力坝在承受水压力的作用下发生倾覆。
抗滑移稳定分析是为了防止重力坝在地基土的滑移力的作用下发生移动。
通过合理选择坝体的高度、坝基的强度和选择合适的岩质材料，可以有效地提高重力坝的稳定性。
在重力坝的优化设计中，可以从以下几个方面进行考虑。
首先是合理选择重力坝的形式，可以是三角形、梯形或者圆弧形等不同形式，根据工程实际情况进行选择。
其次是选择合适的坝基处理措施，包括混凝土垫层、防渗墙等，可以提高坝体的稳定性。
另外，可以考虑采用辅助措施，如设置消能防冲设施、阻水排水系统等，提高重力坝的安全性。
最后，可以进行不同形式的优化设计，如遗传算法、模拟退火算法等，寻找最优设计方案，既能满足工程要求，又能提高工程的经济性和可行性。
综上所述，基于重力坝的应力计算及稳定分析的优化设计是一个综合性的工程问题。
通过合理的应力计算和稳定分析，可以提高重力坝的稳定性。
同时，通过优化设计，可以选择合适的形式和措施，提高工程的安全性和经济性。
因此，在重力坝的设计中，需要综合考虑各种因素，进行全面的分析和优化设计。
基于有限元法的大坝应力、变形及稳定分析
论￣弹性理论数值解法的基本概念和步骤 [９] 如下ꎮ
物ꎬ 主要利用自身重力维持坝身的稳定ꎮ 通常采用
反应坝体及坝基不同部位的应力分布和变位场ꎮ 有
的单元ꎬ 采用位移函数法ꎬ 公式如下:
限元法( ＦＥＡ) 起始于 ２０ 世纪中叶ꎬ 经过 ７０ 年的不
{ ｕｖ} ＝ [ Ｎ] { δ}
断 研 究 发 展ꎬ 并 通 过 ＡＮＳＹＳ、 ＡＤＩＮＡ、 ＭＡＲＫ、
{ δ} ＝ [ Ｋ] －１ { Ｐ}
再利用式(３) 和式(４) 计算单元应力ꎮ
(７)
１ ２ 抗滑稳定计算公式
混凝土重力坝抗滑稳定安全系数是分析其稳定
的重要 指 标ꎮ 混 凝 土 重 力 坝 设 计 规 范 抗 剪 断 公
式 [１０] 采用的是刚体极限平衡法ꎮ 本次以抗剪断公
式为基础ꎬ 将有限元计算出的滑动面各单元 σ ｎｉ 和
计算方法ꎬ 基于有限元法分析后的坝体平面应力
{ σ} ＝ [ Ｄ] { ε}
成果ꎬ 代入调整后抗滑稳定系数公式求出稳定系
{ σ} ＝ [ σ Ｘ σ ｙ τ ｘｙ ] Ｔ
数ꎮ 文中以非溢流坝段为例ꎬ 采用有限单元法计
定性ꎮ
１２８
}
式中ꎬ { σ} —应力矩阵ꎬ [ Ｄ] —弹性矩阵ꎮ
算应 力、 变 形、 抗 滑 稳 定 系 数ꎬ 综 合 评 定 其 稳
２０１９ 年第 １１ 期
设计施工
水利规划与设计
ＤＯＩ: １０ ３９６９ / ｊ ｉｓｓｎ １６７２￣２４６９ ２０１９ １１ ０３２
基于有限元法的大坝应力、 变形及稳定分析
艾子欣ꎬ 陈海霞ꎬ 郭穗丰
( 湖南省常德市水利水电勘测设计院ꎬ 湖南 常德 ４１５０００)
摘要: 文章采用有限元法ꎬ 通过建模、 网格化分、 加载、 求解、 后处理等一系列过程ꎬ 对该坝的应力、 变形及抗
基于三维有限元的某重力坝应力变形特性分析
基于三维有限元的某重力坝应力变形特性分析乌日晗【摘要】重力坝作为主要坝型之一在人类筑坝治水历史长河中发挥着举足轻重的作用,特别是随着我国西部水电能源建设的飞速发展,将有一批高重力坝水电工程投入建设和运营,这些大坝均位于高山峡谷之中,工程规模大,地质环境复杂,其坝与地基稳定安全问题是须长期关注和解决的关键科学技术难题.本文采用三维有限元法分析了某重力坝正常工况下的应力变形特性,对其正常工作状态和应力变形规律进行了深入分析.【期刊名称】《吉林水利》【年(卷),期】2018(000)001【总页数】4页(P32-34,55)【关键词】重力坝;应力;位移;有限元【作者】乌日晗【作者单位】辽宁省蓡窝水库管理局,辽宁辽阳 111000【正文语种】中文【中图分类】TV642.31 引言重力坝的历史非常的悠久，大约在公元前两千九百年，人类就建造了第一座重力坝。
但由于重力坝结构简单，安全稳定，到现在还在被广泛采用[1]。
重力坝是在上游水压力、下游水压力及坝基扬压力等荷载的作用下，主要依靠坝体自身重量产生的抗滑力来维持坝体的稳定，尽管该类坝型体积足够大，本身具有很强的稳定能力[1]，但是由于各种原因，重力坝失事还时有发生。
大坝的失事往往会给国民经济和人民的生活水平造成比较大的影响，美国的圣弗朗西斯大坝，基础沉陷，坝体开裂严重，最终滑移崩溃，伤亡四百多人[2]。
美国奥斯汀大坝，坝基岩体中页岩泥化，坝体抗剪强度剧烈降低，引起坝基渗漏，导致坝体和坝基滑移，使整个坝体遭到毁灭性的破坏，给下游群众带来了灾难性的破坏[3]。
人类一直十分关注和重视重力坝的应力应变状态和坝体的抗滑稳定。
现今，人们对重力坝的抗滑稳定问题认识更深，提出了大量分析解决方法。
但是，由于很多因素都可以影响到大坝的抗滑稳定，而且这些因素中的有些因素还很难确定，到现在仍没有充分的理论依据，很大程度上还是依赖于工程师的工程经验。
因此，分析重力坝在正常运行过程中的应力变形特性、合理评价大坝工作性态对于长期安全运行是非常有必要的。
坝体的有限元建模与应力应变分析1
Project2 坝体的有限元建模与应力应变分析计算分析模型如图2-1 所示, 习题文件名: dam 。
图2－1 坝体的计算分析模型选择单元类型Solid Quad 4node 42 Options… →select K3: Plane Strain定义材料参数EX:2.1e11, PRXY:0.3模型施加约束✓ 分别给下底边和竖直的纵边施加x 和y 方向的约束✓ 给斜边施加x 方向的分布载荷：ANSYS 命令菜单栏: Parameters →Functions →Define/Edit →1) 在下方的下拉列表框内选择x ,作为设置的变量；2) 在Result 窗口中出现{X},写入所施加的载荷函数：1000*{X}；3) File>Save(文件扩展名：func) →返回：Parameters →Functions →Read from file ：将需要的.func 文件打开，任给一个参数名，它表示随之将施加的载荷→OK →ANSYS Main Menu: Solution →Define Loads →Apply →Structural →Pressure →On Lines →拾取斜边；OK →在下拉列表框中，选择：Existing table →OK →选择需要的载荷参数名→OK单元控制 纵边20等分；上下底边15等分结果显示ANSYS Main Menu: General Postproc →Plot Results →Deformed Shape… → select Def + Undeformed →OK (back to Plot Results window)→Contour Plot →Nodal Solu… →select: DOF solution, UX,UY, Def + Undeformed , Stress ,SX,SY,SZ, Def + Undeformed →OK。
有限元等效应力法在重力坝强度分析中的应用
2007年6月水 利 学 报SHUILI XUEBAO第38卷 第6期收稿日期:2006-07-10基金项目:国家自然科学基金资助项目(90510018);新世纪优秀人才支持计划资助项目(NCET-06-0270)作者简介:范书立(1978-),男,山东高唐人,博士生,主要从事水工结构研究。
E -mail:shl fan@文章编号:0559-9350(2007)06-0754-06有限元等效应力法在重力坝强度分析中的应用范书立,陈健云,郭建业(大连理工大学海岸与近海工程国家重点实验室,辽宁大连 116023)摘要:采用有限元法进行重力坝坝体应力分析时,坝踵、坝址的应力计算结果一直无法作为设计坝体断面的依据。
本文将有限单元法所求得的应力合成截面内力,用材料力学公式求出重力坝建基面上任意点对应的线性化等效应力,并分析了影响有限元等效应力的各种因素。
研究结果表明,有限元等效应力法有效地避免了坝踵、坝址的应力集中现象,解决了坝踵应力随网格的变化而变化的问题,计算出的建基面应力具有数值稳定性,得出的应力结果可以按现行的重力坝设计规范规定的应力标准评价重力坝坝体安全度。
因此,结合有限元结果和等效应力法计算的坝踵、坝址应力结果,可为重力坝坝体强度的综合评定提供相应的依据。
关键词:有限元;等效应力;重力坝;设计准则中图分类号:TV6413文献标识码:A1 研究背景有限元法由于使用方法灵活,适用范围广泛,能够有效地解决复杂结构的应力分析,自20世纪60年代末被用于计算混凝土重力坝坝体应力后,已迅速发展成为进行坝工设计和坝工学科定量研究的一个重要手段。
但是,有限元法计算的应力结果与所采用的单元型式和单元网格划分疏密有较大的关系,在坝踵、坝趾等角缘部位存在应力集中现象,网格愈密应力集中程度愈高,而这些部位往往是坝工设计最受关注的控制部位,这给大坝的应力安全评价以及合理确定大坝的应力控制标准带来困难。
因此,有限元法分析中控制应力取值标准一直是坝工界长期关注且尚未得到很好解决的关键问题之一。
基于有限元等效应力的重力坝坝踵应力分析
．
—
—
基于有 限元 等效应力的重力坝坝踵 应力分析
王建伟 徐 ， 礼 姚 江强 ，
（． 昌大 学建 筑 ３程 学 院 ， 西 南 昌 ３ ０３ ； １南 － ． 江 ３０ １
２江 西省 丰城 市公 安 消 防 大队 ， 西 丰城 ３ １ ０ ） ． 江 ３ １ ０
摘 要： 采用有限元进行重力坝坝体应力分析时 ， 坝踵应力的取值一直是坝工界长期关注的问题。基于有限元等效应
Ｙ Ｏ ｉ ｎ － ｉ ｎ ／ Ｃ ｌ ｇ ｆ ｃ ｉｅ ｔ ｒ ｌ ｎ ｉ ｉ Ｅ ｇ ｎ ｅ ｉ ｇ Ｎａ ｃ ａ ｇ Ｕｎｖ ｒ ｉ Ａ Ｊ ａ ｇ ｑａ ｇ ／ ｏｌ ｅｏ ｈ ｔｃ ｕ ａ ｄ Ｃ ｖ ｌ ｎ ｉ ｅ ｒｎ ， ｎ ｈ ｎ ｉｅ ｓｔ ｅ Ａｒ ａ ｙ
Ｂ Ｎ ａ— ｅＡ ａ ｓ ｎ ｔ ｓ ｏ ｇ ｖｙ ａ ｅ ｌａｅ ｎ ｙ ＷＡ ＧＪ ｎ ｗ ｉ ｎｌ ｉｏ ｒ ｓ ｆ ｉ ｍ ｈｅｂｓｄ ｉ ： ｙ ｓ ｓ ｅ ａ ｔｄ ｒ … … ｏ …
—
ｅ ｑ ｉ ａｅ ｔｓ ｒ ｓ ｎ ｌｓｓ ｕｖ ｌ ｎ ｔｅ ｓａ ａｙ ｉ
ｗ ｔ ｏ ｔ ｌ ｆ ｈ ｄ ｗ ｄ ｈ ｉ ｃ ｎｒ ｅｍ ｉｔ． ｈ ｏｏｔ Ｋｅ ｒ ｓ ｓ ｅ ｓｏ ａ ｅ ｌ ｓｒ ｓ ｎ ｅ ； ｒｄ ｃ ｎ ｒ ｌ ｙｗｏ ｄ ： ｔ ｓ ｆ ｍ ｈ ｅ ； ｔ ｓ ｄ ｘ ｇ ｉ ｏ ｔ ｒ ｄ ｅ ｉ ｏ
中图 分 类号 ： ｖ ６ ８１ Ｔ ９ ．
设计 中关注的数值 ， 因此有限元分析中的应力取值 标准一 直是非 常重要但 又存在 许多争议 的问题 。
为 了使有 限元 分析 方 法 能在 工程 上 得 到 广泛
应用 ， 多学者 提 出了各种坝 踵应 力取值 问题 。敖 许
利用有限元法的大坝变形分析研究
利用有限元法的大坝变形分析研究摘要：随着经济社会的发展，水利工程也得到了较快的推进，极大程度丰富了人们的生活。
但是大坝有可能还存在着变形问题，对于其质量造成严重的影响，进而威胁到人们的正常生活。
本文以具体工程为例进行分析，论述了利用有限元法来分析大坝变形。
关键词：有限元；大坝变形；分析引言我国是地震多发的国家，建设重力坝的地区一旦发生地震，将带来灾难性的后果。
在修建大坝设计中将地震荷载纳入考虑范围，考察重力坝在地震荷载下的动力响应特性，研究重力坝在静力、动力荷载下的安全问题具有十分重要的意义。
国内外诸多学者用各种非线性模型分析了混凝土重力坝在地震作用下的反应。
然而，总体来看，大坝安全评价准则目前正处在发展阶段，如何有效合理地利用有限元计算结果受到广大研究者的普遍注。
1、工程概况某水库是一座综合利用的大Ⅱ型水利枢纽工程，2015年12月通过了竣工验收。
总库容10137万立方米，控制流域面积875平方千米，坝型为混凝土面板堆石坝，最大坝高88.3m，坝顶长310m，坝顶高程1138.3m。
大坝变形观测分人工观测和自动化观测两大部分，建立了相对完善、全面的变形观测系统。
大坝变形监测工作范围主要包括：1）坝体外部沉降观测；2）坝体内部沉降观测；3）坝体外部水平位移观测；4）坝体内部水平位移观测；5）溢洪道沉降观测。
由于人工观测环境与测量仪器限制，大坝体外水平位移观测目前未能取得观测数据，本文对此项数据不做分析。
2、监测资料分析变形监测在测量领域内占据着重要的位置，从一个工程的施工到完工，以及后续的运营都需要进行不断地监测，掌握变形的情况，及时解决潜在安全问题，保证工程的正常运营。
在大坝变形监测中，传统的变形监测是采用高精度的监测网对大坝变形要素进行监测，但由于大坝所处地形条件的影响，导致监测网的网形差和监测点的位置精度不精确，影响测量的准确性。
这种方法的劳动强度很大，观测时间较长，没有实现自动化监测。
将大坝历年各期资料汇总成数据库，将地震期间监测数据与震前数据、历年同期数据进行比。
基于振型分解反应谱法的重力坝三维有限元动力分析
104 2. 80
104 3. 00
104 3. 15
0. 167 2 400
0. 167 2 400
0. 167 2 400
0. 167 2 400
用子空间迭代法对大坝进行模态分析 , 再采 用振型分解 反应谱 法进行动力计算 , 分析大 坝的 地震动 力响 应及其 分布 规律 ; 最 后按照 强度最不利应力叠 加原 则 将静 态荷 载和动 态荷 载作 用下的反应进行组合 , 分析了坝体的综合反应。 大坝的自振特 性分析 采用 子空 间迭代 法。为了 保证 计算 精度 , 应 计算尽可能多的振型以进行组合。 但对于实际 工程问 题只须取较低的几个振型就能满足计 算要求 , 本文计算 中取前 10 阶振型进行分析。 地震动力分析 采用振 型分 解反 应谱法。 该水电 站坝 址区 场地的地震基本烈度为 7 度 , 根 据规范 [ 4, 5] 要 求壅水 建筑 物抗 震设防烈度概率水准采用基准期 100 年超越概率 2% 。经地震 危险性分 析 , 相当于 P 100 = 2. 0% 超 越概率水 平的该 水电 站坝 址区的基岩水平向设 计地震 加速度 代表 值为 0. 229 g, 竖 向设 计地震加速度代表值取水平向的 2/ 3, 即 0. 153 g 。该工程场地 为 类场地 , 场地的特征周 期 T g = 0. 20 s, 设 计反应 谱最 大值
3 动力计算结果分析
考虑了以下两种计 算工况 : ( 1) 空库 ( 2) 正 常蓄水 位工 况 , 两种工况计算结果如下。
3. 1 坝体自振特性
坝体前 10 阶振型的自振 频率及 周期见表 2 。由于篇 幅有 限 , 文中省略了各阶振型图。
图 2 坝体有限元网格示意图 表 1 大坝混凝土材料物理力学参数表 ( 静态 ) 阶数 混凝土强度等级 项 目 C20 轴心抗压 ( 设计值 ) / M Pa 轴心抗拉 ( 设计值 ) / M Pa 弹性模量 / M Pa 泊松比 密度 / ( k g m- 3) 10. 0 1. 10 2. 55 C25 12. 5 1. 30 C30 15. 0 1. 50 C35 17. 5 1. 65 104 1 2 3 4 5 6 7 8 9 10 频率 / H z 1. 181 4 2. 201 6 3. 318 4 3. 726 4 4. 423 1 5. 270 1 6. 704 9 6. 785 4 6. 910 8 7. 856 9 周期 / s 0. 846 5 0. 454 2 0. 301 4 0. 268 4 0. 226 1 0. 189 7 0. 149 1 0. 147 4 0. 144 7 0. 127 3 频率 / H z 1. 181 4 1. 860 5 3. 285 7 3. 636 6 3. 916 3 5. 239 1 6. 347 3 6. 500 9 6. 710 4 7. 639 7 周期 / s 0. 846 5 0. 537 5 0. 304 3 0. 275 0 0. 255 3 0. 190 9 0. 157 5 0. 153 8 0. 149 0 0. 130 9 表 2 坝体自振频率及周期 空 库 满 库
基于ADINA的重力坝地震响应分析
文章编号:1000-7709(2008)01-0097-03基于ADINA 的重力坝地震响应分析王伟华1 张燎军2(1.河海大学土木工程学院,江苏南京210098; 2.河海大学水利水电工程学院,江苏南京210098)摘要:针对云南阿海重力坝,利用大型有限元软件A DIN A 分别基于流固耦合理论和Wester gaa rd 附加质量方法建立了动力数值分析模型,并进行了自振特性和地震响应计算。
结果表明,两种模型的自振频率较接近,附加质量模型的地震响应大于流固耦合模型。
在重力坝抗震研究中,利用AD IN A 的流固耦合模型考虑水体的影响更接近实际工程。
关键词:A DI NA ;重力坝;流固耦合;附加质量;频率;地震响应中图分类号:T V 642.3;P 631.4文献标志码:A收稿日期:2007-09-10,修回日期:2007-10-10基金项目:国家自然科学基金资助项目(90510017)作者简介:王伟华(1979-),男,硕士研究生,研究方向为结构抗震和设计,E-ma il:w ang w eihua792005@163.co m通讯作者:张燎军(1962-),男,教授、博导、副所长,研究方向为水工结构的数值仿真分析,E -m ail :ljzhang @hhu .edu .cn 考虑动水压力作用的水工结构动力分析方法分两种: 线性法,把水体当作附加质量作用于水工结构上,与结构一起进行动力分析,如附加质量、Housner 模型。
该方法忽略了流体—结构相互作用的影响,自振频率误差较大; 非线性法,研究液体非线性晃动对结构的影响,如边界元法、ALE 有限元法等,适用于求解流体大幅晃动问题,能充分反映流体域的特征。
本文基于非线性有限元仿真分析软件(ADINA)对云南阿海重力坝采用两种分析方法比较地震响应,结果表明韦斯特伽德(Westerg aard )法求解重力坝地震响应是可行且偏于安全。
1 韦斯特伽德附加质量法对大坝等挡水结构,挡水面前的水体可视为无限水域。
基于有限元的重力坝稳定及动力响应分析
[
８]
谱值βmax＝２．
０ ．
图 ４ 重力坝振型云图
图 ５ 重力坝设计反应谱
通过 Mode
l、
Spe
c
t
rum 分析后,将各阶振型的最大反应组合叠加得到结构的最大反应 ．各阶振型效应可
应用规范规定的均方根法(
SRSS)组合得到 ．
通 过反应谱分析,得到的重力坝各阶位移和应力变化值见表２．其中第１、第６、第１２ 和第１５ 阶位移云图
７４
孙 颖,等:基于有限元的重力坝稳定及动力响应分析
０T(
ïì１＋１
βmax－１),
ï
ïβmax,
β＝ íï
９
æ T０ ö ０．
ï
÷÷
,
ïβmax çç
î
èT ø
０＜T≤０．
１,
０．
１＜T≤T０ ,
T０ ＜T,
(
６)
式中:
T 为体系自振周期;
T０ 为特征周期;
β 为 加 速 度 反 应 谱 标 准 值;
J]．世界地震工程,
２００７(
１):
１１０
Ｇ
１１４．
[
２]周朝,张子琴 ．混凝土重力坝的三维有限元分析[
J]．科技创新与应用,
２０１６(
３６):
２３
Ｇ
２４．
[
３]薛松 ．基于 ANSYS 的高混凝土重力坝有限元静动力分析[
D]．郑州:华北水利水电大学,
２０１８．
[
４]龚勋,辛欣,李萌,等 ．反应谱分析法在重力坝抗震性能分析中的应用研究[
５５６
０．
应用有限元分析软件分析混凝土重力拱坝在静载力作用下的应力应变
应用有限元分析软件分析混凝土重力拱坝在静载力作用下的应力应变摘要：应用有限元分析软件对碾压混凝土重力拱坝进行建模分析，计算得到坝体在静载力下作用下的应力和变形，分析坝体应力、变形及安全评价。
关键词：有限元拱坝应力应变中国分类号：TV311 文献标识码：A随着全国水电事业的大力发展，筑坝技术的不断完善，全国拱坝数量也在不断增加，随之而来的对拱坝静动力学分析及安全评价也十分重要。
本文采用有限元软件对南康河水电站碾压混凝土重力坝进行建模分析，分析坝体的应力情况及变形稳定。
1 有限元分析软件碾压混凝土重力拱坝坝体应力计算采用的是大型商用有限元软件。
ABAQUS是一套功能强大的模拟工程的有限元软件，其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。
2 项目简介南康河水电站位于云南省普洱市西盟县境内，位于南卡江上游的一级支流，以水力发电为主。
坝址控制流域面积1152.5km2，多年平均流量37.8m3/s，水库总库容3095.6万m3，为日调节水库。
工程为Ⅲ等工程，大坝分类为中型工程，主要建筑物为3级建筑物，其他次要建筑物为4级建筑物。
南康河水电站挡水坝为碾压混凝土重力拱坝，按100年一遇洪水（P=1%）设计，1000年一遇洪水（P=0.1%）校核，坝顶高程665.0m，坝底高程659.0m，最大坝高68.0m，坝顶宽度6.0m，拱冠梁底宽22.24m，拱冠梁厚高比0.327，拱顶中心线弧长180.623m，弧高比2.66，顶拱中心角91º。
拱坝基础坐落于弱风化带中下部岩体上，河床基岩主要为中厚层状石英片岩夹云母石英片岩，岩体完整性好，坝址区基岩主要为云母石英片岩、云母长石石英片岩间夹少量绢云片岩。
3 研究内容主要研究南康河水电站碾压混凝土重力拱坝在静载下的应力、变形情况，即研究基本荷载组合以及特殊荷载组合下三维大坝-地基系统的应力和变形，拱坝应力分析采用有限元法，并给出坝体上下游面的应力和位移分布。
基于有限元数值计算的重力坝稳定及应力变形分析
基于有限元数值计算的重力坝稳定及应力变形分析汪小军【摘要】为研究重力坝的稳定性以及其应力分布规律,本文采用有限元数值计算方法,在三种不同工况下,建立两种不同断面的(溢流段与非溢流段)三维实体模型计算,并对其结果进行分析总结.结果表明:坝体在三种工况作用下会出现一定的变形,其中水平位移随着工况的变化出现较大的差异,而竖向位移基本不变;水平和竖向位移出现最大的位置分别为坝顶和坝体坝踵处;应力分布规律基本相同,最大拉应力和最大压应力分别出现在坝踵处和廊道附近;在三种工况下,工况3的坝体模型安全系数最小,而相同工况下溢流段的坝体安全系数比非溢流段的坝体安全系数更小.【期刊名称】《吉林水利》【年(卷),期】2018(000)010【总页数】4页(P25-28)【关键词】重力坝;溢流段;位移;主应力;安全系数【作者】汪小军【作者单位】上饶市科信水利水电勘察设计咨询有限公司,江西上饶 344000【正文语种】中文【中图分类】TV642.31 引言重力坝由具有独特的优势，即简单的结构形式、广泛适用的地形、地质条件等[1，2]，遂使其在我国的水利水电工程建设中被认可，并成为主要的坝体之一[3，4]。
如今的社会已经不能单纯的关注工程形式，更多的关注其自身的安全性，因此针对重力坝更多的是关注坝体的应力、应变以及抗滑稳定性[5，6]。
随着水利工程的不断发展，以及社会各界对安全问题的深入探究可知，对于上述安全问题的主要研究手段包含有限元分析的方法、刚体极限平衡法、可靠度分析法及断裂力学的分析方法等[7，8，9]。
鉴于此，本文基于有限元数值分析方法，对某重力坝工程进行三维有限元数值计算，分析重力坝的位移及应力分布规律，并计算其安全稳定系数，为实际重力坝工程建设提供一定的理论基础。
2 计算模型2.1 模型参数本次模型采用某重力坝工程，将溢流段和非溢流段分开进行建模及分析，其中非溢流段、溢流段的主体材料主要为M10浆砌石和C15混凝土，根据文献资料并结合实际工程选取对应的材料参数，如表1所示，而坝体与坝基接触面的抗剪断摩擦系数与黏聚力分别为0.9和1 MPa。
反应谱分析法在重力坝抗震性能分析中的应用研究
反应谱分析法在重力坝抗震性能分析中的应用研究
龚勋;辛欣;李萌;邓源;李悠然
【期刊名称】《水利天地》
【年(卷),期】2017(003)010
【摘要】本文首先对坝体进行静力分析,采用ANSYS中PLANE182单元模拟某重力坝坝体,探究坝体在自重及水压力作用下的应力与变形;然后运用反应谱分析法对坝体进行地震动力分析,包括采用模态分析法提取模态,确定大坝自振频率及振型,结合反应谱曲线方程得出反应谱谱值,然后进行反应谱分析及模态扩展,最后对模态进行合并.求得坝体在各阶模态下的应力及位移响应.结果表明:在地震设防烈度为8°的情况下,研究坝体的变形及强度均在可控范围内;最大应力、应变主要集中于大坝的中上部、坝踵、坝趾等处,在今后设计和施工中应引起足够的重视.
【总页数】7页(P1-7)
【作者】龚勋;辛欣;李萌;邓源;李悠然
【作者单位】惠州市华禹水利水电工程勘测设计有限公司,广东惠州 516000;三峡大学水利与环境学院,湖北宜昌 443002;三峡大学水利与环境学院,湖北宜昌443002;三峡大学水利与环境学院,湖北宜昌 443002;三峡大学水利与环境学院,湖北宜昌 443002
【正文语种】中文
【中图分类】TV642.3
【相关文献】
1.新旧抗震规范下重力坝抗震安全性能分析 [J], 范书立;赵绍宇;陈健云;童伟
2.Ansys在混凝土重力坝振型分解反应谱分析中的应用 [J], 张传国
3.反应谱分析法在重力坝抗震性能分析中的应用研究 [J], 龚勋;辛欣;李萌;邓源;李悠然;
4.基于ANSYS的石桥水电站重力坝抗震性能分析 [J], 刘东升
5.基于三维黏弹性边界的混凝土重力坝抗震性能分析 [J], 李金友;成卫忠;李同春因版权原因，仅展示原文概要，查看原文内容请购买。
有限元方法与应用谱分析
获得模态解
建模 获得模态解： • 与通常的模态分析步骤相同 • 少量不同之处将在后面讨论
典型命令： /SOLU ANTYPE，MODAL
M5-21
获得模态解
• 模态的提取： – 有效的方法只有Block/PCG Lanczos，supernode或缩减法 – 提取足够的模态以包含频谱的频率范围 – 扩展所有的模态，只有扩展的模态才能用于频谱的求解 • 载荷和边界条件：对于基础激励，一定要约束适当的自由度 • 文件：.mode文件包含有特征向量，并且此文件要用于频谱求解
M5-8
定义及目的
• 谱分析的一种代替方法是进行瞬态分析，但是： 瞬态分析很难应用于随时间无规律变化载荷的分析; 在瞬态分析中，为了捕捉载荷，时间步长必须取得很小，因而费时且昂贵.
M5-9
目录
重力坝简介 谱分析的定义及目的 基本概念和术语的理解
如何进行响应谱分析
M5-10
术语及概念
主题包括： • 频谱的定义 • 响应谱如何用于计算结构对激励的响应： 参与系数 模态系数 模态组合
用的标准
M5-29
定义响应谱
模态组合法： • 在模态组合中，规定的有效门限值能使模态组合时仅仅包含主要模态，有效
门限值是某个模态的模态系数对最大模态系数的比率，为了在组合时包括所 有模态，要采用0值作为门限值 • 输出类型是计算不同的响应量，如位移、速度或加速度
典型命令： MCOMB，...
M5-30
M5-18
建模
模型： • 建模的注意事项与模态分析相同 • 仅考虑线性的单元及材料，忽略各种非线性 • 记住密度的输入，同时如果存在依赖于材料的阻尼，也必须在这一步中定义
M5-19
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
= X2 [ mci ] { Ui } + X2 [ mcc ] { Uc } - [ mci ] { Ui } - [ kcc ] { Uc }
( 7)
由此根据( 7) 式可求得第 j 阶振型对应的上下游边缘正应力。
E E Rj, y u =
Wj 6 M j B + B2
( 8)
E E Rj, y d =
LI Tongchun, HU Jigang, ZHAO Lanhao, ZHANG Liping ( College of Water Conservancy and Hydropower Engineering , Hohai University , Nanjing 210098)
Abstract: This paper presents an equivalent dynamic stress analysis for gravity dams based on the response spectrum method and finite element ( FE) analysis, and proposes a formula to directly calculate the internal forces and equivalent stresses for each vibration mode from the finite element method ( FEM) data of vibration modes. The equivalent stresses of four dams under construct ion or proposed are computed as test cases by the response spectrum method of cantilever model or by FEM, with the equivalent dynamic stresses compared for several typical locations. The results show that the equivalent stress obtained from the FE vibrat ion modes is better than that from cantilever model. The former can reasonably reflect the stress concentration due to the sudden change of dam prof ile. The proposed method of FE equivalent stress analysis is easy to implement . Key words: hydraulic structure; equivalent dynamic stress; response spectrum method; cantilever method; finite element method; gravity dam
第 28 卷 第 5 期 2009 年 10 月
水力发电学报 JOURNAL OF HYDROELECTRIC ENGINEERING
Vol. 28 No. 5 Oct. , 2009
基于反应谱法的重力坝有限元等效动应力响应分析
李同春, 胡继刚, 赵兰浩, 张丽萍
( 河海大学水利水电工程学院, 南京 210098)
文中以四个在建或拟建的重力坝为例, 比较了有限元和悬臂梁两种离散方法在求解重力坝结构自振特性及 反应谱法动等效应力响应方面的差异, 从而说明了采用有限元法求解结构自振特性, 再用反应谱法求解重力坝动 力响应的必要性。
1 反应谱法等效应力计算基本原理
111 悬臂梁法
在重力坝的悬臂梁法地震动力分析中, 通常只考虑坝体在水平地震作用下
和剪切变形产生的水平振动, 由结构动力学知识[ 4] , 柔度矩阵 D可以表达为:
j
E Dij = Dji =
dy k ( yj - yk ) ( y i - y k ) PEI + kPGB
( 1)
k= 1
式中: Dij 为结点 i 作用单位力时在结点j 引起的水平位移; dy k 为结点 k 控制的上
E 移响应 ui = Gi Bikg UiPX2i ; 振型惯性力 Fi = Mai ; 振型应力 R( j ) i =
Bj 2Ij
N
(
k= j+ 1
yk
-
yj ) F ( k) i ; 振型坝底剪力 Si =
N
E F( k) i。
k= 1
根据5水工建筑物抗震设计规范6的规定, 采用振型分解反应谱法计算地震作用效应时, 可由各阶振型的地震
20
水力发电学报
2009 年
的应力沿截面呈非线性分布, 在坝踵和坝趾及断面突变处有明显的应力集中现象。 为此, 本文提出了由有限元分析得到的振型求解振型内力及振型等效应力, 并根据反应谱法理论求解重力坝
等效动应力响应的分析方法。这样使得有限元反应谱和悬臂梁反应谱的应力计算方法是一致的, 其主要差异在 于是否对结构进行简化。
摘 要: 提出了由有限元法得到的振型求解振型内力及振型等效 应力, 并根据反应 谱法理论 求解重力 坝等效 动应力 响
应的 分析方法。本文对四个在建或拟建的重力坝, 用悬臂梁反应谱 法和有限元等效应力法分别进行了计算, 比较了坝 体
几个 典型部位等效动应力响应。结果表明, 在坝体体型有突变的地 方, 有限元等效应力法相较悬臂梁法所 得值有明显 的 增高 , 更敏锐的反映出了应力集中对坝体垂直动应力的影 响。本文 提出的方法简单易行, 比悬臂梁反应谱 法更能反映 坝
( 6)
由于方程( 3) 和( 5) 完全等价, 则由方程( 3) 解出结构振型后, 代入式( 6) 就可直接求解出给定截面上的约束内
力; 其表达式为:
{ f c } = - {f c} = [ kci ] { Ui } + [ kcc ] { Uc } - X2 [ mci ] { Ui } - X2 [ mcc ] { Uc }
Wj 6 M j B - B2
E E 式中: Rj, yu , Rj, yd 分别为上下游边缘正应力; Wj 为计算截面上与第 j 阶振型对应的轴向约束力; Mj 为计算
c Ñ, Ò
m m U ci
cc Ñ , Ò
c Ñ, Ò
f c Ñ, Ò
式中: 下标 c 代表指定截面上的结点; i 表示非指定截面上的结点; { f c } Ñ , { f c } Ò 分别代表指定截面上 Ñ 、Ò两个
子结构间的约束内力, 由作用力与反作用力的关系有:
{fc} Ñ + {fc} Ò = 0
作者认为在计算机内存、计算速度及有限元计算分析软件高速发展的当今, 已无必要再采用简化的悬臂梁法 作为重力坝自振特性的基本分析方法。但是目前的应力控制标准主要是针对材料力学方法制定的, 其主要特点 在于控制应力是根据截面内力按等效线性分布假定由材料力学公式计算得到。而常规的由有限元分析直接得到
收稿日期: 2009-04-21 基金项目: 国家自然科学基金委员会资助项目( 90510017) 作者简介: 李同春( 1963) ) , 男, 博士, 教授. E- mail: ltchhu@ 163. com
由弯曲变形和剪切变形产生的水平振动, 忽略坝体的竖向振动, 也忽略转动惯量 的影响[ 3] 。如图 1 所示, 将重力坝简化为一固结于地基的悬臂梁, 取单位厚度重
力坝坝段为分析对象, 将整个坝段沿高度按任意间距划分成 N - 1 段, 结点号自
下至上为 1 至 N, 结点 1 固定于地基上。取 y 轴向上, 当只考虑坝体由弯曲变形
别为第 i 阶、第 j 阶振型的圆频率。
112 有限元法
有限元反应谱法的分析思路和悬臂梁法基本相同, 其基本方程形式同( 3) , 主要差别在于沿上下游方向有限
元的单元数量可以任意选择。为了和悬臂梁法得到的结果具有可比性, 本文根据有限元内力法[ 5] 的基本原理, 利
用有限元得到的振型, 从有限元的静力平衡原理出发, 直接求解得到不同振型下给定截面上对应的有限元内力。
作用效应按平方和方根组合。当两个振型的频率差的绝对值与其中一个较小的频率之比小于 011 时, 地震作用
效应宜采用完全二次型方根组合:
mm
E E SE =
Qij S iS j
ij
( 4)
Qij =
( 1-
8
Ni Nj ( Ni &C2X ) 2 + 4NiNj CX ( 1+ C2X) + 4( N2i +
首先用该指定截面将整个有限元模型分为坝体( 记为 Ñ) 与地基( 记为 Ò) 两个子结构, 并对这两个子结构在
该指定截面上施加一组大小相等、方向相反的约束内力, 则可建立与方程( 3) 完全等价的另一组方程:
k ii kic
Ui
- X2 mii mic
Ui
0
=
( 5)
k k U ci
cc Ñ , Ò
体体型的影响。
关键词: 水工结构; 等效动应力; 反应谱法; 悬臂梁法; 有限元法; 重力坝
中图分类号: TV312
文献标识码: A
FE equivalent dynamic stress analysis based on response spectrum method for gravity dams
0 前言
水工建筑物抗震设计规范[1] ( 以下简称规范) 规定, 重力坝分析应以同时计算弯曲和剪切变形的动、静材料力 学方法为基本分析方法, 即悬臂梁法[ 2] 。而随着数值计算技术的发展, 有限单元法作为一种比较成熟的计算离散 方法, 在结构数值计算领域已得到广泛的应用, 能够很方便地模拟重力坝的施工过程、材料分区、孔洞、坝体- 库水 - 地基系统的相互作用以及复杂的地质条件和边界条件, 具有材料力学方法无可比拟的优势。因此, 规范也同时 规定, 对于工程抗震设防类别为甲类, 或结构复杂, 或地基条件复杂的重力坝, 宜补充作有限元法动力分析。目前 规范规定重力坝的动力分析方法主要采用振型分解反应谱法。