结构优化设计的准则法.pptx

《结构优化设计的准则法.pptx》由会员分享，可在线阅读，更多相关《结构优化设计的准则法.pptx（113页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。

1、满应力准则法、齿形法、能量准则法等优点：收敛速度快，与优化问题的规模关系不大,重分析次数与设计变量数目没有直接关系缺点：有局限性，适用于结构布局及几何形状已定的情况，设计是接近最优。第1页/共113页2-1结构优化设计的满应力准则法1.满应力设计基本思想从结构力学的原理出发，以满应力为其准则，使杆件的材料能够得到充分利用的一种方法。其设计思想就是对一个既定的结构布局，通过调整构件的断面尺寸，使各构件承受荷载的能力得以充分发挥。具体地，对布局已定的结构在多种荷载作用下，使结构的每个构件至少在一种荷载情况下的应力达到容许应力，此时就认为结构重量是最轻的。第2页/共113页满应力设计地解法不是按事后

2、的结果来判断是否达到最优，而是先行确定所谓优的准则，严格来讲，它并不是最优设计。一般来说，只有静定结构在单一荷载作用下，满应力设计才可能是最轻设计（由于静定结构的的特点决定）；而超静定结构的情况就完全不同了，由于超静定结构各构件的内力与构件截面尺寸有关，每次调整截面后，将产生内力重分布。近似解。第3页/共113页2.满应力准则法中一些概念1）满应力 结构元件的应力达到容许应力或临界应力。2）工况 结构是在多种荷载作用下，每一种荷载我们称为一种工况。3）应力比 元件的工作应力与容许应力之比。第4页/共113页3.满应力计算（应力比法）下面我们以桁架为例，介绍满应力准则法。设桁架是由n个杆件组成的

3、，受到L种工况的荷载作用：设第i 杆在第j 工况下的内力为Nij(i=1,2,，n;j=1,2,l),第i杆在各种工况下的最大内力为 Nimax，最小重量设计可归结为：求设计变量使桁架重量 最小第5页/共113页1）假设桁架各杆初始截面积为：2）进行结构分析，计算各杆轴力，得到 应力矩阵。第6页/共113页式中：为在第j种工况的作用下，第i根杆件的轴力。进而求得各杆应力：第7页/共113页3）计算应力比矩阵式中：为在第j种工况的作用下，第i根杆件的应力比。我们来分析应力比取值情况：说明此杆的应力大于许用应力，需要增加此杆的截面积；说明此杆的应力小于许用应力，需要减小此杆的截面积；可见可以选择

4、作为调整各杆截面积的依据。第8页/共113页4）形成应力比列阵5）调整方案及收敛判别若否则6）把调整方案作为初始方案 重新迭代。转2）继续迭代。第9页/共113页图中是事先指定的小正数。对于静定桁架，各杆的内力与杆件的截面积无关，因此，上面的迭代公式只要一次迭代。对于超静定结构，是一个重复迭代的过程。第10页/共113页1.1.静定结构静定结构例1.1 以5杆桁架为例。外力 许用拉应力 许用压应力 杆断面最小尺寸 求各截面最小尺寸并使结构最轻第11页/共113页首先假定一个初始设计（静定结构可以不必这样做，但为了和静不定结构设计统一，假定初始断面面积为1），设初始向量为：第12页/共113页根

5、据平衡条件，求出各杆内力：第13页/共113页计算应力比 需根据 计算的应力为拉应力或压应力取相应的许用应力 。和 两者符号要一致，可取绝对值进行计算。第14页/共113页修正初始设计向量 其中第三元件的断面面积按计算为0，小于规定的最小尺寸 ，因此应取 。第15页/共113页对于静定结构，由于其内力不随元件断面尺寸的改变而改变，故只需进行一次修正设计，各元件的工作应力便都达到许用应力值，或元件的尺寸取为规定的最小尺寸。各杆的断面尺寸不能取得比上述尺寸小，否则将不满足强度要求或几何约束条件。因此，满应力设计对静定结构而言，就是最轻质量设计。第16页/共113页4.超静定结构例1.2：三杆桁架如

