结构优化敏度分析技术

《结构优化敏度分析技术》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、结构优化的敏度分析技术1敏度分析方法结构敏度分析是为结构优化捉供有关结构约束函数及目标函数的一阶甚至二阶导数信息。结构敏度分析包括有差分法、解析法、解析和差分结合的拟解析法。差分法通用性好,易于实现，但计算最人;解析法实现起来较困难，但计算效率高;拟解析法在计算解析敏度困难时可予采用。这里主要介绍解析法。•(1)拟载荷法考虑线性静力有限元分析的系统方程为畑”(20)如果载荷F与设计变量向量X无关，则由式(20)对X求导可得竺”凤=0SX蹴(21)R=—U令就,则山式(2)可得拟载荷法求结构位移导数的公式为(22)R称为虚拟载荷。冇限元法中，总刚矩阵K为一稀疏带状对称正定矩阵。

2、一般式(20)采用Coleskey三角分解来求解，亦即相当于片1已在结构分析中得到，从而易于获得期畝。这里ao/ar的第£行为码驰1ck

3、砌］arNj从式(22)可知,拟载荷法适介于求解所冇位移对所冇设计变屋或对某一个设计变量的导数。(2)单位载荷法假想仅在位移幻对应的节点和方向上诞加单位载荷马，设其相应的位移响应为少",由式(20)对t求导可得(23)对式(44.5-23)两边前乘单位载荷向量，可得将上式两边转置，并利用K的对称性易得单位载荷法求位移导数的公式为此方法适合于求某些位移自由度对所冇设计变量或对某些设计变量的导数。一般当设计变量较少、位移约束数目较少时用此方法较

4、经济。•(3)性态空间法前两种頌度分析方法是在设计空间中进行的。通常称应力、位移等结构性态响应为性态变虽，因此在性态空间中约束可表示为口T(25)定义伴随变量向量幻，使之满足由式(25)对X求导可得dXarau銘(26)将拟载荷法公式(22)两边左乘必有将式(22)代入式(26)并利用K的对称性可得(28)_哑=0懸語当&一引时，盛，3U，此时式(28)变为訣凱叭9)式(29)为用性态空间法计算位移导数的公式。事实上，单位载荷法可以从这里的性态空间法或前血的拟载荷法导出。性态空间法不但适用于求位移约朿导数，而且也适用于求应力约朿导数。只要对某个约束去求出了伴随变量伽就可以求得

5、其敏度值。•(4)应力敏度分析应力敏度分析用來获取应力导数。考虑应力计算有限元公式式中二——应力；S——应力矩阵;u^移向最。由式(30)对e#求导，并考虑到尺寸优化时S与X无关，故有应力敏度分析公式为3U(31)从式(31)可知，计算应力导数实质上是计算位移导数。山于结构优化中，儿乎所有单元应力均有约束，因此应力导数计算量相当人。考虑到应力约束具有局部特性，故实际应用时可以采用近似方法计算应力导数，其效果是计算量剧减而乂能具有足够的梢度。这种近似计算方法是将结构在某-•迭代步中作暂时静定化处理，冻结结构内力，即近似认为结构内力在某迭代步中与设计变量无关。对于杆，则有9oj3

6、(^)当i=j时，冇当时，有对于其他类型单元，结构应力约束往往采用vonMises当最应力约束。•(5)梁结构位移敏度分析在结构优化中，梁结构是十分复杂而H滩以处理的问题，尤其对于受弯、扭、剪及轴力作用的空间梁结构更是如此。一方而对梁单元来说，单元刚度矩阵十分复杂，它不仅与单元截而积A有关，而几还与抗扭惯性矩厶、抗弯惯性矩△:和以及剪切而积禺、▲等诸多几何尺寸因素有关。另一方血梁单元截面形状种类繁多，不同的截面形状呈现出不同的复杂力学性态。为了减少设计变量数冃及简化问题的复杂程度，许多方法往往假设截面特性参数与截面面积之间有函数关系。例如"=叫小，这里禺、彳根据不同截面选取不

7、同的值。通过这一函数关系的引入，每个梁单元可只设一个设计变量，如面积或抗弯惯性矩厶，使敏度分析及优化过程也变得较为简单。但这种方法的假设与实际工程梁结构差别太人，难以实际应用。这里给出的方法，则可以以单元截面的具体几何尺寸为设计变最。对于空间梁单元，每个节点冇6个自rti度，其单元刚度玄*

8、•分复杂。一般情况下，尸为的函数，而这6个几何特性量与具体截面类型有关，如果直接由单元刚度对具体几何尺寸变最求导，工作最大且求导结果也因截面类型的不同而各异。为了避免各截面类型都对单元刚度分别求导，这里引人了梁结构中间变量即敏度变量的概念。定义梁单元的嫩度变量为单元刚度炬阵对敏度变最的导数

9、为—=—(>=12.3)以单元截面几何尺寸为设计变量时，设计变最集合可用向最"二【巧•耳•…■比r表示，n为设计变虽总数。对敏度变量的导数，可以通过Jacobian变换转变为对设计变虽的导数：低卜碉［百卜学点=h2•…■耐丸签…厠式中Jacobian矩阵;m敏度变量数目。对空间梁结构沪6。当函数f为梁单元刚阵时,则由上式可求得梁单元刚度矩阵对设计变量的导数，进而可再由单位虚载荷法或拟载荷法求得位移导数。2敏度分析的实现在结构分析软件中增加敏度分析系统的原则是，敏度分析系统应作为独立的功能模块,并与原结构