ANSYS APDL-TRANSIENT THERMAL ANALYSIS
Hello Everyone,
I am currently working on ANSYS APDL since a few weeks in order to perform coupled Thermomechanical analysis. At first I have started with the transient thermal analysis on a simple block. Firstly, I have retrieved all the elements and their respective centroid locations(X,Y,Z) to a text file. Secondly, I need to perform transient thermal analysis on all the elements individually in a specific element order using the EKILL and EALIVE technique. The elements which are already saved in the text file are retrieved using *TREAD and are respectively called in a DO loop. I would like to visualize the temperature distributions, fluxes. But I am facing some issues with the SFE command, tables and also calling the element inside the DO loop. So any kind of help or advice is highly appreciated.
ERRORS: 1) Invalid element number = 0. The SFE command is ignored
2) The birth and death capability requires that the NROPT,FULL command be specified before the first SOLVE command and before the first EKILL command. The EKILL command is ignored.
3) Present time 0 is less than or equal to the previous time 0 in a transient analysis.
/clear
/prep7
block,,2,,0.2,,0.1
ET,1,186
MP,ex,1,2.4
MP,PRXY,1,0.4
MP,KXX,1,0.19
MP,C,1,1.62
lsel,all
lesize,all,,,5,,,,1
allsel
esel,s,ELEM,,1,,1
esll,s
vmesh,all
alls
*cfopen,text_file,txt
*GET,elem_no,ELEM,0,COUNT ! number of elements
*GET,elemno,ELEM,0,NUM,MIN ! first element label
*Do,i,1,elem_no,1
*GET,EX,ELEM,elemno,CENT,X
*GET,EY,ELEM,elemno,CENT,Y
*GET,EZ,ELEM,elemno,CENT,Z
*VWRITE,elemno,EX,EY,EZ ! writes the elemno to the text file with the centroid cordinates
(F8.0,',',E10.3,',',E10.3,',',E10.3)
*GET,elemno,ELEM,elemno,NXTH
*ENDDO
*CFCLOSE
alls
*cfopen,Newtext_file,txt
*DIM,TEMP_,TABLE,5
*DIM,ALPHA_,TABLE,6
*TAXIS,TEMP_(1),2,0.0,25.0,50.0,75.0,100.0
TEMP_(0),344.0,388.0,432.0,476.0,520.0
*TAXIS,ALPHA_(1),2,0.0,20.0,40.0,60.0,80.0,100.0
ALPHA_(0),80.0,81.0,82.0,83.0,84.0,85.0
TRNOPT,full
NROPT,full
EKILL,ALL
*Do,j,new_par,,1
EALIVE,new_par
*DIM,new_par,TABLE,125,4
*TREAD,new_par,text_file,txt
/SOLU
ANTYPE,4
AUTOTS,ON
DELTIM,10,2,100
TIME,6300
TUNIF,520
KBC,0
SFE,ALL,1,CONV,1,%ALPHA_%
SFE,ALL,1,CONV,2,%TEMP_%
*GET,HF1,ELEM,0,NXTH
*VWRITE,new_par(j),HF1
(F8.4,',',E10.3)
solve
*ENDDO
*CFCLOSE
Finish
Best Regards,
Srikanth purli,
M.Sc Mechatronics,
University of Siegen,
Germany.
[cid:03b3d0c8-6944-4e0c-851d-dd18ff23efb9]
Srikanth,
Save this to a file and read it in as input. It should work for the first part of your code. It helps me to put data into a matrix before writing to a file, so I can check the output of code for issues.
/clear
/prep7
block,,2,,0.2,,0.1
ET,1,186
MP,ex,1,2.4
MP,PRXY,1,0.4
MP,KXX,1,0.19
MP,C,1,1.62
lsel,all
lesize,all,,,5,,,,1
allsel
esel,s,ELEM,,1,,1
esll,s
vmesh,all
allsel
*GET,elem_no,ELEM,0,COUNT
*dim,myarr,array,elem_no,4
*Do,i,1,elem_no,1
*GET,EX,ELEM,i,CENT,X
*GET,EY,ELEM,i,CENT,Y
*GET,EZ,ELEM,i,CENT,Z
myarr(i,1)=i
myarr(i,2)=EX
myarr(i,3)=EY
myarr(i,4)=EZ
*ENDDO
*cfopen,text_file,txt,'C:'
*VWRITE,myarr(1),myarr(2),myarr(3),myarr(4)
(F8.0,E10.3,E10.3,E10.3)
*CFCLOSE
I hope it helps.
Scott Richardson
(980) 345-1344 (desk)
(980) 201-0522 (mobile)
Husqvarna Group, R&D Center of NA
8825 Statesville Road | Charlotte, NC 28269 | USA
-----Original Message-----
From: Purli, Srikanth [mailto:srikanth.purli@student.uni-siegen.de]
Sent: Monday, May 22, 2017 4:27 AM
To: xansys-temp@xansystest.info
Subject: [Xansys] ANSYS APDL-TRANSIENT THERMAL ANALYSIS
Importance: High
Hello Everyone,
I am currently working on ANSYS APDL since a few weeks in order to perform coupled Thermomechanical analysis. At first I have started with the transient thermal analysis on a simple block. Firstly, I have retrieved all the elements and their respective centroid locations(X,Y,Z) to a text file. Secondly, I need to perform transient thermal analysis on all the elements individually in a specific element order using the EKILL and EALIVE technique. The elements which are already saved in the text file are retrieved using *TREAD and are respectively called in a DO loop. I would like to visualize the temperature distributions, fluxes. But I am facing some issues with the SFE command, tables and also calling the element inside the DO loop. So any kind of help or advice is highly appreciated.
ERRORS: 1) Invalid element number = 0. The SFE command is ignored
2) The birth and death capability requires that the NROPT,FULL command be specified before the first SOLVE command and before the first EKILL command. The EKILL command is ignored.
3) Present time 0 is less than or equal to the previous time 0 in a transient analysis.
/clear
/prep7
block,,2,,0.2,,0.1
ET,1,186
MP,ex,1,2.4
MP,PRXY,1,0.4
MP,KXX,1,0.19
MP,C,1,1.62
lsel,all
lesize,all,,,5,,,,1
allsel
esel,s,ELEM,,1,,1
esll,s
vmesh,all
alls
*cfopen,text_file,txt
*GET,elem_no,ELEM,0,COUNT ! number of elements *GET,elemno,ELEM,0,NUM,MIN ! first element label
*Do,i,1,elem_no,1
*GET,EX,ELEM,elemno,CENT,X
*GET,EY,ELEM,elemno,CENT,Y
*GET,EZ,ELEM,elemno,CENT,Z
*VWRITE,elemno,EX,EY,EZ ! writes the elemno to the text file with the centroid cordinates
(F8.0,',',E10.3,',',E10.3,',',E10.3)
*GET,elemno,ELEM,elemno,NXTH
*ENDDO
*CFCLOSE
alls
*cfopen,Newtext_file,txt
*DIM,TEMP_,TABLE,5
*DIM,ALPHA_,TABLE,6
*TAXIS,TEMP_(1),2,0.0,25.0,50.0,75.0,100.0
TEMP_(0),344.0,388.0,432.0,476.0,520.0
*TAXIS,ALPHA_(1),2,0.0,20.0,40.0,60.0,80.0,100.0
ALPHA_(0),80.0,81.0,82.0,83.0,84.0,85.0
TRNOPT,full
NROPT,full
EKILL,ALL
*Do,j,new_par,,1
EALIVE,new_par
*DIM,new_par,TABLE,125,4
*TREAD,new_par,text_file,txt
/SOLU
ANTYPE,4
AUTOTS,ON
DELTIM,10,2,100
TIME,6300
TUNIF,520
KBC,0
SFE,ALL,1,CONV,1,%ALPHA_%
SFE,ALL,1,CONV,2,%TEMP_%
*GET,HF1,ELEM,0,NXTH
*VWRITE,new_par(j),HF1
(F8.4,',',E10.3)
solve
*ENDDO
*CFCLOSE
Finish
Best Regards,
Srikanth purli,
M.Sc Mechatronics,
University of Siegen,
Germany.
[cid:03b3d0c8-6944-4e0c-851d-dd18ff23efb9]