9-拓扑优化方法.ppt

5　算例与讨论 5　算例与讨论 5　算例与讨论 　　例2　短边简支板减振问题.板尺寸为40cm×20cm,厚度2mm.基体材料与阻尼材料同例1,阻尼层厚1mm.板中心受集中力作用:幅值0.1N,频率1～100Hz.要求结构中心一阶最大加速度a≤3.5m/s2,考虑频率约束30Hz≤f1≤32Hz.计算中初始拓扑基结构划分为4×4个板阻尼胞单元,结构如图5所示.由于对称性,具有相同拓扑敏度值的阻尼胞单元可分为4组:(1,4,7,16)、(2,3,12,15)、(5,6,10,14)和(8,9,11,13).表2给出基体材料结构、阻尼拓扑基结构和结构中分别含指定阻尼胞单元时待控点处加速度、拓扑敏度计算结果.则粗选的阻尼材料拓扑配置为(8+9+11+13),网格细化8倍后得优化的阻尼材料拓扑配置如图6所示. 5　算例与讨论 5　算例与讨论 5　算例与讨论 第六节 结构多学科设计优化 随着计算机技术的进一步发展, 人们开始尝试将多学科的设计综合在一起进行协调优化。事实上, 在汽车、船舶等制造领域, 早已提出了基于全寿命周期的产品开发策略, 也就是并行工程理论。并行工程通过集成产品小组( IPT) 打破传统的学科领域界限,实现多学科设计人员协同设计, 通过设计阶段的快速反馈减少产品开发的迭代设计周期, 加快产品开发过程。而在飞行器设计领域, 也有类似的方法, 一般称为多学科优化设计(MDO) 。 Multidisciplinary Design Optimization MDO 是一种通过充分探索和利用工程系统中的相互作用的协同机制来实现复杂飞行器设计的方法论, 其主要思想是在复杂飞行器设计过程中利用分布的计算机网络来集成多学科的知识, 采用有效的设计优化策略, 组织优化设计过程, 从而能在综合考虑系统整体的情况下, 获得设计最优解。 1. MDO 研究内容包括三大方面: 其中,MDO 算法是MDO 领域内最为重要、也最为活跃的研究课题。本文的目的是对MDO算法及其在飞机设计中应用进行归纳和评述。 ①面向设计的各门学科分析方法和软件的集成; ②探索有效的MDO 算法, 实现多学科(子系统) 并行设计, 获得系统整体最优解; ③MDO 分布式计算机网络环境。 2. 多学科设计优化问题的表述 分析复杂系统的有效方法是按学科(或部件)将复杂系统分解为若干个子系统。根据子系统之间关系, 可将复杂系统划分为两类: 一类是层次系统(Hierarchic System );另一类是非层次系统(Nonhierarchic System )。 层次系统特点是子系统之间信息流程具有顺序性, 子系统之间没有耦合关系, 它是一种“树”结构。非层次系统的特点是子系统之间没有等级关系, 子系统之间信息流是“耦合”在一起, 它是一种“网”结构, 也称为耦合系统。 。现实中的复杂工程系统往往属于非分层系统。非分层系统的设计优化问题, 是目前MDO 研究领域的热点。 多学科设计优化问题, 在数学形式上可简单地表达为: 其中: f 为目标函数; x 为设计变量; y 是状态变量; h i (x, y) 是等式约束; g j (x, y) 是不等式约束。状态变量y, 约束h i 和g j 以及目标函数的计算涉及多门学科。对于非分层系统, 状态变量y, 目标函数f , 约束h i 和g j 的计算, 需多次迭代才能完成; 对于分层系统, 可按一定的顺序进行计算。这一计算步骤称为系统分析。只有当一设计变量x通过系统分析有解时, 才能获得约束和目标函数,这一设计方案被称为一致性设计。 寻找: x 最小化: f = f (x, y) 约束: h i (x, y) = 0 ( i= 1, 2, .,m ) g j (x, y) ≤0 ( j= 1, 2, ., n) 3. 多学科设计优化问题的数学表达 目标函数的灵敏度计算 滤波技术（filtering technique） 为了确保拓扑优化设计解的存在，对求解过程应补充一定的限制条件，其中滤波技术即为常用限制条件。 与网格无关的滤波技术的原理是对目标函数的灵敏度计算公式进行如下修正 式中 式中 dist(e,f) 为单元 e 中心点到单元 f 中心点的距离，rmin为滤波尺寸。 该问题也可用固体各向同性惩罚微结构SIMP或材料性能合理近似RAMP (Rational Approximation of Material Properties) 方法求解 SIMP或RAMP 的区别在于材料的弹性模量插值函数的表达式不同： SIMP或RAMP 的刚度矩阵、柔度矩阵以及柔度矩阵的敏度表达式也不同： 该问题也可用MMA 方法求解 七、渐进结构优化法 渐进结构优化法是近年来兴起的一种解决各类结构优化问题的数值方