主权项 |
1.一种基于KS函数的整体式结构系统设计的组件外形设计方法,其特征在于包括以下步骤:(a)根据组件的CAD模型,利用KS函数构造组件边界的隐函数;将组件区域划分为子区域,对这些子区域进行布尔交、布尔并或布尔非的操作,得到整个组件区域,布尔交、布尔并和布尔非对应的KS函数分别为:<img file="FDA0000366777410000011.GIF" wi="787" he="135" /><img file="FDA0000366777410000012.GIF" wi="841" he="137" /><img file="FDA0000366777410000013.GIF" wi="520" he="82" />其中,p为KS函数的参数;对m个子区域作布尔交,这m个子区域的隐函数分别为Φ<sub>1</sub>,Φ<sub>2</sub>,...,Φ<sub>m</sub>,Φ<sub>max</sub>为Φ<sub>1</sub>,Φ<sub>2</sub>,...,Φ<sub>m</sub>的最大值函数;对n个子区域作布尔并,这n个子区域的隐函数分别为Φ<sub>1</sub>,Φ<sub>2</sub>,...,Φ<sub>n</sub>,Φ<sub>min</sub>为Φ<sub>1</sub>,Φ<sub>2</sub>,...,Φ<sub>n</sub>的最小值函数;KS<sub>con</sub>、KS<sub>dis</sub>和KS<sub>neg</sub>分别为布尔交、布尔并和布尔非对应的KS函数表达式;(b)通过结构的CAD模型建立有限元模型,将设计域划分为规则的固定网格,组件嵌入式分布在设计域的网格中,定义载荷和边界条件;(c)建立组件协同拓扑优化问题的模型为:find η=(η<sub>1</sub>,η<sub>2</sub>,...,η<sub>enum</sub>);S=(s<sub>1</sub>,s<sub>2</sub>,...s<sub>n</sub>),其中s<sub>i</sub>=(x<sub>i</sub>,y<sub>i</sub>,θ<sub>i</sub>)min f(η,S)s.t.KU=F<maths num="0001"><![CDATA[<math><mrow><msub><mi>G</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>η</mi><mo>,</mo><mi>S</mi><mo>)</mo></mrow><mo>≤</mo><mover><msub><mi>G</mi><mi>j</mi></msub><mo>‾</mo></mover><mo>,</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>J</mi></mrow></math>]]></maths><maths num="0002"><![CDATA[<math><mrow><msub><mi>C</mi><mi>k</mi></msub><mrow><mo>(</mo><mi>S</mi><mo>)</mo></mrow><mo>≤</mo><msub><mover><mi>C</mi><mo>‾</mo></mover><mi>k</mi></msub><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mi>K</mi></mrow></math>]]></maths>其中,η为设计域上的单元伪密度向量;enum为设计域网格数目;S为组件的位置设计变量,其中s<sub>i</sub>=(x<sub>i</sub>,y<sub>i</sub>,θ<sub>i</sub>)分别代表第i个组件质心的x坐标、y坐标和方向坐标;n为组件数目;f(η,S)为拓扑优化问题的目标函数;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点整体位移向量;G<sub>j</sub>(η,S)为第j个约束函数,<img file="FDA0000366777410000016.GIF" wi="69" he="89" />为第j个约束函数的上限,J为约束的数目;(d)对上面建立的模型进行一次有限元分析,分别对目标函数和约束函数进行灵敏度分析,求得关于几何设计变量和伪密度设计变量进行灵敏度分析;选取优化算法GCMMA对该问题进行优化设计,得到最优化结果。 |