斜拉桥索力优化的非线性规划模型.docxVIP

斜拉桥索力优化的非线性规划模型 1 最小应变能法及其优化算法 斜拉桥的恒载超载分布是衡量设计优良质量的重要标准之一。一旦斜拉桥结构体系确定,总能找出一组索力,它使结构体系在确定性荷载作用下,某种反映受力性能的目标达到最优。这组索力对应的成桥态就是对应目标下最优的成桥内力状态。求解这组最优索力,并在斜拉桥中加以实施,也就实现了斜拉桥的恒载受力优化。因此,在不改变结构参数的前提下,斜拉桥恒载状态的优化,也就转化为索力的优化问题。 在优化方法中,由于优化目标的不同,又有许多方法,如应变能最小法、弯矩最小法、用索量最小法和最大偏差最小法等。其中最小应变能法和求解过程都是使整个结构的内力分布合理,并以此为目标值。由于其求解出的索力能够使整个结构的恒载内力最均匀和最小,因此具有更强的适应性及更高的实用性。但是,当斜拉索较多时,如果不加任何约束条件,用该法得出的弯矩一般比较小,但两端的索力不均匀;如人为地做索力局部调整,容易将受力状态调乱,难以得到索力和梁塔内力均满意的结果。文献提出有约束的最小应变能法克服了上述不足。调索时既可满足对某些节点位移、截面内力或索力的特殊限制,又可通过目标函数最小化使得整个结构的弯矩及索力均匀,从而满足不同的要求。 几何非线性是大跨度斜拉桥结构体系的重要特征。由于大位移、梁柱效应以及拉索垂度等非线性因素的影响,其目标函数及约束函数很难用拉索索力的显式形式来表达,属于非线性优化范畴。而现有的索力优化算法如最小二乘法、影响矩阵法大多以叠加原理为基础,只适用于线性结构。因此,必须探索一套新的更为有效的最优化计算方法,以用于确定桥梁的合理成桥状态。 针对上述非线性优化问题,本文采用了一种适用于大型结构优化求解的FEGP(Finite E1ement Analysis Based Gradient Programming)过程。它借助于结构有限元分析、结构设计灵敏度分析、可行方向法以及对偶单纯形算法,进行反复迭代计算直至满足结构力学性能的最优化。其中,有限元分析法用于计算结构的位移和应力,可确定每一约束条件是否临界,结构的灵敏度分析目的是确定约束函数的导数;而可行方向法则是优化过程的主体。优化迭代时用户可任意指定收敛精度与迭代次数,不会由于在所要求的约束上下界范围内无可行解而造成的算法失败。采用这种方法对某大跨度预应力混凝土斜拉桥的空间非线性恒载索力进行了分析。 2 数学模型 2.1 合理选择桥梁状态的原则 (1) 桥型长索的索力大,呈增长期 索力既要分布均匀,又要有较大的灵活性。通常短索的索力小,长索的索力大,呈递增趋势,但局部地方应允许索力有突变。如0号索(全漂浮体系的桥型时)和1号索的索力通常用较大的值。在所有的索中,不宜有太大或太小索力的索。 (2) 主梁曲线 主梁弯矩通常是混凝土斜拉桥设计中的难点和重点。在成桥状态下,主梁的恒载弯矩要控制在“可行域”范围内。 (3) 徐变、混凝土后期收缩 在恒载状态下,主塔弯矩应考虑活载和混凝土后期收缩徐变的影响。在活载作用下,塔往江侧的弯曲程度比岸侧大,并且混凝土后期收缩徐变的影响往往使塔往江侧偏。因此,在成桥恒载状态下,塔宜向岸侧有一定的预偏。 (4) 组件的应力 合理的成桥状态应使斜拉桥在各种荷载组合作用下,构件的应力包络合理,充分发挥构件材料的使用价值。 (5) 弯曲应变能计算构件数 边墩和辅助墩支座反力在恒载下要有足够的压力储备,最好在活载下不出现负反力。这种受力通常由配重或设置拉力支座来满足。 上述合理成桥状态的确定原则可表述数学问题为 min.Q(F)s.t.Fˉi≤Fi≤ˉFi(i=1?2???m1)σˉj≤σj(F)≤ˉσj(j=1?2???m2)Μˉk≤Μk(F)≤ˉΜk(k=1?2???m3)Rˉl≤Rl(F)≤ˉRl(l=1?2???m4)}(1) 式中Q(F)=p∑i=1∫sΜ(F)EΙds为主梁和主塔的弯曲应变能,将其设为目标函数;p为主梁和主塔弯曲应变能计算构件数;M(F)为在索力F作用下,计算构件的弯矩;I为构件组合截面的抗弯惯性矩,E为材料的压弯弹性模量。 F=[F1F2…Fm1]为斜拉索的索力,为设计变量,Fi?ˉFi和Fˉi为第i号斜拉索的索力和允许偏差的上下限;σj(F)为在索力F作用下,主梁和主塔第j号控制点的应力,为状态变量,ˉσj和σˉj为该点的容许应力;Mk(F)为在索力F作用下,主梁或主塔第k号控制点的弯矩,为状态变量,ˉΜk和Μˉk为该点的控制弯矩;Rl(F)为在索力F作用下,第l个边墩和辅助墩的支承反力,为状态变量,ˉRl和Rˉl为该支承的控制反力。m1为斜拉索的总数,m2为应力控制点的总数,m3为主梁和主塔的弯矩控制点的总数,m4为边墩和辅助墩的支座数。 因此,在不改变结构参数的前提下,斜拉桥恒载状态的索力优化,也就转化