[发明专利]基于多点位移协调约束的结构拓扑优化方法

说明书

技术领域

本发明涉及一种结构拓扑优化方法。特别涉及一种基于多点位移协调约束的结构拓扑优化方法。

背景技术

参照图1-3。在航空航天、汽车制造等领域，大量零部件如飞机舷窗、装配桁架等，在承受外载荷(集中力、热应力等)或自身重力的情况下，与其连接的构件或局部相关控制点会发生相对位移。如果相对变形过大，就会导致控制点组成的几何构型发生扭曲或翘曲等变形。位移协调就是指各控制点的位移相互协调改变，不会导致控制结构内部的剧烈形变。空间多节点位移协调约束旨在约束相关控制点作为一个整体几何构型在空间内产生刚体位移。如飞机机身1上的未变形的飞机舷窗2若发生整体的转动，整体转动后的飞机舷窗4虽然位移控制点3的位移很大，但是由于位移控制点3之间的相对位移很小，故飞机舷窗不会被破坏；变形后的飞机舷窗5虽然每个位移控制点3的绝对位移都很小，但是位移控制点3之间的相对位移很大，故飞机舷窗会被破坏。

参照图4。文献1“杨德庆，隋允康，刘正兴，孙焕纯，应力和位移约束下连续体结构拓扑优化，应用数学和力学，2000”公开了一种位移约束下连续体结构拓扑优化方法。该方法建立了以重量为目标，考虑应力和位移约束的连续体结构拓扑优化模型，并推导出应力及位移约束与拓扑设计变量间显式关系式。利用对偶规划简化模型，通过对应力拓扑解和位移-应力拓扑解的综合协调，进而对于拓扑协调解采用阈值完成从离散到连续的反演。

文献1中公开的方法虽然能够约束给定节点的位移，但是在多位移约束中，每个节点的位移是独立的，该方法不能保形：优化后的构型不能约束节点原来的相对位移，即该方法不能实现多点位移协调约束。在载荷P1-P3的作用下，该方法单独约束点A、B、C的位移，这三点的在加载后位置变为A'、B'、C'，虽然该方法可以单独约束这三点的位移，但是不能约束这三点的相对位移，所以不具备多点位移协调功能，实用性差。

发明内容

为了克服现有连续体结构拓扑优化方法不能实现位移协调约束而导致实用性差的不足，本发明提供一种基于多点位移协调约束的结构拓扑优化方法。该方法通过建立有限元模型，选取一定的位移约束控制点，将这些位移控制点围成的几何区域覆盖上一层辅助单元；建立拓扑优化模型，通过优化灵敏度分析，求得目标函数和约束条件的灵敏度，采用梯度优化算法进行优化设计，得到优化结果。由于辅助单元与各个控制点相连，通过约束辅助单元的整体柔顺度来达到控制各个节点间相对位移的目的，从而达到了多点位移协调约束的目的，实用性强。

本发明解决其技术问题所采用的技术方案是：一种基于多点位移协调约束的结构拓扑优化方法，其特点是包括以下步骤：

步骤一、建立有限元模型，并对模型施加约束和边界载荷。

步骤二、选取位移约束控制点，将这些控制点围成的几何区域覆盖上一层辅助单元，辅助单元的杨氏模量比结构至少低五个数量级。

步骤三、定义拓扑优化模型：

find X=(x₁，x₂，...，x_n)

min Φ(X)

s.t. KU＝F

C_a≤ε

Gj(X)-G‾j≤0,j=1,...,J]]>

式中，X为设计域上的拓扑优化变量向量；n为形状设计变量个数；Φ(X)为拓扑优化的目标函数；K为有限元模型总体刚度矩阵；F为节点等效载荷向量；U为节点整体位移向量；C_a为辅助单元的柔顺度；ε为辅助单元柔顺度的上限，为极小的正数；G_j(X)为第j个约束函数；为第j个约束函数的上限；J为约束的数量。

步骤四、将模型进行一次有限元分析；通过优化灵敏度分析，求得目标函数和约束条件的灵敏度，选取梯度优化算法进行优化设计，得到优化结果。