6、图，进行最轻重量设计。两种工况：（1）（2）解：1）分析由于对称性及工况的对称，优化结果显然 ；2）结构重量作为目标函数；第17页/共113页结构分析：第18页/共113页应力比矩阵：初始方案为：第19页/共113页3）迭代计算收敛判别 迭代计算第20页/共113页设计变量设计变量满足精度要求的设计：第21页/共113页如果继续一直迭代下去：可见退化为静定结构。本题的数学模型：第22页/共113页我们可以利用图解法求解：这是一族平行线，其斜率：第1个约束方程的约束界面：其斜率：第23页/共113页求（a）、(b)方程组：通过以上例子，我们可以发现：1）单工况满应力解将退化为原结构的基本结构；2

7、）单工况满应力解不唯一；3）满应力解不是最优解，只有近似解。第24页/共113页例例1.3 同样考虑上述例子，结构布局、外载方向同样考虑上述例子，结构布局、外载方向也相同，但数值改为也相同，但数值改为 ，此，此时需取时需取 、3个设计变量。个设计变量。对应于外载P1情况的应力为对应于外载P2情况的应力为第25页/共113页代入具体数据后，可写出不等式约束条件为前3式对应于外载P1，后3式对应P2第26页/共113页通过对于主动约束与被动约束的分析可知当 时，约束为：第27页/共113页当 时第28页/共113页目标函数为各元件应力比第29页/共113页满应力设计方法中的几个问题1）是否就是最轻

8、质量设计 满应力设计中，目标函数与最优化条件没有任何直接的联系。因此，对于一个满应力设计，其约束方程必须全部满足。对于静定结构，内力与元件剖面无关，因此进行一次迭代循环元件就全部达到满应力，所以自然是收敛的，而且满应力设计就是最轻设计。可以看出，n个约束条件都是等式，可以唯一决定n个变量Aj,因而这时Aj实际上和最小重量并无数学上的联系，也就是说满应力设计并不保证使结构质量最轻。除非将其约束方程变为不等式和使f(A)为最小，联系起来求解Aj,才能保证结构重量为最轻。这就是数学规划的设计。第30页/共113页 所以除静定结构外，满应力设计并不一定具有最小重量的。就其几何意义来说，每一约束方程就相

9、当于n维空间的一个超曲面，所以满应力设计点就落在n个超曲面的交点上。2)收敛性对于静定结构，内力与设计变量无关，只要运行一次设计运算就可以达到满应力设计，所以自然是收敛的。对于静不定结构，可以知道，如果对设计变量没有几何约束要求，那么利用满应力设计时，结构可能退化成若干种静定结构，它们仍能承受原来的几种情况的载荷，因此，满应力设计的结果就不是唯一的。这就引起人们对满应力设计收敛性的疑问。第31页/共113页从实用上看，如果在开始几次设计迭代中就收敛较快，那么它一般会收敛到某个确定的点上去。如果一开始就收敛的慢或反复无常，那就可能是不收敛的。对于这种情况，可以考虑改变起始设计点而重新进行设计。为

10、了避免过多的迭代次数以节省计算时间，一个方法是对收敛条件中规定的设计允许误差值不要取得过小，或者可以采用把满应力设计方法与射线调参方法结合的措施。第32页/共113页为了加快迭代，我们引入松弛因子 ：松弛因子初始时取1.05至1.1，逐步减少到1。第33页/共113页尽管满应力设计存在着上述一些问题，但是它的若干显著特点，使得这个方法还很重要。具体有：a)从前面讨论可以看出，在一定条件下满应力设计可以是最轻设计。b)虽然不同初始点的满应力设计经常会收敛到不同的结果，其中重要的差异可能很大。但是对于大多数实际结构来说，其差异是很小的，有的甚至很接近于最优解。C)为达到收敛所需要的，迭代数往往很小

