形状优化的敏度分析.pdfVIP

形状优化的敏度分析 页码，1/4 服务简介 算法原理 优化计算 产品优化 零部件优化 设计软件 客户留言 返回系统首 北京时间 ：20:36:27      今天是： 2012年2月17日 星期五 欢迎访问本站！ 形状优化的敏度分析 形状优化的敏度分析在形状优化问题中占有重要地位。目前大多数形状优化算法都是基于函数的梯度信息的，这靠形 状敏度分析来实现。形状敏度分析主要可分为两类不同的方法:第一类使用离散有限元模型，在有限元系统方程的基础上对 形状设计变量直接求导;第二类使用连续质力学中的物质导数法(Material Derivative Method)，这种方法是一种变分方 法，具有一般性且数值计算精度较高。这两大类方法的差别不在于是否进行结构离散化，而在于何时对结构进行离散。第 一类离散模型是先离散后求导，第二类连续模型是先求导后离散。因此连续模型的物质导数法，结构分析时可以采用有限 元法、边界元法等数值方法。本节介绍第一类方法，有关第二类方法请参见其他文献。 离子模型形状敏度分析又可分为差分法、解析法和拟解析法等。差分法简单，程序易于实现，但计算量太大，一般不 予采用;解析法精度高，但程序实现较困难;拟解析法则将某些导数采用差分法计算而另一些导数则采用解析法计算。本节 介绍常用的位移和应力敏度计算方法。 1  位移敏度计算 对形状优化问题，由有限元系统平衡方法得到 （1） 对上式计算敏度有 （2） 如果形状改变的边界处没有外部载荷，则 ,此时有 （3） 只要已知 ，则可将 作为拟载荷R加在结构上，从而可直接利用有限元分析程序求得位移对形状设计变量 的 导数。实际问题中，当位移约束只在少数几个节点上存在时，可以利用单位虚载荷法来求解，这样效率更高。现在主要问 题在于如何求得精确的 ，这对三维形状优化问题并非易事，要花费大量的数值计算时间。下面给出采用拟载荷法进行 形状敏度分析的方法及步骤。 以三维等参元为例，单元刚阵为 （4） 式中 B——单元应变矩阵; D——弹性矩阵; ——雅可比变换矩阵。 /kuai_su/youhuasheji/gezhongleixingyouhua/6-2.asp 2012/2/17 形状优化的敏度分析 页码，2/4 D与单元形状无关，B和 则与单元形状有关，故将单元刚阵对形状变量 求导可得 （5） 上式计算中一般情况下得不到解析解，需使用数值积分方法: （6） 式中 ——高斯积分权因子 ——高斯积分阶次。 为形函数 对x,y,z的导数。 偏导数 为整体坐标方向 的导数。 类似地有其他偏导数的计算式。式(6)中 （7） /kuai_su/youhuasheji/gezhongleixingyouhua/6-2.asp 2012/2/17 形状优化的敏度分析