CN106529036B - 一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法 - Google Patents
一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法 Download PDFInfo
- Publication number
- CN106529036B CN106529036B CN201610987608.7A CN201610987608A CN106529036B CN 106529036 B CN106529036 B CN 106529036B CN 201610987608 A CN201610987608 A CN 201610987608A CN 106529036 B CN106529036 B CN 106529036B
- Authority
- CN
- China
- Prior art keywords
- contact
- resistance
- thermal
- micro
- bulge
- 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.)
- Active
Links
Images
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
-
- 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/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法,该方法考虑了空气介质热阻对接触热阻的影响。该方法根据微凸体的弹性变形和塑性变形计算结合面实际接触面积和接触载荷,然后分别计算基体热阻、收缩热阻和空气介质热阻,基体热阻和收缩热阻串联再与空气介质热阻并联计算出总接触热阻,最后使用Matlab编写计算程序得到接触热阻和载荷的关系。本发明的特点在于考虑了微凸体空气介质热阻的影响,载荷较小时空气介质热阻影响较大,不能忽略。本发明提供的方法可为电主轴热态分析的边界条件接触热阻的计算提供理论依据。
Description
技术领域
本发明属于电主轴热态特性研究领域,涉及一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法,该方法运用matlab计算并分析了基体热阻、收缩热阻和空气介质热阻的影响。
背景技术
电主轴是数控机床的关键部件之一。电主轴的热特性是影响机床加工精度的重要因素之一。目前,接触热阻的建模方法主要有传统的赫兹接触模型、基于经典力学和统计的G-W模型、基于统计学参数的W-A接触模型、还有基于W-M函数的M-B模型,前三个模型具有尺度依赖性,受仪器分辨率和取样长度的影响,后一个模型具有全面性、确定性和尺度独立性,但没有综合考虑微凸体的基体热阻、收缩热阻和空气介质热阻。本发明考虑了微凸体的基体热阻、收缩热阻和空气介质热阻的影响。
发明内容
本发明的目的是提供一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法,该方法运用分形理论建立接触热阻模型并考虑了微凸体的基体热阻、收缩热阻和空气介质热阻的影响。该方法首先建立结合面实际接触面积和接触载荷方程,然后建立基体热阻、收缩热阻和空气介质热阻从而建立总接触热阻模型,最后使用Matlab编写计算程序得到总接触热阻,弹性阶段接触热阻和间隙介质热阻随结合面外在载荷变化情况。
本发明是采用以下技术手段实现的:
S1、首先用分形理论对整个粗糙表面进行理论描述,单个微凸体进行基体热阻、收缩热阻的理论分析,得到微凸体的接触面积、接触载荷以及临界接触面积。将接触面积和接触载荷积分得到实际接触面积和接触载荷。
S2、考虑空气的热传导建立空气介质热阻模型,两者并联得到接触热阻模型。
S3、按照计算流程编写Matlab程序计算接触热阻随载荷变化曲线图。
本发明的特点在于考虑了微凸体的基体热阻、收缩热阻和空气介质热阻影响,能够对各个热阻进行精确计算。本发明提供的方法可为电主轴的热态特性分析边界条件接触热阻的计算提供指导。
通过下面的描述并结合附图说明,本发明会更加清晰,附图说明用于解释本发明方法及实施例。
附图说明
图1整个粗糙表面的分形表征;
图2基体热阻、收缩热阻和空气间隙热阻的热阻网络;
图3收缩热阻时的单点接触;
图4分形维数D=1.4976分形粗糙度参数G=8×10-10m时,基体热导、收缩热导和无量纲间隙热导随接触载荷的变化关系图;
具体实施方式
本发明实施一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法,下面结合附图,对本发明的实施进行具体说明。
当两个粗糙表面相互接触时,被假设为许多大小不等的圆形微凸体的相互接触,而同一表面微凸体间的变形影响予以忽略。粗糙表面可被看作由大量的、离散的、相互并联的小圆柱形微凸体组成,如图1所示。因此,粗糙表面的彼此接触问题可被简化为接触面积为a′接触高度为d′的圆柱微凸体间的相互接触。接触高度为d′可得:
d′=Z′max-δmax (1)
式中Z′max样本中最高点到最低点的距离,可给出Z′max=L(GL)D-1。L——样本长度。δmax微凸体的最大变形:
当热流通过相互接触的粗糙表面时,接触热阻由基体热阻、收缩热阻和小间隙热阻组成。两个相互接触的粗糙表面的热阻网络如图2所示。基体热阻rbi与收缩热阻rci相互串联形成接触的小热阻单元。间隙热阻rgi与接触的热阻单元并联。对于本接触热阻模型,存在下列假设:(1)热流通道之间无热交换;(2)所有热流通道界面温差相同;(3)通过空隙中的流体介质的传热方式为热传导,不考虑热辐射与热对流。
1.1接触微凸体总接触面积和接触压力
当两个粗糙表面相互接触时,由于相互接触的表面微凸体受到彼此的挤压,从而使微凸体产生弹性或者塑性变形。相互接触的粗糙表面的弹性模量可由等效弹性模量获得,v1,v2,E1,E2分别是两个表面的泊松比和弹性模量。对于接触的一个微凸体的等效半径可表示为R=R1R2/(R1+R2),这里R1和R2分别为相互接触的小微凸体的半径。
真实接触面积和总接触载荷可由积分获得。
微凸体截面积a′的分布函数为:
式中a′微凸体变形后的截面积,当微凸体为弹性变形时a′=2a;a′L最大微凸体的截面积;ψ描述微观接触时微凸体大小分布的域扩张系数,可由超越方程(ψ(2-D)/2-(1+ψ-D/2)-(2-D)/D)/((2-D)/D)=1得。
区分弹性变形状态和塑性变形状态的临界微凸体面积可给出:
式中b=(πq/2)2,q=0.454+0.41ν1;ν1较软材料的泊松比;σy较软材料的屈服强度;K硬度系数,通常K=2.8。
微凸体截面积a′>a′c时为弹性变形阶段,微凸体截面积a′≤a′c时为塑性变形阶段。
对于弹性变形和塑性变形状态下的微凸体,单个微凸体间接触力f与截面积a′的关系可给出:
fp=Kσya′ (6)
式中e、p分别代表弹性变形和塑性变形状态。
整个表面的负载F和真实接触面积Ar可由积分获得:
1.2单个热阻单元的基体热阻
由分形接触模型可知,对于接触面积a的微凸体,在接触高度d′范围内的基体热阻可表达为:
式中k当量导热系数,k=2k1k2/(k1+k2);k1,k2分别为两接触材料的导热系数。
根据弹性理论,微凸体最大变形量可得:
1.3单个热阻单元的收缩热阻
由于接触界面的实际接触面积远小于名义接触面积,热流通过接触界面时,热流线只通过这些离散的接触点会产生热流线收缩现象,即收缩热阻。根据Cooper,Mikic和Yovanovich等人提出的经典截锥体接触模型(如图3所示),单个微凸体接触点处收缩热阻为:
对于接触面积a的微凸体,相应的收缩热阻为:
由于基体热阻与收缩热阻是并联关系,热阻单元的接触热阻可表示为:
1.4间隙中空气介质热阻
对于间隙中空气介质热阻传热情况。对于粗糙界面的小间隙传热,一个无量纲的判别参数克努森数可定义:
d为接触平面平均间隙高度,对于相互接触的两个平面中:
Λ为间隙气体分子平均自由程:
式中T开氏温度;P接触间隙气体压强,可视为标准大气压;δ气体分子直径;kB玻尔兹曼常数,kB=1.38×10-23J/k。
对于常温常压下,间隙为空气时的分子平均自由程Λ=0.064μm。
当kn<0.01时为气体连续传热状态,平板之间热传导发生主要通过气体分子之间的碰撞和能量交换,可适用傅里叶传热定律:
当0.01<kn<10时为温度跳跃状态,平板/气体界面之间的能量交换不完全导致这些界面温度不连续。这时间隙介质的接触热导可以定义为:
式中kg间隙气体介质的导热率。
根据Song and Yovanovich,气体参数M为:
这里α1、α2是气体对两表面间的热适应系数,对于两相同材料微凸体α1=α2;γ是介质气体的比热容比,即气体的定压比热容(Cp)与定容比热容(Cv)的比值,常压常温可取1.40;Pr是气体的普朗特数。
Song and Yovanovich给出了一个气固结合的热适应系数的预测关系式:
当kn>10时为气体自由分子传热状态,分子间的碰撞几乎没有,平板间的传热机制主要是气体分子和平板间的能量交换,这时间隙介质热导为:
1.5粗糙表面的接触热阻
粗糙表面的接触热阻由许多热阻单元与间隙热阻组成。整个粗糙表面的接触热阻数值可由其积分获得。为了获得接触表面的接触热阻这里引入接触热导h,其与接触热阻的关系为h=1/r。单个热阻单元的接触热导由可得:
整个粗糙表面的接触热导由接触单元的热导与间隙热导串联而成:
Hc=Hbc+Hcg (25)
式中,Hcg=1/rgAg=kg/[(M+d)Ag],Ag为间隙传热面积。
上述模型可使用matlab编程来实现。为更具体的说明本方法的有效性,本发明提供了一个计算实例。
以铝合金/铝合金为对象对接触热阻的倒数接触热导进行分析,铝合金的物理性参数。如表1所示。
表1 T=293K时接触材料的特性
根据表1中的参数和建模过程给出的参数代入模型,分形维数D=1.4976分形粗糙度参数G=8×10-10m时的分析结果,微凸体分别发生弹性、塑性变形的收缩热导和间隙热导以及总的收缩热导随接触载荷的变化,可以看出当载荷小时,间隙热阻占很大比例,热流传递受间隙热阻的阻碍,因此不能忽略。
Claims (1)
1.一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法,该方法运用分形理论建立接触热阻模型并考虑了微凸体的基体热阻、收缩热阻和空气介质热阻的影响;该方法首先建立结合面实际接触面积和接触载荷方程,然后建立基体热阻、收缩热阻和空气介质热阻从而建立总接触热阻模型,最后使用Matlab编写计算程序得到总接触热阻,弹性阶段接触热阻和间隙介质热阻随结合面外在载荷变化情况;
其特征在于,
S1、首先用分形理论对整个粗糙表面进行理论描述,单个微凸体进行基体热阻、收缩热阻的理论分析,得到微凸体的接触面积、接触载荷以及临界接触面积;将接触面积和接触载荷积分得到实际接触面积和接触载荷;
S2、考虑空气的热传导建立空气介质热阻模型,两者并联得到接触热阻模型;
S3、按照计算流程编写Matlab程序计算接触热阻随载荷变化曲线图;
本方法在于考虑了微凸体的基体热阻、收缩热阻和空气介质热阻影响,能够对各个热阻进行精确计算;本方法为电主轴的热态特性分析边界条件接触热阻的计算提供指导;
当两个粗糙表面相互接触时,被假设为许多大小不等的圆形微凸体的相互接触,而同一表面微凸体间的变形影响予以忽略;粗糙表面被看作由大量的、离散的、相互并联的小圆柱形微凸体组成;因此,粗糙表面的彼此接触问题被简化为微凸体接触面积为a′接触高度为d′的圆柱微凸体间的相互接触;接触高度为d′得:
d′=Z′max-δmax (1)
式中Z′max样本中最高点到最低点的距离,给出Z′max=L(G/L)D-1;L——样本长度;δmax微凸体的最大变形:
当热流通过相互接触的粗糙表面时,接触热阻由基体热阻、收缩热阻和小间隙热阻组成;基体热阻rbi与收缩热阻rci相互串联形成接触的小热阻单元;间隙热阻rgi与接触的热阻单元并联;对于本接触热阻模型,存在下列假设:(1)热流通道之间无热交换;(2)所有热流通道界面温差相同;(3)通过空隙中的流体介质的传热方式为热传导,不考虑热辐射与热对流;
1.1接触微凸体总接触面积和接触压力
当两个粗糙表面相互接触时,由于相互接触的表面微凸体受到彼此的挤压,从而使微凸体产生弹性或者塑性变形;相互接触的粗糙表面的弹性模量由等效弹性模量获得,v1,v2,E1,E2分别是两个表面的泊松比和弹性模量;对于接触的一个微凸体的等效半径表示为R=R1R2/(R1+R2),这里R1和R2分别为相互接触的小微凸体的半径;
真实接触面积和总接触载荷由积分获得;
微凸体截面积a′的分布函数为:
式中a′微凸体变形后的截面积,当微凸体为弹性变形时a′=2a;a′L最大微凸体的截面积;ψ描述微观接触时微凸体大小分布的域扩张系数,由超越方程(ψ(2-D)/2-(1+ψ-D/2)-(2-D)/D)/((2-D)/D)=1得;
区分弹性变形状态和塑性变形状态的临界微凸体面积给出:
式中b=(πq/2)2,q=0.454+0.41ν1;ν1较软材料的泊松比;σy较软材料的屈服强度;K硬度系数,K=2.8;
微凸体截面积a′>a′c时为弹性变形阶段,微凸体截面积a′≤a′c时为塑性变形阶段;
对于弹性变形和塑性变形状态下的微凸体,单个微凸体间接触力f与截面积a′的关系给出:
fp=Kσya′ (6)
式中e、p分别代表弹性变形和塑性变形状态;
整个表面的负载F和真实接触面积Ar由积分获得:
1.2单个热阻单元的基体热阻
由分形接触模型知,对于接触面积a的微凸体,在接触高度d′范围内的基体热阻表达为:
式中k当量导热系数,k=2k1k2/(k1+k2);k1,k2分别为两接触材料的导热系数;
根据弹性理论,微凸体最大变形量得:
1.3单个热阻单元的收缩热阻
由于接触界面的实际接触面积小于名义接触面积,热流通过接触界面时,热流线只通过这些离散的接触点会产生热流线收缩现象,即收缩热阻;提出的经典截锥体接触模型,单个微凸体接触点处收缩热阻为:
对于接触面积a的微凸体,相应的收缩热阻为:
由于基体热阻与收缩热阻是并联关系,热阻单元的接触热阻表示为:
1.4间隙中空气介质热阻
对于间隙中空气介质热阻传热情况;对于粗糙界面的小间隙传热,一个无量纲的判别参数克努森数定义:
d为接触平面平均间隙高度,对于相互接触的两个平面中:
Λ为间隙气体分子平均自由程:
式中T开氏温度;P接触间隙气体压强,视为标准大气压;δ气体分子直径;kB玻尔兹曼常数,kB=1.38×10-23J/k;
对于常温常压下,间隙为空气时的分子平均自由程Λ=0.064μm;
当kn<0.01时为气体连续传热状态,平板之间热传导发生通过气体分子之间的碰撞和能量交换,适用傅里叶传热定律:
当0.01<kn<10时为温度跳跃状态,平板/气体界面之间的能量交换不完全导致这些界面温度不连续;这时间隙介质的接触热导定义为:
式中,kg间隙气体介质的导热率;
气体参数M为:
这里α1、α2是气体对两表面间的热适应系数,对于两相同材料微凸体α1=α2;γ是介质气体的比热容比,即气体的定压比热容Cp与定容比热容Cv的比值,常压常温取1.40;Pr是气体的普朗特数;
气固结合的热适应系数的预测关系式:
当kn>10时为气体自由分子传热状态,分子间没有碰撞,平板间的传热机制是气体分子和平板间的能量交换,这时间隙介质热导为:
1.5粗糙表面的接触热阻
粗糙表面的接触热阻由热阻单元与间隙热阻组成;整个粗糙表面的接触热阻数值由其积分获得;为了获得接触表面的接触热阻这里引入接触热导h,其与接触热阻的关系为h=1/r;单个热阻单元的接触热导由得:
整个粗糙表面的接触热导由接触单元的热导与间隙热导串联而成:
Hc=Hbc+Hcg (25)
式中,Hcg=1/rgAg=kg/[(M+d)Ag],Ag为间隙传热面积;
弹性接触热导模型、塑性接触热导模型、间隙空气介质热导模型和总接触热导模型使用matlab编程来实现。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610987608.7A CN106529036B (zh) | 2016-11-09 | 2016-11-09 | 一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610987608.7A CN106529036B (zh) | 2016-11-09 | 2016-11-09 | 一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106529036A CN106529036A (zh) | 2017-03-22 |
CN106529036B true CN106529036B (zh) | 2019-10-15 |
Family
ID=58350477
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610987608.7A Active CN106529036B (zh) | 2016-11-09 | 2016-11-09 | 一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106529036B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018229995A1 (ja) * | 2017-06-16 | 2018-12-20 | 三菱電機株式会社 | ブレーキライニング |
CN109271673B (zh) * | 2018-08-22 | 2023-03-24 | 江苏大学 | 一种温差发电片等效热导率的计算方法 |
CN110705147B (zh) * | 2019-09-18 | 2023-08-04 | 北京工业大学 | 一种数控机床主轴热态特性综合性理论建模与分析方法 |
CN112966376B (zh) * | 2021-03-02 | 2023-05-26 | 桂林电子科技大学 | 一种基于分形理论的接触热导建模方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366050A (zh) * | 2013-06-25 | 2013-10-23 | 华中科技大学 | 一种高速电主轴力-热耦合建模方法 |
CN104978465A (zh) * | 2015-07-16 | 2015-10-14 | 北京工业大学 | 一种考虑微凸体的弹塑性变形和空气介质热阻的接触热阻建模方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2012004216A (ja) * | 2010-06-15 | 2012-01-05 | Denso Corp | 電子制御装置 |
-
2016
- 2016-11-09 CN CN201610987608.7A patent/CN106529036B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366050A (zh) * | 2013-06-25 | 2013-10-23 | 华中科技大学 | 一种高速电主轴力-热耦合建模方法 |
CN104978465A (zh) * | 2015-07-16 | 2015-10-14 | 北京工业大学 | 一种考虑微凸体的弹塑性变形和空气介质热阻的接触热阻建模方法 |
Non-Patent Citations (1)
Title |
---|
基于接触热阻的高速精密电主轴热特性分析;刘志峰 等;《北京工业大学学报》;20160131;第42卷(第1期);17-23 * |
Also Published As
Publication number | Publication date |
---|---|
CN106529036A (zh) | 2017-03-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106529036B (zh) | 一种考虑微凸体的基体热阻、收缩热阻和空气介质热阻的接触热阻建模方法 | |
Hopmann et al. | Extrusion dies for plastics and rubber: design and engineering computations | |
Krishnan et al. | Study of the size effect on friction conditions in microextrusion—part I: microextrusion experiments and analysis | |
Han et al. | A new thermal error modeling method for CNC machine tools | |
CN104978465B (zh) | 一种考虑微凸体的弹塑性变形和空气介质热阻的接触热阻建模方法 | |
Yang et al. | Modified Elman network for thermal deformation compensation modeling in machine tools | |
Vera et al. | Factors influencing microinjection molding replication quality | |
Dats et al. | On a multi-physics modelling framework for thermo-elastic-plastic materials processing | |
Kumar et al. | A novel methodology of Combined Compromise Solution and Principal Component Analysis (CoCoSo-PCA) for machinability investigation of graphene nanocomposites | |
Mai-Duy et al. | Dissipative particle dynamics modeling of low Reynolds number incompressible flows | |
Iqbal et al. | Characterization of machining of AISI 1045 steel over a wide range of cutting speeds. Part 2: evaluation of flow stress models and interface friction distribution schemes | |
Li et al. | Calculation method of convective heat transfer coefficients for thermal simulation of a spindle system based on RBF neural network | |
Pandey et al. | Investigations into complete liquefier dynamics and optimization of process parameters for fused deposition modeling | |
Correia et al. | Analysis of friction in the ejection of thermoplastic mouldings | |
Gupta et al. | Investigations of performance parameters in NFMQL assisted turning of titanium alloy using TOPSIS and particle swarm optimisation method | |
Juul et al. | Steady-state numerical modeling of size effects in micron scale wire drawing | |
Huang et al. | Fabrication of microfluidic chip using micro‐hot embossing with micro electrical discharge machining mold | |
Wang et al. | Effect of indirect cryogenic cooling on the machining accuracy and tool vibration in the turning of polysulfone | |
Günel et al. | Investigation of heat sink geometry effect on cooling performance for a fdm 3d printer liquefier | |
Zhu et al. | Research on influencing factors of the thermal contact resistance in thin-wall injection molding | |
Sönmez et al. | Numerical simulation and process optimization of a 3D thin-walled polymeric part using injection compression molding | |
Yardimeden | Investigation of optimum cutting parameters and tool radius in turning glass-fiber-reinforced composite material | |
Muñoz-Sánchez et al. | Hybrid FE/ANN and LPR approach for the inverse identification of material parameters from cutting tests | |
CN111027243B (zh) | 一种基于数值模拟的等温挤压能耗影响因素分析方法 | |
Rodríguez García et al. | Experimental measurement of the thermal conductivity of fused deposition modeling materials with a DTC-25 conductivity meter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 |