CN113343351A - 一种汽车制动盘热-结构耦合分析方法 - Google Patents

一种汽车制动盘热-结构耦合分析方法 Download PDF

Info

Publication number
CN113343351A
CN113343351A CN202110539901.8A CN202110539901A CN113343351A CN 113343351 A CN113343351 A CN 113343351A CN 202110539901 A CN202110539901 A CN 202110539901A CN 113343351 A CN113343351 A CN 113343351A
Authority
CN
China
Prior art keywords
brake disc
analysis
thermal
heat transfer
friction
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
CN202110539901.8A
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.)
Chongqing Changan Automobile Co Ltd
Original Assignee
Chongqing Changan Automobile Co Ltd
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 Chongqing Changan Automobile Co Ltd filed Critical Chongqing Changan Automobile Co Ltd
Priority to CN202110539901.8A priority Critical patent/CN113343351A/zh
Publication of CN113343351A publication Critical patent/CN113343351A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • 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]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/10Numerical modelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Algebra (AREA)
  • Braking Arrangements (AREA)

Abstract

本发明公开了一种汽车制动盘热‑结构耦合分析方法,其包括如下步骤:S1,将整车三维模型导入处理软件中,建立整车外部流场域;S2,对不同车速条件下的整车进行稳态流场分析,得到制动过程中制动盘的特征面上随时间变化的对流换热系数;S3,将制动盘和摩擦片的三维模型导入有限元分析软件中,根据S2得到的对流换热系数设定制动盘与外界的对流换热边界条件;S4,建立多分析步,模拟制动盘的旋转运动、摩擦片在活塞缸推动下贴紧制动盘的瞬时动作以及制动减速过程,进行热‑结构耦合求解及结果后处理。其能够在保证计算效率的同时提高计算精度,使得计算得到的温度场和应力场特征与实际情况相符,利于预测制动盘的热性能。

Description

