基于ANSYS的制动盘热固耦合分析

第5章基于ANSYS的制动盘热固耦合分析ANSYS软件简介ANSYS是美国ANSYS公司设计开发的大型通用有限元计算软件，是一个融结构、热、流体、电、磁。声学于一体的大型通用有限元软件。作为目前最流行的有限元软件之一，它具备功能强大、兼容性好、使用方便、计算速度快等优点，成为了工程师们开发设计的首选，并广泛应用于核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、生物医学、水利、日用家电等一般工业及科学研究领域。其特点如下：1）能实现多场及多场耦合分析2）能实现统一前处理、求解、后处理及多场分析数据库的大型FEA软件3）具有多物理场优化功能的FEA软件4）较为强大的非线性分析功能5）多种求解器分别适用于不同的问题及不同的硬件配置6）支持异种、异构平台的网络浮动，在异种、异构平台上用户界面统一、数据文件全部兼容7）强大的并行计算功能支持分布式并行及共享内存式并行8）多种自动网格划分技术9）良好的用户开发环境制动盘有限元模型的建立本课题设计的制动盘为双摩擦面带中间散热筋板的自通式制动盘。在有限元计算分析时为提高计算效率，在计算过程中对模型做了必要的简化处理，最后选为单摩擦面施加对称约束并进行有限元计算。制动盘和制动片单元类型选用实体单元S0LID5。该问题属于摩擦生热为题，选择C0NTA173接触单元和TARGE170目标单元生成接触对。模型材料特性参数见表5-1。制动盘的几何三维模型见图5-1。对其进行网格划分，划分网格后的有限元模型见图5-2所示。表5-1制动器模型材料的特性参数特性参数制动盘制动片热导率九/(W-m-1-K-1)48.461.212密度p/(kg-m3)72282595比热容c/(J-kg-1-K-1)4191465热膨胀系数a/K-111x10-630x10-4

