2D Interactive Finite Element Code

aspen

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.

普之仁

你好，能否将你05年6月4日公布的这份材料再公布一次，谢谢