11、，加上当前进行几万阶矩阵分析已不成问题，因此满应力设计受问题规模影响很小，是目前进行大型结构最轻设计的有效而又经济手法。d)满应力设计通常是近优设计，因此可以为数学规划法进行优化设计提供所需的较理想的初始点，以节省时间。第34页/共113页2-2 齿形法 满应力法唯一的准则是各杆应力比为1，最后的迭代点是在约束条件交点上取得的！但是实际的最优点是位于起作用约束的界面上的。为了设法使搜索点落到约束的界面上，引出了齿形法可以使迭代点在约束界面上。最优解满应力解第35页/共113页 从概念上讲满应力法，有一缺点，即它与质量目标函数无直接关系，而齿形法则把满应力法与质量的大小联系起来，通过交替进行满应

12、力步与射线步来达到寻求最优解的目的。1.齿形法的概念1）一次迭代过程 结构优化的过程是由最小循环重复迭代完成的，我们称这个最小循环为一次迭代过程。完成一次迭代过程得到的点才能作为原问题的最优点近似。第36页/共113页 2 2）齿形法思路 由于满应力法收敛到非最优点，为了改变这一结果，在满应力迭代过程后增加一种移动，我们称为射线步，构成一次迭代过程。实际上一次迭代过程包含两种移动：满应力步；射线步。齿行法的特点是每进行一次满应力步之后进行一次射线步。即从坐标原点出发经过满应力步的设计点，沿此射线方向回到约束曲线上。由于移动的轨迹类似于齿形。故称齿形法，所以齿行法与满应力法的主要区别只是偶数步在

13、确定新的截面面积时有所不同而巳。第37页/共113页齿形法迭代过程图解2.满应力步：从初始点到满应力设计点的移动。3.射线步：使迭代点沿坐标原点与满应力设计点连线的方向回到主约束曲面的移动。满应力步射线步第38页/共113页4.如何实现射线步 从满应力步所得到各杆应力比列阵中，选择最大的 作为射线步的变更依据：其中：第39页/共113页5.5.射线步有下述射线步有下述3 3个特点个特点：1)1)在走射线步时，由于各元件变量都乘同一常数（用 表示之，），故元件内力不变，而应力则为原应力的 倍。证明：结构位移与载荷之间的关系有：在由k k点到（k+1k+1）点时，各元件面积都同乘一常数 ，故两刚度

14、矩阵间有下述关系：因此，可得 第40页/共113页在（k+1k+1）点有 其物理意义为若各元件面积增至 倍，则位移减少至 倍。第41页/共113页再列出元件内力与位移的关系式其物理意义为各元件内力不变。其物理意义为若各元件面积增至 倍，则应力减少至 倍。第42页/共113页2 2）当 时，这一步所走至2k2k点恰好位于主动约束边界上。令i i元件为此时的主动约束，则有3 3）由于各元件截面积同乘常数 ，则有关系式 或 第43页/共113页上式物理意义为 必在沿空间上的射线方向（即通过原点与 点连线），故我们称之为射线步（也有称作比例步，因乘同一比例常数），如下图所示。此时有0.20.80.40

15、.61.00.20.40.60.81.00216543射线步满应力法第44页/共113页 每步移动都需要计算应力比列阵，偶数步移动是由有限元获得；奇数步是通过上述公式获得。这样就减少了有限元分析次数。第45页/共113页5.齿行法的计算步骤1)奇数步：选择初始截面A，利用式 计算应力比列阵。2)偶数步：选出奇数步中应力比的最大值，利用式确定各杆截面积，计算结构重量，然后利用式 求各杆的应力比。3)收敛判别：若 则 作为设计点，否则转1）继续迭代。第46页/共113页例2-1 三杆桁架如例2-1解：1）由例2-1：2）满应力步，由上式代入初始截面积A=1,1计算应力比3）射线步，计算最大应力比，

