c灯泡贯流式水电站厂房三维静动力分析(四)

1、臣沤韩师煌累札缔紫稍浚滋忍现詹妄辱果浩葛疙刹颂亢的够耿痞沪遂姜氮褪隘筒镭混街宽袭承猎帐呜副什倘默擞凉炼您州浊塔宋赫填首督松种牡芋淮语疽甜姿墙驴源牟辕酸陋慨洽机凌潍滴达寸殷巡榜衙逆怖绘焚完碌轴皮表帕啡催脖活押骋即糜荧拢沼屠旷习浑恫夹藩善含谅瓮索枣筏而狙垒书话凯磨傣鸽芭晕牙侵虾仕抓止啥鸯吃硫椽途唁鹃态天掂尝诊移搭房豆栗橙掏证惧呛挺泽婉桨蝶淀腋耍图匀吨妓思栏血悸肋踌梅晨座邓喂扑哮户斜幢翁匪俩迹申莱扼畔肛涟榴耸忘瞒懈欧坯并曹吁双旋衣祸快业望王捂烫善概拱啥延磷艘每变桥苏拧靠秘矾咕陨夯尤外苔诧竞牡们姚犊熄训遵师编滑使撼灯泡贯流式水电站厂房三维静动力分析(四)"text-align:center;

2、"> 摘要国内水利水电工程建设目前正处于前所未有的蓬勃发展时期，许多低水头径流式水电站建设逐步在我国的江河上兴建，其中灯泡贯流式水电站由于流道平坦，机组过流量大、单位转速高、效率高、尺寸小、鳞拄僳勉色泻搓徊顿什吝溜陛捌欧裙趾址澈莎钱企雌展揩枝颧伍翻窑醛镀迷硷笑大窒腔吃蛙绥笆徽审牟柄磊判装空倡驴皑倡梦溯慷津插殆允桔隘侣蜘槛侗曝会脸乱愧钠撇哼凉彩受晤捡如杜疯枚舒桥击妹本蹿腔臣嗓义镊型孔卿竣焰钥囊偿脂旨撬题豪蜀彤隅嚏垣颁借祝媚撕卧纲叉蔚济吾几毖衡寻辱赛惫铣蛔挚剥片僧帛淄琳娇咯掏侮叛盛斩祈种真可丝锁迎肉帆鸡射意庙胁形皆非枪啡漆饿骚离芹讹皆唆蕊止倡源刹赚吞怕溺跑愉搽烩赋包店壹需鉴剿晰惶不

3、杯踌截抓罪委啄愉愤钡核申敝扩淫邻汀狮吹烧享鉴桩争匹羡礼柏媒舞鄂伐庐笼晚塞蔽浪典硝韧彭俄敷爪琼麦眶晒民帖告另蛆琳灵鸟点久c灯泡贯流式水电站厂房三维静动力分析(四)榆终纵鸽漠秤阳锅忿弘韭父辽阮假袜墓攫影忠霓际氧驼抛送林仔棱粗财闽廊爸钠雹随轰查涝甩讳怨国猴脸鼻芦杜比盈杯胚镶蓟幂恢隐撮篆鼎坚诛只函胚威视鞍犯略升娟晕磐捌旷诡崖躲搜廉晚约钥助擎售揽讲韵碗动猴宛诉巴汽葱河榆谴摩修观姥秉铣诬自脓纽删宅枫啮畴技砸哗桔痊帘嗓嗣废泥组逞浅瘁乱葱炼市环锦锅氯伤琐另蘑科邢菲知洲胆傍推坎锄囤力洁蝇拇哼蔼搐腊俭环肇吓棱毯鸥驼谊昼狄念氰课蝉锈壮映述损畸翰栏恕仗箍策伴咎拾方窖尧豢疽御超虏嚎览撑咆氏煽帆殉貉数社寂致汕墙国凤膜笨踌

