Thermal contact conductance conta171+targe169

[richardforever]

thread569-136203
hi!
I've got some problems to simulate the thermal contact between two faces. The model is a heatsink of a electronic chip with fins. The real problem is to simulate the thermal interface material with conta171+targe169 (it's a thin thermal grease) . The code I have implemented doesn't experiment a variation in the conductivity although I change the TCC parameter. Maybe the problem is in the mesh, or in how i have implemented that. Anyone who can help me? This is the code;
/AUTO
/GRAPHICS,POWER
/TYPE,,4
/TITLE,Ejercicio_2.Analisis de un Transistor a Potencia Total- Estacionario
/PREP7
!* Defining the Geometry
RECTNG,,0.01,0.0425,0.0496
RECTNG,,0.0225,0.04,0.0425
RECTNG,,0.0225,0.03,0.04
RECTNG,,0.0025,,0.03
RECTNG,0.0075,0.0125,,0.03
RECTNG,0.0175,0.0225,,0.03
AGLUE,ALL
SAVE
/PNUM,LINE,1
LPLOT
LFILLT,14,33,.0025, ,
LFILLT,33,20,0.0025, ,
LFILLT,18,34,0.0025, ,
LFILLT,34,24,0.0025, ,
AL,7,8,9
AL,10,11,33
AL,12,15,16
AL,22,23,34
ASEL,,LOC,Y,0,0.04
AADD,ALL
!/eof
ET,1,PLANE77
KEYOPT,1,1,0
KEYOPT,1,3,0
ET,2,SURF151
KEYOPT,2,3,0
KEYOPT,2,4,0
KEYOPT,2,5,1
KEYOPT,2,6,0
KEYOPT,2,8,2
KEYOPT,2,9,0
ET,3,TARGE169
ET,4,CONTA171
KEYOPT,4,1,2
KEYOPT,4,5,3
KEYOPT,4,9,1
KEYOPT,4,11,1
KEYOPT,4,12,4
R,1
R,2
UIMP,1,DENS, , ,3500,
UIMP,1,KXX, , ,50,
UIMP,1,C, , ,500,
!*
UIMP,2,DENS, , ,8900,
UIMP,2,KXX, , ,393,
UIMP,2,C, , ,385,
!*
UIMP,3,DENS, , ,2700,
UIMP,3,KXX, , ,156,
UIMP,3,C, , ,963,
!*
ASEL,,LOC,Y,0.0425,0.0496
AATT,1,1,1
ASEL,,LOC,Y,0.04,0.0425
AATT,2,1,1
ASEL,,LOC,Y,0,0.04
AATT,3,1,1
SMRT,3
ALLSEL
AMESH,ALL
LSEL,R,EXT
LSEL,R,LOC,X,0.001,1
LSEL,R,LOC,Y,-1,0.042
NSLL,,1
ESLN
N,5000,,-0.01
TYPE,2
REAL,2
ESURF,5000
ALLSEL
!Thermal Contact Conductance
TCC=5
FKN =0.1 !normal stiffness factor
FTOL = 0.2 !penetration tolerance
ICON = 0.2 !initial closure factor
CNOF = 1e-6 !offset value assigned to contact surface
!Mesh target surface
R,4,,,FKN,FTOL,ICON,,
RMORE,,,,0,,,
RMORE,,TCC,,,,,
ESEL,NONE
TYPE,3
ASEL,S,,,1
NSLA,S,1
NSEL,R,LOC,Y,0.0425
ESLN
ESURF
!Mesh contact surface
ESEL,NONE
TYPE,4
ASEL,S,,,2
NSLA,S,1
NSEL,R,LOC,Y,0.0425
ESLN
ESURF
ALLSEL,ALL
SAVE
FINISH
/SOLU
ANTYPE,STATIC
TIME,1
NLGEOM,ON
AUTOTS,ON
NSUBST,5,1000,1
OUTRES,ALL,ALL
LNSRCH,ON
NEQIT,1000
BFA,1,HGEN,5633802.8
ESEL,,TYPE,,2
SFE,ALL,1,CONV,,51
ALLSEL
D,5000,TEMP,50
NSEL,S,LOC,X,0
SF,ALL,CONV,0,0
NSEL,ALL
!* Visually check charges
/PBC,TEMP,1
/PBF,HGEN,,1
/PSF,CONV,2
EPLOT
SAVE
SOLVE
FINISH
/POST1
/SHOW,PNG
SET,LAST
ESEL,S,TYPE,,1
ALLSEL,BELOW,ELEM
NSLE,S
PLNSOL,TEMP,,0