框剪结构结构振动分析毕业论文

框剪结构结构振动分析毕业论文二 框剪结构结构振动分析计算模型2.1 运动方程的建立多自由度结构的运动方程为：{})(}]{[}]{[}]{[t F x k x c x M =++ （2.1）式中，[]M －质量矩阵； []c －阻尼矩阵； []k －刚度矩阵；{x}-速度 {x}-加速度 对于梁，柱线刚度比值大于5的强柱弱梁多层框架结构(如图1所示)，在动力分析中采用层间剪切模型。

层间剪切模型是把横梁视为刚性梁，并假定结构质量集中在各层楼面及屋面处。

当横梁发生水平震动时，不考虑柱的轴向变形，横梁只产生水平位移而不产生转角。

一般情况是把每层楼面积屋面简化为一个质点，将剪切型多层框架定一步简化为质点体系（如图2所示），这些质点只能有水平线位移。

下图中i m 为第i 层集中质量，i k 为第i 层柱上，下端错动单位位移所需的水平力，成为第i 层的层间刚度。

它等于该层所有柱的侧位移刚度之和，即∑=)(312i i ci h EI k （2.2）式中，c EI 为柱的抗弯刚度；i h 为第i 层层高。

`图2.1 多层剪切型框架计算图示图2.2 简化计算图示图 2.3 楼层变形图图2.4 简化后计算图示根据刚度系数的定义，容易得到剪切性框架结构的侧位移刚度为：1,11,,1_1_,1=-==-==+=++++iriiii iiiii iiiiikkkkkkkkkk（2.3）于是得到剪切型结构的刚度矩阵为[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----+--+=--nnnnnnkkkkkkkkkkkkkk113322221（2.4）剪切型多层框架的质量矩阵式对角阵：[]⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡=nmmmmM321（2.5）][][][K M C βα+=2.2剪切多层框架结构的自振特性考虑无阻尼系统的自由振动：}0{}]{[}]{[=+xk x M 应用雅克比法，其矩阵特征值问题为：[]{}[]{}020∆=∆M K ω （2.6） 首先要将矩阵化为同阶对称矩阵的标准值问题。

（1）将结构质量阵[M]进行乔列斯基（Cholesky ）分解质量矩阵[M]是对称正定矩阵，根据线性代数理论，一个对称的正定实矩阵总可以分解为一个下三角矩阵[L]及其转置矩阵[]TL 的乘积，即[ M]= [L] ×[]T L （2.7） （2）形成对称矩阵[A]。

将上式带入上上式得，[]{}[][]{}020∆=∆T L L K ω （2.8a ）或 [][]()[]{}[][]{}021∆=∆-TTT L L L L K ω （2.8b ）将[]1-L 左乘上式，并注意到[]()1-TL =[]()TL 1-，则有[][][]()[]{}[]{}02011∆=∆--T T TL L L K L ω （2.9） 令 2ωλ= （2.10）{}[]{}0∆=T L X （2.11）[][][][]()TL K L A 11--= （2.12）则式（a ）变为[]{}{}X X A λ= （2.13）因为[K]为对称矩阵，由上上式可见[A]=[]TA ,即[A]为对称矩阵。

这样就可以利用雅克比法求的上式中对称矩阵[A]的特征值和特征向量。

将结构质量矩阵[M]进行乔列斯基分解。

式（3.21）可以写为：⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡nn n n n n n m m m m m m m m m m m m21332312222111211=⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡nn n n n nn n n n l l l l l l l l l l l l l l l l l l l l 333223221131211321333231222111（2.14）根据矩阵乘法可得1111121111,l l m l m i ==（i=2，3，·，n ） （2.15） jj ij jk j k ik jj ij j i j i ij l l l l l l l l l l m +=+++=∑-=112211 （ 3,2,1.=j i ） （2.16）于是可得计算[L]下三角元素的公式：1111111,i i i l m l m l ==（i=2，3，·n ） （2.17） ()n j l m l j k jk jj jj ,3,2112=-=∑-= （2.18）jjj k jkik ij ij l l l m l ∑-=-=11(n i ,3,2=)(1,3,2-=j j ) (2.19)由于现在是对角阵，处理起来比较简单。

由[]{}[]{}020∆=∆M K ω可得 (2.20) [][][]{}[]{}021∆=∆-M M M K ω（b ）其中[]⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎣⎡=n m m m M0021 (2.30)以[]1-M 前乘（b ）式的两边得[]1-M [][][]{}[]{}021∆=∆-M M MK ω（d ）令 2ωλ={}[]{}0∆=M X（2.31） [][][][]11--=M K MA（2.32）则式（d ）可写成[]{}{}X X A λ= （2.33）⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎢⎣⎡---+--+=--m k m m k m m k m m k m k k m m k m m k mk k A nnn n n n m 113231212122121210][（2.34）2.3剪切型多层框架的水平地震作用单自由度体系的反应谱理论可得，但自由度体系的地震作用为：)()(2t x m t F ω= （2.35）根据反应谱理论，但自由度体系的最大地震作用为：G F α=（a ） （2.36）对于多自由度体系，第j 振型第I 质点的地震作用为：)()(2t x m t F ij j i ji ω= （2.37）由振兴分解法可得，)()()(t X t q X t x j j ji j ji ij ∆==γ （2.38）将其带入式（a ）的)()(2t m X t F j j i ji j ji ∆=ωγ （2.39）利用单自由度反应谱的概念，得到第j 振型第I 质点的最大最大地震作用为：i ji j j j j ji j ji G X G X F γααγ== （2.40）式中，)(t F ji ——j 振型I 质点的水平地震作用；j α——与第j 振型自振周期相应的地震影响系数；i G ——集中于质点I 的重力荷载代表值； ji X ——j 振型I 质点的水平相对位移；j γ——j 振型的参与系数。

由上式即可求得各阶振型下各个质点的最大水平位移，下图为一三质点体系各振型的地震水平作用示意图。

体系简图第一振型第二振型第三振型图2.5三质点体系各振型下的示意图2.4水平地震影响系数地震影响系数α就是单指点弹性体系在重力加速度为单位的质点最大加速度反应。

同时，地震影响系数α也指作用于单质点弹性体系上的水平地震作用与质点重力荷载代表值之比。

在不同烈度下，地震系数k为一具体数值。

则α曲线的形状由β谱决定。

通过地震系数k与动力系数β的乘积，即可得到抗震设计反应谱α－T曲线。

用某一次地震的地面运动加速度计卢梭的反应谱曲线作为依据是不准确是的。

为了满足房屋抗震设计需要，应根据大量强震地面运动加速度纪录算出对应于每一条记录的反应谱曲线，按影响其形状的因素进行分类，再进行统计分析，求出最具代表性的平均曲线作为设计依据，这种曲线称为标准反应谱。

规中采用的抗震设计反应谱曲线就是有上述方法得到的标准反应谱曲线。

图2.6 地震影响系数曲线α为地震影响系数曲线；T 为自振周期（s ）；max α为地震影响系数最大值，按表1确定；g T 为特征周期，与场地条件和设计分组有关，按表2确定；2η为阻尼调整系数，按下式计算，当小于0.55时，取0.5；（2.41）式中 ξ——阻尼比，一般情况下，对钢筋混凝土结构去为0.05，钢结构取0.02；当ξ=0.05时，2η=1.0；γ为曲线的下降段的衰减指数，按下式计算：γ=0.9+ξξ55.005.0+-（2.42）1η为直线下降段的下降斜率调整系数，按下式计算，小于0取0；1η=0.02+()8/05.0ξ-（2.43）表1 水平地震影响系数最大值注：括号中的数值分别用于基本地震加速度为0.15g 和0。

30g 的地区ξξη7.106.005.012+-+=表2 特征周期值注：当计算8度，9度罕遇地震时，特征周期应增加0.05s 。

由上图可以看出，α曲线由4部分组成，及直线上升段()1.00<≤T ；水平段()gT T ≤≤1.0；曲线下降段()g gT T T5≤<；直线下降段0.65≤<T T g 。

同时，上图也可以用下式表示()[]()()()[]()⎪⎪⎩⎪⎪⎨⎧≤<--≤<≤≤<≤-+=0.6552.0)5()(1.01.005.41045.0max 12max 2max2max 2T T T T T T T T T T T T T g g gg g gαηηαηαηηαγγ（2.44）三 程序设计及工程实例3.1程序设计根据前面给出双排桩支护结构力计算，采用Visual Fortrain 6.0 编制可视化的使用计算程序。

计算界面如图3.1所示，该计算程序界面友好，操作灵活，便于用户修改，为工程设计人员提供了设计依据。

计算结果的输出采用列表形式和文档形式两种方式出入，便于计算结果的整理。

3.1 程序总框图3.2工程实例计算五层混凝土框架的水平地震作用。

已知：集中与各楼层出得质量为；6.1081=m t ，6.1082=m t ，4.853=m t ，5.604=m t ，2.405=m t 各层柱的截面尺寸均为400mm*500mm，横梁EI=∞；混凝土的弹性模量E=3.0*1010kN/2m；地震类别为小震，远震；设防烈度为8度；建筑场地为三类场地；拟组合前3个振型的地震作用。

图 3.1 计算图示图 3.2 计算简图由于横梁刚度无限大，所以可以把各楼层出的质量简化为一个质点。

简化以后的模型如图所示3.5。

柱侧位移刚度的计算公式：()∑=i ici hEIk312可得，()3113111212hEIhEIk cc==∑=33105106.12100.312-⨯⨯⨯⨯=36000kN/m()3223221212hEIhEIk cc==∑=33105106.12100.312-⨯⨯⨯⨯=36000kN/m()3333331212hEIhEIk cc==∑=33105106.12100.312-⨯⨯⨯⨯=36000kN/m()3443441212hEIhEIk cc==∑=33104104.8100.312-⨯⨯⨯⨯=47300kN/m()3553551212h EI h EI k c c ==∑=33104104.8100.312-*⨯⨯⨯=47300kN/m 应用所编制的计算程序，进行计算，计算结果见表3.1~表3.8，图3.6～3.13所示。