CN106650178A - 一种基于细分构型的机床龙门结构刚度的优化设计方法 - Google Patents

一种基于细分构型的机床龙门结构刚度的优化设计方法 Download PDF

Info

Publication number
CN106650178A
CN106650178A CN201710041963.XA CN201710041963A CN106650178A CN 106650178 A CN106650178 A CN 106650178A CN 201710041963 A CN201710041963 A CN 201710041963A CN 106650178 A CN106650178 A CN 106650178A
Authority
CN
China
Prior art keywords
gantry structure
particle
floor
lathe
configuration
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.)
Granted
Application number
CN201710041963.XA
Other languages
English (en)
Other versions
CN106650178B (zh
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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201710041963.XA priority Critical patent/CN106650178B/zh
Publication of CN106650178A publication Critical patent/CN106650178A/zh
Application granted granted Critical
Publication of CN106650178B publication Critical patent/CN106650178B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • 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]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • General Factory Administration (AREA)

Abstract

本发明公开了一种基于细分构型的机床龙门结构刚度的优化设计方法。建立基础龙门结构三维模型,通过细分构型的参数驱动修改调整三维模型参数获得一系列的三维模型作为结构模板,并通过有限元分析得到结构性能数据,形成机床龙门结构件样本库;根据机床加工工况的刚度要求,将设计人员输入的机床龙门结构的初始方案结构依据基础构型和细分构型进行优化处理。本发明在机床设计初期预测龙门结构件刚度,为设计高刚度龙门结构提供参考,减少机床龙门结构设计过程中的设计修改及试制实验成本。

Description

