发明名称 基于多点约束的多组件结构系统布局优化设计方法
摘要 本发明公开了一种基于多点约束的多组件结构系统布局优化设计方法,用于解决现有多组件结构系统布局优化设计方法效率低的技术问题。技术方案是建立组件和结构的有限元模型,使用多点约束方法建立组件和结构之间的连接关系,创建优化问题的数学模型,进行结构系统有限元分析并求解灵敏度信息,最终优化得到合理的结果。由于引入多点约束技术建立组件和支撑结构的连接关系,避免了组件附近的结构网格重新划分,不仅减少了有限元模型前处理所需要的时间,而且也减少了整体有限元模型的单元总数目,缩短了有限元分析计算所需的时间。对于二维支撑多组件结构系统布局优化设计问题,整个优化过程所需时间由背景技术的20分钟46秒缩短为5分13秒。
申请公布号 CN103455670B 申请公布日期 2016.02.24
申请号 CN201310376821.0 申请日期 2013.08.27
申请人 西北工业大学 发明人 朱继宏;高欢欢;张卫红;谷小军
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 一种基于多点约束的多组件结构系统布局优化设计方法,其特征在于包括以下步骤:步骤一、通过组件和支撑结构的CAD模型建立有限元模型;为组件和支撑结构的材料属性赋值;定义载荷和边界条件;步骤二、使用多点约束方法,建立组件网格(2)的组件连接节点(3)和结构网格(1)的结构连接节点(4)的位移线性关系式:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>d</mi><mn>3</mn></msubsup><mo>-</mo><msubsup><mi>Nu</mi><mi>d</mi><mn>4</mn></msubsup><mo>=</mo><mn>0</mn><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000836447290000011.GIF" wi="1144" he="78" /></maths>式中,<img file="FDA0000836447290000012.GIF" wi="61" he="78" />表示组件连接节点(3)在d方向上的位移数值;<img file="FDA0000836447290000013.GIF" wi="61" he="70" />表示结构连接节点(4)在d方向上的位移列向量,d为x、y或者z;N是组件连接节点(3)的位置坐标代入到结构网格(1)的型函数得到的系数行向量;当组件连接节点(3)的数目大于等于1时,则建立相应的组件网格(2)上的组件连接节点(3)与结构网格(1)上的节点之间的位移线性关系式,并且将这些位移线性关系式改写为矩阵形式:HU=0   (2)式中,H为上述位移线性关系式中的系数项,U表示由结构网格(1)上所有的节点的位移列向量和组件网格(2)上所有的节点位移列向量组合形成的系统总节点位移列向量;步骤三、建立优化问题的数学模型为:find X=(x<sub>1</sub>,x<sub>2</sub>,...,x<sub>n</sub>),ξ=(ξ<sub>jx</sub>,ξ<sub>jy</sub>,ξ<sub>jθ</sub>),j=1,2,……N<sub>c</sub><maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open = "" close = ""><mtable><mtr><mtd><mi>min</mi></mtd><mtd><mrow><mi>&Phi;</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msup><mi>U</mi><mi>T</mi></msup><mi>K</mi><mi>U</mi></mrow></mtd></mtr></mtable></mfenced>]]></math><img file="FDA0000836447290000014.GIF" wi="413" he="122" /></maths>s.t. KU+H<sup>T</sup>λ=FHU=0   (3)V(X)‑V<sub>0</sub>≤0<img file="FDA0000836447290000015.GIF" wi="1166" he="75" /><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>&Omega;</mi><mi>j</mi></msub><mo>&SubsetEqual;</mo><msub><mi>&Omega;</mi><mi>d</mi></msub><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mo>...</mo><mo>...</mo><msub><mi>N</mi><mi>c</mi></msub></mrow>]]></math><img file="FDA0000836447290000016.GIF" wi="518" he="71" /></maths>式中,X为设计域上的拓扑优化设计变量;n为拓扑优化设计变量个数;ξ为表征组件安装位置和安装角度的几何设计变量,N<sub>c</sub>表示组件个数,ξ<sub>jx</sub>、ξ<sub>jy</sub>和ξ<sub>jθ</sub>分别表示第j个组件在x、y的安装位置坐标和安装角度;Φ为目标函数,表示为结构的总体应变能;K为有限元模型总体刚度矩阵;F为节点等效载荷向量;U为节点整体位移向量;V(X)为支撑结构材料用量,V<sub>0</sub>为材料用量上限;Ω<sub>j</sub>、Ω<sub>j1</sub>和Ω<sub>j2</sub>分别是第j个、第j1和第j2个组件占据的区域,Ω<sub>d</sub>表示拓扑优化设计域;步骤四、将有限元模型进行一次有限元分析;分别对几何设计变量和伪密度设计变量进行灵敏度分析,求得目标函数和约束条件的灵敏度,选取梯度优化算法进行优化设计,得到优化结果。
地址 710072 陕西省西安市友谊西路127号