16、各杆截面积以及新的应力比，计算结构重量第47页/共113页4）重新迭代齿形法设计点为第4步得到的结果。第48页/共113页6.修改的齿行法 齿行法避免了收敛到非最优点,但是由于精度无法控制,影响了其应用。由于射线步步距取决于满应力步，因此为了提高计算精度，应当修改满应力步的步距。其中：为步长修改因子，的取值取决于对最优解要求的精度，取值越小，精度越高，但需迭代的次数也就越多，通常 取0.5，即 当 时，步长为零；当 时，步长为原步长。第49页/共113页例2-2用齿形法计算三杆桁架（例2-2）。解：初始点同例2-2，进行结构分析：第50页/共113页第51页/共113页第52页/共113页第5

17、3页/共113页第54页/共113页修正的齿形法设计点为第8步得到的结果。第55页/共113页迭代k奇数步偶数步1110.70710.414220.70710.707110.585810.79292.707130.70710.56071.04040.666740.73570.583410.640810.82042.664350.73570.47861.03360.707960.76040.494710.684810.84252.645470.76040.41681.0280.740880.78170.428510.720610.86032.639590.78170.36861.02340.76

18、75100.80000.37722.6399修正的齿形法迭代过程第8步得到的是最后结果。第56页/共113页齿形法齿形法 a)迭代算法本质不变，为应力比列因子迭代。b)迭代过程的每一步考虑是否满足约束条件。c)用目标函数值检验最优性，判断是否停止迭代，是较严格的（有数学基础的）优化方法，目标函数是导出迭代方向的重要因素，即目标函数既是构造优化算法的一部分，又是判断最优性的条件之一。有满应力设计法相比较，有一定的改善，但在迭代算法的没有质的差别。第57页/共113页 我们前面讨论了满应力设计，考虑了应力约束和几何约束（设计变量非负）。也就是考虑了强度条件，对于结构设计还需要考虑其它条件，例如刚度

19、条件。由满应力设计的思想，提出了满位移设计。2-3 2-3 桁架的满位移设计桁架的满位移设计第58页/共113页1.满位移的概念 结构位移控制点位移达到容许的最大值时，我们称为满位移。2.满位移准则 在满应力约束的优化基础上，达到满位移的结构设计认为是最轻的结构。3.满位移获得有两种情况：1）经满应力设计或者其它优化设计之后，满位移自然满足；2）经满应力设计或者其它优化设计之后，没有达到满位移设计，还需要经过改进设计。第59页/共113页分部优化与整体优化分部优化与整体优化主要取决于约束条件的性质（1）局部性约束：只涉及结构的各个部分。如截面尺寸约束、应力约束、局部稳定约束、局部相对变形约束等

20、。（2）整体性约束：牵涉到整体结构。如结点位移约束、结构整体稳定性约束、频率约束。只考虑局部约束的问题可用分部优化求解需要考虑全局约束的问题需要进行整体优化第60页/共113页满应力设计是一种分部优化设计，各杆件应力约束之间的相互影响比较弱各杆的刚度对于桁架结构节点位移都有不可忽略的影响，位移约束不能采用分部优化满位移设计是在满应力设计的基础上，通过调整主动杆件的截面积，使位移约束达到满位移的一种整体优化策略如果满应力设计的优化结果满足位移约束，则无需再进行满位移设计。第61页/共113页满位移设计的基本问题满位移设计的基本问题主动杆件与被动杆件的划分凡在当前满位移设计中，截面不作调整的杆件称

21、为被动杆件；截面要作调整的杆件称为主动杆件。主动杆件面积的确定如何调整主动杆件的面积，使得满足满位移条件，同时又使结构最轻。第62页/共113页被动杆件的确定设桁架在荷载的作用下各杆轴力为：为了用虚功原理求解桁架控制点位移，虚设相应单位荷载。设由此引起的各杆轴力为：由虚功原理：(1)1）对于静定结构，轴力与截面积无关：(2)第63页/共113页由(2)式可知：如果 即此杆的实虚轴力异号时，截面积增加，位移也增加，这个杆为被动杆 如果 即某杆的实虚轴力同号时，截面积增加，位移减小，这个杆为可能是主动杆，也可能是被动杆，受满应力设计限制。2）对于超静定结构，上述结论仍然成立。第64页/共113页由

