CN110046437A - 基于果蝇优化算法对栓接结合部间距的优化方法 - Google Patents

基于果蝇优化算法对栓接结合部间距的优化方法 Download PDF

Info

Publication number
CN110046437A
CN110046437A CN201910320600.9A CN201910320600A CN110046437A CN 110046437 A CN110046437 A CN 110046437A CN 201910320600 A CN201910320600 A CN 201910320600A CN 110046437 A CN110046437 A CN 110046437A
Authority
CN
China
Prior art keywords
drosophila
optimization
bolted joint
micro
node
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201910320600.9A
Other languages
English (en)
Inventor
赵永胜
赵海忠
杨聪彬
刘志峰
王建华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN201910320600.9A priority Critical patent/CN110046437A/zh
Publication of CN110046437A publication Critical patent/CN110046437A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06NCOMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00Computing arrangements based on biological models
    • G06N3/004Artificial life, i.e. computing arrangements simulating life
    • G06N3/006Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了基于果蝇优化算法对栓接结合部间距的优化方法,该方法包括如下步骤:对所建立的初始栓接结合部模型进行ANSYS有限元分析,获取结合部各节点的节点编号和节点压力,并且生成ANSYS的命令流文件;基于分形接触理论,使用MATLAB软件,编写运算程序,计算各节点的接触刚度并求和;在isight软件中集成ANSYS和MATLAB组块,并使用其中的optimization组块,引入果蝇优化算法,对栓接结合部螺栓之间的距离long和width进行优化,筛选出致使结合部刚度达到最大值的间距,提供了一种提高高精度机床栓接结合部螺栓间距分布的方法。

Description

