单质点体系地震作用的处理方法

1、单质点体系地震作用的处理方法黄宇 (51116),杨星蕊 (51134),高霰 (52569),冯永艳 (51159)在讨论这个题目之前， 首先我们要明确在我们的计算中， 什么是单质点体系， 哪些结构可以视作单质点体系。根据我们所学的知识，地震作用对结构物产生的影响属于结构动力学部分分析内容。而结构动力分析地基本特征是需要考虑惯性力的作用。所以体系质量的分布及其运动方向（即位置变化特征）是决定结构动力特征的关键因素之一。实际结构都是由变形体构成的，而且质量是连续分布的。所以，确定振动过程中全部质量位置的独立几何参数就需要无限多个，属于无限自由度的体系。因此，在实际计算时，常常是将无限自由度的振

2、动问题近似的转化为有限自由度的方法。而当刚性质体的几何尺寸较小时，又可以忽略因转动引起的惯性力矩，即其视为一个质点。实际工程中，由于横梁的刚度远大于柱子的刚度，分析时，可将横梁视为无限刚性，质量分别集中到两端。因此常将如图所示的排架视为单质点体系。根据动静法列静力平衡方程，质体运动微分方程为：mxkxcxFp t地震波是一种弹性波， 它包括体波和面波， 体波分为纵波和横波， 面波分为瑞雷波和洛夫波。地震波传播速度以纵波最快，横波次之，面波最慢。纵波使工程结构产生上下颠簸，横波使工程结构产生水平摇晃，当体波和面波同时到达时振动最为剧烈。尽管地震运动是三维运动，但是建筑结构通常具备一定的弹性。当其

3、处于弹性状态时，可将三维地面运动对结构的影响分解为三个一维地面运动对结构的影响之和。 地震水平地面运动作用为例进行讨论。当单质点体系受到地震水平地面运动作用时， 将产生相对于地面的水平运动。 质体受到三个力作用：一：惯性力：FIm( xg x) ，其中 xa 为地面加速度，x 为质体相对的加速度；二：阻尼力：FDcx ；三：弹性恢复力：FSkx 。三力平衡，对应的方程式为：mx cx kx mxg方程表示，当结构遇到地震时，相当于将等效外力F= MXa 施加于质量自由度方向上，体系将在动力作用（即地震作用）下产生振动，此时情况为强迫振动即受迫振动，支承处运动对于体系的动力作用就相当于在质量上施

4、加一动力荷载。下面开始求解。.y 2y2 yx g1t .e( t)y t0 x asin ( t )d注： d1 -2因为很小，所以我们取d。一般情况下，结构原始处于静止状态， 所以以上即可作为结构在地震作用下的运动方程。由于 Xg ( t) 不能表示为时间的确定性函数， 而只能用概率统计的方法寻求其统计规律，所以，我们通常的实用方法是地震反应谱法。如图的水塔， 由于质量集中在一个确定的位置，结构可以简化为单质点体系 （如图所示）对其进行地震反映分析。尽管地震的地面运动是三维运动，但因此结构处于弹性状态，可将三维地面运动对结构的影响分解为三个一维地面运动队结构的影响之和。 故以下只讨论单向水

5、平地震对单质点体系的影响。单质点体系在地震水平地面运动作用下，将产生相对于地面的水平运动，如图所示， 此时质点上作用有三种力：式单质点系在单向地面水平运动下的运动方程，在数学上为二阶线性微分方程。在另初始条件（初位移 x (0 ) = 0 , 初速度为 0） 下，方程的解为1t .( t)x td0 x aesin ( t)d该式实际上是用积分形式表达的单质点体系地震位移反映。其中：为无阻尼体系自由振动圆频率； 为阻尼比。 一般工程结构值较小，在 0.010.1 之间； 为有阻尼时体系自由振动圆频率。 一般情况下， 若对式关于时间微分一次和二次，则可得体系地震速度反应和地震加速度，注意到一般值

