应力约束下长悬臂桁架的拓扑优化设计.docxVIP

应力约束下长悬臂桁架的拓扑优化设计 20世纪30年代以来，结构优化问题引起了人们的广泛关注。目前，主要采用数值优化方法，但只有少数简单问题的分析和回答，其中大部分是基于悬臂梁问题的[9、10、11、12、13、14]。从技术的角度来看，高度分散的绿色能源的类型，如电网塔和天线等，在工程中被广泛使用。在长期的工程实践中，人们发现了优化设计的问题。例如，45的应用是合理的结构，但30-60的重量更小。结构是否比30-60的重量小。结构的最小边界是，如果存在比30-60的结构，即这种结构的性质，一般比理论上最优的理论解类杆的连续体积的体积大于主体的体积。结构优化研究中的这个重要问题是，该方法采用分析方法，导出这些问题的答案。 1 腹杆长度n 设桁架结构由n节构成, 结构等宽, 每节的尺寸和形状可以不一样, 如图1所示.其结点和脱离体的受力分析, 如图2所示.两种常见的悬臂梁桁架, 如图3所示.假设腹杆和上下弦杆的内力分别为Ni(i=0, 1, …, 2n) 和Ti(i=1, 2, …,n) .由图2 (a) 中的A0点和图2 (b) 中的C1点的平衡条件得 Ν0cosα0=Ν1cosα1=F/2.(1)N0cosα0=N1cosα1=F/2.(1) 由图2 (c) 中的Ai点和图2 (d) 中的Ci点的平衡条件得 Ν2i-1cosα2i-1=Ν2icosα2i=Ν2i+1cosα2i+1.(2)N2i?1cosα2i?1=N2icosα2i=N2i+1cosα2i+1.(2) 由式 (1) 可得 Νi=F2cosαi,i=0,1,2,?,2n.(3) 由图2 (e) 关于Ai点力矩平衡条件, 可得 Τi=Fh2i-1Σj=0aj,i=1,2,?,n.(4) 由几何关系可知, 腹杆的长度li(i=0, 1, …, 2n) 和每节的长度ai(i=0, 1, …, 2n) 分别为 li=h2cosαi,ai=h2tanαi.}(5) 2 结构的最优设计 单工况应力约束静定优化结构的各杆是满应力的, 所以优化的体积是 V=2nΣi=0Νiσpli+22nΣi=1Τiσp(a2i-1+a2i),(6) 式 (6) 中, σp是许用应力.优化问题可以表述为 Find{n,ai,Νi(i=0,?,2n),Τj(j=1,?,n)},minV.(7) 其约束条件:σi≤σp,2nΣi=1ai=L, 平衡条件为式 (1) , (2) . 将式 (3) ～ (5) 带入式 (6) , 整理得 V=2Fσph[(2n+1)h24+2nΣi=0a2i+nΣi=1(a2i-1+a2i)2i-1Σj=0aj].(8) 悬臂梁长L应满足 L=2nΣi=0ai.(9) 为计算体积V的极值, 利用Lagrange乘子引入式 (9) 的条件, 构造函数为 f=2Fσph[(2n+1)h24+2nΣi=0a2i+nΣi=1(a2i-1+a2i)2i-1Σj=0aj]+λ[L-2nΣi=0ai],(10) 其极值条件为 ?f?ai=2Fσph{[(2-(-1)i]ai+L}-λ=0,i=0,1,2,?,2n.(11) 由式 (11) 可解得 ai=λσph2F-L2-(-1)i,i=0,1,2,?,2n.(12) 带入式 (9) 中, 可解得 λ=4(2n+3)FL(4n+3)σph.(13) 将式 (13) 带入式 (12) 得到各结点的位置, 即 a2i-1=L4n+3,a2i=3L4n+3,i=0,1,2,?,n.(14) 因此, 可得 a2i-1a2i=1/3,i=1,2,?,n.(15) 中间结点位于各节的1/4位置.将式 (14) 带入式 (8) , 并记ˉL=L/h, 可得结构体积为 V(n,ˉL)=Fh2σp(4ˉL22n+34n+3+2n+1).(16) 不同节数时, 体积关于杆长的变化规律如图4所示.式 (16) 取极小值的条件是 ?V?n=phσp[1-12ˉL2(4n+3)2]=0.(17)?2V?n2=96FhˉL2σp(4n+3)30.(18) 式 (18) 显然成立, 所以式 (16) 关于节数n有最小值.由式 (17) 求解得 n=(2√3ˉL-3)/4.(19) 当然, 式 (19) 不一定是整数.对于图3 (b) 结构, 有 α2i-1=30°,α2i=60°,i=0,1,2,?,n.(20) 将式 (5) 带入式 (9) 也可以得到式 (19) .因此, 当n取整数时即为图3 (b) 结构.这也就证明了图3 (b) 结构就是最小重量的拓扑优化结构. 如果式 (19) 不是整数, 需要选择一个适当的整数.为此, 作如下分析.对于一个确定的节数n, 有对应的ˉL范围使得结构体积最小.当长度超过某个临界值ˉLc(n) 时, 节数就应该增加到