基于响应面方法的五心底结构形状优化.docVIP

基于响应面方法的心底结构形状优化* 杜家政 隋允康 龙连春 阳志光 (北京工业大学工程数值模拟中心) 摘 要: 三心底结构是国防、航空、航天等领域经常采用的结构，将响应面方法和模型参数化相结合，建立了以结构重量为目标结构强度为约束的形状优化模型。根据结构形状选取独立的设计变量，推导结构重量的解析表达式和结构重量对设计变量的一阶导数，采用响应面方法结构最大应力拟合表达式，建立形状优化的规划模型。通过PCL（Patran Command language）语言编写函数，建立分析模型并将跟设计变量有关的数值参数化，在优化迭代的过程中不断对结构进行调整，最后得到满足约束条件的最优设计。数值表明，该方法具有很强的灵活性，适用范围较广，优化结果具有一定的稳定性和准确性。 关键词: 响应面法 心底结构 形状优化1 引言 在国防、航空、航天等领域，推进剂贮箱的基本结构均是由圆筒段和箱底等组成，箱底的底形可分为：半球形底、三心底、半椭球形底、锥形底以及由这些底形曲线组合的底。三心底对于壳体结构，常规的形状优化算法[2-]可以使节点在初始设计面内移动以达到改变形状的目的，一般通过二级控制建立自然设计变量跟节点坐标之间的函数关系，从而通过设计变量驱动节点改变边界或者内部圆孔的形状，二级控制的中间变量及其形状函数的确定是至关重要的，而且对于不规则边界无法确定。然而，该算法不能使节点离开结构的初始设计面，无法改变心底结构的深浅以及旋转母线的半径，也就不能得到心底结构的最佳形状。 图1 三心底结构的旋转母线 图2 五心底结构的旋转母线 将响应面方法[-7]和模型参数化相结合，建立以结构重量为目标结构强度为约束的形状优化模型，这样就可以实现心底结构的形状优化。根据结构形状选取独立的设计变量，推导结构重量的解析表达式和结构重量对设计变量的一阶导数，采用响应面方法结构最大应力拟合表达式，就得到了形状优化的规划模型，可以采用规划程序求解。响应面法就是采用一次或二次多项式代替难于显式化的状态函数，通过多次迭代调整，一般都能满足实际工程精度，具有较高的效率，很有使用价值，是一个很有发展前景的计算方法，也是处理形状优化的应力约束时常用的方法之一。在Patran中，分析模型的建立可以通过界面操作完成，也可以通过PCL（Patran Command Language）语言来实现。采用PCL语言编写建立分析模型的程序，将跟设计变量有关的数值参数化，输入不同的参数可以得到不同的结构形状，只要能参数化的量都能作为设计变量。 本文将上述方法应用于心底结构的形状优化，针对心底结构的形状选取独立的设计变量，该设计变量独立于有限元模型，可以对心底结构的半径、高度等参数进行控制，突破了以往形状优化只改变边界和内部圆孔的局限，扩大了形状优化的应用范围。 2 心底结构的形状优化模型 五心底结构一般与其它结构连接在一起，连接部分会对五心底的形状发生影响。为了使问题具有一般性，在结构上增加一个高为h直筒柱，如图所示。五心底母线的形状对结构的性能和承载能力有很大的影响，为了使结构能应用在实际工程中，在保证结构还是五心底的基础上，对结构的参数进行优化，从而降低结构的重量、改善结构性能和提高结构承载能力。θ1, R2, θ2为设计变量。 以θ1, R2, θ2 作为设计变量，以心底的面积（相当于结构重量）最小作为目标，以应力作为约束优化模型如下： 其中是区间因子，防止目标函数和约束的近似造成的迭代震荡取0.。为了将转化为可以求解的模型，必须将模型中的结构面积（目标函数）和结构最大应力显式化。根据模型推导五心底的面积对设计变量的函数，求面积对设计变量的一阶导数，将目标函数做阶泰勒展开。采用响应面方法（RSM）将结构最大应力拟合为设计变量的一次函数，则模型转化为规划形式，用规划求解器求解。 3 目标函数的表达式及二次展开 θ3= π/2-(θ1+θ2)，(R3–R2)cos(θ1+θ2) = 900-R1-(R2 -R1)cosθ1，则：R3= [900-R1-(R2 -R1)cosθ1]/ cos(θ1+θ2)+ R2， 五心底部分面积的计算： （1）球冠部分：（R为球冠的半径，h为球冠的高度） （2）球带部分：如图4所示，任取球带的一小部分，其面积为： 其中，为圆心坐标，为母线的半径。 对于下部的球带： 对于上部的球带： 图心底结构设计变量示意图 其中，， 五心底总面积： 五心底的面积对设计变量的导数： 目标函数一阶泰勒展开，并去掉常数项的表达式如下： 4 用响应面方法（RSM）拟合结构最大应力 响应面方法是将隐式函数显式化的一种有效途径，其基本思想是构造一个具有明确表达形式的多项式来近似隐式函数，然后根据组变量函数值来