一种汽车制动盘热-结构耦合分析方法
技术领域
本发明涉及汽车制动盘性能评估,具体涉及汽车制动盘热-结构耦合分析方法。
背景技术
对盘式制动器而言,在汽车的制动过程中,随着动能不断转换为热量,制动盘的温度会不断升高,由此会带来一些危害。一方面,高温会降低制动盘与摩擦片之间的摩擦系数,产生制动热衰退现象。另一方面,盘面摩擦区域与其他区域的温差过大会使制动盘产生热变形,引发制动抖动、制动热疲劳等现象,影响用户体验感及行驶安全性。因此,在整车开发过程中,有必要采取措施来获取制动盘在制动过程中的温度、热应力及热变形等热性能,从而保障车辆的行驶安全及舒适。
通常,整车厂都会采用台架试验或整车试验的方式来测试制动盘在制动过程中的热性能。由于台架试验的工况难以和真实驾驶工况保持一致(尤其在制动盘外部流场方面),因此往往出现即使台架试验通过,但实车路试仍然出问题的情况。因此,在整车上对制动盘开展热性能试验仍然是必不可少的。整车试验虽然结果可靠,但成本较高,而且只能在整车开发后期进行,只有验证性,不具有预测性,是一种不经济、不高效的手段。相较这两种试验手段而言,采用CAE的方式在整车开发前期对制动盘进行热性能仿真分析,在降低成本、提升效率、规避质量问题等方面具有明显的优势。
当前,采用CAE仿真手段来获取制动盘热性能的主要途径为基于有限元法的热-结构耦合分析。这种方法虽然计算效率高,但往往在计算精度方面难以保障,主要体现在两个方面:①
不能准确考虑制动盘与外部环境的外部对流换热边界条件,致使计算得到的温度场有较大误差;②采用将车辆的动能转换为热量的思维来处理制动热源,与真实的摩擦现象有较大区别,使得计算得到的温度场和应力场特征与实际情况不符。
发明内容
本发明的目的是提供一种汽车制动盘热-结构耦合分析方法,其能够在保证计算效率的同时提高计算精度,使得计算得到的温度场和应力场特征与实际情况相符,利于预测制动盘的热性能。
本发明所述的汽车制动盘热-结构耦合分析方法,其包括如下步骤:
S1,将整车三维模型导入处理软件中,建立整车外部流场域;
S2,对不同车速条件下的整车进行稳态流场分析,得到制动过程中制动盘的特征面上随时间变化的对流换热系数;
S3,将制动盘和摩擦片的三维模型导入有限元分析软件中,根据S2得到的对流换热系数设定制动盘与外界的对流换热边界条件;
S4,建立多分析步,模拟制动盘的旋转运动、摩擦片在活塞缸推动下贴紧制动盘的瞬时动作以及制动减速过程,进行热-结构耦合求解及结果后处理。
进一步,所述S3中将制动盘和摩擦片的三维模型导入有限元分析软件后,设置制动盘和摩擦片的材料属性,划分网格,设定制动盘与摩擦片之间的摩擦系数,然后设定制动盘与外界的对流换热边界条件。
进一步,按照GB5763-2018获取不同温度条件下制动盘与摩擦片之间的摩擦系数。
进一步,所述S3中根据S2得到的制动盘的特征面上的对流换热系数,以传热膜系数的形式设定制动盘与外界的对流换热边界条件,设定制动盘与摩擦片的初始温度与环境温度相同。
进一步,所述S4中根据实际制动工况建立多分析步,具体为:
a.建立制动盘与摩擦片的接触关系,制动盘上的摩擦面为主面,摩擦片上的摩擦面为从面,接触类型为硬接触;
b.对制动盘施加初始角速度,模拟制动盘的旋转运动;
c.加载制动压力载荷,模拟摩擦片在活塞缸推动下贴紧制动盘的瞬时动作;
d.模拟车辆的制动减速过程,时间周期为制动时长;
e.建立分析任务,进行热-结构耦合求解。
进一步,所述S4中的后处理包括对热-结构耦合求解结果进行分析,得到制动盘温度场特征和制动盘应力场特征。
本发明通过联合整车稳态流场分析来获得精确的制动盘与外界的对流换热边界条件,并且对摩擦生热过程进行直接仿真,在此基础上进行热-结构耦合仿真分析和结果后处理,进而在保证计算效率的同时得到了高精度的制动盘温度场特征和制动盘应力场特征,为预测汽车制动盘的热性能提供了一种精准、可靠的新方法。
附图说明
图1是本发明的流程示意图;
图2是本发明制动盘表面划分及编号示意图;
图3是本发明整车流场计算域示意图;
图4是本发明制动盘特征面上的对流换热系数随时间变化曲线图;
图5是本发明制动盘的有限元分析模型示意图;
图6是本发明制动盘与摩擦片之间的摩擦系数随温度变化情况示意图;
图7是本发明制动盘表面温度分布云图;
图8是本发明制动盘表面节点位置示意图;
图9是本发明制动盘表面节点的温升曲线示意图;
图10是本发明制动盘表面应力分布云图。
具体实施方式
下面结合附图对本发明作详细说明。
参见图1,所示的汽车制动盘热-结构耦合分析方法,其包括如下步骤:
S1,将整车三维CAD模型导入处理软件中,该处理软件为STAR-CCM+,在完成几何清理后,根据几何特征对制动盘表面进行划分并编号,参见图2,FI为内侧端面,F2为外侧端面,F3为通风孔内侧端面,F4为通风孔外侧端面,F5为通风孔柱面,F6为外缘面,F7为中心柱面,F8为鹅颈内表面,F9为鹅颈外表面。
建立整车外部流场域,参见图3,L为车长,W为车宽,H为车高,X、Y、Z分别为车长、车宽、车高方向,由于在实际情况下,有气流经发动机舱流向制动盘,会影响制动盘的冷却效果。因此为了提高仿真精度,应考虑发动机舱内部的流场,体现在以多孔介质的方法来处理散热器、冷凝器及中冷器,以MRF的方法来处理风扇。进行网格划分,需要说明的是,对边界层和制动盘周边流场的网格进行加密处理。对整车外部流场赋予空气属性,设置流场入口为速度边界条件,出口为压力边界条件,采用MRF即多重参考系的方法来处理制动盘及车轮的旋转。
S2,设置湍流参数和求解参数,分别在0km/h、10km/h、20km/h、40km/h、60km/h、80km/h、100km/h七个的车速条件下来进行稳态流场求解。对不同车速条件下的整车稳态流场计算结果进行处理,提取制动盘各特征面的对流换热系数。插值处理后,参见图4,得到整个制动过程中,制动盘各特征面上随时间变化的对流换热系数。
S3,将制动盘和摩擦片的三维CAD模型导入ABAQUS有限元分析软件中,设置制动盘和摩擦片的材料属性。参见图5,划分网格,网格类型为C3D8RT。设定制动盘与摩擦片之间的摩擦系数,为了提高仿真精度,参见图6,按照GB5763-2018获取不同温度条件下制动盘与摩擦片之间的摩擦系数。然后根据S2得到的制动盘的特征面上的对流换热系数,以传热膜系数的形式设定制动盘与外界的对流换热边界条件,设定制动盘与摩擦片的初始温度与环境温度相同。
S4,根据实际制动工况建立多分析步,具体为:
a.建立初始分析步,即建立制动盘与摩擦片的接触关系,制动盘上的摩擦面为主面,摩擦片上的摩擦面为从面,接触类型为硬接触。
b.建立分析步1,类型为Static General,对制动盘施加初始角速度,模拟制动盘的旋转运动。
c.建立分析步2,类型为Static General,加载制动压力载荷,模拟摩擦片在活塞缸推动下贴紧制动盘的瞬时动作。
制动压力由公式
Figure BDA0003071256140000041
计算得到,其中,n为夹钳活塞缸数量,Pf为活塞缸内部制动液的压力,Dp为活塞缸内径,A为摩擦片与制动盘的摩擦面积。
d.建立分析步3,模拟车辆的制动减速过程,分析步类型为Coupled temp-displacement,时间周期为制动时长;
e.建立分析任务,进行热-结构耦合求解。
后处理具体是对热-结构耦合求解结果进行分析,得到整个制动过程中制动盘温度场特征和制动盘应力场特征。参见图7,得到整个制动过程中制动盘表面的温度分布特征,参见图8和图9,在划分网格时得到节点1至节点8,对热-结构耦合求解结果进行分析得到制动盘上各个节点的温升情况,参见图10得到制动盘在制动过程中的应力场分布特征以及制动盘变形情况。进而能够根据制动盘温度场特征和制动盘应力场特征对汽车制动盘的热性能进行准确预测,相较于传统的热-结构耦合分析方法,本发明通过联合整车稳态流场分析来获得精确的制动盘与外界的对流换热边界条件,并且对摩擦生热过程进行直接仿真,在此基础上进行热-结构耦合仿真分析和结果后处理,进而在保证计算效率的同时得到了高精度的制动盘温度场特征和制动盘应力场特征,为预测汽车制动盘的热性能提供了一种精准、可靠的新方法。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (6)

