第31卷第2期 2 0 1 0年2月 兵 工 学 报 Vo1.31 No.2 ACTA ARMAMENTARII Feb. 201O 改进遗传算法和有限元计算相结合的优化方法研究 张军挪,,尚利民,李永建 (军械工程学院火炮工程系,河北石家庄050003) 摘要:为了提高优化效率,采用先验法、字节型动态编码方法以及自适应变异技术对遗传算法 进行了改进,并将其与有限元分析计算相结合,提出了一种适用于复杂结构模型的优化分析方法。 利用该优化方法对某型重机架结构参数进行了优化,仿真结果表明,本文提出的优化方法是可 行的,优化效果良好。 关键词:固体力学;遗传算法;有限元;优化;结构参数 中图分类号:TJ201 文献标志码:A 文章编号:1000-1093(2010)02-0135-04 Research on an Optimization Method of Combining Finite Element Calculation with Improved Genetic Algorithm ZHANG Jun—nuo,WANG Rui—lin,SHANG Li-min,LI Yong-jian (Department of Artillery Engineering,Ordnance Engineering College,Shijiazhuang 050003,Hebei,China) Abstract:The genetic algorithm method was modified by a prior method,a dynamic byte coding method and a self-adaptation variation technique in order to improve its optimization efficiency,a new optimiza— tion method combined with the finite element and improved genetic algorithm was brought forward which is suitable for complex structural mode1.The mount structural parameters of a heavy machine gun were optimized by the optimization method.The simulated results show that the method is feasible and effi— cient. Key words:solid mechanics;genetic algorithm;finite element;optimization;structural parameter 对于结构优化方法的研究很多,如遗传算法、矩 阵摄动理论、逆摄动理论、有限元参数优化、神经网 络等,从文献[1—3]看,这些方法所优化的目标主 要是结构的固有特性,而很少以结构的响应特性为 目标进行优化。对于某些物体来说,如武器系统,在 确定结构参数时,不仅要满足武器质量、刚度、强度 等固有特性的要求,还必须要考虑决定武器射击 1遗传算法优化过程 遗传算法与传统的优化方法相比,它不需要求 出目标函数的导数信息;采用概率的而不是确定的 状态转移规则;尤其适用于较难获得目标函数导数 值的优化问题,如含多构件的装配体的优化问题。 采用遗传算法寻优的主要步骤如下 : 1)基因编码;2)初始群体生成;3)群体中个体 精度的结构响应特性。本文提出的改进遗传算法 和有限元计算相结合的优化算法,利用FORTRAN 语言对改进后的遗传算法进行编写,利用ANSYS 软件的APDL进行参数化建模和有限元分析,并将 适应度的计算;4)对群体中个体进行选择、交叉和 变异等遗传操作;5)适应度的再次计算;6)如果满 足停止搜索判据,迭代停止,输出问题的最优解;否 则,转向步骤4). 其作为子程序在编程语言中进行调用。该优化方 法能实现各种优化目标的复杂结构模型的结构优 化分析。 收稿日期:2008—07—25 遗传算法用于结构优化,所处理的问题一般涉 及应力或位移响应,要进行有限元分析,当结构比较 作者简介:张军挪(1978一),男,讲师。E-mail:shiningzjn@163.corn 136 兵 工 学 报 第 卷 复杂时,需要多次迭代运算,所以在设计遗传算法时 要从各角度提高算法的运行效率和收敛速度。 2遗传算法的改进 为了提高遗传算法的全局搜索能力、局部搜索 能力及搜索效率,在保留群体更新换代过程中最优 值的基础上,本文采用了先验法、字节型动态编码方 l按设计要求给定数据输入 l l先验法确定初始个体并编码 l l随机产生种 它个体l 法以及自适应变异技术。 1)先验法:在算法生产初始群体的过程中,并 不全盘使用遗传算法随机产生个体的机制,而同时 使用了先验法来确定满足所有约束条件的个体,这 既加快了遗传算法优化的进程和全局收敛,又给目 标函数中的罚因子和加权函数的选定起到了参考作 用 。 2)字节型动态编码方法:二进制编码方法具有 空间(全局)搜索能力强、效率高、寻优结果鲁棒性 好的优点,而动态编码方法具有局部收敛性好的优 点,字节型动态编码方法是将二者结合在一起的一 种高效的编码方法 。 3)自适应变异:自适应技术使遗传操作过程中 交叉率和变异率随着个体的适应度大小而自适应变 化。算法运行的初期阶段变异概率较大,有利用产 生新的个体;算法运行的后期阶段变异概率较小,有 利用进行局部搜索。 3基于改进遗传算法和有限元计算相结合 的优化算法的实现过程 图1为改进遗传算法和有限元计算相结合的优 化算法优化过程,改进遗传算法是采用FORTRAN 语言编写,有限元计算采用ANSYS软件。 优化过程步骤如下: 1)按设计要求输入确定的数据; 2)生成初始群体,并对其进行适应度计算和统 计; 3)通过选择、交叉和变异等遗传算子操作对群体 更新,给新群体译码,计算适应度并统计群体信息; 4)新旧两代群体个体相比较,保留适应度最大 者; 5)进人世代更新进化循环,判断新群体是否满 足优化终止条件,满足则执行步骤6),否则跳至步 骤3); 6)得到最优方案,退出迭代优化,参数优化结 束,输出优化结果。 从图1可以看出,程序调用和数据交换是需要 解决的核心问题,在FORTRAN中调用ANSYS可以 I初始代群 进行译码 l 程序调用 有限元 兮析计算 l进行适应度计算和统计 I 数据交换 l参数1七_ 建模l ● ————-1进行选择操作 J 自适应交叉和变异产生新群体l l有限元计算l 要 I 程序调用 ● 新 l计算适应 并统计I 数据交换 l结果数据保存l 三 图1优化的实现过程 Fig.1 The process of optimization 通过编写命令语言实现;FORTRAN程序可以通过 创建的.mac文件来告诉ANSYS改变单元属性和载 荷变化,同时,ANSYS把计算结果写在数据文件中 供FORTRAN程序读取,进而实现二者的数据交 换 。 4优化算法的应用 大口径机轻量化以后,武器的动态稳定性成 为影响机射击精度的关键因素,因此,如何优化武 器的结构参数来改善武器的动态稳定性、提高武器 的射击精度已成为武器结构设计中必须要解决的问 题。本文以某型12.7 mm大口径机为例,采用改 进遗传算法和有限元计算相结合的优化方法对该武 器的架结构参数进行了优化。 4.1优化问题的描述 由于本文优化的目标是为了提高武器的射击精 度,而通过试验证明,武器的射击精度与弹出口 瞬间有很好的对应关系,因此以弹出口瞬间偏 离平衡位置的距离为约束条件,同时,基于该武器结 构修改对武器质量的要求,以武器的质量最小为目 标函数,优化结构的质量为m,初始结构的质量为 第2期 改进遗传算法和有限元计算相结合的优化方法研究 137 m .则优化问题可以描述为 minf(X)=m/mn, (1) g ( )=df( )/b一1≤0,i=1,2,…,n (2) d ( )= ̄/《(X)+d2 ( ), (3) XL≤ ≤ u, (4) 式中:g ( )为一系列的位移约束;d (X)为n连发 射击时,每发弹出口瞬间偏离平衡位置的距离; d (X)、d (X)分别为弹出口瞬间在上下和左 右方向偏离平衡位置的距离;b为允许偏离平衡位 置的最大距离;X为设计变量;X .和 ,分别为设 计变量的下限和上限。 由于遗传算法仅靠适应度函数来评估和引导搜 索,因此,优化问题中的约束条件一般用惩罚方法来 处理。惩罚方法的作用是对个体违反约束条件的情 况给予惩罚,并将此惩罚体现在适应度函数设计中。 这样,一个约束优化问题就转化为一个附带惩罚的 无约束问题。本文将位移约束通过罚函数策略在求 解目标函数和适应度函数时进行考虑,其它约束在 计算中单独进行考虑。优化问题可表示为 N mi ( )=m/mo+p×∑(max (g ( )),o), 式中P为正的惩罚因子。 4.2适应度函数的选择及计算 由于优化问题是要求找出满足约束条件,同时 使目标函数达到最小的设计向量解。因此,需要确 定针对结构设计优化问题的个体评价方法,即建立 遗传算法的适应值函数。适应值是对解的质量的一 种度量,在演化过程中适应值是进行选择的唯一依 据,它通常依赖于解的行为与环境(即种群)的关 系。本文的评价准则为:要求适应度函数是非负的, 且适应值越大表明个体的性能越好。利用罚函数法 构造的适应值函数表示为 G ( )= ————————N—— ———————一', 、一,(5) m/m。+p ×∑max (g ( )),0)+c G( )= m0 , (6) m+P× (max (g (X)),0)+c 式中c为足够小的数。对群体中的各个体利用(6) 式计算其适应值函数,并依次进行选择。G(X)表明 在满足位移约束的条件下,原始结构的体积与优化 结构的体积比值越大,得到的结构越大,也就是求解 的目标函数值越小。选择G( )作为适应值函数, 有以下优点:1)由于位移约束与设计变量有关,将 位移约束引入到适应值函数中,避免了单独处理复 杂的隐式约束的困难,从而提高了搜索效率。2)将 原结构的质量与优化结构的质量的比值作为优化指 标,可以直观的反映优化方法的性能,且不受武器结 构的。 4.3优化过程及结果 根据优化目标的结构特点,选择前架杆厚 、 长 、后架杆厚 、长 :四个参数为设计变量,所 建立的有限元分析模型,如图2所示" 。 图2机构的有限元模型 Fig.2 Finite element model of the machine gun 该优化算法选择群体大小m为10,交叉概率 为0.9,变异率为0.05,遗传算法运行300代后得到 的最优结果,如表1所示。 表1 设计变量的优化结果 Tab.1 The optimized results of the design variables mm 墨5 0 I√ 一一优化前 ——优化后 l 6__——1 ——— —— 占——— 广—— 5 时间/s 图3口上下振动位移一时间曲线 Fig.3 Curves of displacement VS.time at up/down direc— tions of the muzzle before and after optimizings 138 兵 工 学 报 第31卷 图3、4分别为架参数优化前后5连发射击时 口上下、左右振动的位移时间曲线对比结果,5连 发射击时,结构参数优化前后弹丸出口瞬间偏离 平衡位置的对比,如表2所示。表2中,方案A、B 分别为结构参数优化前后的响应结果, ( = 1,2,3,4,5)分别为第i发弹出13瞬时在纵向Y和 横向 的位移, 分别为弹出口瞬时在纵向Y 和横向 的位移的平均值。由表2可以看出,结构 修改后口点上下、左右的振动能量分配更为合理, 每发弹出口瞬间更接**衡位置,说明结构修改 后,武器的射击频率与机整体结构的固有频率有 更好的匹配关系,机的射向一致性更好,提高了武 器的射击精度。 图4口左右振动位移一时间曲线 Fig.4 Curves of displacement VS.time at left/right direc— tions of the muzzle before and after optimizing 表2 弹丸出口瞬时位移的比较 Tab.2 Comparisons of instant displacements of the bullets and the muzzle between before and after optimizings for the different rounds mm QIU Zhi—ping,Wang Xiao-jun,Ma Yi.Interval perturbation meth- 5 结论 改进遗传算法和有限元计算相结合的优化算法 不仅具有遗传算法的自适应、全局搜索能力,对目标 函数没有特别要求,无需梯度信息等优良特性,还具 有有限元软件参数化建模容易、模态分析、动力响应 分析等优点,从应用算例中可以看出,优化结果良 好,说明此优化方法可行,能够实现对各种复杂结构 模型的结构优化分析。 参考文献(References) [1] 陈伟,何飞,温卫东.基于结构参数化的有限元分析方法[J]. 0d for computing regions containing complex natural ̄equencies of structures[J].Journal of Beijing University of Aeronautics and Astronautics,2003,29(5):406—409.(in Chinese) [4] 高松竹,蒲家宁,雍歧卫,等.基于遗传算法的导管——气室 系统参数优化研究[J].后勤工程学院学报,2005,24(1):1— 4. GAO Song—zhu,PU Jia-ning,YONG Qi-wei,ct a1.Data optimiza- tion research of pipe and air chamber system based on genetic algo— rithm[J].Journal of Logistics Engineering University,2005,24 (1):1—4.(in Chinese) [5] 朱清华,张呈林,倪先平,等.改进遗传算法的纵列式直升机 总体参数优化设计[J].南京航空航天大学学报,2006,38 (1):l一5. - ZHU Qing—hua,ZHANG Cheng-lin,NI Xian—ping,et a1.General parameter optimization of tandem helicoper based on improved ge- 机械科学与技术,2003,22(11):948—950. CHEN Wei,HE Fei,WEN Wei—dong.Finite element analysis netic algorithm[J].Journal of Nanjing University of Aeronautics and Astronautics,2006,38(1):1—5.(in Chinese) based on structural parametirc design and its application[J].Me- chanical Science and Technology,2003,22(1 1):948—950.(in Chinese) [6] 刘亚营.改进遗传算法及其在神经网络参数优化中的应用 [D].上海:上海海运学院,2003. LIU Ya-ying.Improved genetic algorithm by adding BP algorithm [2] 杨海威,詹永麒,乔俊伟,等.基于结构的神经网络优化结构参 数的应用[J].系统仿真学报,2003,15(8):1116—1119. YANG Hal—wci,ZHAN Yong-qi,QIAO Jun-wei,et a1.Applica- tion in structural parameter optimization by using architecture- into parameter optimization[D].Shanghai:Shanghai Maritime U— niversity,2003.(in Chinese) [7] 张军挪.某型重机结构动力学仿真与参数优化[D].石家 庄:军械工程学院,2007. ZHANG Jan-nuo.Stuctrure dynamics simulation and parameter based neural network[J].Journal of System Simulation,2003,15 (8):1116—1119.(in Chinese) [3]邱志平,王晓军,马一.结构复固有频率区域的区问摄动法 [J].北京航空航天大学学报,2003,29(5):406—409. optimization of a type of heavy machine gun[D].Sh ̄iazhuang: Ordnance Engineering College,2007.(in Chinese)