构形状优化响应面方法改进与应用.pdfVIP

第34卷第1期 北京工业大学学报 Vd．34No．1 2008年1月 JOURNALOFBEIJINGUNIVERSITYOFTECHNOILOGY Jan．2008 结构形状优化响应面方法的改进与应用 隋允康，张 轩，宇慧平 (北京工业大学机械工程与应用电子技术学院工程数值模拟中心，北京 100022) 摘要：为了克服-f--维连续体形状优化过程中铒析敏度求解困难的问题，利用响应面方法将目标函数和约束 函数近似显式化，并建立了位移应力约束下的形状优化序列二次规划模型；为了克服响应面方法不够准确的缺 点，将响应面方法进行改造，使其在设计点上取精确值；提出了改良的试验设计方法，达到了以较少的结构分析 代价构造约束响应面的目的；通过建立二次评价函数，确定了具有自适应能力的运动极限。形成了具有较少计算 量和较高逼近精度的优化策略．算例说明这种策略是有效而稳定的． 关键词：形状优化；响应面方法；敏度分析 中图分类号：TP3l 文献标识码：A 文章编号：0254—0037(2008)01—0031—06 形状优化问题是结构优化中层次较高的一类问题，它可以通过改变结构的边界形状达到改善结构力 学性能、提高经济效益的目的．随着工业水平的发展，形状优化技术已经广泛应用于汽车、航天等领 域【1_2J，越来越受到广大学者和工程人员的重视． 灵敏度分析是建立形状优化模型的前提，也是形状优化的主要难题之一。在传统的结构形状优化中。 一般采用解析法或差分法[3]进行敏度分析．一般而言，解析法进行敏度分析时，理论推导复杂，在有些情 况下甚至是十分困难的；差分法虽然可用于各种问题，但其精度受差分步长影响较大，且计算量较大． 利用响应面方法建立形状优化模型可以方便地解决敏度分析困难的问题，而且响应面方法可以通过 同一组试验点集合得到多组响应值，特别适合于处理带有多个约束或多种约束的结构优化问题．然而响 应面方法也存在不能表达精确点、所需试验点过多等缺陷．另外，其整体逼近效果只能利用统计指标模糊 指出，无法构造合理的运动极限，造成近似模型失真．针对以上问题，本文对响应面方法提出改进，构造通 过精确点的响应面，结合适当的试验设计方法，在减少试验点的同时，保留了约束函数的基本曲率信息，并 利用该信息构造二次评价函数．通过线性响应面对评价函数的相对误差确定具有给定逼近精度的响应面 范围，形成具有自适应能力的理性运动极限．文中算例对改进响应面和理性运动极限进行了分析。说明了 这种优化策略能明显提高优化质量，验证了它的有效性和稳定性． 响应面方法及其改进 响应面方法利用试验数据估计原函数的参数[4】．将此方法用于结构形状优化，通常采用线性函数逼 近优化模型中的约束条件，用二次响应面逼近目标函数．在确定二次响应面系数时，通过对二次项进行变 量代换可将其表达为形式上的线性函数为 ’ l 歹=凤+∑胁f (1) i=1 为了确定参数展，做m次独立试验(m≥z+1)，得到m组实验数据。z≯’表示变量的第i号分量在第 _f次试验中的取值，y(j)表示第．『次试验的函数实验值．利用最小二乘法算出参数鹿． 响应面方法得到的近似函数在整体上最接近试验点集合。但却不通过其中任何一点．将响应面方法 收稿日期：2006．11．28． 基金项目：国家自然科学基金资助项目；北京市自然科学基金资助项目(3042002)． 作者简介：隋允康(1943一)，男，辽宁大连人。教授，博士生导师． 万方数据 32 北京工业大学学报 2008矩 用于结构优化建模，将导致响应面只能无限接近而不能最终收敛到极小值，而结构优化总是希望在每次迭 代中目标函数和约束函数在当前设计点处取值都准确无误． 为了解决以上问题，本文对响应面方法进行改造，基本思想有2点：①响应面函数在点取值恰好等于 该点的试验值y(∞，该