6、很小，则两式可近似简化为.t .e( t)x (t )0 x acos (t)d.t .(t)x g(t )x (t)d0 x aesin (t)d对于工程设计来说，最有用的是结构地震时程反应（内力、变形等）的最大值。由数学原理知，当结构位移反应取最大值时, 结构速度反应为0 。对比上面的表达式得。xgx2 X。m xg x m a x k X m a x。上式左端为质点的最大惯性力，将其定义为地震作用，即F= m xg x max则式 F= k X max 成为一静力方程。按静力分析方法就可得到结构的最大地震位移及相应的内力反应。为计算地震作用，定义如下地震加速度反应谱Sa= xgx max

7、因地震反应谱可预先计算得到，若以确定地震反应谱，则单质点体系的地震作用可根据其子镇周期对应的反应谱值十分简便地按F=mSa(T) 计算得到。又因地震反应谱与地面运动幅值和频谱有关，为分别考虑他们的影响，引进两个参数.xg max地震系数 k=动力系数Sa(T)(T)g.xg max地震反应谱即可表达为：Sa（ T） =kg( T)1. 地震系数的取值可与地震烈度设防标准相联系。设防烈度一般定义为结构设计基准内超越概率为 10的烈度水平， 各地的设防烈度即可取为当地的基本烈度。当结构设计基准期为50 年时，小震烈度比中震烈度约小 1.5度，大烈度比中烈度约大1 度（当基本烈度为 9度时不到1 度

8、）。2. 动力系数的取值与结构的阻尼有关，其图如下， 我国建筑物和构筑物抗震设计规范均取T=0 时结构刚度无限大，结构绝对加速度反应与地面运动加速度完全相同，故=1 ；T=0T 时，一般为直线b ;(T )TgTT T 时，有max对于参数 b,一般取值为 0.651.0 之间。最后得地震作用的计算公式F = ( mg)k( T ) = G( T )G- 单质点体系质量(T ) -地震影响系数，其变化规律同动力系数这就是地震反应谱法的全过程。下面给出反应谱理论一种直接算法的推导：空间结构在X 方向水平地面加速度作用下的动力学方程为：M DC DK DM E x xg其中M、C、 K 分别是质量

9、矩阵、 阻尼矩阵和刚度矩阵；D 是结构位移向量； E x是列向量，对应结构X 方向位移的向量元素为1, 其余向量元素为0。向量元素1和0分别表示质量矩阵中相应的质量元素对于X 方向地面加速度 xg 产生或不产生惯性力。求解结构特征值问题的方程：K, d2 M ,d可以得到结构的n 阶自振频率j j1,2,n 及 n 阶振型 d j j 1,2,n 。将各dj 排列成 nn的振型矩阵 d ，然后用d T 左乘式，得到：*TMU CU KUd M E x xg其中，M*d T M,d ， C *d T C,d* TKdK,d ， Dd U由于， d 关于 M 、 K 的正交性，同时在工程上通常假定

10、d 关于阻尼矩阵C 也具有正交性，所以M*、C*、 K *都是对角矩阵。若设 c j *2 j j m j *，这样式可以表示为在 n 阶振型向量的正交空间下的n 个独立的单变量微分方程问题如下：mj * u jc j * u jk j *d Tj M , E X xg ( j1,2 , n)或 u j2jj u j2j u jtj xg ( j1,2, n)式中， u j是位移向量D 在由振型向量定义的正交空间内所表示的广义位移向量U的第 j 个元素，2*tj 定义为 X 方向地面振动时，结构计入扭转的第jjk j/ m j。这里将阶振型参与系数，表示为：T*TTtjd j M ,E X / mjd j M , E X/ d j M d jnnX ji M i/ X ji2Yji2ri 2ji2 M ii 1i 1比较式与式，根据单质点体系的计算方法，得到广义位移uj 的最大值：ujtjj g /j2，于是，对应于结构