剪力墙位置优化设计研究.docVIP

剪力墙位置优化设计研究 　　摘 要：对剪力墙位置优化设计问题建立了数学模型，针对剪力墙位置分析的离散变量广义0-1规划问题，采用微粒群算法对其优化数学模型进行求解，并由inventor三维建模软件对其进行了模态分析和数字模拟。结果表明，微粒群优化算法对剪力墙位置优化分析和计算是科学和有效的，从而为抗震建筑结构中剪力墙的设计提供有效的方法和依据。 　　关键词：剪力墙；广义0-1规划问题；微粒群；优化设计；数字模拟 　　在抗震结构体系中，剪力墙的布置方式很大程度上决定了剪力墙的刚度和刚度分布，确定剪力墙的原则是剪力墙结构设计和结构优化设计的关键内容[1]。优化剪力墙的配置方式，就是要通过调整和改变剪力墙位置和设置方式，来达到调整和控制抗震的框架剪力墙结构和剪力墙结构体系的每层结构的等效刚度，即在优化过程中不断地修正和改变剪力墙配置的位置和方式，从而达到修改整个抗震结构体系的刚度，使得满足结构体系的地震反应达到安全和经济的最低要求。同时还要尽量保证外部干扰作用下，尽量减少或不出现整个结构体系绕着通过刚心的中心惯性主轴发生扭转振动反应。 　　1 剪力墙优化设计的发展现状 　　剪力墙优化配置问题一直是建筑领域的重要研究课题，且与建筑设计和施工之间有着密切的联系，并对有抗震设防要求的建筑结构的安全和经济有着直接和关键的影响。由于问题的特殊性，剪力墙的优化配置的分析和计算通常涉及离散变量优化问题，实际上剪力墙的优化配置是处理如何用尽量少的剪力墙数量和合理的位置来保证建筑结构体系满足抗震要求的问题。为了寻求有效的离散优化计算方法，分支定界法、隐枚举法、模拟退火技术、蚁群算法和遗传基因算法等都被用来处理离散优化问题的有效算法，特别在大型结构体系的传感器/作动器的优化配置分析方面，取得了长足的进步[2，3]。随着对抗震体系中剪力墙优化设计方法的深入探索和研究，剪力墙的优化布置问题也将成为科研人员研究的热点，因为采用微粒群算法解决大型抗震结构体系中剪力墙优化配置问题有望成为最卓有成效的解决方法之一。 　　2 剪力墙设计中存在的问题 　　由于剪力墙结构体系相当于底部固定顶部自由的悬臂梁体系，剪力墙的数量和布置方式是造成总的混凝土用量和含钢量的根本差异的根本原因[3]。对于结构体系中剪力墙的数量，我国相关建筑设计规范都采用纯剪结构的层间位移角极限值1/1000和结构的基本自振周期0.05n至0.08n（n为结构的总层数）作出严格要求，本文仅考虑剪力墙在满足规范的前提下数量保持不变，只解决剪力墙位置的优化分析和计算问题。 　　3 剪力墙位置优化数学模型 　　3.1 剪力墙结构体系的动力学模型 　　一个n自由度结构在环境干扰和控制力作用下运动方程表示为： 　　（3-1） 　　式中M、C、K分别为n×n阶质量、阻尼和刚度矩阵；q、■和■分别是结构的位移、速度和加速度；q（t）是n维独立位移向量（独立广义坐标向量）；F（t）为p维环境干扰向量；U（t）为r维控制力向量；n×p阶矩阵D和n×r阶矩阵H分别描述了环境干扰和控制力在结构上的分布情况。 　　引入状态向量， 将式（3-1）写成状态方程的形式： 　　（3-2） 　　X（t0）=X0 （初始状态向量） 　　式中， ， ， ， 　　分别是2n×2n阶系统矩阵，2n维状态向量，2n×r阶控制器位置矩阵和2n×p阶环境干扰力位置矩阵，0和I分别为相应维数的零矩阵和单位矩阵。 　　式（3-2）构成剪力墙优化配置问题的动力学模型。其中，U（t）、B和H矩阵的元素均为零。 　　3.2 求解振动体系动态反应的算法 　　剪力墙位置优化设计的模型求解算法是结构体系振动主动控制问题衍化出来的仅有地震干扰作用的状态方程求解的核心内容。从原理上讲，所有现代控制理论的控制算法都可以借鉴过来用于剪力墙优化设计内容的分析方法，本文采用经典线性最优控制法，它的目标函数是剪力墙优化配置算法的理论依据。 　　定义结构体系的二次型性能泛函为 　　（3-3） 　　式子中，Q∈Rn×n为剪力墙结构体系的动态反应的权重矩阵。 　　系统状态最优控制问题就是在无限时间区间[t0，∞]内，寻找最小的Z（t），将系统从初始状态Z0转移到零状态附近，并使式（3-5）定义的性能泛函取极小值。因此，剪力墙最优配置问题的数学描述为 　　求Z（t） 　　Min （3-4） 　　S.t. 　　结构体系状态向量最优问题是一个始端固定、无限时间最优控制问题。 　　本文采用的地震干扰作用是美国的EL_Centro波（NS，1940），其地震加速度时程曲线的数据的最大加速度峰值为3.42m/s2。 　　3.3 微粒群优化算法及其实现 　　微粒群算法（Particle Swarm Optimization，PS