Static structural problem with thermal loads

[Julia123]

Hello,

I am trying to compute the stresses in a turbine shroud under a thermal load. I am doing it in ANSYS Workbench with a static structural box and adding the thermal loads with APDL commands. I also created a fixed support to ground in order to avoid rigid body movements. The APDL code for the thermal loads is at the end of the post.

The result that I obtain for the stresses is more or less the expected one, however, the result of the temperatures is not. I obtain a really hot area in the rear part of the top surface, which is supposed to be colder than the front part.

Does anyone know whether the code or the BC are wrong? I have been thinking about this for already too long and I do not find the mistake.

Thank you in advance and let me know if the information given is not enough.


Julia


------------------------------
APDL CODE

SYS,0
! Temperature distribution in the hot area
*DIM,tab,TABLE,4,2,1,X,TIME
! Time Values
tab(0,1,1) = 0.
tab(0,2,1) = 1.
! X values
tab(1,0,1) = -0.25018
tab(2,0,1) = -0.20692
tab(3,0,1) = -0.20222
tab(4,0,1) = -0.19992
! Load values
tab(1,1,1) = 1060
tab(1,2,1) = 1060
tab(2,1,1) = 1300
tab(2,2,1) = 1300
tab(3,1,1) = 1300
tab(3,2,1) = 1300
tab(4,1,1) = 1255
tab(4,2,1) = 1255
!D,noditos,TEMP,%tab%
BF,noditos,TEMP,%tab%

! Convection in the hot surface
*DIM,convvari35,TABLE,1,2,1,TEMP,,TIME
! Time Values
convvari35(0,1,1) = 0.
convvari35(0,2,1) = 1.
! Time values
convvari35(1,0,1) = 1100
! Load values
convvari35(1,1,1) = 5500
convvari35(1,2,1) = 5500.
*DIM,loadvari35,TABLE,1,1,1,TIME,
! Time values
loadvari35(1,0,1) = 1100
! Load values
loadvari35(1,1,1) = 1100
SF,hot,CONV,%convvari35%,%loadvari35%

! Convection in the sides
*DIM,tab3,TABLE,1,2,1,TEMP,,TIME
! Time Values
tab3(0,1,1) = 0.
tab3(0,2,1) = 1.
! Time values
tab3(1,0,1) = 800
! Load values
tab3(1,1,1) = 1000
tab3(1,2,1) = 1000.
*DIM,tab4,TABLE,1,1,1,TIME,
! Time values
tab4(1,0,1) = 800
! Load values
tab4(1,1,1) = 800
SF,sides,CONV,%tab3%,%tab4%

! Convection in the impingement pocket
*DIM,tab5,TABLE,1,2,1,TEMP,,TIME
! Time Values
tab5(0,1,1) = 0.
tab5(0,2,1) = 1.
! Time values
tab5(1,0,1) = 800
! Load values
tab5(1,1,1) = 4600
tab5(1,2,1) = 4600.
*DIM,tab6,TABLE,1,1,1,TIME,
! Time values
tab6(1,0,1) = 800
! Load values
tab6(1,1,1) = 800
SF,impingement,CONV,%tab5%,%tab6%

! Convection in the tubes
*DIM,tab7,TABLE,1,2,1,TEMP,,TIME
! Time Values
tab7(0,1,1) = 0.
tab7(0,2,1) = 1.
! Time values
tab7(1,0,1) = 800
! Load values
tab7(1,1,1) = 4200
tab7(1,2,1) = 4200.
*DIM,tab8,TABLE,1,1,1,TIME,
! Time values
tab8(1,0,1) = 800
! Load values
tab8(1,1,1) = 800
SF,tubes,CONV,%tab7%,%tab8%
NSEL,ALL

 

[Julia123]

I also tried with this code:

! Convection in the hot surface
SF,hot,CONV,5500,1000

! Convection in the sides
SF,sides,CONV,1000,800

! Convection in the impingement pocket
SF,impingement,CONV,4600,800

! Convection in the tubes
SF,tubes,CONV,4200,800

! Temperature distribution in the hot area
*DIM,tab,TABLE,4,2,1,X
! X values
tab(1,0,1) = -0.25018
tab(2,0,1) = -0.20692
tab(3,0,1) = -0.20222
tab(4,0,1) = -0.19992
! Load values
tab(1,1,1) = 1060
tab(2,1,1) = 1300
tab(3,1,1) = 1300
tab(4,1,1) = 1255
BF,noditos,TEMP,%tab%

NSEL,ALL

But I got exactly the same solution. Does someone have any ideas?