FEMOD2 - 2D Interactive Finite Element Code

FEMOD2 is an interactive program for the setup, solution, and printout of a boundary value problem in a two dimensional region. The finite element method is used.

Reference:

Eric Becker, Graham Carey, J Tinsley Oden,
Finite Elements, An Introduction, Volume I.
Prentice-Hall, Englewood Cliffs, New Jersey, 1981.

Files you may copy include:

• femod2.f90, the source code;

The list of routines includes:

• FEMOD2 is a two dimensional interactive or batch finite element code.
• ADDMAT adds a given value to the band matrix A in position (I,J).
• C_CAP capitalizes a single character.
• C_EQI is a case insensitive comparison of two characters for equality.
• C_TO_DIGIT returns the integer value of a base 10 digit.
• COMRPN can translate formulas you type in and evaluate them.
• DIGIT_TO_C returns the character representation of a decimal digit.
• EBC_APPLY enforces essential boundary conditions by modifying the system.
• EBC_READ reads the essential boundary conditions.
• EBC_WRITE prints out the essential boundary conditions.
• EKF_SET sets up an element stiffness matrix and right hand side.
• EKF_WRITE writes an element matrix and right hand side.
• ELEMENT_READ reads the element order.
• ELEMENT_WRITE prints the element definitions.
• ENORMS computes the RMS and energy norm error in the approximate solution.
• EXACT_READ reads a formula for the exact solution.
• FLUX computes flux terms.
• FUNSCL evaluates a scalar function of one or more scalar arguments.
• FUNVAL evaluates a scalar, vector or matrix valued function.
• GETMAT returns ENTRY = A(I,J).
• GK_BANDWIDTH computes the bandwidth of the global stiffness matrix.
• GKF_WRITE prints out a matrix.
• I_GCD finds the greatest common divisor of I and J.
• I_TO_S_ZERO converts an integer to a string, with zero padding.
• INDATA stores and retrieves data.
• INICOM initializes data for COMRPN.
• I_LCM computes the least common multiple of two integers.
• MALPHA returns .TRUE. if STRING contains only alphabetic characters.
• MATWRT prints out a scalar, vector or array.
• MATERIAL_READ reads coefficients K, B, and F for each material.
• NBC_APPLY enforces natural boundary conditions by modifying the system.
• NBC_READ reads the natural boundary conditions.
• NBC_WRITE prints out the natural boundary conditions.
• NODE_READ reads the node data.
• NODE_WRITE prints the coordinates of the nodes.
• PCHECK performs some simple checks on the parentheses in a formula.
• POINT_APPLY modifies the linear system by applying point loads.
• POINT_READ reads point loads.
• POINT_WRITE prints the point load data.
• POST post-processes the solution.
• PREP preprocesses the problem.
• PROC processes the problem, that is, it computes the solution.
• RPNCHK examines IRPN, looking for a complete RPN expression.
• RPNSET converts the infix formula into an RPN formula.
• RPNVAL evaluates the symbolic functions in an RPN formula.
• S_BLANK_DELETE removes blanks from a string, left justifying the remainder.
• S_BLANKS_DELETE replaces consecutive blanks by one blank.
• S_EQI is a case insensitive comparison of two strings for equality.
• S_INDEXI is a case-insensitive INDEX function.
• S_TO_R reads a real number from a string.
• SCOPY copies a vector, X, to a vector, Y.
• SEIDEL carries out the Gauss-Seidel iterative linear system solution.
• SELEM sets the node numbers and material types for each element.
• SETGKF sets up the global stiffness matrix and right hand side.
• SETINT sets the quadrature points and weights.
• SETMAT sets A(I,J) = ENTRY.
• SETNBC sets natural boundary conditions.
• SGB_CHECK checks the dimensions of a general band matrix.
• SGB_FA factors a matrix stored in LINPACK general band storage.
• SGB_SL solves a system factored by SGB_FA.
• SHAPE evaluates the shape functions as functions of the normalized
• SYMADD adds a symbol name to the list of symbolic names.
• SYMLST prints out information about a variable.
• SYMVAL sets, evaluates, or deletes a variable.
• TOKENS parses a character string into recognized symbols and constants.
• XSITOX converts coordinates in the master XSI system to the X, Y system.

Return to the FORTRAN software page.