一种基于细分构型的机床龙门结构刚度的优化设计方法
技术领域
本发明属于先进设计领域,特别是涉及了一种基于细分构型的机床龙门结构刚度的优化设计方法。
技术背景
随着社会经济的发展,机床作为工作母机的需求量与日俱增,如何提高机床的性能是一个研究的热点。设计机床所要考虑的性能指标,主要包括了整机重量、结构刚度、动态响应特性、机床精度等。龙门结构是龙门式机床的关键结构,龙门结构的刚度是评价龙门式机床性能的重要指标,提高龙门结构的刚度,对于减轻机床重量、提升机床动态响应特性、提高机床精度都有重要的意义。目前国内外对于机床龙门结构的刚度设计方法,主要是通过结构设计、实验测试或仿真分析、结构修改,三个阶段循环往复的方式进行龙门结构刚度性能的逐步的改进,设计过程繁琐、设计成本高。
发明内容
为了解决背景技术中存在的问题,本发明的目的在于提供一种基于细分构型的机床龙门结构刚度的优化设计方法,简化龙门结构的设计流程和降低龙门结构的设计成本。
为达到上述目的,本发明采用了以下技术方案:
(1)建立带有结构性能数据及可调参数的机床龙门结构件样本库
建立典型机床龙门结构形式的基础龙门结构三维模型,通过参数驱动修改调整三维模型参数获得一系列的三维模型作为结构模板,并通过有限元分析得到结构性能数据,形成机床龙门结构件样本库;
(2)根据机床龙门结构件样本库通过细分构型对机床龙门结构进行刚度优化
根据机床加工工况的刚度要求,将设计人员输入的机床龙门结构的初始方案结构依据基础构型和细分构型进行优化处理。
细分构型包括重心位置、上下圆角、顶部肋板、侧面肋板、底部肋板几个方面。
所述的结构性能数据包括空载下静力分析后的最大变形量、最大应力值和最大应力发生位置,加载下静力分析后的最大变形量、最大应力值和最大应力发生位置,以及空载下模态分析后的一阶模态、二阶模态、三阶模态、四阶模态、五阶模态和六阶模态。
所述的可调参数是指三维模型在细分构型中的参数能够用来调整结构性能数据的变形参数。可调参数具体包括立柱的长宽高、横梁的长宽高、上下圆角倒角半径、顶部肋板类型、侧面肋板类型、底部肋板类型、顶部肋板参数、侧面肋板参数和底部肋板参数。
所述的典型的机床龙门结构形式包括定梁动柱、定梁定柱、动梁动柱和动梁定柱的四种基础构型,具体是以横梁和立柱是否可以运动为标准将机床龙门结构划分为定梁动柱、定梁定柱、动梁动柱和动梁定柱的四种基础龙门结构形式。
参数驱动修改是指针对机床龙门结构中上圆角、下圆角、侧面肋板、顶部肋板和底部肋板的五种细分构型的参数修改。
所述的有限元分析包括空载下的静力分析、加载下的静力分析和空载下的模态分析。
所述步骤(1)中具体包括以下步骤:
(1.1)建立定梁动柱、定梁定柱、动梁动柱和动梁定柱四种基础构型的四种基础龙门结构三维模型,通过对基础龙门结构三维模型中的上圆角、下圆角、侧面肋板、顶部肋板和底部肋板的五种参数调整初步采样来获得不同参数组合下的不同三维模型,然后针对每个三维模型通过空载下的静力分析、加载下的静力分析和空载下的模态分析获得每个三维模型对应的结构性能数据,由此获得初步三维模型及其结构性能数据;
(1.2)在可调参数相近的两个初步三维模型之间构建插入新模型,采用以下公式计算新模型的结构性能数据L(x):
式中,(u0,v0)与(u1,v1)为可调参数相近的两个初步三维模型,其中u表示初步三维模型的参数,v表示初步三维模型的结构性能数据;
然后将新模型和初步三维模型及其各自的结构性能数据组成机床龙门结构样本库,从而通过新模型对机床龙门结构样本库进行扩展。
具体实施中初步三维模型的参数u即是以机床龙门结构作为粒子,以机床龙门结构中细分构型的可调参数作为粒子的编码,可调参数的种类和粒子编码的维数相同均为N,粒子编码的维数(即可调参数的种类)包括针对肋板编码的肋板类型、肋板的长宽高和肋板间距的五维以及针对圆角编码的圆角位置和弧度的两维。
设计人员输入的机床龙门结构的初始方案结构包括龙门结构的三维模型、材料、立柱的长宽高、横梁的长宽高、上下圆角倒角半径、顶部肋板类型、侧面肋板类型、底部肋板类型、顶部肋板参数、侧面肋板参数和底部肋板参数。
机床加工工况包括机床加工方式、龙门受力大小和受力方向。
所述步骤(2)具体是:判断需刚度优化的龙门结构归属的基础构型,依据基础构型确定龙门结构可优化的细分构型,利用机床龙门结构件样本库,调整细分构型的可调参数对机床龙门结构的刚度进行优化。
细分构型的可调参数包括龙门结构的重心位置、上下圆角、顶部肋板、侧面肋板和底部肋板的几个方面。
针对所述龙门结构的重心位置的优化是是指通过对龙门结构的三维模型进行分析,依据龙门立柱的基础构型计算自重,调整龙门结构中立柱的长宽高和横梁的长宽高使龙门结构重心在自身的几何中心,即使得龙门结构重心位置在竖直中心面(横梁宽度方向二分之一处)和水平中心面(立柱高度方向二分之一处)。
所述的上下圆角、顶部肋板、侧面肋板和底部肋板的优化具体是:
针对所述上下圆角的优化是调整上下圆角倒角半径以加强龙门立柱的刚度。
针对所述顶部肋板、侧面肋板和底部肋板的优化均是调整肋板类型和肋板参数以加强龙门立柱的刚度。
肋板类型是指肋板形状,包括一字型、米字型、V字型、M字型、X字型,编码值分别为0、1、2、3、4;肋板参数包括肋板的长宽高和肋板间距。
所述的上下圆角、顶部肋板、侧面肋板和底部肋板的优化具体是:
1)利用动态粒子群算法,以机床龙门结构作为粒子,以机床龙门结构中细分构型的可调参数作为粒子的编码,可调参数的种类和粒子编码的维数相同均为N,粒子编码的维数(即可调参数的种类)包括针对肋板编码的肋板类型、肋板的长宽高和肋板间距的五维以及针对圆角编码的圆角位置和弧度的二维;
2)将初始方案结构作为初始粒子,采用以下公式进行粒子进化:
式中,为粒子x第t次迭代后的的编码值;为粒子x第t+1次迭代后的编码值;为粒子x第t次迭代后编码值的增量;为粒子x第t+1次迭代后编码值的增量;为粒子x第t次迭代后的最优位置,即该粒子在迭代过程中结构性能最优的编码值;为第t次迭代后所有粒子中的最优位置,即所有粒子在迭代过程中结构性能最优的编码值;w为惯性权重,即编码增量在迭代过程中的权重,取值区间为[0,1];r1、r2分别表示粒子在进化过程中对于自身最优位置、群体最优位置的增量权重,具体为[0,1]区间随机数。
3)利用机床龙门结构件样本库计算进化后的粒子所表示的机床龙门结构的结构性能数据;
4)进行粒子的优劣的判断,将进化后的粒子所表示机床龙门结构的结构性能数据和机床加工工况的刚度要求进行对比,得到粒子所代表的龙门结构刚度匹配度是否满足要求的结果,若匹配度满足要求则以进化后的粒子中编码所代表的可调参数作为优化设计后的机床龙门结构的参数,若匹配度不满足要求则迭代重复上述步骤2)~3)直到匹配度满足要求。
所述步骤4)中,进化后的粒子所表示机床龙门结构的结构性能数据高于等于机床加工工况的刚度要求1%,则匹配度满足要求,否则匹配度不满足要求。
所述步骤3)具体为:采用以下公式的欧氏距离计算进化后的粒子与机床龙门结构件样本库中各个三维模型对应的粒子之间的相似度:
式中,d(x,y)为粒子间的欧氏距离,α为布尔值,编码标识位相同时取1,不相同时取0;xi为粒子编码第i位的值;yi为机床龙门结构件样本库中各个三维模型对应的粒子第i位的编码值;i表示粒子编码的第i位,n表示粒子的总维度,即粒子编码的位数。
再取相似度最高的若干个粒子(具体实施中可取4-5),用该若干粒子对应的三维模型的结构性能数据采用以下加权公式计算获得进化后粒子的结构性能数据:
式中,d(x,ai)为进化后粒子与相似度最高若干个粒子中的第i个粒子之间的相似度,f(ai)为相似度最高若干个粒子中的第i个粒子的结构性能数据。
本发明的有益效果是:
本发明方法在机床龙门结构设计的初期进行,简化龙门结构的设计流程和降低龙门结构的设计成本,为设计符合刚度需求的龙门结构件提供设计依据,解决了现有技术中龙门结构刚度性能改进繁琐、成本高的问题。
附图说明
图1是本发明的方法流程图。
图2是本发明细分构型的示意图。
图3是本发明机床龙门结构件样本库构建流程图。
图4是本发明动态粒子编码方式示意图。
具体实施方式
下面结合附图和实施例对本发明做进一步说明。
如图1所示,本发明的实施例如下:
(1)建立带有结构性能数据及可调参数的机床龙门结构件样本库
建立典型机床龙门结构形式的基础龙门结构三维模型,通过参数驱动修改和改变模型结构调整模型参数获得一系列的三维模型作为结构模板,并通过有限元分析得到结构性能数据,形成机床龙门结构件样本库;
步骤(1)中具体包括以下步骤:
(1.1)建立定梁动柱、定梁定柱、动梁动柱和动梁定柱四种基础构型的四种基础龙门结构三维模型;
(1.2)以定梁定柱式龙门结构为例,先针对基础龙门结构三维模型通过空载下的静力分析、加载下的静力分析和空载下的模态分析,分析的参数设置为:杨氏模量1.5*1011,泊松比0.28,密度7200kg/m3
获得每个三维模型对应的结构性能数据,结构性能数据包括空载下静力分析后的最大变形量、最大应力值和最大应力发生位置,加载下静力分析后的最大变形量、最大应力值和最大应力发生位置,以及空载下模态分析后的一阶模态、二阶模态、三阶模态、四阶模态、五阶模态和六阶模态。其中
(a)空载下的静力分析的最大变形量、最大应力值和最大应力发生位置,得到空载下的静力分析结果,立柱最大变形发生在立柱顶端面上,最大变形量为19um,而应力最大则发生在立柱内侧拐角处,应力最大值为2.547MPa
(b)加载下的静力分析的最大变形量、最大应力值和最大应力发生位置,得到加载下的静力分析结果,加载力为11000N,作用于龙门结构横梁中间位置,立柱最大变形发生在立柱外侧的滑轨上,最大变形量为76.4um,而应力最大则发生在立柱内侧拐角处,应力最大值为7.882MPa
(c)空载下的模态分析的一阶模态、二阶模态、三阶模态、四阶模态、五阶模态和六阶模态,得到空载下的一阶模态为54.007Hz,立柱沿着X轴摆动;二阶模态为67.149Hz,立柱沿着Z摆动;三阶模态为95.77Hz,立柱沿着Y轴扭转;四阶模态为105.95Hz,立柱外侧板的变形变大;五阶模态为106.64Hz,立柱外侧板的变形变化更加剧烈;六阶模态为114.58Hz。
由此获得初步三维模型及其结构性能数据。
(1.3)在通过对基础龙门结构三维模型中的上圆角、下圆角、侧面肋板、顶部肋板和底部肋板的五种参数进行参数驱动修改,初步采样来获得不同参数组合下的不同三维模型。参数驱动修改可以是单个参数的修改,也可以是几个参数的组合修改,所修改参数及需记录的数据如图3所示,对于修改后的模型重复上述(1.2)步骤,完成整个系列结构的有限元分析,记录数据。
(1.4)需要运用拉格朗日插值公式,在可调参数相近的两个初步三维模型之间构建插入新模型,并计算新模型的结构性能数据。
通过新模型的结构性能数据为参数相近的两个样本库中三维模型的结构性能数据的均值,如顶部肋板厚度通过有限元分析得到了30mm、40mm的结构性能数据,通过拉格朗日插值公式得到35mm的结构性能数据为30mm、40mm两组数据的平均值。最后将新模型和初步三维模型及其各自的结构性能数据组成机床龙门结构样本库,从而通过新模型对机床龙门结构样本库进行扩展。
(2)根据机床龙门结构件样本库通过细分构型对机床龙门结构刚度优化
设计人员输入机床龙门结构的初始方案结构,初始方案结构包括龙门结构的三维模型,材料为铸铁,立柱的长1350mm、宽700mm,高3020mm,横梁的长3400mm、宽780mm、高445mm,上下圆角倒角半径均为100mm、顶部肋板类型、侧面肋板类型、底部肋板类型、顶部肋板参数、侧面肋板参数和底部肋板参数。实施实例的初始方案结构在加工过程中所受切削抗力为11000N,作用于横梁中间位置,原有构型在空载下最大形变为19.6μm,加载下的最大变形量为76.4um,一阶模态为54.007Hz。
机床加工工况的刚度要求为空载下最大变形小于18.6μm,加载下的最大变形量小于72um,一阶模态大于70HZ。
因此接下来依据基础构型和细分构型进行优化处理:
需刚度优化的龙门结构归属的基础构型为定梁定柱式,然后从龙门结构的重心位置、上下圆角、顶部肋板、侧面肋板和底部肋板的几个方面的细分构型对机床龙门结构的刚度进行优化。如图2所示,包括上圆角1、下圆角2、顶部肋板3、侧面肋板4、底部肋板5。
针对龙门结构的重心位置的优化是通过对龙门结构的三维模型进行分析,依据龙门立柱的基础构型计算自重,调整龙门结构中立柱的长宽高和横梁的长宽高使龙门结构重心在自身的几何中心,即使得龙门结构重心位置在竖直中心面(横梁宽度方向二分之一处)和水平中心面(立柱高度方向二分之一处)。初始方案的结构重心位置在自身的几何中心,不需要进行调整。
针对上下圆角的优化是调整上下圆角倒角半径,针对顶部肋板、侧面肋板和底部肋板的优化均是调整肋板类型和肋板参数,优化具体是:
(2.1)以机床龙门结构作为粒子,以机床龙门结构中细分构型的可调参数作为粒子的编码,可调参数的种类和粒子编码的维数相同均为N,粒子编码的维数(即可调参数的种类)包括针对肋板编码的肋板类型、肋板的长宽高和肋板间距的五维以及针对圆角编码的圆角位置和弧度的两维。
具体实施中,初始方案结构的龙门结构为定梁定柱式,在这五个细分构型均可进行优化,编码共19位排列方式如图4所示。
(2.2)提取初始方案结构的结构参数进行编码作为初始粒子,对初始粒子编码的各个数值位的值进行随机的增减,创建20个新的编码,形成初始的粒子种群,进行粒子的进化;
(2.3)利用机床龙门结构件样本库计算进化后的粒子所表示的机床龙门结构的结构性能数据;
(2.4)将进化后的粒子所表示机床龙门结构的结构性能数据和机床加工工况的刚度要求进行对比,得到粒子所代表的龙门结构刚度匹配度是否满足要求的结果,
若进化后的粒子所表示机床龙门结构的结构性能数据高于等于机床加工工况的刚度要求1%,则匹配度满足要求,以进化后的粒子中编码所代表的可调参数作为优化设计后的机床龙门结构的参数。否则匹配度不满足要求则迭代重复上述步骤直到匹配度满足要求。
具体实施经多次迭代进化后,得到满足匹配度的粒子,考虑制造加工,对编码的各项参数进行取整修正,得到当机床龙门结构编码在顶部肋板为一字型,长度为3200mm,宽度为80mm,高度为440mm,间距为40mm;侧面肋板为一字型,长度为1080mm,宽度为80mm,高度为220mm,间距为60mm;底部肋板为一字型,长度为3200mm,宽度为60mm,高度为560mm,间距为40mm;立柱上端内侧倒角半径取R=200mm,下端内侧倒角半径取R=300mm时达到了设计要求。以上述参数构建模型,进行仿真分析后得到该机床龙门结构的结构性能数据是空载下最大变形量为18.2μm,加载下的最大变形量为71.3um一阶模态为72HZ,符合设计需求。