基于果蝇优化算法对栓接结合部间距的优化方法
技术领域
本发明涉及机床基础制造领域,特别是涉及基于果蝇优化算法优化机床结合部螺栓之间间距的研究。
背景技术
栓接结合部作为组成机床的基本部件,由于栓接结构在机床的组装和拆卸过程中的方便、简易,而被广泛用在机床设备中,在高精度机床中,栓接结合部是最常见的结合部形式,这些结合部不但破坏结构的连续性,而且会增加整体的复杂性,栓接结合部的质量直接影响着机床的加工精度,作为提高机床加工精度的重要影响因素,栓接结合部质量的优化显得尤为的重要,影响结合部特性的参数很多,主要有:结合部外部载荷、材料、加工方法及表面粗糙度、结合部结构类型与尺寸等[1],
而螺栓间距作为影响栓接结合部质量关键因素之一,对结合部的质量有着重要的影响,。
在生产实践中,机床上螺栓间距的布置多是采用国家标准来进行设计,间距在一个较大的区间范围之内变化,对于普通机床,因为没有过高的要求,国家标准完全适用,但是在高精度机床中,结合部的质量作为影响加工质量的一个关键因素,螺栓的粗略布置就对结合部的质量产生很大影响,进而影响机床的加工质量,所以有必要对高精度机床的螺栓间距进行优化,提高结合部质量,进而提高加工质量。传统分析中,主要采用有限元分析计算不同螺栓间距对工件的影响,没有采用有效的优化算法对螺栓间距进行优化。
本文采用果蝇优化算法,基于ANSYS有限元分析和MATLAB数值计算,并采用联合优化软件isight,对建立所建立模型进行螺栓间距的优化,使其在设定的参数范围内达到结合部刚度最大的预期设定的目标。
发明内容
本发明目的是:为了能够通过优化高精度机床螺栓之间的距离,进而提高栓接结合部刚度,最终影响机床的加工精度,通过使用果蝇优化算法,以有限元分析软件ANSYS和计算软件MATLAB为优化工具,采用isight集合功能,对所建立的栓接结合部模型进行螺栓间距优化,使结合部刚度达到最大。因此本发明提供了一种基于果蝇优化算法模拟栓接结合部间距优化的计算方法。
本发明所采取的技术方案是:
基于果蝇优化算法对栓接结合部间距的优化方法,包括如下步骤:
步骤一、在ANSYS中对所建立的初始栓接结合部模型进行有限元分析,获取结合部各节点的节点编号和节点压力,同时生成ANSYS的命令流文件;
步骤二、以分形接触理论为基础,运用MATLAB软件,编写计算程序,计算各节点的接触刚度;
步骤三、基于果蝇优化算法,采用isight集合功能,将ANSYS和MATLAB软件集合到isight中,确定优化间距长度(long)和宽度(width)的范围、粒子群大小、迭代次数和最终的优化目标等参数,不同的参数设定直接影响优化的能力,进行循环迭代,寻找最大的结合部刚度值。
所述步骤一具体为:
建立栓接结合部模型,进行有限元划分,施加约束和模拟螺栓预紧力,进行有限元静力分析,提取各节点的节点压力;
所述步骤二具体为:
利用结构函数法计算出结合部的分形参数,结合部栓接表面法的轮廓曲线利用Weierstrass-Mandelbrot函数表示
式中z—粗糙表面轮廓高度;
x—表面采样长度坐标;
D—轮廓曲线的分形维数;
G—粗糙表面的特征尺度系数;
γ—谱密度的尺寸参数;
φn—随机相位;
n—频率指数;
对于独立的微凸体形貌表述,根据单频率(γn=1/l)的W-M函数来获得,独立微凸体Z0(x)数学公式为:
对于一个给定接触点截面面积为a′的微凸体,其变形量δ由z0(x)在x=0处求得:
在微观尺度下,认为a′=l2,微凸体的顶点曲率半径为:
由Hertz理论可知单个微凸体的弹性载荷为Fe=4Er3/3R,式中r为真实接触面积半径
将其带入Fe=4Er3/3R,可得
在MB模型中,微凸体接触点数与截面面积关系为
其中a′L为最大微接触点平截面积
根据刚度定义可得到单个微凸体的法向接触刚度为:
结合微凸体截面积分布函数积分得到结合部的总刚度为:
根据以上公式,在MATLAB中编写出节点刚度计算公式,并对各节点刚度进行加和计算。
所述步骤三具体为:
果蝇优化算法作为常用优化方法,又被称之为FOA算法,可以总结为以下几个重要的步骤:
1)随机初始果蝇群体位置X,Y;
2)配置独立果蝇使用嗅觉寻找食物的随机搜索方向与距离,rand()是寻优过程中随机生成的一个0-1之间的一个值:
Xi=X+rand(),Yi=Y+rand()
3)估算果蝇与搜索原点之间的距离D,之后计算味道浓度的判定值S,S为距离参数D的倒数值:
D=(Xi2+Yi2)0.5,S=1/D
4)将味道浓度判定数值引用到味道浓度判定函数来求求独立果蝇个体的位置味道浓度smell:
smell=Function(S)
5)寻到得到果蝇群体中所有独立果蝇味道浓度最大的果蝇:
[bestsmellbestindex]=max(smell)
6)保留味道浓度最大的果蝇的最佳味道浓度值与X、Y的坐标,之后果蝇种群会利用视觉朝向该最大值位置飞去:
smellbest=bestsmell
X=X(bestindex)
Y=Y(bestindex)
7)进入果蝇算法迭代寻优,重复之上的步骤2)~5),并且判断味道的浓度是否优于之前迭代的味道浓度,若是则进行步骤6);
所述步骤四具体为:
进入isight软件中optimization组块,在此插入果蝇优化算法,设置总迭代次数、粒子群数量等参数,选择变量为long和width,并设置变化区间,设置约束和优化目标参数。MATLAB输出总刚度结果计入optimization,optimization变换long和width的数值,等于变化了一个果蝇的飞行方向,循环迭代进入ANSYS输入,optimization会根据总刚度值的结果变换long和width的变化方向,使long和width的变化向着总刚度最大的目标变化,等于变化了果蝇群体的飞行方向。
与现有技术相比较,基于isight、ANSYS和MATLAB软件和果蝇优化算法的螺栓间距的研究,可以有效的对给定的参数在一定区间内进行寻优,筛选出最优的间距距离,使结合部刚度达到最大值。
附图说明
图1结合部正面分布图。
图2单个微凸体形貌图。
图3 isight、ANSYS和MATLAB优化流程图。
具体实施方式
以下结合附图和实施例对本方法进行详细说明。
基于果蝇优化算法模拟栓接结合部间距优化的计算方法,包括如下步骤:
步骤一、在ANSYS中对所建立的初始栓接结合部模型进行有限元分析,获取结合部各节点的节点编号和节点压力,同时生成ANSYS的命令流文件;
步骤二、以分形接触理论为基础,运用MATLAB软件,编写计算程序,计算各节点的接触刚度;
步骤三、基于果蝇优化算法,采用isight集合功能,将ANSYS和MATLAB软件集合到isight中,确定优化间距长度(long)和宽度(width)的范围、粒子群大小、迭代次数和最终的优化目标等参数,不同的参数设定直接影响优化的能力,进行循环迭代,寻找最大的结合部刚度值。
所述步骤一具体为:
建立栓接结合部模型,进行有限元划分,施加约束和模拟螺栓预紧力,进行有限元静力分析,提取各节点的节点压力;
所述步骤二具体为:
利用结构函数法计算出结合部的分形参数,结合部栓接表面法的轮廓曲线利用Weierstrass-Mandelbrot函数表示
式中z—粗糙表面轮廓高度;
x—表面采样长度坐标;
D—轮廓曲线的分形维数;
G—粗糙表面的特征尺度系数;
γ—谱密度的尺寸参数;
φn—随机相位;
n—频率指数;
对于独立的微凸体形貌表述,可以根据单频率(γn=1/l)的W-M函数来获得,独立微凸体的形貌如图2所示:
数学公式为:
对于一个给定接触点截面面积为a′的微凸体,其变形量δ由z0(x)在x=0处求得:
在微观尺度下,可以认为a′=l2,微凸体的顶点曲率半径为:
由Hertz理论可知单个微凸体的弹性载荷为Fe=4Er3/3R,式中r为真实接触面积半径
将其带入Fe=4Er3/3R,可得
在MB模型中,微凸体接触点数与截面面积关系为
其中a′L为最大微接触点平截面积
根据刚度定义可得到单个微凸体的法向接触刚度为:
结合微凸体截面积分布函数积分得到结合部的总刚度为:
根据以上公式,在MATLAB中编写出节点刚度计算公式,并对各节点刚度进行加和计算。
所述步骤三具体为:
果蝇优化算法作为常用优化方法,又被称之为FOA算法,可以总结为以下几个重要的步骤:
1)随机初始果蝇群体位置X,Y;
2)配置独立果蝇使用嗅觉寻找食物的随机搜索方向与距离,rand()是寻优过程中随机生成的一个0-1之间的一个值:
Xi=X+rand(),Yi=Y+rand()
3)估算果蝇与搜索原点之间的距离D,之后计算味道浓度的判定值S,S为距离参数D的倒数值:
D=(Xi2+Yi2)0.5,S=1/D
4)将味道浓度判定数值引用到味道浓度判定函数来求求独立果蝇个体的位置味道浓度smell:
smell=Function(S)
5)寻到得到果蝇群体中所有独立果蝇味道浓度最大的果蝇:
[bestsmellbestindex]=max(smell)
6)保留味道浓度最大的果蝇的最佳味道浓度值与X、Y的坐标,之后果蝇种群会利用视觉朝向该最大值位置飞去:
smellbest=bestsmell
X=X(bestindex)
Y=Y(bestindex)
7)进入果蝇算法迭代寻优,重复之上的步骤2)~5),并且判断味道的浓度是否优于之前迭代的味道浓度,若是则进行步骤6);
所述步骤四具体为:
进入isight软件中optimization组块,在此插入果蝇优化算法,设置总迭代次数、粒子群数量等参数,选择变量为long和width,并设置变化区间,设置约束和优化目标参数。MATLAB输出总刚度结果计入optimization,optimization变换long和width的数值,等于变化了一个果蝇的飞行方向,循环迭代进入ANSYS输入,optimization会根据总刚度值的结果变换long和width的变化方向,使long和width的变化向着总刚度最大的目标变化,等于变化了果蝇群体的飞行方向。Isight整体设置流程如图3所示:
基于isight、ANSYS和MATLAB软件和果蝇优化算法的螺栓间距的研究,可以有效的对给定的参数在一定区间内进行寻优,筛选出最优的间距距离,使结合部刚度达到最大值。

