ansys的优化设计在斜拉桥设计中的应用.docxVIP

ansys的优化设计在斜拉桥设计中的应用 ans是集结构、加热、液体、电、噪声等分析于一体的大型通用有限分析软件。它拥有丰富和完善的单元库、材料模型库和求解器,具有完善的前后处理和强大的接口,能高效地求解各类结构的静力、动力、振动、线性和非线性、模态分析、谐波响应分析、瞬态动力分析、断裂力学等问题。同时,ANSYS还有结构优化设计的功能,可以方便地解决工程结构的优化设计问题。本文讨论利用ANSYS的优化设计来确定斜拉桥的成桥初始最优恒载索力。 1 anasas索力优化设计 优化设计是一种寻找确定最优设计方案的设计。所谓最优化是指在满足某种限制的条件下达到给定目标的最佳效果。设计方案的任何方面都是可以优化的,比如说: 尺寸、形状、支撑位置、制造费用、材料特性等。实际上所有可以参数化的ANSYS选项都可作优化设计。 在ANSYS的优化模块中,有3大变量: 设计变量、状态变量、目标函数。设计变量为自变量,优化结果的取得就是通过改变设计变量的数值来实现的。每个设计变量都有上下限,它定义了设计变量的变化范围。在斜拉桥的索力优化中,初始索力自然就作为设计变量了。载荷和边界条件也可以作为设计变量。状态变量是约束设计的数值,用来体现结构设计应该满足功能上或性能上的要求以及其他一些要求。状态变量可有上下限,或只有单方面的限制。在本问题中,可以取梁上下缘的应力或索力均匀性约束、控制截面内力约束等作为状态变量。目标函数就是要尽量减小的数值,它必须是设计变量的函数。在斜拉桥合理成桥受力状态确定中,可以取弯曲应变能、斜拉索用量和弯矩平方和等作为目标函数。设计变量、状态变量、目标函数统称为优化变量。 优化设计通常包括以下几个步骤: ① 生成分析文件。ANSYS程序运用分析文件构造循环文件,进行循环分析。在分析文件中,模型的建立必须参数化,通常是使用优化变量作为参数。模型建立完成后进行求解,定义分析类型和分析选项,施加载荷,指定载荷步,完成有限元计算。然后进入后处理器中进行参数化提取结果,并赋值给相应的参数,即状态变量和目标函数。② 建立优化过程的参数,进入优化处理器,指定分析文件。声明优化变量,即指定哪些参数是设计变量,哪些参数是状态变量,哪个参数是目标函数。然后选择优化工具或优化方法。ANSYS提供了2种优化方法: 零阶方法(直接法)和一阶方法(间接法)。零阶方法使用所有因变量(状态变量和目标函数)的逼近,该方法是通用的方法,可以有效地处理绝大多数的工程问题。一阶方法使用因变量的一阶偏导数,此方法精度很高,但消耗机时较多。ANSYS还提供了用户优化接口,以方便用户编写适合于自己问题的优化方法来使用。之后就是指定优化循环控制方式,进行优化分析。 2 基于优化技术的初始恒载维护能力 2.1 u3000子优化问题的求解 目前国内外斜拉桥确定恒载索力的方法大致有:刚性支承连续梁法、零位移法、最小弯曲能量法、内力平衡法、影响矩阵法、考虑活载效应的分步算法、用索量最小法等。从节省材料、降低造价的角度出发,静载作用下的弯曲能量最小的优化方法是多种优化方法中较为理想的方法,本文即采用本方法来计算。这种方法就是以结构的弯曲应变能为目标函数来进行索力的优化。 U=∫sΜ2(s)2EΙds(1)U=∫sM2(s)2EIds(1) 设主梁和索塔所积蓄的能量分别为: 建立目标函数如下: U=Ug+?Ut 式中: ?表示索塔与主梁的能量代价之比。然后以索力为设计变量,以主梁的应力为状态变量即可求解。其数学表述为: 最小值: minU({x}) 约束条件: {索力上下限∶{X}l≤{Xj}≤{X}u(j=1?2???n)截面允许应力∶{S}l≤{Sj}≤{S}u(j=1?2???m) 用惩罚函数法将其转化为无约束的单目标优化问题,则罚函数为: Ρ(X?q)=f(x)+n∑j=1Ρx(Xj)+qm∑j=1Ρs(Sj)(2) 式中:Px、Ps为受约束的设计变量和状态变量的惩罚因子。转化为无约束优化问题后,可利用梯度法(最速下降法)、牛顿法、变尺度法等进行求解。以梯度法为例,其迭代公式为: X(j+1)=X(j)+sjd(j) 式中:X(j+1)和X(j)分别为设计变量第j+1次和第j次迭代的结果,sj为最优步长因子,d(j)为第j次搜索方向。迭代的收敛条件为: |f(j)-f(j-1)|≤τ?|f(j)-f(b)|≤τ 式中:f(j)和f(j-1)分别为目标函数第j次和第j-1次迭代的结果,f(b)为最优目标函数,τ为目标函数的公差。 2.2 基于apdl的计算方法 以需要优化的索力作为设计变量,主梁的最大应力为状态变量,以弯曲应变能为目标函数来进行求解。由于ANSYS未提供直接输入初始索力的方法,因此,调整斜拉桥初始索力的过程必须通过改变结构的初应变来实现。在建模和加载时可以利用胡