Claims (10)

1.一种基于细分构型的机床龙门结构刚度的优化设计方法,包括以下步骤:
(1)建立带有结构性能数据及可调参数的机床龙门结构件样本库:建立基础龙门结构三维模型,通过细分构型的参数驱动修改调整三维模型参数获得一系列的三维模型作为结构模板,并通过有限元分析得到结构性能数据,形成机床龙门结构件样本库;
(2)根据机床龙门结构件样本库通过细分构型对机床龙门结构刚度进行优化:根据机床加工工况的刚度要求,将设计人员输入的机床龙门结构的初始方案结构依据基础构型和细分构型进行优化处理。
2.根据权利要求1所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述的结构性能数据包括空载下静力分析后的最大变形量、最大应力值和最大应力发生位置,加载下静力分析后的最大变形量、最大应力值和最大应力发生位置,以及空载下模态分析后的一阶模态、二阶模态、三阶模态、四阶模态、五阶模态和六阶模态。
3.根据权利要求1所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述的典型的机床龙门结构形式包括定梁动柱、定梁定柱、动梁动柱和动梁定柱的四种基础构型。
4.根据权利要求1所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述的细分构型包括龙门结构的重心位置、上下圆角、顶部肋板、侧面肋板和底部肋板的几个方面,所述的参数驱动修改是指针对机床龙门结构中上圆角、下圆角、侧面肋板、顶部肋板和底部肋板的五种细分构型的参数修改。
5.根据权利要求1所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述步骤(1)中具体包括以下步骤:
(1.1)建立定梁动柱、定梁定柱、动梁动柱和动梁定柱四种基础构型的四种基础龙门结构三维模型,通过对基础龙门结构三维模型中的上圆角、下圆角、侧面肋板、顶部肋板和底部肋板的五种参数调整初步采样来获得不同参数组合下的不同三维模型,然后针对每个三维模型通过空载下的静力分析、加载下的静力分析和空载下的模态分析获得每个三维模型对应的结构性能数据,由此获得初步三维模型及其结构性能数据;
(1.2)在可调参数相近的两个初步三维模型之间构建插入新模型,采用以下公式计算新模型的结构性能数据L(x):
L ( u ) = v 0 + v 1 - v 0 u 1 - u 0 ( u - u 0 )
式中,(u0,v0)与(u1,v1)为可调参数相近的两个初步三维模型,其中u表示初步三维模型的参数,v表示初步三维模型的结构性能数据;
然后将新模型和初步三维模型及其各自的结构性能数据组成机床龙门结构样本库。
具体实施中初步三维模型的参数u即是以机床龙门结构作为粒子,以机床龙门结构中细分构型的可调参数作为粒子的编码,可调参数的种类和粒子编码的维数相同均为N,粒子编码的维数(即可调参数的种类)包括针对肋板编码的肋板类型、肋板的长宽高和肋板间距的五维以及针对圆角编码的圆角位置和弧度的两维。
6.根据权利要求1所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述步骤(2)具体是:判断需刚度优化的龙门结构归属的基础构型,依据基础构型确定龙门结构可优化的细分构型,利用机床龙门结构件样本库,调整细分构型的可调参数对机床龙门结构的刚度进行优化,包括重心位置的优化以及上下圆角、顶部肋板、侧面肋板和底部肋板的优化。
7.根据权利要求6所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:针对所述龙门结构的重心位置的优化是是指通过对龙门结构的三维模型进行分析,依据龙门立柱的基础构型计算自重,调整龙门结构中立柱的长宽高和横梁的长宽高使龙门结构重心在自身的几何中心。
8.根据权利要求6所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述的上下圆角、顶部肋板、侧面肋板和底部肋板的优化具体是:针对所述上下圆角的优化是调整上下圆角倒角半径以加强龙门立柱的刚度,针对所述顶部肋板、侧面肋板和底部肋板的优化均是调整肋板类型和肋板参数以加强龙门立柱的刚度。
9.根据权利要求8所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述的上下圆角、顶部肋板、侧面肋板和底部肋板的优化具体是:
1)以机床龙门结构作为粒子,以机床龙门结构中细分构型的可调参数作为粒子的编码,可调参数的种类和粒子编码的维数相同均为N,粒子编码的维数(即可调参数的种类)包括针对肋板编码的肋板类型、肋板的长宽高和肋板间距的五维以及针对圆角编码的圆角位置和弧度的两维;
2)将初始方案结构作为初始粒子,采用以下公式进行粒子进化:
V q t + 1 = wV q t + r ( P q t - X q t ) + r 2 ( P G t - X q t )
X q t + 1 = X q t + V q t + 1
式中,为粒子x第t次迭代后的的编码值;为粒子x第t+1次迭代后的编码值;为粒子x第t次迭代后编码值的增量;为粒子x第t+1次迭代后编码值的增量;为粒子x第t次迭代后的最优位置,即该粒子在迭代过程中结构性能最优的编码值;为第t次迭代后所有粒子中的最优位置,即所有粒子在迭代过程中结构性能最优的编码值;w为惯性权重,即编码增量在迭代过程中的权重,取值区间为[0,1];;r1、r2分别表示粒子在进化过程中相对于自身最优位置、群体最优位置的增量权重,具体为[0,1]区间随机数。
3)利用机床龙门结构件样本库计算进化后的粒子所表示的机床龙门结构的结构性能数据;
4)将进化后的粒子所表示机床龙门结构的结构性能数据和机床加工工况的刚度要求进行对比,得到粒子所代表的龙门结构刚度匹配度是否满足要求的结果,若匹配度满足要求则以进化后的粒子中编码所代表的可调参数作为优化设计后的机床龙门结构的参数,若匹配度不满足要求则迭代重复上述步骤2)~3)直到匹配度满足要求。
10.根据权利要求8所述的一种基于细分构型的机床龙门结构刚度的优化设计方法,其特征在于:所述步骤3)具体为:采用以下公式计算进化后的粒子与机床龙门结构件样本库中各个三维模型对应的粒子之间的相似度:
d ( x , y ) = α ( x 1 - y 1 ) 2 + α ( x 2 - y 2 ) 2 + ... + α ( x n - y n ) 2 = Σ i = 1 n α ( x i - y i ) 2
式中,d(x,y)为粒子间的欧氏距离,α为布尔值,编码标识位相同时取1,不相同时取0;xi为粒子编码第i位的值;yi为机床龙门结构件样本库中各个三维模型对应的粒子第i位的编码值;i表示粒子编码的第i位,n表示粒子的总维度,即粒子编码的位数;
再取相似度最高的若干个粒子,用该若干粒子对应的三维模型的结构性能数据采用以下加权公式计算获得进化后粒子的结构性能数据:
f ( x ) = Σ i = 1 n d ( x , a i ) f ( a i ) Σ i = 1 n d ( x , a i )
式中,d(x,ai)为进化后粒子与相似度最高若干个粒子中的第i个粒子之间的相似度,f(ai)为相似度最高若干个粒子中的第i个粒子的结构性能数据。
CN201710041963.XA 2017-01-20 2017-01-20 一种基于细分构型的机床龙门结构刚度的优化设计方法 Active CN106650178B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710041963.XA CN106650178B (zh) 2017-01-20 2017-01-20 一种基于细分构型的机床龙门结构刚度的优化设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710041963.XA CN106650178B (zh) 2017-01-20 2017-01-20 一种基于细分构型的机床龙门结构刚度的优化设计方法