Claims (5)

1.基于果蝇优化算法对栓接结合部间距的优化方法,其特征在于:包括如下步骤,
步骤一、在ANSYS中对所建立的初始栓接结合部模型进行有限元分析,获取结合部各节点的节点编号和节点压力,同时生成ANSYS的命令流文件;
步骤二、以分形接触理论为基础,运用MATLAB软件,编写计算程序,计算各节点的接触刚度;
步骤三、基于果蝇优化算法,采用isight集合功能,将ANSYS和MATLAB软件集合到isight中,确定优化间距长度(long)和宽度(width)的范围、粒子群大小、迭代次数和最终的优化目标等参数,不同的参数设定直接影响优化的能力,进行循环迭代,寻找最大的结合部刚度值。
2.根据权利要求1所述的基于果蝇优化算法对栓接结合部间距的优化方法,其特征在于:所述步骤一具体为,
建立栓接结合部模型,进行有限元划分,施加约束和模拟螺栓预紧力,进行有限元静力分析,提取各节点的节点压力。
3.根据权利要求1所述的基于果蝇优化算法对栓接结合部间距的优化方法,其特征在于:所述步骤二具体为,
利用结构函数法计算出结合部的分形参数,结合部栓接表面法的轮廓曲线利用Weierstrass-Mandelbrot函数表示
式中z—粗糙表面轮廓高度;
x—表面采样长度坐标;
D—轮廓曲线的分形维数;
G—粗糙表面的特征尺度系数;
γ—谱密度的尺寸参数;
φn—随机相位;
n—频率指数;
对于独立的微凸体形貌表述,根据单频率(γn=1/l)的W-M函数来获得,独立微凸体Z0(x)数学公式为:
对于一个给定接触点截面面积为a′的微凸体,其变形量δ由z0(x)在x=0处求得:
在微观尺度下,认为a′=l2,微凸体的顶点曲率半径为:
由Hertz理论可知单个微凸体的弹性载荷为Fe=4Er3/3R,式中r为真实接触面积半径
将其带入Fe=4Er3/3R,可得
在MB模型中,微凸体接触点数与截面面积关系为
其中a′L为最大微接触点平截面积
根据刚度定义可得到单个微凸体的法向接触刚度为:
结合微凸体截面积分布函数积分得到结合部的总刚度为:
根据以上公式,在MATLAB中编写出节点刚度计算公式,并对各节点刚度进行加和计算。
4.根据权利要求1所述的基于果蝇优化算法对栓接结合部间距的优化方法,其特征在于:所述步骤三具体为,
果蝇优化算法通过以下几个步骤实现:
1)随机初始果蝇群体位置X,Y;
2)配置独立果蝇使用嗅觉寻找食物的随机搜索方向与距离,rand()是寻优过程中随机生成的一个0-1之间的一个值:
Xi=X+rand(),Yi=Y+rand()
3)估算果蝇与搜索原点之间的距离D,之后计算味道浓度的判定值S,S为距离参数D的倒数值:
D=(Xi2+Yi2)0.5,S=1/D
4)将味道浓度判定数值引用到味道浓度判定函数来求求独立果蝇个体的位置味道浓度smell:
smell=Function(S)
5)寻到得到果蝇群体中所有独立果蝇味道浓度最大的果蝇:
[bestsmellbestindex]=max(smell)
6)保留味道浓度最大的果蝇的最佳味道浓度值与X、Y的坐标,之后果蝇种群会利用视觉朝向该最大值位置飞去:
smellbest=bestsmell
X=X(bestindex)
Y=Y(bestindex)
7)进入果蝇算法迭代寻优,重复之上的步骤2)~5),并且判断味道的浓度是否优于之前迭代的味道浓度,若是则进行步骤6)。
5.根据权利要求1所述的基于果蝇优化算法对栓接结合部间距的优化方法,其特征在于:所述步骤四具体过程如下,
进入isight软件中optimization组块,在此插入果蝇优化算法,设置总迭代次数、粒子群数量等参数,选择变量为long和width,并设置变化区间,设置约束和优化目标参数;MATLAB输出总刚度结果计入optimization,optimization变换long和width的数值,等于变化了一个果蝇的飞行方向,循环迭代进入ANSYS输入,optimization会根据总刚度值的结果变换long和width的变化方向,使long和width的变化向着总刚度最大的目标变化,等于变化了果蝇群体的飞行方向。
CN201910320600.9A 2019-04-20 2019-04-20 基于果蝇优化算法对栓接结合部间距的优化方法 Pending CN110046437A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910320600.9A CN110046437A (zh) 2019-04-20 2019-04-20 基于果蝇优化算法对栓接结合部间距的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910320600.9A CN110046437A (zh) 2019-04-20 2019-04-20 基于果蝇优化算法对栓接结合部间距的优化方法

