水平地震作用下基础滑移隔震MDOF结构地震反应振型分解时程分析法研究及计算机仿真.docVIP

水平地震作用下基础滑移隔震MDOF结构地震反应振型分解时程分析法研究及计算机仿真 　　【摘 要】目的为基础滑移隔震MDOF结构地震反应时程分析计算提供精确解析解算法.方法通过对结构滑移、啮合两种状态边界条件的判断及振型分解法运动方程的建立，并通过对两种运动状态下解的正则坐标的坐标变换实现了两种运动状态的衔接.结果通过计算机编程及具体的工程实例进行仿真计算，并与同一实例的数值近似算法结论进行对比，发现本算法精确度高且性能稳定。结论振型分解法同样适用于基础滑移隔震MDOF结构地震反应时程分析计算并为其它数值近似算法提供精确解析解。 　　【关键词】基础滑移隔震；时程分析法；振型分解法；解析解 　　建筑减震[1]的优势表现在地震历程的中后期，建筑隔震[1]的优势表现在地震历程的前期，地震发生时破坏力最强的阶段往往发生在地震历程的前几秒[2]，所以建筑基础滑移隔震逐步成为热门研究课题，基础滑移隔震结构地震反应时程分析算法常用的是等数值算法[3]，其时间积分步长的合理选择对计算结果的精度的影响很大，并直接影响到计算结果的稳定性，决定时间积分步长的主要因素有两个：（1）结构的最小自振周期，（2）输入地震波谱中所包含的最高频率。实践证明数值算法会有较大的时间损耗，而且计算精度及稳定性也不好保证，本文就是在这种情况下寻求精确解析解且计算在整个历时过程无时间衰减的基础滑移隔震多自由度（MDOF）结构地震反应时程算法的研究。它通过对结构滑移、啮合两种状态边界条件的判断及振型分解法运动方程的建立，并通过对两种运动状态下解的正则坐标的坐标变换实现了两种运动状态的衔接，从而实现了计算机的仿真，成功找到了本问题的精确解析解算法。 　　1.无阻尼多自由度结构在滑移状态下的运动方程及求解 　　设有一层的剪切型具有滑移隔震基础的MDOF结构，基础的质量为，向上每层的质量依次为：结构在时刻处于滑动状态，为基础在本时间段在总整体时间坐标系里的绝对位移，为第层的绝对位移函数；为第层柱子的总刚度，为结构的总重量，为基础和地面滑动平台之间摩擦材料的动滑动摩擦系数，为基础和地面滑动平台之间摩擦材料的静滑动摩擦系数，为已知的地震动函数，为提高计算精度，本文采用的是摩擦力模型，则基础下表面所受的滑动摩擦力为，则有运动方程： 　　（1） 　　其中； 　　； 　　，它们皆为阶矩阵，方程（1）解的形式为： 　　（2） 　　由于结构在滑动时会出现一个平动振型，与其对应的结构自频是零，那么上式中的正则坐标中将出现分母为零的状态，这样一来计算机运算时将会逸出并致使计算中断，所以必须对上式中的与对应的正则坐标进行化简！ 　　================ 　　+ 　　+ 　　====================+ 　　+。 　　其中： ；。 　　这里是结构的第个固有频率；是结构与对应的第个固有振型。 　　分析：若当时，与其对应的平动振型为：，若在时刻； 　　，则有： 　　，这里，由此可以发现结构的整体滑动的原因就是由于本平动振型的出现使结构会发生大幅度的平移的可能。 　　2.无阻尼多自由结构在啮合状态下的运动方程及求解 　　若在时刻基础滑移结构由于静滑动摩擦力的过大会导致基础和地面平台之间咬合在一起运动，这时结构的运动状态将和上节所说的大不一样：这时基础的运动状态是明确的，它的速度和加速度都和地震动函数对应的一样，这时候的结构实质上是一个个自由度的结构，我们可以将一层楼盖的质量改为，向上依次为，结构的各层位移函数依次为。则结构的运动方程如下： 　　（3） 　　和式（1）相比这里的刚度矩阵和质量矩阵已发生了变化：阶数皆为，且这时的 　　= 　　=。 　　这里的K11为（1）式中刚度矩阵左上角的第一个元素，为时刻基础的位移函数，它相当于（1）式中的。式（3）的解的形式为： 　　（4） 　　注意（4）式中的与（2）中的是不一样的，它们的数据的具体取值和（3）式中的、有关。同时注意（3）式所对应的结构的固有频率不会有零的情况存在，所以（4）式的计算机运算不会溢出。 　　3. 结构啮合、滑移两种运动状态切换临界边界条件的判定方法 　　当基础的速度和地面平台的速度相同时刻将是运动状态转换的关键时刻，这时结构是继续滑动还是和地面啮合运动？又当二者咬合到什么时候又开始滑动？若是继续滑动的话滑动摩擦力的方向是保持不变还是相反？这两个问题实质上就是边界条件的判定问题。 　　当=0 时，结构的基础和地面不再发生相互运动，这时候静滑动摩擦力将会发生作用，而且通常它比动滑动摩擦力要大得多，这样一来结构的减震效果将会变差，所以要重视静滑动摩擦力对结构的影响。 　　3.1 结构滑动后是否啮合的判定方法 　　显然当=0 时，如果，结构继续滑