Publications (2)

Publication Number Publication Date
CN106650178A true CN106650178A (zh) 2017-05-10
CN106650178B CN106650178B (zh) 2019-11-19

Family

ID=58842071

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710041963.XA Active CN106650178B (zh) 2017-01-20 2017-01-20 一种基于细分构型的机床龙门结构刚度的优化设计方法

Country Status (1)

Country Link
CN (1) CN106650178B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090252A (zh) * 2017-11-27 2018-05-29 华侨大学 一种数控铣床横梁的优化设计方法及优化设计的横梁
CN108363827A (zh) * 2017-12-28 2018-08-03 清华大学 一种基于相似理论的混联机床静刚度分析方法
CN110059382A (zh) * 2019-04-04 2019-07-26 北京工业大学 一种非均匀加筋层的重型机床复合基础设计方法
CN113312701A (zh) * 2021-04-30 2021-08-27 中铝材料应用研究院有限公司 基于拓扑及尺寸优化的全铝客车车身门立柱结构设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004171227A (ja) * 2002-11-19 2004-06-17 Fujitsu Ltd 設計支援装置、設計支援方法、およびプログラム
CN103927411A (zh) * 2014-04-01 2014-07-16 深圳市大族激光科技股份有限公司 一种龙门式激光切割机床的设计方法
CN104156501A (zh) * 2014-07-15 2014-11-19 天津大学 一种机床整机静刚度的优化设计方法
CN105528503A (zh) * 2016-02-17 2016-04-27 中国科学院沈阳自动化研究所 一种基于结构分解的大型构件动态优化设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004171227A (ja) * 2002-11-19 2004-06-17 Fujitsu Ltd 設計支援装置、設計支援方法、およびプログラム
CN103927411A (zh) * 2014-04-01 2014-07-16 深圳市大族激光科技股份有限公司 一种龙门式激光切割机床的设计方法
CN104156501A (zh) * 2014-07-15 2014-11-19 天津大学 一种机床整机静刚度的优化设计方法
CN105528503A (zh) * 2016-02-17 2016-04-27 中国科学院沈阳自动化研究所 一种基于结构分解的大型构件动态优化设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
朱金波等: "《基于多目标优化的龙门加工中心动横梁筋板结构优化设计方法》", 《组合机床与自动化加工技术》 *
李博: "《粒子群算法及其在桁架结构优化设计中的应用研究》", 《价值工程》 *
王一江,张晓艳,庞学慧: "《基于ANSYS的龙门铣床龙门结构模态分析》", 《计算机应用技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108090252A (zh) * 2017-11-27 2018-05-29 华侨大学 一种数控铣床横梁的优化设计方法及优化设计的横梁
CN108363827A (zh) * 2017-12-28 2018-08-03 清华大学 一种基于相似理论的混联机床静刚度分析方法
CN110059382A (zh) * 2019-04-04 2019-07-26 北京工业大学 一种非均匀加筋层的重型机床复合基础设计方法
CN110059382B (zh) * 2019-04-04 2022-11-29 北京工业大学 一种非均匀加筋层的重型机床复合基础设计方法
CN113312701A (zh) * 2021-04-30 2021-08-27 中铝材料应用研究院有限公司 基于拓扑及尺寸优化的全铝客车车身门立柱结构设计方法
CN113312701B (zh) * 2021-04-30 2024-03-19 中铝材料应用研究院有限公司 基于拓扑及尺寸优化的全铝客车车身门立柱结构设计方法