1.一种汽车制动盘热-结构耦合分析方法,其特征在于,包括如下步骤:
S1,将整车三维模型导入处理软件中,建立整车外部流场域;
S2,对不同车速条件下的整车进行稳态流场分析,得到制动过程中制动盘的特征面上随时间变化的对流换热系数;
S3,将制动盘和摩擦片的三维模型导入有限元分析软件中,根据S2得到的对流换热系数设定制动盘与外界的对流换热边界条件;
S4,建立多分析步,模拟制动盘的旋转运动、摩擦片在活塞缸推动下贴紧制动盘的瞬时动作以及制动减速过程,进行热-结构耦合求解及结果后处理。
2.根据权利要求1所述的汽车制动盘热-结构耦合分析方法,其特征在于:所述S3中将制动盘和摩擦片的三维模型导入有限元分析软件后,设置制动盘和摩擦片的材料属性,划分网格,设定制动盘与摩擦片之间的摩擦系数,然后设定制动盘与外界的对流换热边界条件。
3.根据权利要求2所述的汽车制动盘热-结构耦合分析方法,其特征在于:按照GB5763-2018获取不同温度条件下制动盘与摩擦片之间的摩擦系数。
4.根据权利要求1或2所述的汽车制动盘热-结构耦合分析方法,其特征在于:所述S3中根据S2得到的制动盘的特征面上的对流换热系数,以传热膜系数的形式设定制动盘与外界的对流换热边界条件,设定制动盘与摩擦片的初始温度与环境温度相同。
5.根据权利要求1或2所述的汽车制动盘热-结构耦合分析方法,其特征在于:所述S4中根据实际制动工况建立多分析步,具体为:
a.建立制动盘与摩擦片的接触关系,制动盘上的摩擦面为主面,摩擦片上的摩擦面为从面,接触类型为硬接触;
b.对制动盘施加初始角速度,模拟制动盘的旋转运动;
c.加载制动压力载荷,模拟摩擦片在活塞缸推动下贴紧制动盘的瞬时动作;
d.模拟车辆的制动减速过程,时间周期为制动时长;
e.建立分析任务,进行热-结构耦合求解。
6.根据权利要求1或2所述的汽车制动盘热-结构耦合分析方法,其特征在于:所述S4中的后处理包括对热-结构耦合求解结果进行分析,得到制动盘温度场特征和制动盘应力场特征。
CN202110539901.8A 2021-05-18 2021-05-18 一种汽车制动盘热-结构耦合分析方法 Pending CN113343351A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110539901.8A CN113343351A (zh) 2021-05-18 2021-05-18 一种汽车制动盘热-结构耦合分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110539901.8A CN113343351A (zh) 2021-05-18 2021-05-18 一种汽车制动盘热-结构耦合分析方法

