CN106951610A - 一种基于近似模型的插秧机秧箱结构优化方法 - Google Patents
一种基于近似模型的插秧机秧箱结构优化方法 Download PDFInfo
- Publication number
- CN106951610A CN106951610A CN201710128019.8A CN201710128019A CN106951610A CN 106951610 A CN106951610 A CN 106951610A CN 201710128019 A CN201710128019 A CN 201710128019A CN 106951610 A CN106951610 A CN 106951610A
- Authority
- CN
- China
- Prior art keywords
- seedling
- optimization
- design
- variable
- rigidity
- 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
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
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)
- Transplanting Machines (AREA)
Abstract
本发明公开了一种插秧机秧箱结构优化方法,对秧箱结构进行初始性能的有限元分析,确定刚度富裕点;对影响秧箱动静态性能的关键因素进行灵敏度分析,确定设计变量,以秧箱的模态、刚度作为优化设计的约束条件,质量作为目标函数,构建秧箱结构优化数学模型;采用优化拉丁超立方试验设计产生样本点,选用精度较高的Kriging模型,构建质量、模态和刚度与设计变量关系的近似模型,对优化数学模型进行简化;采用智能优化算法对优化数学模型进行求解,得出质量最小值和对应的设计变量取值,并对优化后的秧箱静动态性能进行校验分析。本发明通过构建刚度、模态和质量与设计变量关系的近似模型,能够减少计算复杂性,因此可加快优化进程。
Description
技术领域:
本发明涉及插秧机结构设计技术领域,主要涉及一种基于近似模型的插秧机秧箱结构优化方法。
背景技术:
插秧机秧箱的轻量化设计属于典型的高维、非线性多学科优化问题,其基本原则是在满足静刚度和振动约束的前提下使秧箱质量最小。利用传统的计算机辅助设计方法,容易出现计算复杂或性能设计过剩等问题。本发明通过结合灵敏度分析与近似模型方法,确定了优化目标、约束条件和设计变量,将多目标优化问题转化为单目标优化问题,解决了仿真精度与计算复杂性的矛盾。
发明内容:
本发明目的就是为了弥补已有技术的缺陷,提供一种基于近似模型的插秧机秧箱结构优化方法,实现结构的轻量化,在满足结构静动态性能的前提下,使其质量达到最小,通过近似模型的应用进行有约束的优化,实现秧箱的轻量化,加快了计算进程,缩短了优化周期。解决了仿真精度与计算复杂性的矛盾;在牺牲一定精度的基础上,有效的减少了计算的复杂性。
本发明是通过以下技术方案实现的:
一种基于近似模型的插秧机秧箱结构优化方法,其特征在于:首先通过对箱结构进行初始性能的有限元分析,确定刚度富裕点;对影响秧箱静动态性能的关键因素进行灵敏度分析,确定设计变量,以秧箱的模态、刚度作为优化设计的约束条件,质量作为目标函数,构建秧箱结构优化数学模型;然后采用优化拉丁超立方试验设计产生样本点,选用精度较高的Kriging模型,构建质量、模态和刚度与设计变量关系的近似模型,对优化数学模型进行简化;最后采用智能优化算法对优化数学模型进行求解,得出质量最小值和对应的设计变量取值,并对优化后的秧箱静动态性能进行校验分析。
其具体包括以下步骤:
步骤1:秧箱结构初始性能有限元分析:包括构建秧箱结构有限元模型,进行刚度分析与模态分析,确定初始条件下秧箱的最大硬点刚度与模态第一阶固有频率,并确定刚度富裕点;
步骤2:秧箱结构优化数学模型的构建:包括确定优化目标、约束条件和设计变量,优化目标为秧箱质量最小;约束条件为满足秧箱静动态性能要求,即满足静刚度和振动约束;变量包括形状变量和尺寸变量;对变量进行质量灵敏度分析,灵敏度较大的变量作为设计变量;根据加工工艺及插秧机秧箱性能要求确定各设计变量的取值范围;
步骤3:近似模型的构建:包括确定试验设计方法、初始样本点、近似模型类型及其精度评估准则;根据步骤2的分析结果,构建以刚度、模态和质量与设计变量关系的近似模型,并进行精度评估,如不满足条件,则增加样本点重新构建近似模型,直至满足精度要求;
步骤4:秧箱结构优化:根据步骤3所构建的近似模型,通过有约束的单目标优化方法得到满足静动态性能条件下,秧箱质量的最小值及对应设计变量取值;
步骤5:秧箱静动态性能校验分析:校验秧箱设计变量最优值,如满足设计目标则设计完成,如不满足设计目标,则修改设计变量的初始取值进入步骤2,按照上述步骤构建、优化和分析,直至满足设计目标则设计完成。
所述步骤2中的秧箱形状变量主要有各个经轴架的截面宽度和高度;尺寸变量主要有各个经轴架及秧箱单元的厚度。
所述步骤3中的近似模型包括秧箱质量、刚度、模态三个响应量与设计变量关系的近似模型,刚度、模态与设计变量关系的近似模型分别拟合最大硬点刚度和秧箱的第一阶固有频率;近似模型精度评估准则为确定性系数法,取R2≥0.9。
所述步骤5中的校核验证秧箱变量最优值,包括验证秧箱刚度和模态性能。
本发明的优点是:
(1)本发明在满足结构静动态性能要求的前提下进行插秧机秧箱结构优化,摆脱了经验设计带来的不确定性及盲目性;
(2)本发明通过近似模型的应用,解决了仿真精度与计算复杂性之间的矛盾,在牺牲一定精度的基础上,可在设计变量给定范围内更快地找到最优值,大大缩短了计算时间。
附图说明:
图1为本发明的设计流程图。
图2为本发明实施案例中插秧机秧箱有限元模型。
图3为本发明实施案例中经轴架1的截面设计变量定义图。
图4为本发明实施案例中经轴架2的截面设计变量定义图。
图5为本发明实施案例中经轴架3的截面设计变量定义图。
图6为本发明实施案例中经轴架4的截面设计变量定义图。
具体实施方式:
参见附图。
下面结合附图和实施例对本发明进一步说明。
机械的轻量化有益于推进节能减排,本发明提供一种结构优化方法,实现结构的轻量化,在满足结构的动静态性能的前提下,使其质量达到最小,通过近似模型的应用进行有约束的优化,实现秧箱的轻量化,解决了仿真精度与计算复杂性之间的矛盾,加快了计算进程,缩短了优化周期。
步骤1:秧箱结构有限元分析。包括构建秧箱结构有限元模型,进行刚度分析与模态分析,确定初始条件下秧箱的最大硬点刚度与模态第一阶固有频率,并确定刚度富裕点。
步骤2:秧箱结构优化数学模型的构建。包括确定优化目标、约束条件和设计变量。优化目标为秧箱质量最小;约束条件为满足秧箱静动态性能要求,即满足静刚度和振动约束;变量包括形状变量和尺寸变量。对变量进行质量灵敏度分析,灵敏度较大的变量作为设计变量。根据加工工艺及插秧机秧箱性能要求确定各设计变量的取值范围。
步骤3:近似模型的构建。包括确定试验设计方法、初始样本点、近似模型类型及其精度评估准则。根据步骤2的分析结果,构建以刚度、模态和质量与设计变量关系的近似模型,并进行精度评估,如不满足条件,则增加样本点重新构建近似模型,直至满足精度要求。
步骤4:秧箱结构优化。根据步骤3所构建的近似模型,通过有约束的单目标优化方法得到满足静动态性能条件下,秧箱质量的最小值及对应设计变量取值。
步骤5:秧箱静动态性能校验分析。校验秧箱设计变量最优值,如满足设计目标则设计完成,如不满足设计目标,则修改设计变量的初始取值进入步骤2。
作为本发明的进一步改进,所述步骤2中的秧箱形状变量主要有:各个经轴架的截面宽度和高度;尺寸变量主要有:各个经轴架及秧箱单元的厚度。
作为本发明的进一步改进,所述步骤3中的近似模型包括秧箱质量、刚度、模态三个响应量的近似模型,刚度、模态近似模型分别拟合最大硬点刚度和秧箱的一阶固有频率;近似模型精度评估准则为确定性系数法,取R2≥0.9。
作为本发明的进一步改进,所述步骤5中的秧箱静动态性能校验,包括验证秧箱刚度和模态性能。
实施例:现以某型插秧机秧箱的结构优化设计为例进行说明,秧箱的初始设计有限元模型如图2所示,图中1为秧箱单元,2为经轴架1,3为经轴架2,4为经轴架3,5为经轴架4。
步骤1:秧箱结构有限元分析。秧箱有限元模型如图2所示。进行刚度和模态分析,分析结果为秧箱初始质量为10.631Kg,最大硬点刚度为22.712N·mm-1,第一阶固有频率为30.271Hz。
步骤2:秧箱结构优化数学模型的构建。包括确定优化目标、约束条件和设计变量,设计变量灵敏度分析。秧箱的设计变量包括形状变量和尺寸变量,形状变量为各个经轴架的截面宽度和高度,尺寸变量为经轴架及秧箱单元的厚度。通过CAE计算,秧箱设计变量的灵敏度排序为厚度(秧箱单元、经轴架2、3、4、1)、截面高度(经轴架2、1、3)、截面宽度(经轴架4、3、1)。
步骤3:近似模型的构建。采用优化拉丁超立方试验设计产生样本点,选用精度较高的Kriging模型,得到确定性系数R2均大于0.9。
步骤4:秧箱结构优化。将秧箱的模态、刚度作为优化设计的约束条件,质量作为目标函数进行优化,得到质量最小时秧箱的变量值,并进行校核,优化前后的秧箱设计变量及质量、刚度、模态值如表1所示。
表1秧箱优化前后对比
步骤5:秧箱静动态性能校验分析。优化后秧箱的质量为9.458Kg,第一阶固有频率为30.967Hz,刚度为23.205N·mm-1,秧箱的静动态性能均满足设计目标要求,质量减少了1.173Kg,实现了秧箱的轻量化。
本专利与现有技术相比:
(1)对设计变量进行灵敏度分析,选出对优化目标影响较大的变量,在一定程度上降低了优化问题的难度。
(2)考虑了秧箱实际工况下的模态与刚度,在满足实际性能需要的前提下,实现了秧箱的轻量化设计。
(3)通过构建结构响应的近似模型替代真实的结构响应,加快了优化进程,缩短了优化周期。
(4)将多目标优化问题转化为单目标优化问题,降低了优化问题的难度。以上内容是结合具体的优选实施方式对本发明所做的进一步说明,不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干简单推演或替换,都应当视为属于本发明的保护范围。
Claims (5)
1.一种基于近似模型的插秧机秧箱结构优化方法,其特征在于:首先通过对箱结构进行初始性能的有限元分析,确定刚度富裕点;对影响秧箱静动态性能的关键因素进行灵敏度分析,确定设计变量,以秧箱的模态、刚度作为优化设计的约束条件,质量作为目标函数,构建秧箱结构优化数学模型;然后采用优化拉丁超立方试验设计产生样本点,选用精度较高的Kriging模型,构建质量、模态和刚度与设计变量关系的近似模型,对优化数学模型进行简化;最后采用智能优化算法对优化数学模型进行求解,得出质量最小值和对应的设计变量取值,并对优化后的秧箱静动态性能进行校验分析。
2.根据权利要求1所述的一种基于近似模型的插秧机秧箱结构优化方法,其特征在于:具体包括以下步骤:
步骤1:秧箱结构初始性能有限元分析:包括构建秧箱结构有限元模型,进行刚度分析与模态分析,确定初始条件下秧箱的最大硬点刚度与模态第一阶固有频率,并确定刚度富裕点;
步骤2:秧箱结构优化数学模型的构建:包括确定优化目标、约束条件和设计变量,优化目标为秧箱质量最小;约束条件为满足秧箱静动态性能要求,即满足静刚度和振动约束;变量包括形状变量和尺寸变量;对变量进行质量灵敏度分析,灵敏度较大的变量作为设计变量;根据加工工艺及插秧机秧箱性能要求确定各设计变量的取值范围;
步骤3:近似模型的构建:包括确定试验设计方法、初始样本点、近似模型类型及其精度评估准则;根据步骤2的分析结果,构建以刚度、模态和质量与设计变量关系的近似模型,并进行精度评估,如不满足条件,则增加样本点重新构建近似模型,直至满足精度要求;
步骤4:秧箱结构优化:根据步骤3所构建的近似模型,通过有约束的单目标优化方法得到满足静动态性能条件下,秧箱质量的最小值及对应设计变量取值;
步骤5:秧箱静动态性能校验分析:校验秧箱设计变量最优值,如满足设计目标则设计完成,如不满足设计目标,则修改设计变量的初始取值进入步骤2,按照上述步骤构建、优化和分析,直至满足设计目标则设计完成。
3.根据权利要求2所述的插秧机秧箱结构优化,其特征在于:所述步骤2中的秧箱形状变量主要有各个经轴架的截面宽度和高度;尺寸变量主要有各个经轴架及秧箱单元的厚度。
4.根据权利要求2所述的插秧机秧箱结构优化,其特征在于:所述步骤3中的近似模型包括秧箱质量、刚度、模态三个响应量与设计变量关系的近似模型,刚度、模态与设计变量关系的近似模型分别拟合最大硬点刚度和秧箱的第一阶固有频率;近似模型精度评估准则为确定性系数法,取R2≥0.9。
5.根据权利要求2所述的插秧机秧箱结构优化,其特征在于:所述步骤5中的校核验证秧箱变量最优值,包括验证秧箱刚度和模态性能。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710128019.8A CN106951610B (zh) | 2017-03-06 | 2017-03-06 | 一种基于近似模型的插秧机秧箱结构优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710128019.8A CN106951610B (zh) | 2017-03-06 | 2017-03-06 | 一种基于近似模型的插秧机秧箱结构优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106951610A true CN106951610A (zh) | 2017-07-14 |
CN106951610B CN106951610B (zh) | 2020-12-01 |
Family
ID=59467906
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710128019.8A Active CN106951610B (zh) | 2017-03-06 | 2017-03-06 | 一种基于近似模型的插秧机秧箱结构优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106951610B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984951A (zh) * | 2018-08-08 | 2018-12-11 | 上海理工大学 | 一种优化机床整机动态特性的方法 |
CN109710957A (zh) * | 2017-10-26 | 2019-05-03 | 北京长城华冠汽车科技股份有限公司 | 多连杆铰接刚度优化方法及系统 |
CN110990939A (zh) * | 2019-10-10 | 2020-04-10 | 西北工业大学 | 防冰腔结构可靠性设计方法 |
CN111062094A (zh) * | 2019-10-08 | 2020-04-24 | 珠海格力电器股份有限公司 | 一种加工中心主轴箱的拓扑优化设计方法 |
CN116127651A (zh) * | 2023-03-09 | 2023-05-16 | 上海理工大学 | 一种大型减速器组合式箱体结构的静刚度匹配设计方法 |
CN116541975A (zh) * | 2023-04-11 | 2023-08-04 | 哈尔滨工业大学 | 一种基于代理模型的核电用塔机结构系统动态优化设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1760877A (zh) * | 2005-11-03 | 2006-04-19 | 上海交通大学 | 多项结构性能交互式的轿车车身结构轻量化方法 |
CN105608286A (zh) * | 2016-01-12 | 2016-05-25 | 上海交通大学 | 基于粒子群优化算法的车辆保险杠结构优化方法 |
CN105653768A (zh) * | 2015-12-25 | 2016-06-08 | 上海交通大学 | 基于粒子群优化算法的轿车车身轻量化结构的实现方法 |
-
2017
- 2017-03-06 CN CN201710128019.8A patent/CN106951610B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1760877A (zh) * | 2005-11-03 | 2006-04-19 | 上海交通大学 | 多项结构性能交互式的轿车车身结构轻量化方法 |
CN105653768A (zh) * | 2015-12-25 | 2016-06-08 | 上海交通大学 | 基于粒子群优化算法的轿车车身轻量化结构的实现方法 |
CN105608286A (zh) * | 2016-01-12 | 2016-05-25 | 上海交通大学 | 基于粒子群优化算法的车辆保险杠结构优化方法 |
Non-Patent Citations (5)
Title |
---|
LI DAI1 .ET AL: "The software development on kinematics Parameters Optimization of rice transplanter separating-planting mechanism", 《IEEE XPLORE》 * |
张娜娜: "高速水稻插秧机车架的轻量化设计", 《农业工程学报》 * |
朱德泉 等: "2ZGK-6 型可调宽窄行高速水稻插秧机设计与试验", 《农业工程学报》 * |
朱茂桃: "基于Kriging模型的车门刚度和模态优化", 《汽车工程》 * |
罗云峰: "基于刚度和模态的重型卡车车架轻量化设计", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109710957A (zh) * | 2017-10-26 | 2019-05-03 | 北京长城华冠汽车科技股份有限公司 | 多连杆铰接刚度优化方法及系统 |
CN109710957B (zh) * | 2017-10-26 | 2023-11-24 | 北京长城华冠汽车科技股份有限公司 | 多连杆铰接刚度优化方法及系统 |
CN108984951B (zh) * | 2018-08-08 | 2022-12-27 | 上海理工大学 | 一种优化机床整机动态特性的方法 |
CN108984951A (zh) * | 2018-08-08 | 2018-12-11 | 上海理工大学 | 一种优化机床整机动态特性的方法 |
CN111062094A (zh) * | 2019-10-08 | 2020-04-24 | 珠海格力电器股份有限公司 | 一种加工中心主轴箱的拓扑优化设计方法 |
CN111062094B (zh) * | 2019-10-08 | 2021-06-18 | 珠海格力电器股份有限公司 | 一种加工中心主轴箱的拓扑优化设计方法 |
CN110990939A (zh) * | 2019-10-10 | 2020-04-10 | 西北工业大学 | 防冰腔结构可靠性设计方法 |
CN110990939B (zh) * | 2019-10-10 | 2024-04-19 | 西北工业大学 | 防冰腔结构可靠性设计方法 |
CN116127651A (zh) * | 2023-03-09 | 2023-05-16 | 上海理工大学 | 一种大型减速器组合式箱体结构的静刚度匹配设计方法 |
CN116127651B (zh) * | 2023-03-09 | 2023-08-11 | 上海理工大学 | 一种大型减速器组合式箱体结构的静刚度匹配设计方法 |
US11995385B1 (en) | 2023-03-09 | 2024-05-28 | University Of Shanghai For Science And Technology | Matching design method of static stiffness for combined box body structure of large reducer |
CN116541975A (zh) * | 2023-04-11 | 2023-08-04 | 哈尔滨工业大学 | 一种基于代理模型的核电用塔机结构系统动态优化设计方法 |
CN116541975B (zh) * | 2023-04-11 | 2024-01-30 | 哈尔滨工业大学 | 一种基于代理模型的核电用塔机结构系统动态优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106951610B (zh) | 2020-12-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106951610A (zh) | 一种基于近似模型的插秧机秧箱结构优化方法 | |
CN111191316B (zh) | 一种基于响应面的建筑自然通风性能优化模型及优化方法 | |
CN109344524B (zh) | 一种薄板结构加强筋分布优化方法 | |
CN107330516A (zh) | 模型参数训练方法、装置及系统 | |
CN107169191A (zh) | 一种风机叶片建模方法 | |
CN107391891A (zh) | 一种基于模型融合方法的大展弦比机翼优化设计方法 | |
CN111597631B (zh) | 基于自适应代理模型的汽车风阻系数优化方法 | |
CN104750948B (zh) | 一种处理飞行器设计中多极值多约束问题的优化方法 | |
CN105893669A (zh) | 一种基于数据挖掘的全局仿真性能预测方法 | |
CN103678762A (zh) | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 | |
CN114792037B (zh) | 一种超材料隔振器的序贯稳健性优化设计方法 | |
CN105653799B (zh) | 金属隔膜贮箱一体化设计方法 | |
CN115115796A (zh) | 一种基于语义Web技术的建筑数字孪生框架及建模方法 | |
CN114676522B (zh) | 融合gan和迁移学习的气动形状优化设计方法及系统及设备 | |
CN104484527B (zh) | 一种离散结构拓扑优化过程中均布载荷自动动态修改方法 | |
CN106295015B (zh) | 一种渐开线直齿圆柱齿轮副的齿廓修形方法及与其配套的专用参数化cad系统 | |
CN105956318A (zh) | 基于改进分裂 h-k 聚类方法的风电场机群划分方法 | |
CN113221278B (zh) | 一种车载光电平台照准架轻量化方法 | |
CN109948253B (zh) | 薄板无网格Galerkin结构模态分析的GPU加速方法 | |
CN105608267A (zh) | 一种多变量全局优化算法 | |
CN110826182A (zh) | 一种基于顶点法和序贯优化策略的飞行器结构气动弹性鲁棒设计方法 | |
Yang et al. | Research on Optimization of Injection Molding Process Parameters of Automobile Plastic Front‐End Frame | |
CN106354913A (zh) | 一种用于斜齿轮传动系统参数化设计平台及使用方法 | |
Wang et al. | Variable-thickness optimization method for shell structures based on a regional evolutionary control strategy | |
US8271518B2 (en) | Nearest neighbor search method |
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 |