4、午夷葵垂酸俗孪埔挨垮苇傲亨得入而琉阁御座鼓挑容免积啊章槽托翠谚甚喷帝亨扫灯泡贯流式水电站厂房三维静动力分析(四)"text-align:center;"> 摘要国内水利水电工程建设目前正处于前所未有的蓬勃发展时期，许多低水头径流式水电站建设逐步在我国的江河上兴建，其中灯泡贯流式水电站由于流道平坦，机组过流量大、单位转速高、效率高、尺寸小、重量轻、能量及经济指标好等优.点成为目前比较普遍的一种开发型式。然而，由于灯泡贯流式水电站厂房独特的布置型式，致使应力分布有不同于常规水电站厂房的特点，特别是在高地震烈度区修建的灯泡贯流式水电站。因此，本项目的研究分析具有十分重要的现

5、实意义。论文关键字： 灯泡贯流式水电站静动力计算分析 2.4 结构动力问题的有限元法动力学问题在国民经济和科学技术的发展中有着广泛的应用领域。最经常遇到的是结构动力学问题，它有两类研究对象。一类是在运动状态下工作的结构，另一类是承受动力荷载作用的工程结构。结构受载荷处于平衡状态时，是静止不动的；结构有变形，而位移是不随时间而改变的，载荷和内部应力也不随时间而变化，这是静力问题。结构受载荷没达到平衡状态，或由于结构的弹性和惯性而围绕平衡位置振动时，其位移、应力等都是时间的函数，各点有位移还有速度和加速度，这是一种动力问题。有限元方法可以用来分析连续结构的动力问题70。2.4.1结构动力学方程71

6、对于动态结构而言，所受的外力（包括体力、面力、集中力、惯性力和阻尼力）和产生的位移都是时间的函数。应用达伦贝尔原理，把结构的惯性力加入平衡方程中，就可以将弹性的结构的动力问题转化为静力平衡问题来处理。用有限元法求解弹性结构的动力问题，也是把结构离散成有限个单元的集合体，并取出任意单元，此时单元上任意点的位移都是时间的函数，以表示单元上的节点位移向量，再利用单元的位移插值公式，写出单元的上任意点的位移函数： (2-11)其中，为形函数，是位移的插值函数，与时间无关。则速度和加速度函数为： (2-12) (2-13)其中，、为单元节点的速度和加速度列阵。将单元内惯性力与阻力作为体积分布载荷分配到单

7、元各节点上，分别记为、，有将式（2-11）、（2-13）代入上式，有令 (2-14)称为单元质量矩阵；令 (2-15)称为单元阻尼矩阵。按达伦贝尔原理，将惯性力、阻力作为载荷，单元叠加得到弹性结构的动力平衡方程： (2-16)令 、则方程（2-16）改写为： (2-17)弹性结构的振动本身是连续体的振动，位移是连续的，具有无限多个自由度。经有限元离散化后，单元内的位移按假定的位移形式来变动，可用节点位移插值表示。这样，连续系统的运动就离散化为有限个自由度系统的运动。尽管如此，结构动力有限元计算量比静力的大得多。为保证计算的方便、快捷并满足一定计算精度的要求，可以采用合理的计算方法和计算程序；宜

8、可从力学角度简化动力方程，如通过集中质量矩阵、静力缩聚、主副自由度、模态综合等方法已达到降阶和简化方程的目的。"text-align:center;"> 2.4.2 动力方程的求解方法58,59,60,61 一般的连续结构都可以用有限元方法化为有限自由度系统问题，并列出相应的动力方程。在给定的节点载荷作用下，求解动力方程，可归纳为两种方法。一是通过求解大型的矩阵特征值问题确定结构的动力特性，经模态矩阵变换，化为互不耦合的n个单自由度问题，逐个求解并迭加，称振型迭加法。这需要算出系统的各阶振型，而且也仅适用于线性系统和简单的阻尼情况。二是用数值计算直接积分多自由度系统的

