**RSICC CODE PACKAGE PSR-478**

**1. NAME AND TITLE**

GRFPAK: Code System to Plot CORTES FEM Results.

**2. CONTRIBUTOR**

Oak Ridge National Laboratory, Tennessee through the Energy Science and Technology Software Center, Oak Ridge, Tennessee.

**3. CODING LANGUAGE AND COMPUTER**

FORTRAN IV; IBM360 (P00478I36000).

**4. NATURE OF PROBLEM SOLVED**

GRFPAK is a graphics package written for the CORTES finite-element programs. It includes three plotting routines to assist in analyzing, interpreting, and presenting CORTES results. One plotting routine displays stresses or stress indices by means of contour curves drawn within either the outside or inside surface outline of a quarter section of the tee-joint. Using this routine, one can also obtain plots of the finite-element mesh as viewed from any point in space. A second plotting routine gives a stress versus distance plot along any specified line of nodes. A third routine displays cross-sectional views of the finite-element mesh for both the undeformed (original) and deformed configuration. The deformed configurations are drawn using an exaggerated scale specified by the user.

**5. METHOD OF SOLUTION**

All of the plotter software calls are channeled through one of the several multiple entries in subroutine PLTINT. The segment of GRFPAK designed for contour, mesh, and displacement plotting can produce the finite-element mesh for the original geometry and for each load case of the displaced geometry, and contour plots for any surface in any quadrant for any stress (or temperature). These contours may be displayed on the mesh or on the outline of the tee. The node line plot segment produces the plots of the longitudinal (along the node line) stresses and transverse (normal to the node line) stresses along the 0 degree node line (x-y plane) and 90 degree node line (y-z plane) for a given surface. The global (x,y,z) coordinates for the nodes along the 0 degree node line are given a 45 degree rotation to a new x'-y' coordinate system. The abscissa values in the plots represent x' along this line. For the 90 degree node line the abscissa values for the plots represent the y global coordinate values. For THFA, the temperature values along these node lines are plotted. If more than one time-step is being plotted for THFA or TSA, all of the time-steps will appear on one graph. The cross-section segment of GRFPAK produces plots of the x=0 and z=0 planes of the tee for SA and EP. The entire cross section or only a specified number of nodes along the axes from the control node may be represented.

**6. RESTRICTIONS OR LIMITATIONS**

Requests for plots relating to displacements are invalid for the CORTES programs TSA, THFA, and SHFA. The maximum value for the variable MTOT in the MAIN program is 22,000.

**7. TYPICAL RUNNING TIME**

Execution time is approximately 12 seconds for sample data which produces nine plots. NESC executed the sample problem in 35 seconds on an IBM370/195.

**8. COMPUTER HARDWARE REQUIREMENTS**

GRFPAK was developed on IBM 360 and requires card reader, tape unit, line printer, and CalComp plotting equipment. The sample problem which consists of 675 elements required 290K bytes of memory.

**9. COMPUTER SOFTWARE REQUIREMENTS**

The software was developed under OS/360 and requires a Fortran IV compiler.

**10. REFERENCE**

P.G. Fowler and J.W. Bryson, "User's Manual for the CORTES Graphics Package GRFPAK," ORNL/NUREG/TM-127 (August 1977).

**11. CONTENTS OF CODE PACKAGE**

Included are the reference document and software on one 3.5" diskette which includes source and sample input with IBM 360 JCL sample file. No executable is included with this package.

**12. DATE OF ABSTRACT**

April 2000.

** KEYWORDS:** FINITE ELEMENT METHOD; PLOTTING