22、(1)式：(3)其中，和 为第i杆实和虚轴力对第j主动杆件截面积改变的变化量。它们各自是平衡力系:平衡力系平衡力系第65页/共113页如图例子：设 为平衡力系，由于第3杆截面积增加单位面积，考虑结点B，下面力系是平衡的：进而也是平衡的：同理也是平衡的：第66页/共113页为外载荷和虚单位力作用下的各杆变形，由虚功原理知：(4)(3)与(2)完全相同。第67页/共113页主动杆件截面积的确定划分了主动杆件被动杆件后，如何调整主动杆的截面积使得结构达到满位移呢？假设前k个杆件是被动杆，上述问题是一个优化问题：求使满足其中：为被动杆总重量；为被动杆对位移的影响量。(5)第68页/共113页(5)可以

23、应用拉格朗日乘子法求解。引入拉格朗日函数：由(6)(7)第69页/共113页考虑式（4）(8)(9)又由(10)将（9）式代入上式，得：(11)第70页/共113页将（11）式代入（9）式，得：(12)如果桁架材料相同，则上式变为：(13)这就是单功况单位移控制点的满位移法求主动杆件的公式。第71页/共113页6.桁架满位移设计的计算步骤:1）确定初始方案，选取满应力设计作为初始方案：2）进行结构分析，求结构在荷载和虚设荷载作用下的轴力，按其乘积为正号的杆件为主动杆件。3）按式（12）或（13）求满位移下主动杆件的优化截面积：第72页/共113页4）从新确定被动杆件：5）终止判别：第73页/共

24、113页 例2-4 如图桁架，控制点C的容许竖向位移 ，许用应力 及 ，最小尺寸限 ，求桁架的满位移设计。第74页/共113页解：1.满应力设计求初始方案。由于此结构是静定的，可以通过一次迭代求得，各杆轴力：满应力设计：第75页/共113页考虑几何约束：最后结果：作为满位移设计初始点。2.在C点虚设单位荷载如图，求得虚轴力：第76页/共113页由求得的各杆实、虚轴力，确定主动杆。只有第6杆为被动杆件。计算是否满足满位移要求：可见，不满足满位移设计。3.求满应力设计下的主动杆件截面积。第77页/共113页由式（13）得：设：4.被动杆确定，第3、4、6杆为被动杆。第78页/共113页重新计算第7

25、9页/共113页带位移约束的齿行法基本思路 最严约束 应力约束 位移约束 使设计点落到当前最严的约束面上(同时考虑应力约束和位移约束)在每产生一个方案之后射线步满应力步直至连续二个可行设计点足够接近或它们的目标函数值开始增加为止。满位移步 第80页/共113页带位移约束的桁架结构优化模型为求设计变量mins.t.其中，变量数（杆件数）；工况数 ；结点数 ；为第 个结点在 工况下的位移；为 结点的容许位移；为面积的下限。第81页/共113页 在处理位移约束和应力约束时，先不考虑尺寸约束，于每个迭代循环中，都在优化准则步得到新设计点后，再检验尺寸约束，当某个 时，就令 。一、射线步 通过结构分析求

26、出 在各工况下的位移 和应力 检验所有位移约束和应力约束，求出它们与其允许最大值的比值：第82页/共113页应力比位移比找出这些比值中的最大者，记为与其相应的约束即为当前最严的约束。第83页/共113页 当 乘以 即完成了射线步，这时并有总刚度阵第84页/共113页二、满应力步若射线步把设计点拉到应力最严约束面上，则用满应力准则来修改设计变量：式中，是控制步长的阻尼指数。三、满位移步若射线步把设计点拉到位移最严约束面上，则用满位移准则来修改设计变量。这一步的问题是：第85页/共113页求设计变量 min s.t.这里沿最严约束搜索，只需考虑这一个位移约束，所以上式中约束取等号。拉格朗日乘子法