Publications (1)

Publication Number Publication Date
CN110046437A true CN110046437A (zh) 2019-07-23

Family

ID=67278230

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910320600.9A Pending CN110046437A (zh) 2019-04-20 2019-04-20 基于果蝇优化算法对栓接结合部间距的优化方法

Country Status (1)

Country Link
CN (1) CN110046437A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111600204A (zh) * 2020-06-05 2020-08-28 杭州交联电力设计股份有限公司 一种单组电气设备的自动排布方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080300831A1 (en) * 2006-12-19 2008-12-04 Board Of Governors For Higher Education, State Of Rhode Island And Providence System and method for finite element based on topology optimization
CN105972185A (zh) * 2016-07-10 2016-09-28 北京工业大学 基于果蝇优化算法的摆线轮工作齿廓计算方法
CN106777566A (zh) * 2016-11-29 2017-05-31 湖南科技大学 基于自适应果蝇算法的开关磁阻电机结构参数优化方法
CN107367937A (zh) * 2017-08-07 2017-11-21 陕西科技大学 一种基于自适应果蝇优化算法的pid参数优化方法
CN107703491A (zh) * 2017-11-13 2018-02-16 浙江大学 基于改进的果蝇优化算法优化rbf神经网络的海杂波最优软测量仪表及方法
CN108664739A (zh) * 2018-05-14 2018-10-16 北京工业大学 基于改进粒子群算法的栓接结合部螺栓间距的优化方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080300831A1 (en) * 2006-12-19 2008-12-04 Board Of Governors For Higher Education, State Of Rhode Island And Providence System and method for finite element based on topology optimization
CN105972185A (zh) * 2016-07-10 2016-09-28 北京工业大学 基于果蝇优化算法的摆线轮工作齿廓计算方法
CN106777566A (zh) * 2016-11-29 2017-05-31 湖南科技大学 基于自适应果蝇算法的开关磁阻电机结构参数优化方法
CN107367937A (zh) * 2017-08-07 2017-11-21 陕西科技大学 一种基于自适应果蝇优化算法的pid参数优化方法
CN107703491A (zh) * 2017-11-13 2018-02-16 浙江大学 基于改进的果蝇优化算法优化rbf神经网络的海杂波最优软测量仪表及方法
CN108664739A (zh) * 2018-05-14 2018-10-16 北京工业大学 基于改进粒子群算法的栓接结合部螺栓间距的优化方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111600204A (zh) * 2020-06-05 2020-08-28 杭州交联电力设计股份有限公司 一种单组电气设备的自动排布方法
CN111600204B (zh) * 2020-06-05 2022-03-18 杭州交联电力设计股份有限公司 一种单组电气设备的自动排布方法