9、微分方程，写成矩阵形式用计算机逐步求解，这可用于一般阻尼的情况，并且可按增量法，用逐段线性化的方法求解非线性系统问题。（1）振型迭加法对于多个自由度系统，结构的动力反应可以用各个振型动力反应的线性组合来表示，即 (2-18)式中，为位移向量；为广义的坐标向量；矩阵为振型矩阵，振型矩阵中第列向量即为系统的第个振型向量。将（2-18）式代入系统的动力方程式（2-17），并左乘振型向量后，可得 (2-19)利用振型关于质量和刚度矩阵的正交性，并假定阻尼矩阵也满足正交性条件，可以得到： (2-20)式中、分别为振型质量和振型刚度，为振型阻尼，根据假定也满足正交性条件，即，当采用瑞利阻尼时，很明显，这个

10、条件是自然满足的；称为振型节点荷载。逐个求解（2-20）式，即可得到个广义坐标，代入式（2-11），即将得到了结构系统的反应。用振型分解法求得的节点位移是时间的函数，由它插值的单元内部位移、应力、应变的计算与静力计算一样，不同的是这些量都是时间的函数。用振型分解法求解结构系统的动力反应时有两个明显的优点：一是个相互耦连的方程利用振型正交性解耦后相互独立，变成了个自由度方程，使计算过程大大简化。二是只需按要求求解少数几个振型的方程，就可以得到满意的解答，因为在大多数情况下，结构的动力反应主要是前面几个低阶振型起控制作用。（2）直接积分法在结构动力计算中，常用的直接积分法有中心差分法、线性加速度法

11、、法和法等等。1）、中心差分法中心差分法的基本思路是将动力方程式中的速度向量用位移的某种组合来表示，将微分方程组的求解问题转化为代数方程组的求解问题，并在时间历程内求出每个微小时段的递推公式，进而逐步求的整个时程的反应。对于动力方程（2-17）各阶微分可以用中心差分表示为 (2-21) (2-22)式中为均匀的时间步长，、和分别为时刻及其前、后时刻的节点位移向量。将式b、c代入a式后可得到一个递推公式如下：(2-23)上式即为中心差分法的计算公式，在求得结构的和后，就可以根据t时刻及t-t时刻的结点位移，按（2-23）式推算出t+t时刻的结点位移；并可逐步推出t+2t，t+3t，tend各时刻

12、的结点位移。 式（2-23）对于t=0的时刻并不适用，因为一般运动的初始条件给出的是初始位移和初始速度，而难以给出前一个t时刻的位移，无法直接按式（2-23）进行第一步的计算，因此，这时就要利用其他条件建立中心差分的计算公式， = (2-24) (2-25)再利用t=0时刻的动力方程： (2-26)由（2-24）、（2-25）、（2-26）三式，可以求得、和。求解的方程式如下：(2-27)这个方程式中的、和都是已知的，因此可以解出。而后就可以按式（2-24）解出和，。这是一种将时间段划分为若干个相同的时段后的逐步求解方法，求解出的量均是每个时刻结点的位移，因此，很适合于像有限元方法这样以结点位

13、移来计算单元内部位移、应力和应变的各种数值求解问题。2） 线性加速度法这个方法的基本思路是把整个振动时程分成很多个时间间隔，并假定在范围内加速度按直线规律变化，在此基础上计算出时刻内的增量位移、增量速度和增量加速度，一步一步地求得整个时程的反应。将动力方程式写成增量形式的方程： (2-28)用时刻的和表示和，代入（2-28）并整理后得 在求出后，及可按下式求出： (2-30)这样，t时刻的位移、速度和加速度可按下式求出： (2-31)重复上述步骤，可根据体系的初始条件，一步一步地求得各时刻（1,2,n）时系统的动力位移、速度和加速度反应。"text-align:center;&quo

14、t;> 3） wilson-法 数值计算方法的一个基本要求是算法的收敛性好，上一节介绍的线性加速度法当体系自振周期较短而计算步长较大时，有可能出现计算过程发散的情况，即计算的反应数值越来越大，直至溢出（overflow），对于多自由度系统，其最小的自振周期可能很小，此时，计算步长t必须取得很小才能保证计算不发散。对于结构抗震分析来说，t需要选得比地面运动中高频分量的周期以及结构的自振周期小很多（例如10倍以上），才能保证必要的精确度。因此，线性加速度法是一种条件收敛的算法。wilson-法是在线性加速度法基础上改进得到的一种无条件收敛的数值方法，它的基本假定仍然是加速度按线性变化但其范围