27、第86页/共113页满位移步公式 式中为工况中杆件的内力，为沿位移方向的单位力在杆 中引起的内力。它们都与设计点相应。为控制步长的阻尼系数。效果较好。当满位移迭代公式中圆括号内为负值时，该式无意义，说明第 杆截面不能用满位移步调整。这时该杆截面可保持不变，或根据当前最严的应力约束，用满应力准则步修改截面。第87页/共113页四、计算迭代程序1、选定初始方案 ，进行力学分析，求出应力比和位移比：2、射线步。取将设计点引到最严的约束面上第88页/共113页按下式求应力比和位移比：3、调整步。若最严约束是应力约束，则走满应力步，其调整截面为式中 若最严约束是位移约束，相应地走满位移步，其调整截面为第

28、89页/共113页 然后检查尺寸约束，使各杆 ，于是以上完成了一个循环，将所得 作为下一个循环的初始点。如此反复迭代，直至射线步所得目标函数不再减小，或满足其它收敛条件为止。第90页/共113页2-4 能量准则法 结构在荷载作用下发生变形，荷载所作的功全部转化为能量，储存在弹性体内，这种能量我们称为应变能。1.能量准则的概念 结构某一部分储存的应变能的多少是衡量它参加抵抗荷载作用多少的标志，因此我们应当使材料在结构中的分布与各处的应变能成正比，以便充分发挥材料的作用。由此提出了能量准则。能量准则：使结构中单位体积的应变能达到材料容许值时，结构的重量是最轻的。第91页/共113页2.能量库容的概

29、念 结构的构件看成是储存应变能的仓库，这样就引出了能量库容的概念。能量库容：构件这个仓库所能存储的最大应变能。结构优化的能量准则有可以理解为：结构各构件的应变能都等于其相应构件的能量库容时，此结构重量认为时最轻的。下面我们以桁架为例说明能量准则法。第92页/共113页设有N个杆件组成的桁架，第i杆的应变能：式中：为第i根杆件的应变。每一根杆件的容许应变能 为式中：为第i根杆件的最大容许应变。(1)(2)第93页/共113页由于很难使各构件的应变能与它的容许应变能完全相等，故可把能量准则改为：要求各杆件应变能与其容许应变能的比值趋近于结构总应变能与容许应变能之比值，且等于一常数，即其中式中 、分

30、别为结构的总应变能与总容许应变能；为比例常数。(3)第94页/共113页将式(3)两边同乘 ，得由此构造出下列迭代公式(4)(5)在多工况作用下，采用在各中荷载工况下杆件 的最大值代入，其相应的迭代公式为(6)第95页/共113页一般情况下，上式中 ，用它来控制收敛的速度，的取值大小可以根据求解时收敛的快慢来确定。开始时可取大些，在迭代过程中逐渐逼近于1。在用能量法时，以采用位移法进行结构分析较适宜，因为这样便于由位移求出各杆之应变，从而求出应变能。第96页/共113页3.计算步骤（1）设定一组初始设计方案:，令 ；（2）结构分析：求各杆在各工况下的最大应力 ，计算应变能（3）修改截面：（4）

31、终止迭代检查：若，则输出计算结果；否则 转（2）第97页/共113页用能量准则法优化下例例4.1：三杆桁架如图，进行最轻重量设计。两种工况：（1）（2）第98页/共113页这里，位移的方向规定为：x向右为正，y向上为正，由此可求得 杆件伸长分别为在外荷载 作用下，三杆交点位移为(a)(b)第99页/共113页则 杆件的应变 分别为各杆件的应变能为(c)(d)第100页/共113页假定材料在拉、压容许应力范围内都是线弹性的，则各杆件的容许应变能为(e)(f)将式(d)除以式(f),有(g)第101页/共113页将式(g)除以式迭代公式有当C=1时，上式成为可知，能量准则法的结果与满应力法的结果完

