CN102855359A - 一种汽车车轮变厚度轮辋的优化设计方法 - Google Patents
一种汽车车轮变厚度轮辋的优化设计方法 Download PDFInfo
- Publication number
- CN102855359A CN102855359A CN2012103357098A CN201210335709A CN102855359A CN 102855359 A CN102855359 A CN 102855359A CN 2012103357098 A CN2012103357098 A CN 2012103357098A CN 201210335709 A CN201210335709 A CN 201210335709A CN 102855359 A CN102855359 A CN 102855359A
- Authority
- CN
- China
- Prior art keywords
- wheel rim
- software
- rim
- wheel
- stress
- 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
Images
Landscapes
- Tires In General (AREA)
Abstract
一种汽车车轮变厚度轮辋的优化设计方法,它有四大步骤:一:建立等截面轮辋有限元模型,计算其在径向载荷作用下的应力,确定强度约束条件;二:通过建立轮辋参数化CAD模型,设定轮辋典型尺寸,更新轮辋几何;三:将轮辐、轮辋CAD模型导入CAE软件中,模拟其在径向实验下的应力水平,设定轮辐、轮辋的连接关系、边界条件、载荷工况、划分网格,进行求解分析,得到轮辋最大应力;四:在优化平台软件中集成CAD软件和CAE软件,选择轮辋典型设计尺寸为设计变量,轮辋最大应力为约束条件,轮辋质量最小为优化目标,选择优化算法进行优化,直到得到最优结果。本发明提高了验证校核的准确性,“设计—计算—修改”优化过程自动运行,缩短了产品开发时间。
Description
技术领域
本发明涉及一种汽车车轮变截面轮辋的尺寸优化设计方法,尤其涉及一种汽车车轮变厚度轮辋的优化设计方法,属于汽车及机械工程技术领域。
背景技术
无内胎钢制车轮中的轮辐和轮辋两部分分别制造,然后焊接在一起。目前,无内胎钢制车轮的轮辐普遍采用的等强度结构,而无内胎钢制车轮的轮辋则还是按照等厚度结构设计,设计流程是通过反复的“设计—计算(或实验)—修改”来实现的。此种轮辋设计方法有以下不足之处:1、轮辋的截面厚度是按最危险载荷截面设计的,载荷小和载荷大的部位采取相同的厚度。这样,不仅浪费了钢材,也增加了车轮的质量;2、对轮辋尺寸设计的依据主要靠经验,CAE仅用来检验设计结果,而不是用来驱动产品设计,无法得到最优的产品结构,严重增加了设计时间和开发成本。
发明内容
本发明的目的,是为了解决上述问题,提供一种汽车车轮变厚度轮辋的优化设计方法。该方法不但提高了验证校核的准确性,而且设计优化过程“设计—计算—修改”可以自动运行,大大缩短产品开发时间。
本发明所提供的是一种变截面轮辋设计方法,它是集成优化方法,即将优化算法和CAD建模、CAE分析集成,首先建立CAD参数化模型、CAE分析模型,选择具有代表性的截面设计参数进行建模分析,优化器对这些模型进行优化,寻找最优的设计参数取值。
本发明通过以下技术方案予以实现:
本发明一种汽车车轮变厚度轮辋的优化设计方法,其特征在于,它包括以下步骤:
步骤一:建立等截面轮辋有限元模型,计算其在径向载荷作用下的应力;在径向载荷作用下轮辋与轮胎接触部位的应力分布近似服从余弦波状,波形中心夹角对称于压力方向,载荷作用的最大偏转角θ0的范围为30度至40度之间,θ0的含义见图2,它是径向载荷作用的最大偏转角。在进行仿真计算时,设应力分布为余弦波状且θ0为36度,将钢圈平均划分为10等份。作用力的间隔恰好为36度,依次施加该径向作用力,共计进行分析10次完成一完整作用力周期。
车轮径向分布力与最大径向分布力间的关系:
公式1中Wr为角度为θ时,对应的等效车轮径向分布力;
W0为等效的最大径向分布力;
对公式1进行积分得:
即: (公式2)
公式2中W:径向集中力
W0为等效最大径向分布力
b为胎圈座受力宽度
rb为胎圈座半径
θ0为径向分布载荷作用的最大偏转角
由于在车轮轮胎座上的分布力是作用在轮辋上的两侧,所以,公式2中所用的加载径向集中力(W)为试验加载力的一半。
从CAD软件导入轮辐、轮辋模型到CAE软件中,将径向载荷加在轮辋胎圈座上,设定轮辐、轮辋的连接关系为绑定关系、边界条件为约束轮辐内侧面的所有自由度、划分网格,保存有限元模型,进行求解分析,得到轮辋最大应力,作为步骤四优化过程的约束条件。
步骤二:通过CAD软件建立轮辋参数化模型,等截面轮辋截面由圆弧和直线组成,变截面轮辋的设计方法是采用样条曲线代替等截面轮辋中的圆弧线,添加尺寸约束实现参数化,这样可以通过控制参数的变化来变化轮辋截面形状。取截面控制参数D1、D2、D3等作为设计参数,导出轮辋几何STP格式文件供步骤三的CAE软件使用;
步骤三:将轮辋CAD模型STP格式文件导入CAE软件中,模拟其在径向实验下的应力水平,进行下列操作:
a.打开步骤一的有限元分析模型,删除原轮辋几何,导入更新后的轮辋几何;
b.设定轮辐、轮辋的连接关系为绑定关系、边界条件为约束轮辐内侧面的所有自由度、划分网格,进行求解分析,得到轮辋最大应力。
步骤四:在优化平台软件中集成CAD软件和CAE软件,选择轮辋典型设计尺寸为设计变量,等截面轮辋最大应力为约束条件,轮辋质量最小为优化目标,选择优化算法进行优化,不断重复步骤二、三,直到得到最优结果。
其中,步骤一中所述的CAD软件为SOLIDWORKS软件。
其中,步骤一中所述的CAE软件为ABAQUS软件。
其中,步骤四中所述的优化平台软件为ISIGHT软件,优化算法为多岛遗传算法。
本发明的优点在于:(1)轮辋为从中间向两边逐渐变薄的等强度变截面结构,按轮辋截面所受载荷的大小而变化的,载荷大的部位截面厚,载荷小的部位截面薄。这样既降低了质量,又保证了强度。(2)在轮辋设计初期采用CAE分析取代了传统的实验验证校核强度的方法,不但提高了验证校核的准确性,而且设计优化过程“设计—计算—修改”可以自动运行,大大缩短产品开发时间。
附图说明
图1是本发明实施例提供的等强度轮辋优化设计方法的流程图
图2是本发明实施例提供的车轮径向载荷分布示意图
图3是本发明实施例提供的轮辋参数化CAD模型的示意图
图4是本发明实施例提供的轮辋参数化CAD模型截面示意图
图5是本发明实施例提供的轮辐轮辋装配体的网格单元示意图
图6是本发明实施例提供的轮辐轮辋装配体的载荷和边界条件示意图
图7是本发明实施例提供的轮辐轮辋装配体在径向载荷工况下计算所得的应力云图
图8是本发明实施例提供的轮辋优化原理图
图9是本发明实施例提供的轮辋优化后的CAD模型截面示意图
图中符号说明如下:
图2中Bead Seat为载荷作用的位置,即胎圈座。Wr:车轮径向分布力;W0:最大径向分布力;b:为轮胎座受力宽度;θ0:径向载荷作用的最大偏转角。
图8中EXCEL表示EXCEL软件,SOLIDWORKS表示三维建模软件SOLIDWORKS软件,ABAQUS表示有限元分析软件ABAQUS软件,ISIGHT表示优化平台软件ISIGHT软件,abaqus.rpt表示存储ABAQUS计算结果最大应力值的文件。
具体实施方式
下面将结合附图和实施例对本发明实施例中的技术方案进行清楚、完整的描述。请参阅图1,图1是本发明实施例提供的等强度轮辋优化设计方法的流程图。
本发明一种汽车车轮变厚度轮辋的优化设计方法,它包括下列步骤:
步骤一:对于本发明实施例提供的轮辋,径向载荷分布如图2所示,具体加载参数如下:
W=88200/2=44100N,b=28mm,rb=280mm,θ0=36度,
则:
从CAD软件导入轮辐、轮辋模型到CAE软件中,将径向载荷加在轮辋胎圈座上,设定轮辐、轮辋的连接关系为绑定关系、边界条件为约束轮辐内侧面的所有自由度、划分网格,保存有限元模型,进行求解分析,得到等截面轮辋在径向载荷作用下的最大应力为349.49Mpa;
步骤二:如图3、图4所示,通过SOLIDWORKS软件建立轮辋参数化模型,添加轮辋截面特征尺寸D1,D7,D6,D2,D10,D3,D14,D15,D5,D17,D4,D19和D18,输出轮辋几何LW.step文件;
步骤三:将轮辋CAD模型导入ABAQUS软件中,模拟其在径向实验下的应力水平,进行下列操作:
a.打开步骤一建立的有限元分析模型,删除原轮辋几何,导入更新后的轮辋几何LW.step文件;
b.如图5、图6所示,设定轮辐、轮辋的连接关系为绑定关系、边界条件为约束轮辐内侧面的所有自由度、划分网格,进行求解分析,得到轮辋应力分布,应力云图如图7所示,输出轮辋最大应力到abaqus.rpt文件中。
步骤四:在ISIGHT中集成SOLIDWOKS软件和ABAQUS软件,以轮辋截面特征尺寸D1,D7,D6,D2,D10,D3,D14,D15,D5,D17,D4,D19和D18为设计变量,以等截面轮辋最大应力为约束条件,以轮辋质量最小为优化目标,选择多岛遗传算法进行优化,得到最优结果。图8为优化原理图,集成方法如下:
a.将设计变量写到EXCEL中,通过Visual Basic编程实现EXCEL对轮辋SOLIDWORKS模型特征尺寸的控制,输出轮辋质量到EXCEL中并输出轮辋几何LW.step文件;
b.通过ABAQUS的批处理文件自动进行原轮辋几何的删除,LW.step的导入,边界条件、载荷等的加载和网格划分,并计算输出轮辋最大应力到abaqus.rpt文件中;
优化结果如下表所示,优化后轮辋质量由31.22kg下降到28.97kg,减重7.2%,达到了节约成本的目的,优化后的变截面轮辋如图9所示。
Claims (4)
1.一种汽车车轮变厚度轮辋的优化设计方法,其特征在于:它包括以下步骤:
步骤一:建立等截面轮辋有限元模型,计算其在径向载荷作用下的应力;在径向载荷作用下轮辋与轮胎接触部位的应力分布近似服从余弦波状,波形中心夹角对称于压力方向,载荷作用的最大偏转角θ0的范围为30度至40度之间,θ0是径向载荷作用的最大偏转角,在进行仿真计算时,设应力分布为余弦波状且θ0为36度,将钢圈平均划分为10等份,作用力的间隔恰好为36度,依次施加该径向作用力,共计进行分析10次完成一完整作用力周期;
车轮径向分布力与最大径向分布力间的关系:
公式1中,Wr为角度为θ时,对应的等效车轮径向分布力;
W0为等效的最大径向分布力;
对公式1进行积分得:
即: (公式2)
公式2中,W:径向集中力
W0为等效最大径向分布力
b为胎圈座受力宽度
rb为胎圈座半径
θ0为径向分布载荷作用的最大偏转角
由于在车轮轮胎座上的分布力是作用在轮辋上的两侧,所以,公式2中所用的加载径向集中力W为试验加载力的一半;
从CAD软件导入轮辐、轮辋模型到CAE软件中,将径向载荷加在轮辋胎圈座上,设定轮辐、轮辋的连接关系为绑定关系,边界条件为约束轮辐内侧面的所有自由度、划分网格,保存有限元模型,进行求解分析,得到轮辋最大应力,作为步骤四优化过程的约束条件;
步骤二:通过CAD软件建立轮辋参数化模型,等截面轮辋截面由圆弧和直线组成,变截面轮辋的设计方法是采用样条曲线代替等截面轮辋中的圆弧线,添加尺寸约束实现参数化,这样可以通过控制参数的变化来变化轮辋截面形状;取截面控制参数D1、D2、D3作为设计参数,导出轮辋几何STP格式文件供步骤三的CAE软件使用;
步骤三:将轮辋CAD模型STP格式文件导入CAE软件中,模拟其在径向实验下的应力水平,进行下列操作:
a.打开步骤一的有限元分析模型,删除原轮辋几何,导入更新后的轮辋几何;
b.设定轮辐、轮辋的连接关系为绑定关系,边界条件为约束轮辐内侧面的所有自由度、划分网格,进行求解分析,得到轮辋最大应力;
步骤四:在优化平台软件中集成CAD软件和CAE软件,选择轮辋典型设计尺寸为设计变量,等截面轮辋最大应力为约束条件,轮辋质量最小为优化目标,选择优化算法进行优化,不断重复步骤二、三,直到得到最优结果。
2.根据权利要求1所述的一种汽车车轮变厚度轮辋的优化设计方法,其特征在于:步骤一所述的CAD软件为SOLIDWORKS软件。
3.根据权利要求1所述的一种汽车车轮变厚度轮辋的优化设计方法,其特征在于:步骤一所述的CAE软件为ABAQUS软件。
4.根据权利要求1所述的一种汽车车轮变厚度轮辋的优化设计方法,其特征在于:步骤四中的优化平台软件为ISIGHT软件,优化算法为多岛遗传算法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210335709.8A CN102855359B (zh) | 2012-09-11 | 2012-09-11 | 一种汽车车轮变厚度轮辋的优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210335709.8A CN102855359B (zh) | 2012-09-11 | 2012-09-11 | 一种汽车车轮变厚度轮辋的优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102855359A true CN102855359A (zh) | 2013-01-02 |
CN102855359B CN102855359B (zh) | 2015-06-03 |
Family
ID=47401945
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210335709.8A Expired - Fee Related CN102855359B (zh) | 2012-09-11 | 2012-09-11 | 一种汽车车轮变厚度轮辋的优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102855359B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103818184A (zh) * | 2014-03-18 | 2014-05-28 | 江铃汽车股份有限公司 | 一种提高汽车车轮强度的成型工艺 |
CN104281727A (zh) * | 2014-06-04 | 2015-01-14 | 河南金牌山模具有限公司 | 一种基于Excel及CAD平台的轮胎模具设计方法 |
CN104778326A (zh) * | 2015-04-23 | 2015-07-15 | 华东交通大学 | 一种基于优先数系的起重机车轮组适应性设计方法 |
CN105385839A (zh) * | 2014-09-09 | 2016-03-09 | 中国科学院沈阳自动化研究所 | 一种激光冲击强化自动化控制系统和方法 |
CN106326592A (zh) * | 2016-08-31 | 2017-01-11 | 浙江金固股份有限公司 | 一种轻量化钢制车轮的设计方法 |
CN107066715A (zh) * | 2017-03-31 | 2017-08-18 | 柳州顺联科技有限公司 | 工程机械中轮胎优化设计方法 |
CN107256304A (zh) * | 2017-06-02 | 2017-10-17 | 中北大学 | 一种枪弹弹头发射强度可靠性评估方法 |
CN108256183A (zh) * | 2018-01-02 | 2018-07-06 | 北京汽车股份有限公司 | 轮辋模型的确定方法和装置 |
CN109033635A (zh) * | 2018-07-26 | 2018-12-18 | 中车青岛四方车辆研究所有限公司 | 一种s型辐板车轮优化设计方法 |
CN109598023A (zh) * | 2018-11-02 | 2019-04-09 | 东风汽车车轮随州有限公司 | 一种高强度车轮轮辐结构设计方法 |
CN110489889A (zh) * | 2019-08-23 | 2019-11-22 | 中国电建集团成都勘测设计研究院有限公司 | 导入式结构有限元分析处理方法 |
CN112199880A (zh) * | 2020-10-27 | 2021-01-08 | 北京机电工程总体设计部 | 一种天线罩参数化及强度快速分析仿真系统 |
CN112199771A (zh) * | 2020-09-24 | 2021-01-08 | 燕山大学 | 车轮轮辋形状优化方法 |
CN112560185A (zh) * | 2020-12-22 | 2021-03-26 | 柳州市智甲金属科技有限公司 | 一种提高有限元分析中轮辋径向加载计算精度的方法 |
CN114896953A (zh) * | 2022-04-26 | 2022-08-12 | 东风汽车集团股份有限公司 | 车身侧倾刚度性能指标正向分解方法及分解装置 |
-
2012
- 2012-09-11 CN CN201210335709.8A patent/CN102855359B/zh not_active Expired - Fee Related
Non-Patent Citations (2)
Title |
---|
张红梅等: "基于有限元分析技术的等强度轮辋设计研究", 《现代制造技术与装备》 * |
李青海: "工程机械轮辋的计算机辅助设计与结构分析优化", 《中国优秀硕士学位论文全文数据库(电子期刊)工程科技II辑》 * |
Cited By (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103818184A (zh) * | 2014-03-18 | 2014-05-28 | 江铃汽车股份有限公司 | 一种提高汽车车轮强度的成型工艺 |
CN104281727A (zh) * | 2014-06-04 | 2015-01-14 | 河南金牌山模具有限公司 | 一种基于Excel及CAD平台的轮胎模具设计方法 |
CN105385839A (zh) * | 2014-09-09 | 2016-03-09 | 中国科学院沈阳自动化研究所 | 一种激光冲击强化自动化控制系统和方法 |
CN104778326A (zh) * | 2015-04-23 | 2015-07-15 | 华东交通大学 | 一种基于优先数系的起重机车轮组适应性设计方法 |
CN106326592B (zh) * | 2016-08-31 | 2019-11-05 | 浙江金固股份有限公司 | 一种轻量化钢制车轮的设计方法 |
CN106326592A (zh) * | 2016-08-31 | 2017-01-11 | 浙江金固股份有限公司 | 一种轻量化钢制车轮的设计方法 |
CN107066715A (zh) * | 2017-03-31 | 2017-08-18 | 柳州顺联科技有限公司 | 工程机械中轮胎优化设计方法 |
CN107256304A (zh) * | 2017-06-02 | 2017-10-17 | 中北大学 | 一种枪弹弹头发射强度可靠性评估方法 |
CN107256304B (zh) * | 2017-06-02 | 2020-07-07 | 中北大学 | 一种枪弹弹头发射强度可靠性评估方法 |
CN108256183B (zh) * | 2018-01-02 | 2021-09-17 | 北京汽车股份有限公司 | 轮辋模型的确定方法和装置 |
CN108256183A (zh) * | 2018-01-02 | 2018-07-06 | 北京汽车股份有限公司 | 轮辋模型的确定方法和装置 |
CN109033635A (zh) * | 2018-07-26 | 2018-12-18 | 中车青岛四方车辆研究所有限公司 | 一种s型辐板车轮优化设计方法 |
CN109033635B (zh) * | 2018-07-26 | 2022-10-21 | 中车青岛四方车辆研究所有限公司 | 一种s型辐板车轮优化设计方法 |
CN109598023A (zh) * | 2018-11-02 | 2019-04-09 | 东风汽车车轮随州有限公司 | 一种高强度车轮轮辐结构设计方法 |
CN110489889A (zh) * | 2019-08-23 | 2019-11-22 | 中国电建集团成都勘测设计研究院有限公司 | 导入式结构有限元分析处理方法 |
CN112199771A (zh) * | 2020-09-24 | 2021-01-08 | 燕山大学 | 车轮轮辋形状优化方法 |
CN112199771B (zh) * | 2020-09-24 | 2022-03-15 | 燕山大学 | 车轮轮辋形状优化方法 |
CN112199880A (zh) * | 2020-10-27 | 2021-01-08 | 北京机电工程总体设计部 | 一种天线罩参数化及强度快速分析仿真系统 |
CN112560185A (zh) * | 2020-12-22 | 2021-03-26 | 柳州市智甲金属科技有限公司 | 一种提高有限元分析中轮辋径向加载计算精度的方法 |
CN114896953A (zh) * | 2022-04-26 | 2022-08-12 | 东风汽车集团股份有限公司 | 车身侧倾刚度性能指标正向分解方法及分解装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102855359B (zh) | 2015-06-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102855359B (zh) | 一种汽车车轮变厚度轮辋的优化设计方法 | |
CN103942392B (zh) | 一种基于全生命周期的汽车底盘技术参数稳健设计方法 | |
CN102087670A (zh) | 基于多体动力学的汽车悬架与减振系统仿真试验及设计平台 | |
CN106202732B (zh) | 一种渐开线直齿圆柱齿轮副的齿向修形方法及与其配套的专用参数化cad系统 | |
CN105138796A (zh) | 一种用于整车振动噪声仿真的模态轮胎建模方法 | |
CN109977460A (zh) | 一种基于车身断面参数化的多目标优化设计方法 | |
CN111737816A (zh) | 一种非充气防爆车轮的轻量化设计方法 | |
CN105046012B (zh) | 一种考虑车轮侧倾的汽车车轮双轴疲劳试验仿真方法 | |
CN102567581A (zh) | 一种汽车控制臂的设计方法 | |
CN112084585A (zh) | 一种造型钢轮轻量化设计方法及装置 | |
CN109284582A (zh) | 一种车身沥青阻尼胶轻量化布置优化方法 | |
CN108491612A (zh) | 为复合管液压胀形工艺提供选材方案的有限元模拟方法 | |
CN1818914A (zh) | 一种计算机辅助的轿车仿真分析方法 | |
CN104915490A (zh) | 动车组头型气动反设计方法及装置 | |
CN116432426A (zh) | 车载动力电池评估方法、可读存储介质及计算机设备 | |
CN109657393A (zh) | 用于轮胎与底盘电控系统匹配的仿真平台及匹配仿真方法 | |
Messana et al. | Design, optimization and manufacturing of an aluminum wheel rim for the IDRAkronos vehicle prototype | |
CN106326592B (zh) | 一种轻量化钢制车轮的设计方法 | |
CN115017762A (zh) | 一种基于Bezier曲线的非充气轮胎的优化设计方法 | |
CN115455754A (zh) | 一种基于数字孪生的矿山液压支架设计方法 | |
CN115544746A (zh) | 一种多属性目标驱动的铝制副车架优化设计方法及系统 | |
CN115391916A (zh) | 一种车轮双轴疲劳仿真分析方法、装置、设备及介质 | |
CN112084572B (zh) | 一种车身造型阶段进行车身断面结构优化方法 | |
CN113408043B (zh) | 一种整车轮胎匹配方法及装置 | |
CN116522510B (zh) | 一种用于车轮性能仿真的轮胎材料参数的反演方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150603 Termination date: 20190911 |