NEA-1094 EURDYN.

1. NAME OR DESIGNATION OF PROGRAM:  EURDYN.

3. DESCRIPTION OF PROGRAM OR FUNCTION

The EURDYN computer codes are under development at JRC-Ispra since 1973 for the simulation of non- linear dynamic response of fast-reactor components submitted to impulsive loading due to abnormal working conditions. They are thus  mainly used in reactor safety analysis but can apply to other fields. Indeed the codes compute the elasto-plastic transient response of 2-D and thin 3-D structures submitted to fast dynamic loading generated by explosions, impacts... and represented by time  dependent pressures, concentrated loads and prescribed displacements, or by initial speeds.
Two releases of the structural computer codes EURDYN 01 (2-D beams and triangles and axisymmetric conical shells and triangular tores), 02 (axisymmetric and 2-D quadratic isoparametric elements) and 03 (triangular plate elements) have already been produced in 1976(1) and 1980(2). They include material (elastoplasticity using the classical flow theory approach) and geometrical (large displacements and rotations treated by a corotational technique) nonlinearities.
The present version (Release 3) has been completed mid-1982 and is documented in EUR 8357 EN.
The new features of Release 3, as compared to the former ones, roughly consist in:
- full large strain capability for 9-node isoparametric elements   (EURDYN 02),
- generalized array dimensions,
- introduction of the radial return algorithm for elasto-plastic   material modelling,
- extension of the energy check facility to the case of prescribed   displacements,
- possible interface to a post-processing package including time   plot facilities (TPLOT).
The theoretical aspects can be found in refs. 2,4,5,6,7,8.

4. METHOD OF SOLUTION

- Finite element space discretization.
- Explicit time integration.
- Lumped masses.
- EURDYN 01: 2-D corotational formulation including constant strain   triangles (plane or axisymmetric), beams and conical shells, this  last element being particularly useful for the study of thin   axisymmetric structures.
- EURDYN 02: 2-D full large strain formulation; discretization of  the structure is obtained through 9-node biquadratic isoparametric    elements (plane or axisymmetric).
- EURDYN 03: 3-D corotational formulation; discretization of the  (thin) structure is achieved by means of triangular plate   elements.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

No 3-D continuum element available (only thin 3-D structures can be treated).

6. TYPICAL RUNNING TIME

AMDAHL 470/V8 under IBM/VMS: from 200 to more than 1000 elements * time steps/CPU second, depending on problem. NEA 1094/01: The test case described on page 37 of the report has been executed by NEA-DB on an IBM 3084 computer in 0.42 seconds of CPU time.

NEA-DB executed the test case included in this package  on an IBM 370/084 computer in 0.21 seconds of CPU time.

7. UNUSUAL FEATURES OF THE PROGRAM:

8. RELATED AND AUXILIARY PROGRAMS

TPLOT/3 for post processing. TPLOT  is an interactive data management system for transient problems. It  includes TEKTRONIX PLOT-10 graphical facilities.

10. REFERENCES

- J.P. Halleux, F. Casadei:
  "TPLOT: an Interactive Data Management for Transient Problems,
  2nd edition"
  EUR- report to be published.

- J.P. Halleux:
  EURDYN: Computer Programs for the Nonlinear Transient Analysis of
  Structures Submitted to Dynamic Loading
  EUR 8357 EN  (1982)
- J. Donea and S. Giuliani:
  "EURDYN: Computer Programs for the Transient Dynamic Analysis of
  Large Displacement, Small Strain Problems with Material Non
  linearities (User's Manual)"
  EUR 5474 e (1976)
- J. Donea, S. Giuliani, J.P. Halleux:
  "Recent Improvements of the Nonlinear Transient Dynamic Structural
  Computer Programs EURDYN"
  EUR 6694 EN (1980).
- J. Donea, S. Giuliani, J.P. Halleux:
  "Theoretical Aspects of the EURDYN Computer Programs for Non
  Linear Transient Dynamic Analysis of Structural Components"
  EUR 5473 e (1976)
- J.P. Halleux:
  Modelling of Elasto-Plastic Material Behaviour.
  EUR 7329 EN  (1981)
- J. Donea, S. Giuliani, J.P. Halleux:
  "Prediction of the Nonlinear Dynamic Response of Structural
  Components Using F.E."
  Reprint from "Nuclear Engineering and Design" 37(1976),95-114
- J.P. Halleux:
  "Nonlinear Transient Analysis of Thin Axisymmetric Structures".
  Proceedings of the International Conference on Numerical Methods
  for Nonlinear Problems, University College Swansea,
  2-5 September 1980, Vol. 1, Pineridge Press  (1980).
- J.P. Halleux:
  "Transient Nonlinear Finite Element Analysis of Axisymmetric
  Structures by Lumped-Explicit Methods",
  Paper M 6/8, SMIRT 6, Paris  (August 1981).
- J.P. Halleux and F. Casadei:
  "Transient Large Strain Analysis of Plane and Axisymmetric
  Structures by Means of Biquadratic Finite Elements".
Eng. Comput. 1984, Vol.1 December

11. MACHINE REQUIREMENTS

200Kbytes + data (semi-variable dimension) disk files: maximum 7.

Main storage requirement to run the test case is 1836KB

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:  IBM-VMS.

VS2 03.8 (IBM 3084).

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS:

15. NAME AND ESTABLISHMENT OF AUTHORS

     J. Donea, S. Giuliani, J.P. Halleux
     Commission of the European Communities
     Joint Research Centre
     ISPRA Establishment, Italy

File name	File description	Records
NEA1094_06.001	EURDYN-3.3 INFORMATION FILE	68
NEA1094_06.002	EURDYN-3.3 SOURCE PROGRAM	2731
NEA1094_06.003	EURDYN-3.3 JCL	15
NEA1094_06.004	EURDYN-3.3 Sample problem input data	33
NEA1094_06.005	EURDYN-3.3 Sample problem printed output	267