CN109598023A - 一种高强度车轮轮辐结构设计方法 - Google Patents
一种高强度车轮轮辐结构设计方法 Download PDFInfo
- Publication number
- CN109598023A CN109598023A CN201811298065.3A CN201811298065A CN109598023A CN 109598023 A CN109598023 A CN 109598023A CN 201811298065 A CN201811298065 A CN 201811298065A CN 109598023 A CN109598023 A CN 109598023A
- Authority
- CN
- China
- Prior art keywords
- spoke
- model
- stress
- load
- analysis
- 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
Links
- 238000013461 design Methods 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 238000013401 experimental design Methods 0.000 claims abstract description 9
- 238000011161 development Methods 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims description 23
- 238000005452 bending Methods 0.000 claims description 10
- 229910052739 hydrogen Inorganic materials 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000009826 distribution Methods 0.000 claims description 3
- 230000002068 genetic effect Effects 0.000 claims description 3
- 230000008676 import Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000005070 sampling Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 238000012356 Product development Methods 0.000 description 2
- 230000008878 coupling Effects 0.000 description 2
- 238000010168 coupling process Methods 0.000 description 2
- 238000005859 coupling reaction Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000009916 joint effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012797 qualification Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
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/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种高强度车轮轮辐结构设计方法，包括以下步骤：1）建立轮辐参数化模型；2）建立轮辐有限元模型；3）对步骤2）得到的有限元模型进行实验设计；4）建立轮辐近似模型；5）对步骤4）得到的近似模型进行优化，直到得到最优结果。本发明通过对9个关键参数的建模，可以更好地贴近现实，能从全局范围内完成所有组合方案的对比，提高了产品设计可靠性；本设计方法在有限元模型的基础上增加了近似模型的构建与优化，大大缩短了产品的开发时间，也使得产品的设计精度有的很大的提升。
Description
技术领域
本发明涉及汽车及机械工程技术领域，尤其涉及一种高强度车轮轮辐结构设计方法。
背景技术
车轮是汽车的关键部件，其重量直接影响整车油耗，近年来，客户对车轮轻量化要求越来越高，产品开发人员除了研究高强度材料的使用，更重要的是根据车轮设计边界条件，设计出一种最优的高强度结构，通过减少材料厚度来降低重量；
轮辐结构强度主要由9个关键参数共同影响，每个参数都有一个变化范围，这些参数有成百上千种组合，如何确定每一个参数值，找到同等料厚下最优结构组合方案，是车轮研发人员一直关注的技术难题，
传统设计方法主要采用经验设计，通过人工修改参数，反复建模，然后进行CAE对比分析，花大量的时间也很难找到最优结构，虽然后来采用优化软件对有限元模型进行优化，进一步降低了设计周期，但由于有限元模型还是比较复杂，动诸几天的计算周期对于本来就紧张的开发周期来说还是太长了，这样的设计方法会严重影响产品开发进度，甚至丢掉产品开发资格。
发明内容
本发明的目的是克服现有技术中存在的缺陷和不足，提供一种设计周期短、计算精度高的高强度车轮轮辐结构设计方法。
为实现以上目的，本发明的技术解决方案是：
一种高强度车轮轮辐结构设计方法，其特征在于，包括以下步骤：
1）建立轮辐参数化模型：依据开发技术要求及制动钳轮廓确定轮辐初始形状，利用CAD软件对车轮轮辐进行参数化建模，对轮辐9个关键尺寸R1、R2、R3、R4、R5、R6、D、H、A 进行参数化，导出关键参数表文件.txt；进行宏程序的录制以及脚本运行程序的编写，能通过脚本运行宏程序自动实现轮辐模型的更新；
2）建立轮辐有限元模型：将步骤1）得到的参数化模型导入CAE软件中，进行车轮弯曲加载有限元模型的搭建，具体操作如下：
①将轮辐3D参数化模型和制动钳3D模型导入CAE软件中进行接触应力分析；
②将轮辐3D参数化模型导入CAE软件中进行弯曲加载最大应力分析；
3）对步骤2）得到的有限元模型进行实验设计，具体操作如下：
①对步骤2）中的分析操作指令进行批处理操作，使得CAE建模与分析能够自动运行；
②在优化软件中对轮辐9个关键参数进行实验设计与计算，采用优化拉丁方对输入采样点的空间分布进行优化，保证设计点在空间尽可能的均匀分布，采样数目≥55个，提交计算得到实验设计结果，并对结果进行评价，导出结果文件；
4）建立轮辐近似模型：在优化软件中，将步骤3）生成的结果文件导入，采用创建响应面模型方法进行近似模型的建立，定义R1、R2、R3、R4、R5、R6、D、H、A 共9个参数为输入因子，接触应力和最大应力为响应因子，拟合成近似模型，对近似模型进行误差分析，控制误差＜5%；
5）对步骤4）得到的近似模型进行优化，选取轮辐9个关键参数R1、R2、R3、R4、R5、R6、D、H、A 为设计变量，以轮辐最大接触应力为0为约束条件，轮辐最大应力最小化为优化目标，直到得到最优结果。
所述步骤1）中①的具体操作如下：赋予材料属性，装配、创建分析步、提取轮辐内表面与卡钳外表面，创建接触关系、网格划分、提交计算，得到最大接触应力；输出轮辐最大接触应力到touch_stress.rpt文件，并将本次操作的指令记录导入到touch.py文件中。
所述步骤1）中②的具体操作如下：创建加载参考点、赋予材料属性，装配、创建分析步、参考点与安装面耦合约束、载荷加载和边界约束、网格划分、提交计算，得到最大弯曲应力；输出轮辐最大接触应力到load_stress.rpt文件，并将本次操作的指令记录导入到load.py文件中。
所述步骤3）中①的具体操作包括：编写运行touch.py的批处理运行程序run_touch.bat、编写运行load.py的批处理运行程序run_load.bat。
所述步骤3）、4）、5）中的优化软件为集成了CATIA和ABAQUS软件的Isight平台，所述步骤5）中的优化算法为多岛遗传算法。
与现有技术相比，本发明的有益效果为：
本发明通过对9个关键参数的建模，可以更好地贴近现实，能从全局范围内完成所有组合方案的对比，提高了产品设计可靠性；本设计方法在有限元模型的基础上增加了近似模型的构建与优化，大大缩短了产品的开发时间，也使得产品的设计精度有了很大的提升。
附图说明
图1是本发明的设计流程图。
图2是本发明实施提供的车轮轮辐参数化CAD模型示意图。
图3是本发明实施提供的制动钳CAD模型示意图。
图4是本发明实施提供的计算车轮弯曲加载轮辐与制动钳接触CAE模型示意图。
图5是本发明实施提供的车轮轮辐弯曲强度CAE模型图。
图6是本发明实施提供的是采用优化拉丁方得到的DOE样本点。
图7是本发明实施提供的是采用创建响应面模型，得到的近似模型误差分析示意图。
图8是本发明优化前后轮辐结构曲线对比。
图9是本发明优化前的轮辐最大应力云图。
图10是本发明优化后的轮辐最大应力云图。
具体实施方式
下面将结合附图和实施例对本发明实施例中的技术方案进行清楚、完整的描述。本发明的设计流程图参见图1。
一种高强度车轮轮辐结构设计方法，包括以下步骤：
1）建立轮辐参数化模型：依据开发技术要求及制动钳轮廓确定轮辐初始形状，利用CAD软件对车轮轮辐进行参数化建模，对轮辐9个关键尺寸R1、R2、R3、R4、R5、R6、D、H、A 进行参数化，导出关键参数表文件.txt；进行宏程序的录制以及脚本运行程序的编写，能通过脚本运行宏程序自动实现轮辐模型的更新，轮辐模型示意图和制动钳模型示意图参见图2-图3；
2）建立轮辐有限元模型：将步骤1）得到的参数化模型导入CAE软件中，进行车轮弯曲加载有限元模型的搭建，具体操作如下：
①将轮辐3D参数化模型和制动钳3D模型导入CAE软件中进行接触应力分析，赋予材料属性，装配、创建分析步、提取轮辐内表面与卡钳外表面，创建接触关系、网格划分、提交计算，得到最大接触应力；输出轮辐最大接触应力到touch_stress.rpt文件，并将本次操作的指令记录导入到touch.py文件中；
②将轮辐3D参数化模型导入CAE软件中进行弯曲加载最大应力分析，创建加载参考点、赋予材料属性，装配、创建分析步、参考点与安装面耦合约束、载荷加载和边界约束、网格划分、提交计算，得到最大弯曲应力；输出轮辐最大接触应力到load_stress.rpt文件，并将本次操作的指令记录导入到load.py文件中；车轮弯曲加载轮辐与制动钳接触CAE模型示意图参见图4，车轮轮辐弯曲强度CAE模型图参见图5，车轮弯曲加载应力云图参见图6；
3）对步骤2）得到的有限元模型进行实验设计，具体操作如下：
①编写运行touch.py的批处理运行程序run_touch.bat、编写运行load.py的批处理运行程序run_load.bat，通过运行批处理程序，使得CAE建模与分析能够自动运行；
②在Isight中对轮辐9个关键参数进行实验设计与计算，采用优化拉丁方对输入采样点的空间分布进行优化，保证设计点在空间尽可能的均匀分布，采样数目最小满足（9+1）*（9+2）/2=55个，提交计算得到DOE结果，并对结果进行评价，导出结果文件，DOE样本文件参见图7；
4）建立轮辐近似模型：在Isight中，将步骤3）生成的结果文件导入，采用创建响应面模型方法进行近似模型的建立，定义R1、R2、R3、R4、R5、R6、D、H、A 共9个参数为输入因子，接触应力和最大应力为响应因子，采用三阶多项式拟合成近似模型，对近似模型进行误差分析，控制误差＜5%，近似模型误差分析示意图参见图8；
5）对步骤4）得到的近似模型进行优化，选取轮辐9个关键参数R1、R2、R3、R4、R5、R6、D、H、A 为设计变量，以轮辐最大接触应力为0为约束条件，轮辐最大应力最小化为优化目标，选用多岛遗传算法进行优化，直到得到最优结果。
优化前的轮辐最大应力云图参见图9，轮辐最大应力364MPa；优化后的轮辐最大应力云图参见图10，轮辐应力为269.07MPa，应力降低26.1%，优化结果参见下表
Claims (5)
1.一种高强度车轮轮辐结构设计方法，其特征在于，包括以下步骤：
1）建立轮辐参数化模型：依据开发技术要求及制动钳轮廓确定轮辐初始形状，利用CAD软件对车轮轮辐进行参数化建模，对轮辐9个关键尺寸R1、R2、R3、R4、R5、R6、D、H、A 进行参数化，导出关键参数表文件.txt；进行宏程序的录制以及脚本运行程序的编写，能通过脚本运行宏程序自动实现轮辐模型的更新；
2）建立轮辐有限元模型：将步骤1）得到的参数化模型导入CAE软件中，进行车轮弯曲加载有限元模型的搭建，具体操作如下：
①将轮辐3D参数化模型和制动钳3D模型导入CAE软件中进行接触应力分析；
②将轮辐3D参数化模型导入CAE软件中进行弯曲加载最大应力分析；
3）对步骤2）得到的有限元模型进行实验设计，具体操作如下：
①对步骤2）中的分析操作指令进行批处理操作，使得CAE建模与分析能够自动运行；
②在优化软件中对轮辐9个关键参数进行实验设计与计算，采用优化拉丁方对输入采样点的空间分布进行优化，保证设计点在空间尽可能的均匀分布，采样数目≥55个，提交计算得到实验设计结果，并对结果进行评价，导出结果文件；
4）建立轮辐近似模型：在优化软件中，将步骤3）生成的结果文件导入，采用创建响应面模型方法进行近似模型的建立，定义R1、R2、R3、R4、R5、R6、D、H、A 共9个参数为输入因子，接触应力和最大应力为响应因子，拟合成近似模型，对近似模型进行误差分析，控制误差＜5%；
5）对步骤4）得到的近似模型进行优化，选取轮辐9个关键参数R1、R2、R3、R4、R5、R6、D、H、A 为设计变量，以轮辐最大接触应力为0为约束条件，轮辐最大应力最小化为优化目标，直到得到最优结果。
2.根据权利要求1所述的高强度车轮轮辐结构设计方法，其特征在于，所述步骤1）中①的具体操作如下：赋予材料属性，装配、创建分析步、提取轮辐内表面与卡钳外表面，创建接触关系、网格划分、提交计算，得到最大接触应力；输出轮辐最大接触应力到touch_stress.rpt文件，并将本次操作的指令记录导入到touch.py文件中。
3.根据权利要求1所述的高强度车轮轮辐结构设计方法，其特征在于，所述步骤1）中②的具体操作如下：创建加载参考点、赋予材料属性，装配、创建分析步、参考点与安装面耦合约束、载荷加载和边界约束、网格划分、提交计算，得到最大弯曲应力；输出轮辐最大接触应力到load_stress.rpt文件，并将本次操作的指令记录导入到load.py文件中。
4.根据权利要求1所述的高强度车轮轮辐结构设计方法，其特征在于，所述步骤3）中①的具体操作包括：编写运行touch.py的批处理运行程序run_touch.bat、 编写运行load.py的批处理运行程序run_load.bat。
5.根据权利要求1所述的高强度车轮轮辐结构设计方法，其特征在于：所述步骤3）、4）、5）中的优化软件为集成了CATIA和ABAQUS软件的Isight平台，所述步骤5）中的优化算法为多岛遗传算法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811298065.3A CN109598023A (zh) | 2018-11-02 | 2018-11-02 | 一种高强度车轮轮辐结构设计方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811298065.3A CN109598023A (zh) | 2018-11-02 | 2018-11-02 | 一种高强度车轮轮辐结构设计方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109598023A true CN109598023A (zh) | 2019-04-09 |
Family
ID=65957950
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811298065.3A Pending CN109598023A (zh) | 2018-11-02 | 2018-11-02 | 一种高强度车轮轮辐结构设计方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109598023A (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112084585A (zh) * | 2020-07-31 | 2020-12-15 | 东风汽车车轮随州有限公司 | 一种造型钢轮轻量化设计方法及装置 |
| CN114537035A (zh) * | 2022-03-18 | 2022-05-27 | 中新(重庆)超高强材料研究院有限公司 | 一种轻量化汽车车轮 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202080052U (zh) * | 2011-04-14 | 2011-12-21 | 东风汽车车轮有限公司随州车轮厂 | 一种高强度车轮 |
| CN102855359A (zh) * | 2012-09-11 | 2013-01-02 | 北京航空航天大学 | 一种汽车车轮变厚度轮辋的优化设计方法 |
| CN103942390A (zh) * | 2014-04-18 | 2014-07-23 | 吉林万丰奥威汽轮有限公司 | 一种汽车轮毂的优化和设计方法 |
| CN104077439A (zh) * | 2014-06-19 | 2014-10-01 | 西北工业大学 | 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法 |
| CN104239655A (zh) * | 2014-10-14 | 2014-12-24 | 大连大学 | 一种汽车轮毂的轻量化设计方法 |
-
2018
- 2018-11-02 CN CN201811298065.3A patent/CN109598023A/zh active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202080052U (zh) * | 2011-04-14 | 2011-12-21 | 东风汽车车轮有限公司随州车轮厂 | 一种高强度车轮 |
| CN102855359A (zh) * | 2012-09-11 | 2013-01-02 | 北京航空航天大学 | 一种汽车车轮变厚度轮辋的优化设计方法 |
| CN103942390A (zh) * | 2014-04-18 | 2014-07-23 | 吉林万丰奥威汽轮有限公司 | 一种汽车轮毂的优化和设计方法 |
| CN104077439A (zh) * | 2014-06-19 | 2014-10-01 | 西北工业大学 | 一种新型高强钢轮辐拉深冲孔复合工序的数值模拟方法 |
| CN104239655A (zh) * | 2014-10-14 | 2014-12-24 | 大连大学 | 一种汽车轮毂的轻量化设计方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112084585A (zh) * | 2020-07-31 | 2020-12-15 | 东风汽车车轮随州有限公司 | 一种造型钢轮轻量化设计方法及装置 |
| CN114537035A (zh) * | 2022-03-18 | 2022-05-27 | 中新(重庆)超高强材料研究院有限公司 | 一种轻量化汽车车轮 |
| CN114537035B (zh) * | 2022-03-18 | 2023-09-29 | 中新(重庆)超高强材料研究院有限公司 | 一种轻量化汽车车轮 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109977460B (zh) | 一种基于车身断面参数化的多目标优化设计方法 | |
| CN106934074B (zh) | 全局最优涡扇发动机进气道减噪设计方法 | |
| CN111597631B (zh) | 基于自适应代理模型的汽车风阻系数优化方法 | |
| CN104063550A (zh) | 一种基于云计算平台的多物理场cae系统 | |
| CN106951610B (zh) | 一种基于近似模型的插秧机秧箱结构优化方法 | |
| CN109598023A (zh) | 一种高强度车轮轮辐结构设计方法 | |
| CN107423512A (zh) | 一种大型压铸机尾板的轻量化设计方法 | |
| CN106919763A (zh) | 一种产品结构尺寸优化方法 | |
| CN104484511A (zh) | 一种基于仿真分析的机器人结构动态特性设计方法 | |
| CN108595800B (zh) | 一种网壳结构建模方法 | |
| CN103434585A (zh) | 客车车身结构减重方法 | |
| Sun et al. | Multi-objective optimization design of ladle refractory lining based on genetic algorithm | |
| CN106156534B (zh) | 一种特种集装箱基于Creo的钣金批量排样方法 | |
| CN115964882A (zh) | 工艺参数设计优化方法、处理器及制造设备 | |
| CN108959793B (zh) | 一种伞齿轮双频分段循环感应加热数值模拟方法 | |
| CN113094769B (zh) | 一种模块单元化叶片设计方法、装置、设备及存储介质 | |
| CN111444619B (zh) | 一种注塑模具冷却系统在线分析方法及设备 | |
| Zhao et al. | CNC thermal compensation based on mind evolutionary algorithm optimized BP neural network | |
| CN114925578B (zh) | 基于粒子群算法和abaqus的圆钢管混凝土短柱自动优化方法 | |
| CN110096742A (zh) | 一种预测模型及粒子群算法的预锻成形模具设计方法 | |
| CN113297758A (zh) | 一种大型复杂筋板件预锻成形初始毛坯优化设计方法 | |
| CN113761633A (zh) | 一种剪力墙结构优化设计方法 | |
| CN108256271B (zh) | 一种led灯具的轻量化设计方法 | |
| CN112464396B (zh) | 一种基于疏密度和局部复杂度的自适应代理模型在码垛机器人小臂驱动连杆优化中的应用 | |
| CN111967198B (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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190409 |
|
| RJ01 | Rejection of invention patent application after publication |