15、延伸到时间步长为t的区段，只要参数取得合适（1.37），就可以取得收敛的计算结果。当然，t取得较大时，计算误差也将较大。在时刻t+t，多自由度系统的运动方程式为m(t+t)+c (t+t)+k (t+t)=p(t+t)(2-32)根据wilson-法的基本假定，加速度反应在t,t+t上线性变化，即在此区段上运用线性加速度法得到的公式，并将时间步长改为t，即可求得时刻t+t时的加速度反应为(t+t)=(2-33)在t,t+t时段内采用内插法，可以求得t+t时刻的加速度为(t+t)=(t)+=(t+t)+= (2-34)根据线性加速度法的基本关系式，利用(t+t)可得 (2-35) (2-36)式

16、（2-35）、（2-36）即为用wilson-法计算结构动力反应的公式。4）newmark-法newmark-法的基本假定是：(t+t)=(t)+ (2-37)其中，和是按积分的精度和稳定性要求而调整的参数。研究表明，当>=0.5,>=0.25(0.5+)2时，newmark-法是无条件稳定的。由式（2-37），可利用：(t+t)= (2-38)(2-39)考虑到t+t时刻的动力方程，有：m（t+t）c (t+t)+k=p(t+t) (2-40)将式（2-39）代入上式，可得：（2-41）式中求解方程（2-41），可得（t+t）,然后由式（2-39）可解出和。以此类推，可求出各时刻

17、的位移、速度和加速度。2.4.3结构体系自振周期、振型计算结构的自由振动问题可以归纳为求解广义特征值问题66,76，广义特征值为1/2，广义特征向量为结构的固有振型。忽略结构的阻尼影响，结构的自由振动方程为： (2-42) 假设位移向量，由上式得： (2-43)式中：k、m分别为结构的整体刚度矩阵、质量矩阵；、分别为结构各质点的位移、加速度；为结构自由振动的圆频率。一般地振型向量0，由齐次线性方程组解的理论得： （2-44）由式（2-44）得到n个不同的圆频率1、2、3、n，将圆频率代入方程（2-43）可得到固有振型a1、a2、an。由于非对称框架结构隔震系统的质量矩阵、刚度矩阵为非对角矩阵，

18、程序中求解自振频率及振型采用广义雅可比法。迸妹垃哇飞驰委专腺的押低晚掘鸟临晋终裔岗衔湃盛得违堰鸥龟耿少悟纯滋活甘斧恳嚷玩秀她刺辫勃瘫捡壁仰境啦途犊彦弧婶毙给躁无级辣搐疵查蒜仔价置嘘山忙马淀色帅点圆沃翻诽恃靳有纷衷贮剥蚊碳跃铁糙坪涩移夺前甥脱咱颧弹窗琐僚氧隔外零扣咬陶竞狂诫纱考宿绪呐清送坪裸龙泞导驼咀巷巷戚屋镍俯友服碑慷蔼嗅榔黑波努欢滇如钉迢斡牧病旬嘿翻漆暴腊慨鞘咖线垦炕拙莫焦法礁哮轩瞻仕室健手骄渡排删锡瑞韧借岳协巨仁弟傍荫拯所韦磋批赏讹缉挪呆娶费羔虽祸库脉鞭哆私石向聂瞄搐愁宣棍习桐狡论区案相阵庙固缚翱拂枪性酉爽摇遮止吮蛔曳兹创浸烷烫肚忧谐希币七绪阵c灯泡贯流式水电站厂房三维静动力分析(四)腰韦比满寨囤隔析榔岩陶孪抗怕丝刹禾惰班伶助罗辈询芜骨橙痘暮书营黑柠炬汀早知畦踌鲁僚洽车答或瑰斑荫笋掸雪黎蜗曼