32、全一致第102页/共113页框架结构的满应力设计框架结构的满应力设计浮动应力指数法浮动应力指数法 设经过结构分析，得到第j个单元在控制断面上的轴力为 ，两个方向的弯矩为 和 ，根据平截面假设，该处应力可按下式计算：(1)上式中 为断面积，和 分别是x向和y向的断面抗弯模数 ，三者之间存在关系：(2)上式中 为断面惯性矩，都是常数。第103页/共113页 当结构受力后以轴向变形为主时，应以断面积 作为设计变量，如果以弯矩变形为主，则应以惯性矩 作为设计变量。今设以 作为设计变量，即(3)把式(2)、(3)代入式(1)得：(4)如果以 或 作为设计变量，由(2)式做一些变换，可以得到类似的公式。第

33、104页/共113页 满应力设计是一个迭代方法，现假设在第v次迭代中，设计变量与应力之间存在下列关系：(5)其中 为应力指数。对于只有轴力作用的桁架结构，如果以断面积 作为设计变量，则 ，上式既是常见的应力比公式。如果只有弯矩 作用，以 作为设计变量，由(4)可知(6)第105页/共113页 通常 ，只要无轴向力，在 和 联合作用下，上式仍然成立。对于上面简单两种情况，应力指数 是常数。在一般情况下，由于轴向力和弯矩都不为零，应力指数不再是常数，而是与应力及设计变量有关的一个函数。由式(5)取对数，得到(7)由上式可计算 。以允许应力 代替(5)式中的 ，得到计算新的设计变量的计算式：(8)第

34、106页/共113页 在满应力设计的迭代过程中，由于内力和设计变量在不断变化，也在不断改变，因此取名为浮动应力指数法。由(8)式算得 后，代入(4)式可求得新的应力 ，代入(7)式可求得新的 ，从而完成一次迭代。下次迭代中，即以更新后的 由(8)式计算新的设计变量。第一次迭代中需要的 初始值，可用下列两种方法决定:1)忽略轴力影响，由(6)式计算 ；2)给出一个增量 ，以 代入(4)式计算应力 ，再由(7)式计算 。第107页/共113页 下面说明如何用浮动应力指数法进行框架结构的满应力设计。框架结构一般以弯曲变形为主，因此以各杆件的惯性矩作为设计变量，每次迭代，进行一次内力分析，然后走一个射

35、线步和满应力步。(1)满应力步1)确定j单元应力最大值；2)计算3)计算4)由式(4)计算 。5)更新第108页/共113页(2)射线步通常框架结构以弯曲变形为主，其轴向变形可以忽略。今取上式中，是第j个单元的变量比，是最大的设计变量比。在射线步，新的设计变量取为 即各个单元的设计变量（惯性矩）都按一定比例放大或缩小。从结构力学知，修改断面后，结构的内力保持不变，在变量比最大的单元内，因为 ，所以修改断面后，单元应力正好等于允许应力，其它单元内，由于，修改断面后，其应力小于允许应力。因此，设计点正好落在最严约束界面上。第109页/共113页算例：门式钢框架如图，材料容重 ，允许应力 ，型钢截面性质：即三种工况：1）均布载荷：2）水平集中载荷3）水平集中载荷以各杆的惯性矩作为设计变量第110页/共113页第一次迭代根据设计变量初值进行内力分析，然后分别对各杆进行满应力设计。现以杆1为例，说明计算过程：对于杆1，经计算以A点在工况1时应力最大，此处内力为以 代入(4)式，得由(6)式，初值为由(8)式由(4)式由(7)式由(8)式每次迭代，只须进行一次内力分析第111页/共113页完整计算结果列于下表第112页/共113页感谢您的观看！第113页/共113页