Similar Documents

Publication Publication Date Title
CN109726484B (zh) 基于独立连续映射法的多材料连续体结构拓扑优化设计方法
US7191161B1 (en) Method for constructing composite response surfaces by combining neural networks with polynominal interpolation or estimation techniques
CN109359353A (zh) 一种基于拓扑优化的机床立柱结构优化设计方法
Knust et al. Preform optimization for hot forging processes using genetic algorithms
CN102968524A (zh) 一种二维变曲率型材零件工艺模型的建模方法
CN104281751B (zh) 一种基于特征的涡轮冷却叶片参数化构建系统及方法
CN102682172B (zh) 基于参数分类的超临界机翼极多参数优化设计方法
Rai et al. Application of artificial neural networks to the design of turbomachinery airfoils
CN105781626B (zh) 一种大子午扩张涡轮的三维正交非定常设计方法
CN110727995A (zh) 叶片形状的构造方法、叶片以及计算机设备
CN102262692A (zh) 飞机翼面蒙皮亚音速颤振优化方法
CN108664739A (zh) 基于改进粒子群算法的栓接结合部螺栓间距的优化方法
Zhao et al. Studies on optimization of metal forming processes using sensitivity analysis methods
CN110046437A (zh) 基于果蝇优化算法对栓接结合部间距的优化方法
CN109145393A (zh) 一种结构仿生轻量化设计方法
CN110852000B (zh) 一种车身结构优化方法
CN104750948A (zh) 一种处理飞行器设计中多极值多约束问题的优化方法
CN107169171A (zh) 一种动车组头型设计方法
CN115788598B (zh) 一种涡轮叶片气膜孔参数化控制与设计方法
CN110928326A (zh) 一种面向飞行器外形的测量点差异性规划方法
CN105930574A (zh) 基于遗传及神经网络算法的内燃机进气道模型设计方法
CN106874526A (zh) 叶轮机叶片的生产坐标的生成方法和装置
CN109359360B (zh) 一种基于局部特征的结构应力处理方法
US10325413B2 (en) Generating smooth optimized parts
CN106001933B (zh) 激光切割修边线的优化方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190723