一种基于改进PGSA的轮辐式双层索网预应力优化方法与流程

一种基于改进pgsa的轮辐式双层索网预应力优化方法
技术领域
1.本发明涉及建筑工程领域，特别是一种基于改进pgsa的轮辐式双层索网预应力优化方法。


背景技术：

2.随着近几十年来优化算法的发展，面对错综复杂的实际问题，优化不再局限于传统算法的求解，智能优化算法应运而生，如遗传算法、蚁群算法、粒子群算法、模拟退火算法等。但轮辐式双层索网预应力优化，用这些算法得到的结果往往出现目标值过大，位移控制不合理等问题。
3.申请人将模拟植物生长算法(pgsa)用于轮辐式双层索网结构的预应力优化分析。因模拟植物生长算法在生长的过程中，生长步长的取值是固定，申请人将该算法用于轮辐式双层索网结构的找力分析，因索力为连续变量，可行域较广，过大的生长步长会导致优化效果不良，而过小的生长步长会导致优化效率低下，pgsa的搜索速度需要加强。


技术实现要素：

4.本发明提供一种基于改进pgsa的轮辐式双层索网预应力优化方法，以至少解决现有技术中优化分析方法结果目标值过大，位移控制不合理、搜索速度需要加强的问题。
5.本发明目的在于提供一种基于改进pgsa的轮辐式双层索网预应力优化方法，包括以下步骤：
6.s101设置pgsa的计算参数；
7.s102随机生成初始生长点；
8.s103将生长点输出为设计变量数据；
9.s104读取s103输出的设计变量数据，建立轮辐式双层索网模型；
10.s105对轮辐式双层索网模型进行分析计算；
11.s106获得分析计算的结果，所述结果至少包括杆件内力、节点位移中的一种或多种；
12.s107读取分析结果，按约束条件和目标函数对各个生长点进行评价；
13.s108读取s107评价；若不满足终止条件，继续计算生长点的历史最优位置与全局最优位置更新生长点的位置与生长步长，并返回s103；若满足终止条件，输出最优解和最优值。
14.进一步地，所述优化方法还包括初始生长步长设定过程，所述初始生长步长设定过程包括：
15.设定最大生长次数为m，设置初始pgsa生长步长为w1，设置终态pgsa的生长步长w
m
，建立生长步长的初始放大系数为
16.更进一步地，所述基于改进pgsa的轮辐式双层索网预应力优化方法还包括生长步
长更新过程，所述生长步长更新过程包括：
17.s201设定第i
‑
1次生长步长为w
i
‑1，第i次生长步长为w
i
，其中中
18.s202导入获得更新后的生长步长。
19.进一步地，所述设计变量数据为轮辐式双层索网径向索的预应力数据。
20.更进一步地，所述设计变量数据设定方法为：
21.根据平衡状态下轮辐式双层索网的径向索预应力大小差异，对径向索进行分组，并根据分组后的径向索预应力设定为设计变量数据。
22.进一步地，所述目标函数为撑杆上节点竖向位移的径向位移d。
23.更进一步地，所述目标函数具体为：
24.d＝min{|d1|+|d2|+|d3|+|d4|}，其中，d为优化目标值；d1，d2，d3分别为内、中、外圈撑杆上节点竖向位移；d4为中部受拉环上节点径向位移。
25.更进一步地，所述约束条件包括：
26.单层网壳最大竖向位移限值其中l为网壳跨度。
27.更进一步地，所述约束条件还包括：
28.容许长细比，所述容许长细比为
29.强度，所述强度f满足
30.更进一步地，所述约束条件还包括：
31.稳定性，所述稳定性满足公式
32.本发明相对于现有技术，将经过改进的pgsa算法应用到轮辐式双层索网结构优化，能根据实际问题所设置的优化目标，在大范围可行域空间内寻找满足优化目标的较优解，优化目标值远小于现有优化算法得到的结果，使得轮辐式双层索网结构的位移控制在了更加合理的范围内，具有优化效果好，效率高的特点。
附图说明
33.图1为本发明实施例轮辐式双层索网模型图；
34.图2为本发明实施例pgsa具体流程图；
35.图3为本发明实施例生长点放大系数变化情况对照表；
36.图4为本发明实施例matlab和ansys协同工作流程图；
37.图5为本发明实施例流程图。
具体实施方式
38.为了使本技术领域的人员更好地理解本发明方案，下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。
39.本发明实施例提供一种基于改进pgsa的轮辐式双层索网预应力优化方法，包括以下步骤：
40.s101设置pgsa的计算参数；
41.s102随机生成初始生长点；
42.s103将生长点输出为设计变量数据；
43.s104读取s103输出的设计变量数据，建立轮辐式双层索网模型；
44.s105对轮辐式双层索网模型进行分析计算；
45.s106获得分析计算的结果，所述结果至少包括杆件内力、节点位移中的一种或多种；
46.s107读取分析结果，按约束条件和目标函数对各个生长点进行评价；
47.s108读取s107评价；若不满足终止条件，继续计算生长点的历史最优位置与全局最优位置更新生长点的位置与生长步长，并返回s103；若满足终止条件，输出最优解和最优值。
48.其中，本发明实施例采用matlab编写pgsa与ansys协同工作运行文件，matlab与ansys的协同工作流程如图1、4所示。本发明实施例节点连接形式和结构边界约束条件包括，节点连接形式：钢构之间刚接，拉索与钢构之间铰接，撑杆与钢构铰接；结构边界约束条件：v型柱采用铰接。
49.本发明实施例采用ansys通用有限元分析软件，采用几何非线性和材料非线性分析，利用牛顿
‑
拉斐逊迭代求解，通过考虑应力刚化效应，先找出轮辐式索网结构径向索的索力分布，再根据优化目标包括顶部钢结构节点的竖向位移、撑杆顶部的位移、支座反力，可快速确定预应力水平。如图1所示，分析模型包含了钢构件、拉索，荷载只考虑上部恒载。
50.本发明相对于现有技术，将经过改进的pgsa算法应用到轮辐式双层索网结构优化，能根据实际问题所设置的优化目标，在大范围可行域空间内寻找满足优化目标的较优解，优化目标值远小于现有优化算法得到的结果，使得轮辐式双层索网结构的位移控制在了更加合理的范围内，具有优化效果好，效率高的特点。
51.可选的，所述优化方法还包括初始生长步长设定过程，所述初始生长步长设定过程包括：
52.设定最大生长次数为m，设置初始pgsa生长步长为w1，设置终态pgsa的生长步长w
m
，建立生长步长的初始放大系数为
53.其中，在步骤s102过程中，随机生成初始生长点，并进行初始生长步长设定过程，设置初始pgsa生长步长为w1，设置终态pgsa的生长步长w
m
。
54.特别的，所述基于改进pgsa的轮辐式双层索网预应力优化方法还包括生长步长更
新过程，所述生长步长更新过程包括：
55.s201设定第i
‑
1次生长步长为w
i
‑1，第i次生长步长为w
i
，其中中
56.s202导入获得更新后的生长步长。
57.其中，s108中读取s107评价；若不满足终止条件，继续计算生长点的历史最优位置与全局最优位置更新生长点的位置与生长步长，并返回s103，在执行s103过程中同时进行生长步长更新过程，以更新后的生长步长进行执行s104
‑
s107过程，实现再次生长更新，获得新的分析结果。
58.传统的pgsa模拟植物生长算法，在生长的过程中生长步长的取值是固定，该算法用于轮辐式双层索网结构的找力分析，因索力为连续变量，可行域较广，过大的生长步长会导致优化效果不良，而过小的生长步长会导致有效效率的低下。本发明实施例为加强pgsa的搜索速度，对pgsa每一次生长后的生长步长进行更新。在本发明实施例开始时，生长点处在可行域的随机位置，各可生长点间的距离较大，在算法优选的初期，侧重于利用算法的全局搜索能力，快速找出位于全局最优解附近的生长点。越到后期，搜索更为精细，充分利用算法的局部搜索能力，精确地找到全局最优解。
59.本发明实施例在生长点搜索过程中，利用改变生长步长的取值，调整该算法的优化性能，算法初期搜索时，采用大步长，过大的生长步长使得算法拥有优秀的全局搜索能力；而过小的生长步长使得算法表现出卓越的局部搜索能力，在搜索进行中，不断调低其值，从而达到该算法在不同时期拥有不同的优化性能。本发明实施例在最后一次生长时，生长步长恰好为w
m
。
60.假设最大生长次数为100次，如图3所示，生长点的初始放大系数依次为10、20、10，生长点放大系数在算法运行过程中的变化情况，从图中可以看出，生长点的放大倍数的变化趋势大致为先急剧下降，后缓慢减小，符合算法刚开始优选较为粗劣，而后期优选较为精细的目的。
61.可选的，所述设计变量数据为轮辐式双层索网径向索的预应力数据。
62.特别的，所述设计变量数据设定方法为：
63.根据平衡状态下轮辐式双层索网的径向索预应力大小差异，对径向索进行分组，并根据分组后的径向索预应力设定为设计变量数据。
64.其中，本发明实施例以轮辐式双层索网的径向索的预应力作为设计变量，根据平衡状态下径向索预应力大小的差异，将径向索分为4组，因此设计变量为4个。
65.可选的，所述目标函数为撑杆上节点竖向位移的径向位移d。
66.其中，如图1所示，本发明实施例为创造良好的建筑造型，轮辐式双层索网结构在施工张拉过程中，严格控制上网壳的竖向变形，而撑杆的上节点竖向位移在整体位移中起到重要作用，因此将撑杆的上节点竖向位移考虑到目标函数中。
67.特别的，所述目标函数具体为：
68.d＝min{|d1|+|d2|+|d3|+|d4|}，其中，d为优化目标值；d1，d2，d3分别为内、中、外圈撑杆上节点竖向位移；d4为中部受拉环上节点径向位移。
69.其中，本发明实施例因轮辐式双层索网结构的受力理念为自平衡体系，轮辐式双层索网结构支座出的径向位移一定程度上表现了该结构的自平衡效果，即支座处径向位移越小，说明因预应力引起的自平衡效果更优秀，该设计预应力更加合理。
70.特别的，所述约束条件包括：
71.单层网壳最大竖向位移限值其中l为网壳跨度。
72.特别的，所述约束条件还包括：
73.容许长细比，所述容许长细比为
74.其中，压弯构件[λ]＝150，拉弯构件[λ]＝300；
[0075]
强度，所述强度f满足
[0076]
特别的，所述约束条件还包括：
[0077]
稳定性，所述稳定性满足公式
[0078]
其中，本发明实施例稳定性指压弯构件的稳定性，式中：本例中f＝310n/mm2；所述约束条件还包括拉索的应力限值，拉索在各个工况下的最大应力值不超过668n/mm2。
[0079]
在ansys中，轮辐式双层索网结构的撑杆采用link8单元，拉索采用link10单元，用施加初始应变法给索单元加预应力。构件力学参数如下表所示：
[0080]
材料弹性模量(mpa)温度膨胀系数(/℃)密度(kg/m3)钢构件、支架、撑杆2.06
×
1051.2
×
10
‑57.85
×
103拉索1.6
×
1051.2
×
10
‑57.85
×
103[0081]
由于顶部网壳和拉索连接处的节点构造受力情况特殊，在钢构件、支架、撑杆、拉索的密度计算时引入1.1的增大系数，密度为7.85
×
103×
1.1＝8.635
×
103(kg/m3)。
[0082]
采用分布更新法、ansys优化方法与本发明实施例改进后的pgsa分别对轮辐式双层索网结构的设计预应力进行优化，对比结果如下表所示：
[0083][0084]
由上表可见，改进的pgsa经过约56次生长更新后，最优值趋向于定值19.3mm(而分布更新法和ansys优化法在经过56次迭代计算后的目标函数值分别是46.7和34.6)。在此之后，生长求得的最优解不再变化，在经过少次的生长后，生长空间不再产生新的生长点且逐渐变为空集。
[0085]
在经过改进的pgsa优化之后，本发明实施例轮辐式双层索网结构优化目标值远小于用分布更新法以及ansys优化算法得到的结果，说明改进的pgsa的优化结果使得上轮辐式双层索网结构的位移控制在了更加合理的范围内。
[0086]
在采用改进的pgsa对轮辐式双层索网结构进行预应力优化分析时，其能根据实际问题所设置的优化目标，在大范围可行域空间内寻找满足优化目标的较优解，且其优化效果非常好，优化效率极高，适合用于轮辐式双层索网结构的优化分析。
[0087]
最后应当说明的是，以上实施例仅用以说明本发明的技术方案而非对其限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解，技术人员阅读本技术说明书后依然可以对本发明的具体实施方式进行修改或者等同替换，但这些修改或变更均未脱离本发明申请待批权利要求保护范围之内。