发明名称 基于KS函数的整体式结构系统设计的组件外形设计方法
摘要 本发明公开了一种基于KS函数的整体式结构系统设计的组件外形设计方法,用于解决现有方法子区域间存在重合边界导致组件产生多余边界的技术问题。技术方案是用KS函数构造组件边界的隐函数,将组件区域划分为多个子区域,通过对这些子区域进行布尔交、布尔并或布尔非操作,得到整体的组件区域,对这些子区域的隐函数分别作与布尔交、布尔并和布尔非相对应的KS函数操作,得到描述该组件边界的隐函数。对于组件的子区域存在重合边界的情况,由于采用KS函数描述方法得到的组件边界的隐函数在重合边界处不会产生多余边界。本发明方法可以正确判断单元相对于组件边界的位置,建立与实际模型相符的材料插值模型,使得优化顺利进行。
申请公布号 CN103412998A 申请公布日期 2013.11.27
申请号 CN201310354898.8 申请日期 2013.08.14
申请人 西北工业大学 发明人 张卫红;周莹;朱继宏;谷小军
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 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>&eta;</mi><mo>,</mo><mi>S</mi><mo>)</mo></mrow><mo>&le;</mo><mover><msub><mi>G</mi><mi>j</mi></msub><mo>&OverBar;</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>&le;</mo><msub><mover><mi>C</mi><mo>&OverBar;</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对该问题进行优化设计,得到最优化结果。
地址 710072 陕西省西安市友谊西路127号