Publications (1)

Publication Number Publication Date
CN113343351A true CN113343351A (zh) 2021-09-03

Family

ID=77469105

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110539901.8A Pending CN113343351A (zh) 2021-05-18 2021-05-18 一种汽车制动盘热-结构耦合分析方法

Country Status (1)

Country Link
CN (1) CN113343351A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114239132A (zh) * 2021-11-25 2022-03-25 合众新能源汽车有限公司 一种复合材料汽车制动盘的温度场仿真方法及装置
CN117648772A (zh) * 2023-12-07 2024-03-05 青岛青特众力车桥有限公司 一种鼓式制动器热结构直接耦合分析方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663186A (zh) * 2012-04-06 2012-09-12 重庆长安汽车股份有限公司 一种制动器热衰退参数化分析方法
WO2016169249A1 (zh) * 2015-04-21 2016-10-27 中车青岛四方机车车辆股份有限公司 一种高速行驶交通工具仿真模型前处理模板及处理方法
CN107704644A (zh) * 2017-07-20 2018-02-16 同济大学 一种考虑热机耦合效应的制动尖叫预测方法
CN108279257A (zh) * 2018-02-28 2018-07-13 苏州科技大学 测量强迫对流加热棒簇努塞尔数的装置及方法
CN111125858A (zh) * 2018-11-01 2020-05-08 广州汽车集团股份有限公司 一种整车冷却系统的仿真计算方法及系统
CN112464373A (zh) * 2020-11-25 2021-03-09 中汽研(天津)汽车工程研究院有限公司 多循环制动工况车辆盘式制动器温度变化高精度仿真方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102663186A (zh) * 2012-04-06 2012-09-12 重庆长安汽车股份有限公司 一种制动器热衰退参数化分析方法
WO2016169249A1 (zh) * 2015-04-21 2016-10-27 中车青岛四方机车车辆股份有限公司 一种高速行驶交通工具仿真模型前处理模板及处理方法
CN107704644A (zh) * 2017-07-20 2018-02-16 同济大学 一种考虑热机耦合效应的制动尖叫预测方法
CN108279257A (zh) * 2018-02-28 2018-07-13 苏州科技大学 测量强迫对流加热棒簇努塞尔数的装置及方法
CN111125858A (zh) * 2018-11-01 2020-05-08 广州汽车集团股份有限公司 一种整车冷却系统的仿真计算方法及系统
CN112464373A (zh) * 2020-11-25 2021-03-09 中汽研(天津)汽车工程研究院有限公司 多循环制动工况车辆盘式制动器温度变化高精度仿真方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
夏毅敏等: "准高速机车制动盘三维温度场及热应力场分析", 《现代制造工程》 *
张森等: "汽车通风盘式制动器的流固热多物理场耦合建模与分析", 《机械工程学报》 *
朱晴等: "汽车制动盘温度场瞬态分析方法的研究", 《汽车技术》 *
杨莺等: "矿井提升机制动盘温度场与应力场的数值模拟", 《机械强度》 *
葛振亮等: "盘式制动器热弹性耦合分析", 《烟台大学学报(自然科学与工程版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114239132A (zh) * 2021-11-25 2022-03-25 合众新能源汽车有限公司 一种复合材料汽车制动盘的温度场仿真方法及装置
CN117648772A (zh) * 2023-12-07 2024-03-05 青岛青特众力车桥有限公司 一种鼓式制动器热结构直接耦合分析方法

Similar Documents

Publication Publication Date Title
CN107145663B (zh) 车轮多目标优化设计方法
CN113343351A (zh) 一种汽车制动盘热-结构耦合分析方法
CN107704644A (zh) 一种考虑热机耦合效应的制动尖叫预测方法
CN102854010B (zh) 一种基于道路循环工况的发动机零件疲劳寿命计算方法
CN111814261A (zh) 整车冷却特性数据确定方法、装置、设备及可存储介质
Raj et al. Numerical investigation of fluid flow and heat transfer characteristics on the aerodynamics of ventilated disc brake rotor using CFD
Simmonds et al. Complete body aerodynamic study of three vehicles
Fukano et al. Development of disc-brake design method using computer simulation of heat phenomena
Thuresson CFD and Design Analysis of Brake Disc
CN104615835B (zh) 一种发动机中冷器分析方法
Belhocine et al. A predictive tool to evaluate braking system performance using a fully coupled thermo-mechanical finite element model
CN102663186A (zh) 一种制动器热衰退参数化分析方法
Bhambare et al. Brake duty cycle simulation for thermal design of vehicle braking system
Cakmak et al. Estimation of friction surface temperature of a dry clutch
CN115964891A (zh) 一种汽车风扇的仿真分析方法、装置及存储介质
Hunt et al. A coupled simulation approach to race track brake cooling for a GT3 race car
CN110826278B (zh) 基于有限元的硅油风扇离合器散热性能分析方法
KR102070335B1 (ko) 타이어 회전저항 예측방법
Patel et al. Optimization and performance analysis of an automobile radiator using CFD—a review
Thorat et al. CFD investigation of airflow through fan and cooler system of an off road vehicle
Hobeika Wheel modelling and cooling flow effects on car aerodynamics
More et al. Cfd analysis of automotive ventilated disc brake rotor
D'Hondt et al. Modeling and Identification of an Electric Vehicle Braking System: Thermal and Tribology Phenomena Assessment
Zhang et al. On an efficient simulation approach for estimating the cooling performance of automotive vented disc brakes under the scenario of emergency braking
Xiao Transient simulation of heat transfers for vehicle exhaust system

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210903