弹性模量E/GPa1751.5泊松比0.30.25图5-1制动盘的几何三维模型图5-2图5-1制动盘的几何三维模型图5-2划分网格后的有限元模型图5.3施加边界约束及制动盘温度场计算分析热分析模型边界条件的确定是温度场分析的关键所在，本文采用理论推导与实际分析相结合的方法，制动盘施加轴向转动，制动片施加轴向的移动。制动前，制动盘温度等于环境温度，计算制动盘所在的工作环境温度确定，这里取工作环境平均温度为25°C。制动初速度为100Km/h,施加到制动片上的力为89ON,依次进行温度场-应力场计算分析。图5-3分别给出了制动过程中10s、16.4s（最高温度时刻）、30s（制动结束时刻）及60s（停车后30s）时制动盘温度场的分布云图,图5-4为制动盘表面最高温度随时间变化曲线。a）升温过程中，t=10sNODALSOLUTIONSTBP»1SUB=4TIHE=.100E-03TBHPPawenGrachicsErACET=l3MN=-.0634883V=.6826*DI3T».154418-.06348820.32240.70Sa）升温过程中，t=10sNODALSOLUTIONSTBP»1SUB=4TIHE=.100E-03TBHPPawenGrachicsErACET=l3MN=-.0634883V=.6826*DI3T».154418-.06348820.32240.70S81.479101.8€S122.251163.022.100E-03.006123-.022116-.005661b）最高温度时刻，t=16.4sc）制动结束时刻，t=30sd）制动结束后30s,t=60sc）制动结束时刻，t=30sd）制动结束后30s,t=60sNODALSOLUTION8TEF»1SUB=1TIME=.500E-04TEMP(AVG»PoireEGzaphxcaEFACET=1AVRE3"Ma?DMX=.105B-043MN=-.01S883KX■91.704XV—.133535YV--.71B488ZV=.6826*DI3T=.154418*X?s.006123・YF=-.022116亠ZF—.005661A-28-8.€540158810.17520.36630.55740.74950.94€1.13171.32281.51391.704ZEJO图5-3制动开始后几个典型时刻制动盘温度场分布云图图5-4摩擦面最高温度随时间变化曲线分析可知，制动开始后，制动盘摩擦面温度急剧上升，在制动开始10s时，盘面温度已达到128°C，且温度集中分布在盘体摩擦面一侧，散热筋板根部的最高温度约为43C。这是因为，在这一时间段内，由于热量尚未传到散热筋上，制动盘的热量输入远大于制动盘的散热能力，热量在盘内积聚，使盘体温度急剧增加。同时，从温度变化曲线可以看出，由于热流密度随时间的增加而降低，温度变化梯度呈下降趋势。随着制动过程的继续，温度变化梯度进一步降低，热能逐渐向盘体内、外径向方向及散热筋上传递，制动盘的散热能力增强，到制动开始16.4s时，制动盘的热量输入与热量散失达到平衡，盘面温度达到最高为183C,此时散热筋板根部温度约为61C。此后，随着制动时间的延续，热能继续向盘面内、外径向方向和散热筋板传递，制动盘散热能力增强，盘面温度开始下降，温度变化梯度也随热流密度的降低进一步减小，到制动结束时，制动盘摩擦面中部和散热筋板根部的温度相差不到70C。制动结束后，辐射成为热量散失的主要形式。随着制动时间的进一步推移，

盘体与散热筋板滞见得温度差逐渐减小，到制动结束后30s时，散热筋板的温度已接近61°C，与盘体的最高温度106°C仅差30°C。5.3制动盘热应力场计算及分析由制动盘的温度场瞬态分布计算结果，再利用ANSYS软件可以求得制动盘的应力场。运用中的制动盘除受热膨胀引起的应力外，还存在闸片的压力和摩擦力，及由于制动盘旋转旋转引起的离心力、振动载荷以及压装载荷等作用产生的应力，但这些因素对制动盘总应力的影响一般不高于20MPa（参考文献［19］）,远不及热应力的影响，因此，在结构分析中忽略这些因素，仅考虑热应力的影响。经ANSYS求解，同样可以得到任意时刻制动盘任意位置的应力值和制动盘的变形情况。图5-5为制动过程中第10s（升温过程10s），11.5s（整个制动盘的最高应力时刻），30s（制动结束时刻），60s（制动结束30s）时制动盘上应力分布云图，图5-6为制动盘上的合位移分布云图。NODALSOLUTIONSTEP-1SUB-4TIME-.05NODALSOLUTIONSTEP-1SUB-4TIME-.05SEQV(AVG)PowerGraphlcsEFACET=1AVRES>MatDMX147E-04SMN>253384SMX=.228E+09XV■.212852YV—.482523ZV849627*DIST=.1527C2•XF=.008765*YF--.01336S■ZF-.011846A-ZS-1.294Z・BUFFER253384.256E+08IJ.509E+081J.762E+08.101E+09■■.127E409.152E+09.177E+091.203E+09.228E+09SUB-5TIME*.06SEQV(AVG)PowerGraphicsEFACET-1AVRES-MatDMX«.184E-04SMN=305964SMX».271E+09XV212852YV»-.4B2523ZV=.849627*DIST».152762•XF-.000765•YF—.013365♦ZF=.011846A-ZS=1.294Z-BUFFER■吕口30S964.304E+08.6O5E+OB.9O7E+O8.121E+09.1S1E+09.181E+09.211E+09.241E*09•271E+09a）升温过程中，t=10sb）最大应力时刻，t=11.5sNODALSOLUTIONSTEP=1SUBTIME-.04SEQV(AVG)PowerGraphicsEFACET=1AVRES-MatDMX113E-04SMN»200047SMX=.184E+09XV>.21285230S964.304E+08.6O5E+OB.9O7E+O8.121E+09.1S1E+09.181E+09.211E+09.241E*09•271E+09a）升温过程中，t=10sb）最大应力时刻，t=11.5sNODALSOLUTIONSTEP=1SUBTIME-.04SEQV(AVG)PowerGraphicsEFACET=1AVRES-MatDMX113E-04SMN»200047SMX=.184E+09XV>.212852YV--.482523ZV-.849627*DIST=.152762♦XF=.008765•YF--.013365•ZF>.011846A-ZS-1.294Z-BUFFER■L_J吕200047.206E+08.411E+O8.615E+O0.019E+O8.102E+09.123E+09•143E+09.164E+09.184E+09NODALSOLUTIONSTEP=1SUBTIME-.03SEQV(AVG)PowerGraphicsEFACET=lAVRES»d4atDMX-.80€E-05SMN»150146SMX=.139E+09XV-.212852YV—.482523ZV-.849627♦DI3T=.152762•XF=.008765•YF—.013365・ZF>.01184(A-ZS-L.294Z-BUFFER■nL_Jg150146.1S6E*O8.3L1E+O8.465E+08.620E+08.775E+08.929E+08.108E+09.124E+09.139E+09c）制动结束时刻，t=30sd）制动结束后30s，t=60s图5-5制动开始后几个典型时刻制动盘应力分布云图NODALSOLUTIONSTEP=1SUB»3TIB4E-.04USUM(AVG)RSYS=0PowerGraphicsEFACETAVRE3-MatDMX=.113B-04SMX=.113E・04XV063364YV--.478584ZV».B75751*DIST=.152762*XF=.008765•YF>-.013365•ZF•.011846A-ZS»6.042Z-BUFFER05■0505050505OS0404NODALSOLUTIONSTEP=1SUB»3TIB4E-.04USUM(AVG)RSYS=0PowerGraphicsEFACETAVRE3-MatDMX=.113B-04SMX=.113E・04XV063364YV--.478584ZV».B75751*DIST=.152762*XF=.008765•YF>-.013365•ZF•.011846A-ZS»6.042Z-BUFFER05■0505050505OS0404125E251E376E501E627E7S2E077E100E113ENODALSOLUTIONSTE?=1SUB-4TIME-.05USUM(AVO)RSYS=OPowerGraphicsEFACET-1AVRES-MatDMX=.147E-04SMX=.147E-04XV-.063384YV—.478584ZV-.875751•DIST=.152762•XF-.008765•YF—.013365•2F>.011846A-ZS-6.042.491E-05.13IE-04.147E-04a）升温过程中,t=10sb）最大合位移时刻,t=11.5sZF>.011846A-ZS=6.042NODALSOLUTIONSTEP=1SUB-2TIME-.03USUM(AVG)RSYS=0PowerGraphicsEFACET-1AVRES-MatDMXN.806E-05SMX=.806E-05XV».0633S4YV—.478584ZV-.875751♦DIST=.152762♦XF=.008765•YF—.0133€SBUFFER448E537E621ENODALSOLUTIONSTEP-1SUB«1TIME-.02USUM(AVO)RSYS=0PowerGraphicsEFACET-1AVRES=MatDMX=.510E-05SMX=.510E-05XV-.063384YV—.478584ZV=.875751•DIST».152762•XF-.00876S*YF—・013365*ZF=.011846A-ZS=6.042Z-BUFFERs口昌吕■.567E-06.113E-05,1706-05.227E-05.283E-05.340E-05.397E-05.454E-OS.510E-05c）制动结束时刻,t=30sd）制动结束后30s,t=60s图5-6制动开始后几个典型时刻制动盘合位移分布云图从上面分析图可以看出,不同时刻制动盘的最高应力值及出现的位置均有所不同。在制动开始初期,由于制动盘摩擦面的温度始终是最高区域,制动盘摩擦面的热膨胀也最大，此处的应力也就最大，且在制动开始约10s时，制动盘面的应力达到228MPa，合位移为0.11mm。而整个制动过程中制动盘的最大应力值为271MPa，最大位移值为0.15mm，发生于制动开始后11.5s，出现在制动盘长散热筋板。这是因为，热量传导至散热筋板需要一定的时间，制动开始时散热筋板的温度还比低，此时散热筋板的形状基本保持不变，而随着盘体的热膨胀增加，使得盘体热变形呈鼓形，且由于长散热筋板的径向跨度较大，及此处存在有应力集中现象，所以此处应力水平也最高。随着制动过程进一步延续，盘体和散热筋板之间的温度差逐渐减小，散热筋板根部的应力也逐渐降低，而由于制动盘辐板与盘体距离较远，辐板上的温度始终较低，使得辐板在整个制动过程中的形状基本保持不变，随着盘体的膨胀，温度较低的辐板约束着盘体的径向及轴向变形，使得此处逐渐成为高应力区，在制