RSICC CODE PACKAGE PSR-487

1. NAME AND TITLE

SAMCR: 2-D Code System for Elastodynamic Fracture Analysis.

2. CONTRIBUTOR

University of Maryland, College Park, Maryland, through the Energy Science and Technology Software Center, Oak Ridge, Tennessee.

3. CODING LANGUAGE AND COMPUTER

Fortran V; UNIVAC1100 (P00487U110000).

4. NATURE OF PROBLEM SOLVED

SAMCR is a two-dimensional, elastodynamic, finite element code for the stress analysis of moving cracks. The code can be used to provide useful information about the Mode I fracture behavior of bodies under a combination of mechanical, thermal, and pressure loadings for materials that permit a linear-elastic fracture mechanics assumption. Major features of SAMCR are: efficient explicit finite element formulation based on four-node isoparametric quadrilateral elements; crack advance calculations based on a crack tip restraining force model; J-integral formulation for calculation of the dynamic stress intensity factor K; flexibility in the specification of the crack tip velocity versus K fracture constitutive relationship; thermal effects, both for thermal strains and for the crack tip velocity versus K constitutive relationship; energy calculations (strain, fracture, kinetic, and dissipated/damping) as a monitor of solution accuracy; automatic time step selection; static analysis using dynamic relaxation; and an analysis restart option. The code has been shown to perform well in modelling dynamic behavior of both uncracked and cracked structures and has been demonstrated to provide useful information regarding run-arrest events in polymeric laboratory samples and large thermally shocked steel cylinders. A graphics postprocessor, TEKPLOTS, is included which using data from SAMCR output files can generate the following plots: stress intensity factor versus time or crack tip location, crack tip location versus time, crack tip velocity versus time or crack tip location, energy components versus time, nodal displacements versus time, nodal accelerations versus time, and element strains versus time or normalized crack tip location.

5. METHOD OF SOLUTION

The mathematical formulation is based on the conventional variational formulation of finite element theory using the principle of virtual works. Four-noded quadrilateral isoparametric elements are used together with explicit time integration and a restraining nodal force model of incremental crack advance. The code operates in a predictive fracture application mode, i.e., the increment of crack extension at each time step is automatically computed based on the calculated stress intensity value K at the crack tip and the user-specified crack tip velocity versus K constitutive relationship for the fracturing material.

6. RESTRICTIONS OR LIMITATIONS

None noted.

7. TYPICAL RUNNING TIME

NESC executed the sample problem in less than 4 CP minutes on the UNIVAC1100/62.

8. COMPUTER HARDWARE REQUIREMENTS

68,000 words of memory are required for execution on the Univac 1100/62. TEKPLOTS requires a Tektronix 401x graphics terminal for generating the graphics displays.

9. COMPUTER SOFTWARE REQUIREMENTS

SAMCR ran under the EXEC 1100 operating system. A Fortran V compiler is required. SAMCR was tested when it was initially released by NESC in August 1982. It was not tested when it was transferred to RSICC and released in July 2001.

10. REFERENCES

C.W. Schwartz, R. Chona, W.L. Fourney, and G.R. Irwin, "SAMCR: A Two-Dimensional Dynamic Finite Element Code for the Stress Analysis of Moving Cracks," NUREG/CR-3891, ORNL/Sub/79-7778/3 (November 1984).

Usage Notes for Program "TEKPLOTS."

11. CONTENTS OF CODE PACKAGE

Included are the reference document and software on one-3.5" DOS HD diskette containing a PC WINDOWS ® self-extracting executable in PkWare ® format.

12. DATE OF ABSTRACT

July 2001.

KEYWORDS: FRACTURE MECHANICS; REACTOR SAFETY