结构优化设计的准则法

P1
2 A12
2 A2 2 A1A2
P1
应力比矩阵：
μ
A2 2 A12
2 A1 2 A1 A2
P1
[
]
2 A1
P1
2 A12
2 A1 A2
[ ]
A2
P1
2 A12
2 A1 A2
[
]
A2 2 A1
1
2 A12
2 A1 A2
2 A1
4
3 A2
12
22
32
A2
2 A1A3
P1 A1 A2
A2 A3
2A1
A2
2
A1
代入具体数据后，可写出不等式约束条件为
A3 2A1 A2 2 2 2 A2 A1A2 (a)
A3 2A1 A2 2 2 A1A2
(b)
A3
2 A1 A2
4 1.5
A2
A1 A2
(c)
2)
收敛判别 1 i 0.1 迭代计算
A(1)
[ 1
A(0) 1
2
A(0) 2
]T
1[
2 2
2 1 ]T
设计变量 cm2
1
W (l)
A(0)
1 1
0.707 0.414
A(1)
0.707 0.414
1.094 0.774
A(2)
0.7735 0.3204
1.0541 0.8153
例1.3 同样考虑上述例子，结构布局、外载方向也 相同，但数值改为 P1 40000N P2 20000N ，此时需
取 A1 、A2 、A3 3个设计变量。
对应于外载P1情况的应力为
11 21
31
2 A1A3
P1 A1 A2
A2
A3
A2 2 A3
A2
2
A3
对应于外载P2情况的应力为
……
……
0.707 0.414
3.8284
0.094 0.227
0.0541 0.1847
……
2.414
2.5081
……
A(96)
0.9898 0.0144
1.0001 0.9898
0.0001 0.0101
2.8141
满足精度要求的设计：
A(96)
0.9898 0.0144
W 2.8141 l
1.静定结构
例1.1 以5杆桁架为例。 外力 P 10000
许用拉应力 7000 许用压应力 3500
杆断面最小尺寸 Amin 0.8
求各截面最小尺寸并使 结构最轻
a④
a ①
⑤
③
② P 10000
• 首先假定一个初始设计（静定结构可以
不必这样做，但为了和静不定结构设计
统一，假定初始断面面积为1），设初
始向量为：
A(0)
A(0) 1
A(0) 2
A3(
0)
1 11 ;
A(0) 4
1
A5(0) 1
• 根据平衡条件，求出各杆内力：
N1 10000
N
N2 N3
100000
;
N
4
10000
N5 14140
N1
N2
N3
A(0) 1
A(0) 2
A(0) 3
i max ij 1 jL
若 1 i
否则 Ai i Ai
A* A
6）把调整方案作为初始方案A [ A1 A2 L An ]T 重新 迭代。转2）继续迭代。
图中ε是事先指定 的小正数。
对于静定桁架， 各杆的内力与杆 件的截面积无关， 因此，上面的迭 代公式只要一次 迭代。
对于超静定结构， 是一个重复迭代 的过程。
/
A2,L , N2L / LL
A2
[ij ]
N n1
/
An
,
Nn2
/
An
,L
,
N nL
/
An
3）计算应力比矩阵
11 /[1 ],12 /[1 ],L ,1L /[1 ]
μ
[1, 2,L
, L]
21 /[ 2 ], 22
L
L
/[2 ],L ,2L /[2 ]
L L
[ij ]
n1 /[n ],n2 /[n ],L ,nL /[ n ]
2 A2 A1A2 2 A1 A2 2 2
A3
4 1.5
A2
A1 A2
2 A1 A2
A3
2 A1 A2 A1A2 2 A1 A2
当 2A1 A2 2 2 0 时
A3
2 A2 A1A2 2 A1 A2 2 2
且需
A1A2 2 A2 0
A3
4 1.5
A2
A1
A2
2 A1 A2
2）工况
结构是在多种荷载作用下，每一种荷载我们称为 一种工况。
3）应力比
元件的工作应力与容许应力之比。
i i /[i ]
3.满应力计算（应力比法） 下面我们以桁架为例，介绍满应力准则法。
设桁架是由n个杆件组成的，受到L种工况 的荷载作用：
设第i 杆在第j 工况下的内力为 Nij(i=1,2,…，n; j=1,2,…,l),第i杆在各种工 况下的最大内力为 Nimax，最小重量设计可归结
A3
2 A1 A2 A1A2 2 A1 A2
目标函数为 W 2A1 A2 2A3 各元件应力比
11
21
31 12
22
32
2
A2
2 A3
2 2 A3
1
2.667
A1
2 A1A3 A1A2 A2 A3 1.333A1
2 A1
A2
2 A1
满应力设计方法中的几个问题
所以除静定结构外，满应力设计并不一定具有最 小重量的。就其几何意义来说，每一约束方程就 相当于n维空间的一个超曲面，所以满应力设计 点就落在n个超曲面的交点上。
• 2)收敛性
• 对于静定结构，内力与设计变量无关，只要运行 一次设计运算就可以达到满应力设计，所以自然 是收敛的。
• 对于静不定结构，可以知道，如果对设计变量没 有几何约束要求，那么利用满应力设计时，结构 可能退化成若干种静定结构，它们仍能承受原来 的几种情况的载荷，因此，满应力设计的结果就 不是唯一的。这就引起人们对满应力设计收敛性 的疑问。
为：求设计变量
A A1 A2L An T
n
使桁架重量 W i Ai li 最小 i 1
1）假设桁架各杆初始截面积为：
X [x1, x2,L , xn ]T [ A1, A2,L , An ]T
2）进行结构分析，计算各杆轴力，得到 应力矩阵。
N11, N12,L , N1L
N
[N1, N2,L
式中：ij为在第j种工况的作用下，第i根杆件的应力比。
我们来分析应力比取值情况：
ij 1
说明此杆的应力大于许用应力，需要增加 此杆的截面积；
ij 1 说明此杆的应力小于许用应力，需要减小 此杆的截面积；
可见可以选择 ij作为调整各杆截面积的 依据。
4）形成应力比列阵
1
2
M
n
5）调整方案及收敛判别
• 优点：收敛速度快，与优化问题的规模关 系不大,重分析次数与设计变量数目没有直 接关系
• 缺点：有局限性，适用于结构布局及几何 形状已定的情况，设计是接近最优。
2-1结构优化设计的满应力准则法
1.满应力设计基本思想
从结构力学的原理出发，以满应力为其准则，使杆件 的材料能够得到充分利用的一种方法。
初始方案为：
A [1,1]T
0.1
A2 2A1
μ
1
2 A12
2 A1A2
2A1
4 3
A2
3）迭代计算
μ(0)
A(0) 2
1
2[
A(0) 1
]2
2 A1(0)
A(0) 2
2A(0) 1Fra bibliotek4 3A(0) 2
2
A(0) 1
2 2 2 1
2(2 3
• 对于静定结构，由于其内力不随元件断 面尺寸的改变而改变，故只需进行一次 修正设计 ，各元件的工作应力便都达到
许用应力值，或元件的尺寸取为规定的 最小尺寸。各杆的断面尺寸不能取得比 上述尺寸小，否则将不满足强度要求或 几何约束条件。因此，满应力设计对静 定结构而言，就是最轻质量设计。
4.超静定结构
结构分析：
W l(2 2A1 A2)
A1A2 2 A12 2 A12 2 A1A2
P1
N
2 A1A2 2 A12 2 A1A2
P1
A1A2 2 A12 2 A1A2
P1
σ
A2 2 A12
2 A1 2 A1A2
P1
2 A12
A1 2 A1A2
例1.2：三杆桁架如图，进行最轻重量设计。两种工况：
（1） （2）
P1 20000N P2 0N P1 0N P2 20000N
[ ] 15000Pa [ ] 20000Pa
解：1）分析由于对称性及工况的
对称，优化结果显然 A1 A3 ；
2）结构重量作为目标函数；
A [ A1, A2 ]T
A3
2 A1 A2
2 1.5
A2
A1 A2
(d )
A3 2A1 A2 2A1 A1A2
(e)
A3 2 A1 A2 2A1 A2 A1A2 ( f )
前3式对应于外载P1，后3式对应P2