*GET CMSEL NDNEXT
CMSEL ,S ,F1
*GET ,N_COUNT ,NODE ,0 ,COUNT
*GET ,N_MIN ,NODE , ,NUM ,MIN
*DIM ,UX_DATA ,ARRAY ,N_COUNT,2
NU = N_MIN
*DO ,I ,1 ,N_COUNT ,1
UX_DATA(I,1) = NU
*GET ,UX_DATA(I,2) ,NODE ,NU ,U ,X
NU = NDNEXT(NU)
*ENDDO
UX_SUM = 0
*DO ,I ,1 ,N_COUNT ,1
UX_SUM = UX_SUM + UX_DATA(I,2)
*ENDDO
MY_N_COUNT = N_COUNT
MY_UX_SUM = UX_SUM
MY_AUX_SUM = UX_SUM/N_COUNT
*DIM ,UY_DATA ,ARRAY ,N_COUNT,2
NU = N_MIN
*DO ,I ,1 ,N_COUNT ,1
UY_DATA(I,1) = NU
*GET ,UY_DATA(I,2) ,NODE ,NU ,U ,Y
NU = NDNEXT(NU)
*ENDDO
UY_SUM = 0
*DO ,I ,1 ,N_COUNT ,1
UY_SUM = UY_SUM + UY_DATA(I,2)
*ENDDO
MY_UY_SUM = UY_SUM
MY_AUY_SUM = UY_SUM/N_COUNT
_____________________________________________________________________________________
*DIM ,D_DATA ,ARRAY ,4 ,3
INI_TEMP = 20 ! Initial Temperature 20 C
INI_HF = 1 * INI_TEMP ! Initial Heat Flux 1W at
D_DATA(1,1) = 1 ! FIRST LOAD STEP
D_DATA(1,2) = INI_TEMP ! FIRST TEMP
D_DATA(1,3) = INI_HF ! FIRST HF
F ,HF_POINT ,HF ,INI_HF
ALLSEL
SOLVE
*DO ,I ,2 ,4 ,1 ! CURRENT TIME 2, END TIME 4, INC TIME 1
/POST1
CMSEL ,S ,RT_POINT
*GET ,N_NUM ,NODE , ,NUM ,MIN
*GET ,R_TEMP ,NODE ,N_NUM ,TEMP ! Result Temperature at RT_POINT(Named Selection list)
ALLSEL
/SOLU
TIME ,I
E_HF = 1*R_TEMP ! EACH Heat Flux 1W at
F ,HF_POINT ,HF ,E_HF
D_DATA(I,1) = I ! EACH LOAD STEP
D_DATA(I,2) = R_TEMP ! EACH TEMP
D_DATA(I,3) = E_HF ! EACH HF
SOLVE
FDELE,ALL,ALL
*ENDDO
/OUTPUT,D_DATA,TXT
*VWRITE,D_DATA(1,1),D_DATA(1,2),D_DATA(1,3)
(f8.0,2x ,f8.0,2x ,f8.0)
/OUTPUT
____________________________________________________________________________________________
/PREP7
*DIM ,D_DATA , ,5 ,4
*DIM ,T_EMIS ,TABLE ,3 ,1 ,1
T_EMIS(1,0) = 22 ,500 ,1000 ! Temperature Range 22 to 1000 <- Linear interpolation
T_EMIS(1,1) = 0.3 ,0.8 ,1 ! Emissivity Range 0.3 to 1 <- Linear interpolation
T_EMIS(0,0) = 0
T_EMIS(0,1) = 1.0
!------------------ ! T_EMIS TABLE DATA
! | 0 ,1
!---|--------------
! 0 | 0 ,1
! 1 | 22 ,0.3
! 2 | 500 ,0.8
! 3 | 1000 ,1
!------------------
I_TEMP = 22 ! Initial Temperature
I1_EMIS = T_EMIS(I_TEMP,1) ! Initial Emissivity at I_TEMP FOR R1(Named Selection R1)
I2_EMIS = T_EMIS(I_TEMP,1) ! Initial Emissivity at I_TEMP FOR R2(Named Selection R2)
D_DATA(1,1) = I_TEMP
D_DATA(1,2) = I_TEMP
D_DATA(1,3) = I1_EMIS
D_DATA(1,4) = I2_EMIS
/SOLU
CMSEL ,S ,R1
SF ,ALL ,RDSF ,I1_EMIS ,1
CMSEL ,S ,R2
SF ,ALL ,RDSF ,I2_EMIS ,1
SPCTEMP ,1 ,22
/AUX12
STEF ,5.669e-008 ! Specifies Stefan-Boltzmann radiation constant.
TOFFST ,273.15 ! Temperature offset from absolute zero
RADOPT , ,0.0001 ,0 ,1000 ,0.1 ,0.1 ! Specifies Radiosity Solver options.
HEMIOPT,10 ! Hemicube Resolution for View Factor Calculation
CMSEL ,S ,R1
ESLN ,ALL
RSURF,,,4
CMSEL ,S ,R2
ESLN ,ALL
RSURF,,,5
VFSM , ,1 ,2
VFCALC ,F_FAC
VFOPT ,OFF
ALLSEL
FINI
/SOLU
TIME ,1
SOLV
*DO ,I ,2 ,5 ,1 ! CURRENT TIME 2, END TIME 5, INC TIME 1
CMSEL ,S ,R1
*GET ,N_N1 ,NODE , ,NUM ,MIN
*GET ,R_T1 ,NODE ,N_N1 ,TEMP ! Result Temperature at RT_POINT(Named Selection list)
CMSEL ,S ,R2
*GET ,N_N2 ,NODE , ,NUM ,MIN
*GET ,R_T2 ,NODE ,N_N2 ,TEMP ! Result Temperature at RT_POINT(Named Selection list)
ALLSEL
I1_EMIS = T_EMIS(R_T1,1) ! Initial Emissivity at I_TEMP FOR R1(Named Selection R1)
I2_EMIS = T_EMIS(R_T2,1) ! Initial Emissivity at I_TEMP FOR R2(Named Selection R2)
/SOLU
TIME ,I
CMSEL ,S ,R1
SF ,ALL ,RDSF ,I1_EMIS ,1
CMSEL ,S ,R2
SF ,ALL ,RDSF ,I2_EMIS ,1
D_DATA(I,1) = R_T1
D_DATA(I,2) = R_T2
D_DATA(I,3) = I1_EMIS
D_DATA(I,4) = I2_EMIS
SOLVE
SFDELE,ALL,ALL
*ENDDO
/OUTPUT,D_DATA,TXT
*VWRITE,D_DATA(1,1),D_DATA(1,2),D_DATA(1,3),D_DATA(1,4)
(f8.0,2x ,f8.0,2x ,f8.0,2x ,f8.0)
/OUTPUT
NSORT,U,Y ! Sorts nodes based on UY deflection *GET,DMAX,SORT,,MAX ! Parameter DMAX = maximum deflection ! Derived data for line elements are accessed through ETABLE: ETABLE,VOLU,VOLU ! VOLU = volume of each element SSUM ! Sums the data in each column of the element ! table *GET,VOLUME,SSUM,,ITEM,VOLU ! Parameter VOLUME = total volume ETABLE,Stress,LS,1 ESORT,ETAB,Stress,,1 *GET,Stress,SORT,,MAX ! Parameter SMAXJ = max. value of SMAX_J