框架结构动力性能拓扑优化中若干基本问题研究.docVIP

　　框架结构动力性能拓扑优化中若干基本问题研究 1绪 论 1.1研究背景 在我国大力倡导创新驱动发展战略的今天，越来越需要借助先进的设计手段来实现高端复杂产品的自主研发和创新设计。传统结构设计主要依赖工程师的灵感和经验，通过反复试错(trial and error)实现。这种产品设计模式的缺点是设计周期长、效率低且一般不能保证所得到的设计具有优化的特性。随着计算机硬件水平、计算力学以及数学规划理论的快速发展，旨在从力学(物理学)基本原理出发，借助优化方法，通过求解按照设计需求设定的优化问题，以理性方式得到最优设计方案的的结构优化技术应运而生。结构优化技术的出现，为工程师提供了一种高效、可靠、系统的设计工具⑴。它改变了传统的设计流程，有助于缩短设计周期、节省研发费用以及提高设计产品的性能和竞争力。目前结构优化的理论与方法己经在航空航天、机械运载、土木工程等领域获得了许多成功的应用。 1.2结构优化的理论与方法简介 1.2.1结构拓扑优化 根据设计变量的不同，结构优化设计问题大致可分为三种类型(参见图1.1): (1)尺寸优化：其设计变量为结构的几何尺寸(如杆件的横截面积等)；（2)形状优化：其设计变量为结构或构件的几何形貌(如杆系结构的节点位置、结构中的孔洞轮廓等)；（3)拓扑优化：其设计变量为描述材料空间分布方式的特征函数或构件连接关系的指示函数。相对于尺寸优化和形状优化而言，拓扑优化主要面向结构的概念设计阶段，潜在的优化收益最大。但由于本质上是一个离散0-1优化问题，其求解难度也是最大的。自上世纪80年代以来，拓扑优化就一直是结构优化领域中的一个研究热点。结构拓扑优化包含连续体结构拓扑优化和离散结构拓扑优化(主要以析架、框架结构为主)。离散结构拓扑优化研究的历史可以追溯到19世纪末20世纪初Max等[5]提出了离散结构拓扑优化的基结构(ground structure)方法，标志着离散结构拓扑优化的研究进入了一新阶段。以析架结构拓扑优化为例，所谓基结构方法是在给定的设计域内预先布置足够多的杆件，然后通过尺寸优化方法求解相应的优化问题。最后将优化结果中那些面积达到了下限的杆件删除以获得结构的最优拓扑。Dom等人的工作之后，在基结构框架下涌现出了大量针对析架以及框架结构的拓扑优化工作，基结构方法逐渐成为了离散结构拓扑优化的主流方法。人们普遍认为：从基结构出发，只要给与杆件横截面积一个很小的下限，就可以采用尺寸优化的方法来求解拓扑优化问题。直到Sved和Ginos[6]给出了一个反例，人们才发现基结构方法并不是对所有拓扑优化问题都是有效的。 . 2框架结构频率分析旳动力刚度法及s - oulU-Euler梁的动力刚度法简介 结构的自由振动分析能够求出结构的自振频率及振型，而自振频率及振型是结构动力特性的反映，在工程实际中应用广泛，也是抗震结构设计的重要依据。杆系结构包括了析架、框架和混合结构，作为一种工程中常见的结构形式，其振动特性的精确求解也显得尤为必要，在传统结构力学中，杆系结构自由振动分析常用的有三种方法[126]:集中质量法：将结构质量集中在有限个点上，将无限自由度问题简化为有限自由度问题；Ritz法(能量法)：将振型设成有限项试函数(也称为形函数)的线性组合，用能量原理确定这些组合系数，从而求得频率和振型的近似解；普通的有限元法：形函数釆用低阶多项式近似的方法。本文采用一种特殊的有限元方法一一动力刚度法，这种方法是一种适用于计算机分析的精确方法。该方法是建立在杆件的振动微分方程上导出的精确单元刚度阵。下面简单介绍动力刚度阵中单元刚度阵的形成。这里需要申明，本文中除特殊说明外，所有讨论的杆系结构都是由等截面杆件组成的。 2.2 s和s -算法，简称oulli-Euler梁的动力刚度。进一步研究了轴力作用下的Bemoulli-Euler梁动力刚度。由于Bernoulli-Euler梁理论忽略了剪切变形和转动惯量的影响，所以Cheng等研究了基于Timoshenko梁理论下的动力刚度法。Hoosheako梁理论下轴力对刚度阵的影响。Soshenko梁的动力刚度阵。在实际的工程结构中，许多杆系结构受力不完全是单纯的弯曲和拉伸，因此Baneijee等__[81]对考虑弯扭稱合及轴力影响的Bemoulli-Euler梁的动力刚度进行了解析推导，给出了其具体的形式。李俊等研究了考虑翘曲效应的动力刚度阵。Baneijee等[83]采用符号运算工具推导了 Timoshenko梁的动力刚度阵的解析表达式，之后进一步考虑了轴力[84]和翅曲效应[85]的影响。此外，李俊等[86]还对受轴力作用的复合材料Timoshenko梁的动力刚度阵进行了解析推导。BaneIjee[87]推导了复合夹层梁自由振动的动力刚度阵，之后将其应用到考虑扭转的Tim