发明名称 考虑拔模制造约束的拓扑优化设计方法
摘要 本发明公开了一种考虑拔模制造约束的拓扑优化设计方法,用于解决在不规则网格模型的拓扑优化设计中引入拔模制造约束的技术问题。技术方案是采用拓扑优化变量定义方式,将拓扑优化设计变量定义为处于拔模方向上同一列相邻两个单元密度的比。通过该方式,将拔模约束引入到设计变量定义中。解决了在拓扑优化中引入大量约束的问题,并且将拔模约束条件引入到设计变量中,达到设计结构能够通过拔模制造的设计要求。
申请公布号 CN102622486B 申请公布日期 2013.10.09
申请号 CN201210071248.8 申请日期 2012.03.17
申请人 西北工业大学 发明人 谷小军;朱继宏;侯杰;张卫红;王丹
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种考虑拔模制造约束的拓扑优化设计方法,其特征在于包括以下步骤:(a)通过结构的CAD模型建立有限元模型,定义载荷和边界条件;(b)建立拓扑优化模型:find  X=(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>n</sub>)min  Φ(X)s.t.  KU=F                          (1)<maths num="0001"><![CDATA[<math><mrow><msub><mi>G</mi><mi>j</mi></msub><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>-</mo><msub><mover><mi>G</mi><mo>&OverBar;</mo></mover><mi>j</mi></msub><mo>&le;</mo><mn>0</mn><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>J</mi></mrow></math>]]></maths>0<δ≤x<sub>i</sub>≤1,i=1,...,n其中,X为设计变量;n为设计变量个数;Φ(X)为目标函数;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点整体位移向量;G<sub>j</sub>(X)为第j个约束函数;<img file="FDA00003100629700012.GIF" wi="62" he="87" />为第j个约束函数的上限;J为约束的数量;δ=10<sup>-3</sup>以避免刚度矩阵奇异;对拔模方向上处于同一列的单元,其伪密度在拔模方向上的顺序依次为ρ<sub>i</sub>,ρ<sub>i+1</sub>,...,ρ<sub>i+m</sub>,m=1,...,M;拔模方向为模型的分模/型面(4)的法向方向,单元i距离模型的分模/型面(4)最近,单元i+m距离模型的分模/型面(4)最远,则:<img file="FDA00003100629700013.GIF" wi="1488" he="236" />对处于拔模方向上同一列的单元ρ<sub>i</sub>,ρ<sub>i+1</sub>,...,ρ<sub>i+m</sub>,m=1,...,M,0<δ≤x<sub>i</sub>≤1则:<maths num="0002"><![CDATA[<math><mfenced open='{' close=''><mtable><mtr><mtd><mn>0</mn><mo>&lt;</mo><msub><mi>&rho;</mi><mi>i</mi></msub><mo>&le;</mo><mn>1</mn></mtd></mtr><mtr><mtd><mn>0</mn><mo>&lt;</mo><mfrac><msub><mi>&rho;</mi><mrow><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><msub><mi>&rho;</mi><mi>i</mi></msub></mfrac><mo>&le;</mo><mn>1</mn></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mo>&CenterDot;</mo></mtd></mtr><mtr><mtd><mn>0</mn><mo>&lt;</mo><mfrac><msub><mi>&rho;</mi><mrow><mi>i</mi><mo>+</mo><mi>m</mi></mrow></msub><msub><mi>&rho;</mi><mrow><mi>i</mi><mo>+</mo><mi>m</mi><mo>-</mo><mn>1</mn></mrow></msub></mfrac><mo>&le;</mo><mn>1</mn></mtd></mtr></mtable></mfenced></math>]]></maths>由上式得:1>ρ<sub>i</sub>≥ρ<sub>i+1</sub>≥,...,≥ρ<sub>i+m</sub>≥0,m=1,...,M,即满足拔模约束条件;(c)将模型进行一次有限元分析;通过优化灵敏度分析,求得目标函数和约束条件的灵敏度;进行优化设计,得到优化结果。
地址 710072 陕西省西安市友谊西路127号