ANSYS计算的接触热阻数量级有问题

wxwieandy

基于接触热阻的定义：R=(T1-T2)/q , 使用 ANSYS APDL计算了两个接触界面的接触热阻（沿厚度方向一维导热，不考虑热对流和热辐射），计算的界面接触热阻的数量级：10-7，但是参阅的几篇文献计算的接触热阻的的数量级都是10-4，我建模时用的单位是微米，导热系数的单位也已经转换为微米，但是计算的接触热阻转换为国际制单位以后，对应的数量级为10-7，而不是10-4。不知道自己的问题出在什么地方，请大家帮忙分析一下，谢谢！

计算代码如下：
/PREP7
Ln=127.2
Lr=10.6
RECTNG,0,890.4,0,1000
WPOFFS,,1006
CSYS,4
ARSYM,Y,1
WPOFFS,,-6
CSYS,4
*DO,i,1,7,1
a_i1=1/2*(2*i-1)*Ln-1/2*Lr
a_i2=1/2*(2*i-1)*Ln+1/2*Lr
K,2*i+7,a_i1,0,0
K,2*i+8,a_i2,0,0
*ENDDO
KSEL,S,LOC,Y,0
KSEL,U,KP,,3
KSEL,U,KP,,4
KGEN,2,9,22,1,,12
*DO,j,9,22,1
L,j,j+14
*ENDDO
LSEL,S,LENGTH,,12
CM,L1,LINE
LSEL,S,,,3
LSEL,A,,,L1
LSBL,3,L1,,DELETE,KEEP
ALLS
LSEL,S,,,7
LSEL,A,,,L1
LSBL,7,L1,,DELETE,KEEP
ALLS
*DO,i,1,7,1
A,2*i+7,2*i+8,2*i+22,2*i+21
*ENDDO
ET,1,55     !定义热分析单元
UIMP,1,KXX,,,43.2E-6       !定义材料的热导率,注意单位换算为微米
UIMP,2,KXX,,,37.04E-6       !定义粗糙峰的热导率
UIMP,3,KXX,,,236E-6        !定义铝的热导率
LSEL,S,LENGTH,,10.6
LESIZE,ALL,,,3
LSEL,S,LENGTH,,12
LESIZE,ALL,,,4
LSEL,S,LENGTH,,116.6
LESIZE,ALL,,,38
LSEL,S,LENGTH,,58.3
LESIZE,ALL,,,19
ALLS
SMRT,1
MSHAPE,0,2D
MSHKEY,0
ALLS
AMESH,ALL
ALLS
LSEL,S,LENGTH,,12
ASLL,S
ESLA,S
EMODIF,ALL,MAT,2
ALLS
ASEL,S,,,1
ESLA,S
EMODIF,ALL,MAT,3
ALLS
NUMMRG,ALL
NUMCMP,ALL
ALLS
/SOLU
LSEL,S,,,1
NSLL,S,1
D,ALL,TEMP,300
LSEL,S,,,5
NSLL,S,1
D,ALL,TEMP,350
ALLS
SOLVE
!--------------------------------提取节点的温度和热流密度（高温侧）---------------
ALLS
KSEL,S,LOC,Y,12
LSLK,S,1
NSLL,S,1
*GET,NUM_NODE_H,NODE,,COUNT !提取高温侧的节点数目
*GET,NUM_MIN_H,NODE,,NUM,MIN !提取节点的最小编号
*GET,NUM_MAX_H,NODE,,NUM,MAX !提取节点的最大编号



j=1
*DO,i,NUM_MIN_H,NUM_MAX_H
*IF,NSEL(i),EQ,1,THEN   !如果节点i在当前选择集中，则
TEMP_H_j=TEMP(i)          !提取节点i的温度
*GET,FLUX_H_j,NODE,i,TF,SUM    !提取节点i的热流密度
j=j+1
*ELSE
*ENDIF
*ENDDO


TEMP_H_k=0
FLUX_H_k=0
*DO,j,1,NUM_NODE_H
TEMP_H_k=TEMP_H_j+TEMP_H_k  !计算高温侧节点的温度总和
FLUX_H_k=FLUX_H_j+FLUX_H_k  !计算高温侧节点的热流密度总和
*ENDDO



!--------------------------------提取节点的温度和热流密度（低温侧）---------------
ALLS
CSYS,4
KSEL,S,LOC,Y,0
LSLK,S,1
NSLL,S,1
*GET,NUM_NODE_L,NODE,,COUNT !提取低温侧的节点数目
*GET,NUM_MIN_L,NODE,,NUM,MIN !提取节点的最小编号
*GET,NUM_MAX_L,NODE,,NUM,MAX !提取节点的最大编号



j=1
*DO,i,NUM_MIN_L,NUM_MAX_L
*IF,NSEL(i),EQ,1,THEN   !如果节点i在当前选择集中，则
TEMP_L_j=TEMP(i)          !提取节点i的温度
*GET,FLUX_L_j,NODE,i,TF,SUM    !提取节点i的热流密度
j=j+1
*ELSE
*ENDIF
*ENDDO


TEMP_L_k=0
FLUX_L_k=0
*DO,j,1,NUM_NODE_L
TEMP_L_k=TEMP_L_j+TEMP_L_k  !计算低温侧节点的温度总和
FLUX_L_k=FLUX_L_j+FLUX_L_k  !计算低温侧节点的热流密度总和
*ENDDO
T1=TEMP_H_K-TEMP_L_K
F1=FLUX_H_K+FLUX_L_K
THER_R1=TEMP_H_K/(NUM_NODE_H)-TEMP_L_K/(NUM_NODE_L)
THER_R2=FLUX_H_K/(2*NUM_NODE_H)+FLUX_L_K/(2*NUM_NODE_L)
THER_R=THER_R1/THER_R2        !界面接触热阻