拱坝等效应力计算pdf

　　第五章 有限元法的拱坝应力控制标准的探讨 41 第五章 有限元法的拱坝应力控制标准的探讨 5.1 概述 目前，有限元法已经广泛应用于拱坝的应力分析，它在计算理论上比多拱梁法先进， 而且适合各种复杂的体形和结构布置形式，可以比较合理的考虑拱坝的整体作用和复杂 地基的影响，同时也可以考虑大孔口、复杂基础、重力支墩、不规则外形等多种因素的 影响，并可以进行仿真计算。但是，在弹性理论的假定下，有限元应力成果往往与所采 用的单元类型和单元网格划分的尺寸有着密切的关系，特别是在建基面附近存在着严重 的应力集中现象，使得有限元应力成果的数值稳定性较差，难以建立相应的应力控制标 准。实际工程中，由于岩体内存在着大小不等的各种裂隙，应力集中现象将有所缓和， 不一定像计算结果那么严重，计算中如何考虑应力集中，是有限元法应用于拱坝的“瓶 颈” [29][30][68] 。我国一些学者提出了有限元等效应力法克服了这一“瓶颈”，新编的《混 [25] 已经正式规定在拱坝的设计中采用有限元等效应 凝土拱坝设计规范》（SL282—2003 ） 力法，为有限元法在拱坝设计中的应用开拓了良好的前景，由于有限元法的强大计算功 能，它将逐步取代多拱梁法，成为拱坝设计的主要方法。 为了解决拱坝的弹性有限元应力成果集中现象，人们采取了许多措施，如规定应力 点的位置、采用厚壳单元、利用有限元计算出坝体应力沿面积积分，得到截面内力，然 后用材料力学方法计算结点应力，即所谓的有限元等效应力法。 傅作新教授[30]提出的有限元等效应力法是基于有限元法的分析结果，将有限元所求 得的应力合成为截面内力，然后求出对应的线性化应力。该方法为有限元分析结果规范 化提出了一种新思路，但是该法截面内力是拟合出来的，不能精确满足内力平衡条件。 此外，该方法要求沿径向网格必须布置3层以上的单元格才能求得结果。 李同春教授[69]提出的基于有限元内力法的等效应力算法和改进的拱坝等效应力分 析方法[70]，该方法假定拱和梁向的正应力与剪应力在拱和梁的截面上的面单元内呈双线 性分布，根据静力平衡原理导出以结点应力为未知量，截面约束内力为右端顶的平衡方 程。在此基础上导出根据截面约束内力求解等效结点应力的公式，然后根据截面约束内 力求解等效结点应力的公式，并在此基础上导出拱和梁向应力为直线分布的上下游面等 效应力的求解方程和主应力的求解公式，但是计算模型复杂。 针对目前拱坝内力计算流行着以二次曲线表示应力分布，即σ a +bx +cx 2 ，朱伯 芳院士在文献[68]指出：从有限元应力计算内力时，以直接采用数值积分方法为宜，利用 三点应力决定二次曲线再求内力的方法是不妥当的，会带来不必要的误差。与此同时， 新修订的水利行业《混凝土拱坝设计规范》（SL282—2003 ）在有限元等效应力法的相关 条款中尽管明确了等效应力的概念和原理，但是还存在以下的问题[71] ：（1）没有明确有 限元计算的模型和内力计算公式，不同的模型和算法可能导致很大的差异；（ 2 ）没有分 42 基于有限单元法的高拱坝体形优化设计 清需要作等效线性化处理的主要应力，笼统地规定由截面内力按材料力学公式计算应 力；（ 3 ）拉应力控制标准的制定原则不够科学，规定的容许拉应力仍无法适合200m 以上 的高拱坝。钱向东教授引入压力容器“分析设计法”中应力分类的概念，进一步阐述有 限元等效应力法的原理和有关计算方法，通过对若干国内已建拱坝的分析，提出基于有 限元等效应力法的拱坝强度设计准则。 5.2 有限元等效应力的计算 5.2.1 截面内力的计算 从力学的观点分析，双曲拱坝是一种具有不规则边界的变厚度壳体，在弹性壳体理 论中，一般取距上下游面等距离的中间曲面（简称中面）为研究对象。用有限元法计算 ( ′ ′ ′) ′ x ,y ,z { } 拱坝，得到整体坐标系 中的应力 σ ，通过拱轴线上任意结点 取局部坐标系 i (x, y, z)，如图 5-1 所示，其中 x 轴平行于拱轴的切线方向，y 轴平行于半径方向，z 轴 { } 为铅直方向，局部坐标系中的应力 σ 。 图 5-1 应力计算坐标系 Fig.5-1 Stress calculated coordinates 则局部坐标系下的应力由下式计算： {} [ ]{ }′ σ T σ （5-1 ） 其中：[T]为坐标转化矩阵。 ⎡ l 2 m2 n2 2l m 2m n 2l n ⎤ ⎢ 1 1 1 1 1 1 1 1 1 ⎥ l 2 m2 n2 2l m 2m n 2l n ⎢ 2 2 2 2 2 2 2 2 2 ⎥ ⎢ l 2 m2 n2 2l m 2m n 2l n ⎥ [ ] 3 3 3 3 3 3 3 3 3 （5-2 ） T ⎢ ⎥ ⎢l l m m n n l m +l m m n +m n l n +l n ⎥ 1 2 1 2 1 2 1 2 2 1 1 2 2 1 1 2 2 1 ⎢l l m m n n l m +l m m n +m n l n +l n ⎥ ⎢ 2 3 2 3 2 3 2 3 3 2 2 3 3 2 2 3 3 2 ⎥ ⎢l l m m n n l m +l m m n +m n l n +l n ⎥ ⎣ 1 3 1 3 1 3 1 3 3 1 1 3 3 1 1 3 3 1 ⎦ ( ) 式中， l ,m ,n 为坐标轴 x 、y 、z 的方向余弦。 i i i ( ) 假设在局部坐标系中的三个主应力分量 σ ,σ ,τ 沿坝体厚度线性分布。根据拱坝 z x zx 第五章 有限元法的拱坝应力控制标准的探讨 43 应力分析的假设，在用弹性有限元法求的拱坝应力后，便可以用等效应力求出拱坝上下 游坝面的三个应力分量，然后由该处的微分体的平衡条件求得另外三个平行于坝面的应 力分量，就能进一步取出拱坝上下游坝面的主应力。 梁的水平截面在拱中心线上取单位宽度，在y点的宽度为1+y ，r 为中心线] 径，沿厚度方向对梁的应力及其矩进行数值积分，得到梁的截面内力如下 ： 梁的水平截面竖向力为： t 2 ⎛ y ⎞ W − σ ⎜1+ dy （5-3 ） b ∫−t x ⎟ 2 ⎝ r ⎠ 梁的弯矩为： t M − 2 (y −y )σ ⎛⎜1+y ⎞dy （5-4 ） b ∫−t 0 z ⎟ 2 ⎝ r ⎠ 梁的切向剪力为： t 2 ⎛ y ⎞ Q − τ ⎜1+ dy （5-5 ） b ∫−t zx ⎟ 2 ⎝ r ⎠ 梁的径向剪力为： t 2 ⎛ y ⎞ V − τ 1+ dy （5-6 ） b ∫−t zx ⎜ ⎟ 2 ⎝ r ⎠ 梁的扭矩为： t − 2 ( − )⎛ +y ⎞ M τ y y ⎜1 dy （5-7 ） b ∫−t zx 0 ⎟ 2 ⎝ r ⎠ 式中， y0 为梁的截面形心坐标； 单位高度水平拱圈的径向截面宽度为 1，沿厚度方向对拱应力及距积分得到拱的内 力如下： 拱的水平推力为： t 2 H − σ dy （5-8 ） a ∫−t x 2 拱的弯矩为： t 2 M − σ ydy （5-9 ） a ∫−t x 2 拱的径向剪力为： t V 2 τ dy （5-10 ） a ∫−t xy 2 5.2.2 坝体上下游坝面应力的计算 通过数值积分计算得到拱坝的内力后，就可以通过材料力学法[2][72][73]计算坝体上下 游三个应力分量，然后通过微分体平衡条件计算上下游的其它三个应力分量； 44 基于有限单元法的高拱坝体形优化设计 （1） 坝体下游面的应力计算 水平面上悬臂梁的铅直正应力为： W M σzd b − b (t −lg) （5-11 ） A I b b 径向平面上拱的水平正应力为： H M t a a σ − × （5-12 ） zd A I 2 a a 作用在水平面上悬臂梁的切向水平剪力为： Qb M b ( ) τzxd τxzd − + t −lg （5-13 ） A I b b 其他三个应力分量为： τ τ (σ −P )tan η −τ tan φ （5-14 ） xyd yxd xd d d xzd d τ τ (σ −P )tan φ −τ tan η （5-15 ） yzd zyd zd d d xzd d σ P +τ tan η +τ tan φ （5-16 ） yd d xyd d yzd d 式中，ϕ 为径向铅直平面内坝面与铅直线的夹角，η 为水平面内坝面与拱中心线 d d 切线的夹角。 平行于坝面的三个应力为： ′ 2 ′ 2 ′ − σ σ sec φ P tan φd zd zd d （5-17 ） 2 ′ 2 ′ + − σ (1 tan φ ) P tan φ zd d d d ′ 2 2 + + σ σ cos η σ sin η 2τ sin η cos η xd xd d yd d xyd d d ( + + 2 ) 1 （5-18 ） σ 2τ tan η σ tan η xd xyd d yd d 2 1+tan ηd τ′ τ′ (τ cos η −τ sin η sec φ) ′ （5-19 ） xzd zxd zcd d zyd d d 通过以上三个应力分量就可以确定下游面的等效主应力为： 2 ′ ′ ′ ′ + − σ σ ⎛σ σ ⎞ 2 zd xd zd xd ( )′ σpd ± ⎜ ⎟ + τxzd （5-20 ） 2 ⎝ 2 ⎠ 式中 tanφ tanφ cosη′ （5-21 ） d d d ⎛ R ⎛R ⎞2 ⎞ R t 3 ⎜⎜1+4 Rd +⎜⎜Rd ⎟⎟ ⎟⎟ u u ⎝ u ⎠ I b ⎜ ⎟ （5-22 ） 36r ⎜ 1+Rd ⎟ ⎜ R ⎟ ⎝ u ⎠ t R +2R lg × u d （5-23 ） 3 R +R u d 第五章 有限元法的拱坝应力控制标准的探讨 45 3 t I （5-24 ） a 12 t ⎛R −R ⎞ y 0 ⎜⎜ u d ⎟⎟ （5-25 ） 6 R +R ⎝ u d ⎠ t t 其中：Ab t ；Aa t ；Ru r + ；Rd r − 2 2 ( ′ ′ ) 当 σ −σ 0 时，上式根式用正号；反之用负号。 zd xd （2 ） 上游坝面的应力计算为 上游坝面的应力计算和下游一样，主要计算公式如下： W M b b σ + lg （5-26 ） zu A I b b H M t a a σ + × （5-27 ） zu A I 2 a a Qb M b τ τ − − lg （5-28 ） zxu xzu A I b b ( ) τ τ −σ −P tan η −τ tan φ （5-29 ） xyu yxu xu u u xzu u ( ) τ τ −σ −P tan φ −τ tan η （5-30 ） yzu zyu zu u u xzu u σ P −τ tan η −τ tan φ （5-31 ） yu u xyu u yzu u ′ 2 ′ 2 ′ − σ σ sec φ P tan φ zu zu u u （5-32 ） 2 ′ 2 ′ + − σ (1 tan φ ) P tan φ zu

　　2、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。

　　3、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。

　　4、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

　　2014年高考新课标Ⅱ理科数学试题及答案(精校版,解析版,word版).doc

　　T∕CECS G-D66-01-2019 公路波纹钢涵洞技术规程.pdf

　　（高清版）DB3301∕T 0209-2018 农村生活垃圾处理技术管理规范.pdf

　　原创力文档创建于2008年，本站为文档C2C交易模式，即用户上传的文档直接分享给其他用户（可下载、阅读），本站只是中间服务平台，本站所有文档下载所得的收益归上传人所有。原创力文档是网络服务平台方，若您的权利被侵害，请发链接和相关诉求至 电线) ，上传者