Also Published As

Publication number Publication date
CN106650178B (zh) 2019-11-19

Similar Documents

Publication Publication Date Title
CN106650178A (zh) 一种基于细分构型的机床龙门结构刚度的优化设计方法
CN107862141B (zh) 一种轨道车辆低噪声正向设计方法
CN109359353B (zh) 一种基于拓扑优化的机床立柱结构优化设计方法
CN101950319B (zh) 一种高速立式加工中心大跨度横梁拓扑优化设计方法
CN103995914B (zh) 基于动态特性分析的磨齿机立柱结构优化设计的方法
Zhang et al. Optimal topology design of internal stiffeners for machine pedestal structures using biological branching phenomena
CN103034766B (zh) 一种确定层合板复合材料的铺层角方向与厚度的方法
CN106202841B (zh) 一种机床床身三点支撑位置的优化设计方法
CN108416158A (zh) 一种复合材料机床床身内部结构优化设计方法
CN104765912A (zh) 一种铝板冲压工艺的稳健性优化方法
CN107423512A (zh) 一种大型压铸机尾板的轻量化设计方法
CN102622486A (zh) 考虑拔模制造约束的拓扑优化设计方法
CN110580363B (zh) 一种搅拌摩擦焊机器人底座结构的拓扑优化设计方法
CN108038308A (zh) 一种铝合金压铸减震塔的结构设计方法
CN103310064B (zh) 一种采用极端尺寸调整的数控机床横梁结构优化设计方法
CN102430712B (zh) 一种空间复杂曲面芯盒模具的制造方法
CN108256215A (zh) 一种基于结构拓扑优化的龙门加工中心滑枕及其优化设计方法
CN107025354A (zh) 一种基于极差分析的车窗升降板成形工艺优化方法
CN103279614B (zh) 一种机床结构大件筋板布局及设计方法
CN112884198B (zh) 结合门限回归和改进支持向量机面板坝坝顶沉降预测方法
CN106202686A (zh) 一种涡轮盘等温模锻预成形坯料的多目标设计方法
CN103207938B (zh) 一种重型龙门机床超跨距横梁垂直导轨面起拱曲线设计方法
CN112100747A (zh) 车身骨架拓扑优化方法、装置、设备及介质
CN107092745A (zh) 一种基于方差分析的车窗升降板成形工艺优化方法
CN208467737U (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
GR01 Patent grant
GR01 Patent grant