Learning Forum

[askari]

Dear all,

 

Following my question in the previous forum in this link, we are modeling a pipeline under dynamic loading on a supercomputer, and need to extract the element number with maximum absolute total mechanical strain. For now, we are only getting the value of maximum and minimum strain using contour display in each time step via *GET, par, PLNSOL, 0, MAX and *GET, par, PLNSOL, 0, MIN methods.

Below, is our code for post-processing:

*** Code block begins***

! Select the pipe elements (soil springs are not selected)

ESEL,S,TYPE,,0.9,1.1

! Produce an element display

EPLOT 




! Display elements with shapes determined from the real constants etc.

/ESHAPE,1  

! Activates a coordinate system for printout or display of element and nodal results

RSYS,SOLU 




! Specifies the device and other parameters for graphics displays

/SHOW,WIN32C

! Read the first data set (subsetp and time are ignored).

SET,FIRST




! Control the graphics options on subsequent displays

/PLOPTS,INFO,3 ! Legend on

/CONTOUR,ALL,18 ! Specifies the uniform contour values on the stress displays

/PNUM,MAT,1 ! Controls entity numbering

/NUMBER,1  ! Specifies whether numbers are used for display

/REPLOT,RESIZE 




! Read the last data set (substep and time are ignored)

SET,LAST




*get,totalStep,active,,set,LSTP




*dim,ResultData,array,totalStep,3

/graph,power

/show,figure_temp.png

*DO,i,1,totalStep,1

SET,i,

plns,EPTO,X

*get,ResultData(i,1),active,,set,time

*get,ResultData(i,2),plnsol,0,max

*get,ResultData(i,3),plnsol,0,min

*enddo

Finish

*MWRITE,ResultData,Output,txt,,jik

(F14.9, F14.9, F14.9)

/show,close

/dele,figure_temp,png

*** Code block ends***

 

We appreciate your help and support in advance,