Explicit Analysis RADIOSS Ebook

Download as pdf or txt
Download as pdf or txt
You are on page 1of 438

A Study Guide

1
2
Table of Contents

1 About This Study Guide ........................................................... 9

2 HyperWorks For Teaching .......................................................15

3 Explicit Analysis - Some Fundamentals ...................................17


3.1 Introduction ...................................................................................................... 17

3.2 Full Dynamic / Impact Vs Quasistatic Simulations ............................................ 21

3.3 What Do We Solve In Structural Crashworthiness? .......................................... 22

3.4 Typical Application of Crashworthiness Simulations In Various Industries ...... 23

3.5 Numerical Modeling Approaches...................................................................... 26

3.6 Time Step Calculation in RADIOSS..................................................................... 28

3.7 Hourglass ........................................................................................................... 30

3.8 Some Points To Note While Running A Dynamic Simulation ............................ 31

4 RADIOSS Files .........................................................................35


The Starter Input Deck ............................................................................................ 37

Engine File ............................................................................................................... 39

5 Frequently Asked Questions (RADIOSS Files) ..........................45


5.1 How To Set The Number Of Processors For A Run ........................................... 45

3
5.2 Parallelism ......................................................................................................... 47

6 Element Formulations ............................................................51


6.1 1D Elements (Line elements) ............................................................................ 51

6.2 Overview of Element Library – 1D Elements .................................................... 53

6.3 2D Elements (Shell Elements) ........................................................................... 61

6.4 3D Elements (Solid Elements) ........................................................................... 73

6.5 Frequently Asked Questions (Elements) ........................................................... 82

6.6 Frequently Asked Questions (Property) ............................................................ 87

7 Example: S - Beam Crash.........................................................89


Physical Problem Description.................................................................................. 90

Analysis, Assumptions and Modeling Description .................................................. 91

8 Materials .............................................................................. 101


8.1 Isotropic Elasticity ........................................................................................... 102

8.2 Isotropic Elasto-Plastic .................................................................................... 103

8.3 Composite and Anisotropic ............................................................................. 116

8.4 Viscous ............................................................................................................ 117

8.5 Hydrodynamic ................................................................................................. 120

8.6 Failure Models................................................................................................. 121

8.7 Explosives ........................................................................................................ 121

8.8 Generating Material Data ............................................................................... 122

8.9 Frequently Asked Questions (Materials/Failure) ......................................... 123

4
9 Example: Law Characterization - Tensile Test........................ 125

10 Interface / Contact Modeling .............................................. 149


10.1 Frequently Asked Questions (Contact Interfaces) ..................................... 166

11 Kinematic Conditions .......................................................... 169


11.1 Frequently Asked Questions (Kinematic Conditions).................................... 186

12 Load Definition ................................................................... 193


12.1 Frequently Asked Questions (Initial Stresses) ............................................... 199

13 Output Requests................................................................. 201


13.1 Frequently Asked Questions (Post-processing) ............................................ 209

14 Example: Three Point Bending (HyperMesh) ....................... 213

15 Example: Three Point Bending (HyperCrash) ....................... 227

16 Time Step Control ............................................................... 261


16.1 Frequently Asked Questions (Time Step) ...................................................... 272

5
17 Helpful Tools ...................................................................... 275

18 Checks During & After Analysis ........................................... 283

19 Best Practice ....................................................................... 301

20 Postprocessing - Crashworthiness....................................... 309


20.1 Measures – Distance Between And Position ................................................ 309

20.2 Section Cuts................................................................................................... 313

20.3 Vector Plots ................................................................................................... 319

20.4 Tracking Systems ........................................................................................... 324

20.5 Tracing ........................................................................................................... 329

20.6 Exploded View ............................................................................................... 331

20.7 Synchronizing Data And Windows ................................................................ 333

20.8 Overlaying Images and Videos ...................................................................... 339

20.9 Exporting A Deformed Shape ........................................................................ 351

20.10 Crash Tools For Plotting .............................................................................. 352

21 Exercise: Bird Strike Analysis on An Aircraft Leading Edge .. 359

22 Example: Crash Tube .......................................................... 371

6
23 Example: Bumper Analysis .................................................. 381

24 Example: Honeycomb Crash-Box ........................................ 383

25 Capstone Project: Cell Phone Drop ..................................... 409

26 Capstone Project: Roof Crush ............................................. 413

27 Appendix: Recommended Videos ....................................... 417

28 Appendix: Conversion Tables And Constants ...................... 433

7
8
1 About This Study Guide
This study guide aims to provide a basic introduction into the exciting and challenging
world of explicit Finite Element Analysis. As with our other eBooks we deliberately
kept the theoretical aspects as short as possible. Instead, we concentrated on
providing information a new RADIOSS user typically needs (or is interested in) to get
started.
The chapter entitled: Explicit Analysis - Some Fundamentals is meant to provide a
rather high-level introduction into the topic. More detailed information about e.g.
element formulations in RADIOSS, time step control, or interfaces/ contact
formulations is then provided in subsequent chapters.

In this new edition two additional comprehensive examples have been added: “S-
Beam Crash” and “Law Characterization Tensile Test”. Also new is the example “Three
Point Bending (HyperMesh)” - the story presented therein is the same as in “Three
Point Bending (HyperCrash)”. However, in the latter the model setup is based on
HyperCrash.
We moved the exercises related to the modeling and analysis of a “Crash Tube”, a
“Crash Bumper”, a “Honeycomb Crash -Box” (of a student racing car) to the end of
the book. You will see that these examples are a kind of repetition of the information
provided in earlier chapters. Hence, working your way through these examples should
be rather straight forward.

9
Certainly, more demanding are the Capstone Projects “Cell Phone Drop” and “Roof
Crush”. These projects will definitely test your RADIOSS skills. All the model files for
the tutorials in this book can be downloaded from here.

As stated above the focus of this book is on the analysis. A basic understanding of the
preprocessor HyperMesh (or HyperCrash) makes your life much easier. The same
accounts for postprocessing of results with HyperView and HyperGraph. May be the
appendix “Postprocessing - Crashworthiness” is all you need to know. If not, then
please have a look at the following free eLearning courses available in our Learning &
Certification Program (https://certification.altairuniversity.com/).

• Learn Post Processing with HyperView (topics addressed: HyperWorks Desktop


Environment, Animation and View Controls, Strength Analysis, HyperGraph Plotting,
Modal & Frequency Response (NVH) Analysis, Crashworthiness, Publishing Results)
• Learn Preprocessing with HyperMesh (topics addressed: Interface and Basic
Interaction, Geometry, 2D Meshing, Solids and Solid Meshing, 1D Meshing and
Connectors, Analysis Set Up)
• In addition, the free study guide: “Practical Aspects of Finite Element Simulation”
(www.altairuniversity.com/academic-program/) may also give you a good
introduction.
Of course, there is also an eLearning course about RADIOSS available
• Learn Explicit FEA with RADIOSS (topics addressed: Introduction to RADIOSS Elements,
Summarizing Explicit FEA Process with a Crash Box; Introduction to Material Library
and Contact Modeling, Introduction to Rigid Walls and Boundary Conditions, Time
Step Control, Model Set Up and Analysis, Debugging Techniques, Fundamentals of
Crash Analysis).

Additionally, we offer:
• Highly discounted seminars & workshops at colleges and/or at Altair facilities
• An Altair moderated (http://forum.altairhyperworks.com/) and much more ...

10
Now, before you get started, please recall that regardless of the project complexity
you are working on, the reliability of your results is only as good as the person who
uses the software: The solver is only capable of working with the problem that you
define; it cannot design the part for us and does not have engineering judgment, so
the ultimate design decisions lie with YOU. A very popular statement which
summarizes the above (and which is mentioned in many documents) simply says:
Garbage in, Garbage out
Eventually, as this study guide is “Work in Progress” the more advanced topics and
technologies such as composites and composite failure / damage etc. will be added in
subsequent editions. So stay tuned ...

And now - enjoy this study guide and let us know whether it helped you to successfully
solve your explicit FEA projects.
Best regards
Dr. Matthias Goelke
On behalf of “The HyperWorks University Team”

11
Acknowledgement
“If everyone is moving forward together, then success takes care itself”
Henry Ford (1863 -1947)
A very special Thank You goes to:
• Nitin Gokhale from “Finite to Infinite” for all his passion about CAE, the inspirational
collaboration, and the friendship
• Sanjay Nainani for reformatting and updating the book
• George P. Johnson (Altair India) for the thorough review, for many suggestions and
edits which significantly improved the reading and learning experience!
• Rahul Rajan (Altair India) who recorded close to 5 hours of additional eLearning
material (all embedded in this eBook) and Rahul Ponginan, Prakash Pagadala and
Premanand Suryavanshi for their very valuable support (Altair India)
• Andy Barr, John Brink, Paul Sharp and Erwan Mestres for the “Intro to RADIOSS for
Impact” training material which is included to a large extend in this eBook (Altair
USA)
• Sean Putman and Elizabeth White for stimulating discussions (Altair USA)
• Hossein Shakourzadeh, Dario Mendolicchio, Jean Michel Terrier, Lionel Morancay,
Jean-Baptiste Mouillet, Erwan Mestres (Altair France)
• Marian Bulla & Annemarie Bulla, Irina Komarova, Vladimir Trickov, Jan
Grasmannsdorf, Kristian Holm, Juergen Kranzeder, Jacob Tremmel, Moritz Günther,
Nicolas Torgau and Marco Wenzel for discussions and advice (Altair Germany)
• The entire HyperWorks Documentation Team (Altair USA) for putting together 1000’s
of pages of documentation - Thank you very much.
• Lastly, the entire RADIOSS development team deserves huge credit for your passion
& dedication!

12
Disclaimer
Every effort has been made to keep the book free from technical as well as other
mistakes. However, publishers and authors will not be responsible for loss, damage in
any form and consequences arising directly or indirectly from the use of this book.
© 2018 Altair Engineering, Inc. All rights reserved. No part of this publication may be
reproduced, transmitted, transcribed, or traslated to another language without the
written permission of Altair Engineering, Inc. To obtain this permission, write to the
attention Altair Engineering legal department at:
1820 E. Big Beaver, Troy, Michigan, USA, or call +1-248-614-2400.

13
14
2 HyperWorks For Teaching
Leading universities across the globe are using HyperWorks computer aided
engineering (CAE) simulation software for teaching and research in the fields of:
• Structural Analysis
• Computational Fluid Dynamics (CFD)
• Optimization
• Multi-Body Simulation (MBS)
• Electro-Magnetic Field Simulation
• Numerical Methods & System Simulation
Altair has commercial expertise to share with the academic community. By including
real life scenarios in your teaching material, Altair can help you add value to your
engineering design courses.
Our unique licensing system allows universities to use the entire (full version)
HyperWorks suite in a very flexible and cost-efficient way as sketched out below:

Since teaching differs from campus to campus, and from region to region we are
very interested to discuss your needs with you personally.
Please let us know your requirements by sending an email notification to
[email protected]
We are more than happy helping and assisting you with your teaching activities.

15
16
3 Explicit Analysis - Some
Fundamentals
This chapter is based on the fabulous book:
Practical Finite Element Analysis by Nitin Gokhale et al. from Finite To Infinite. We -
the Altair University Team - are very thankful for the courtesy to re-use some of his
material in this Study Guide!
We also included material from our RADIOSS Standard Seminar Classes. More
technical information is available in the HyperWorks Help Documentation:
HyperWorks Solvers >> RADIOSS >> Theory Manual.

3.1 Introduction

Past few decades have seen an increasing application of CAE for simulation of crash
phenomenon particularly due to the development of high computing machines and
parallel computing techniques. The increase in safety standards can be attributed to
the improvement of structural crashworthiness performance through Finite Element
Analysis.
The effect of crash and impacts on structures is one problem and the second one
which is of prime importance is the safety of occupants. We find that occupant safety
simulation offers reasonably accurate results today which can save a lot of testing
time and overall design cycle time. Today’s dummies have closer dynamic properties
and also include the correct load carrying capacity to allow interaction with the
structure as compared to earlier dummy models which didn’t take into account loads
entering into the body. The CAE development for these applications were delayed due
to unavailability of high end computing power and it can be said that such simulations
are barely 20 years old. Although there has been a tremendous increase in application
of software to solve problems related to automotive industry, aerospace industry and

17
drop test of components, still the best is yet to come and research has been very
active in all aspects of this application ranging from basic physical phenomena
understanding to development of efficient numerical algorithms and finally the
development of a general-purpose software which can solve almost all kinds of
problems. We now find the simulation of tests prescribed by the ECE and FMVSS
standards and regulations quite common application. Other simulations widely used
are the bird strike impact analysis and jet engine blade containment analysis
pertaining to FAR in aerospace industry.
The main aim of this chapter is to introduce basics of structural crashworthiness
simulations in terms of the solution methodology and discuss some of the practical
simulations performed in the industry.

Types Of Nonlinearity

There are three basic sources of nonlinearity: Geometric nonlinearity, Material


nonlinearity and Contact nonlinearity.

Geometric nonlinearities can be related to large strain and large displacement.


Geometric nonlinearity takes into account that the geometry cross section may
change as a result of large deformation, where as in for linear static analysis the cross
section is considered as constant.

Material nonlinearity is associated with the inelastic behavior of a component or


system. Material nonlinearities occur when the stress exceeds the yield stress or
force-displacement law is not linear, or when material properties change with the
applied loads. The material will no longer behave in a linear fashion. The causes for
material non-linearity which can be permanent deformations and cracking.

Contact nonlinearity arises when boundary conditions in a FE model change during


the course of the analysis. The boundary conditions could be added to or removed
from the model due to boundary nonlinearity as the analysis progresses.

18
Implicit And Explicit Analysis

Implicit and Explicit methods are approaches used in numerical analysis for obtaining
numerical approximations to the solutions of time-dependent ordinary and partial
differential equations. The implicit approach is useful in problems where time
dependency of the solution is not an important factor [like static, modal analysis etc.]
whereas for Explicit Dynamics approach is most helpful in solving high deformation
time dependent problems [Crash, Impact etc].

When An Explicit Analysis…

An explicit dynamics analysis is used to determine the dynamic response of a structure


due to stress wave propagation, impact or rapidly changing time-dependent loads.
Momentum exchange between moving bodies and inertial effects are usually
important aspects of the type of analysis being conducted. This type of analysis can
also be used to model mechanical phenomena that are highly nonlinear.
Nonlinearities may stem from the materials, (hyperelasticity, plastic flows, failure),
from contact (high speed collisions and impact) and from the geometric deformation
(buckling and collapse). Events with time scales of less than 1 second (usually of order
milliseconds) are efficiently simulated with this type of analysis. Based on all these
conditions the user will select an explicit solver to solve the analysis. Typical examples
of explicit analysis are crash analysis, drop test, bird strike, fan blade out…etc.

19
Advantages/Disadvantages

20
3.2 Full Dynamic / Impact Vs Quasistatic
Simulations
The distinction between these two is quite obvious when we see the load application
in following figure, a suddenly applied load (mathematically a heavy side or step input)
will result in a displacement two times that of the static displacement of the system.
The term Quasi Static simulation is used when we use a dynamic code to produce a
static result. It must be remembered that we will get a static or close to static result
only if the ramp up time is sufficiently large (i.e. the ramp up slope is not too steep so
as to get a close result as same as of that as Full dynamic simulation).
The term Quasi Static is used to represent a slow dynamic process as opposed to a
fast-dynamic process.

Full Dynamic or Impact Simulation Quasi Static Simulation

A dynamic code always produces oscillations in the result and the remedy is to reduce
this by using a larger ramping.

21
3.3 What Do We Solve In Structural
Crashworthiness?
These are basically the equilibrium equations of transient dynamics. These are written
as:

(1)

Where, M is the mass matrix, C is the damping matrix and K is the stiffness matrix.
All that we do is just determine the evolution of the basic quantities such as
displacement, velocity and acceleration given initial conditions on displacement and
velocity with respect to time.
All other quantities can be derived from these and most important are the element
stresses, plastic strains, contact forces and the energies such as kinetic, potential
energy and overall energy absorption characteristics popularly called as energy
management.
Most software would commonly solve the dynamic equilibrium equation in an implicit
way but the most popular way that should be used for highly nonlinear problems is to
use explicit time integration scheme such as a central difference scheme. There are
several advantages of such a procedure and the most important is that it leads to an
algorithm which can be easily programmed, does not require any matrix inversion
procedure and further is extremely suitable for a fast-parallel computing
Methodology.

22
Why & When To Run An Explicit Analysis (14 minutes)

(Video by Rahul Rajan, Altair India)

3.4 Typical Application Of Crashworthiness


Simulations In Various Industries
Automotive
Automotive industry has probably the widest application of crash simulation.
Simulating the crashworthiness of the vehicle in terms of very simple models based
on the spring mass damper systems was the focus when the computers were very
slow. Nowadays software such as RADIOSS and others have very wide practical
aspects such as use of special seat belt elements and development of dummies for
occupant safety. Several standards have originated in various countries and it is
impossible to give a detailed account of these. In the following we discuss some of the
most essential applications related to the automotive industry. There are very
detailed procedures such as five-star ratings and no attempt will be made here to take
a detailed account of these.

23
This is an image displaying the von
Mises stress contour plot of an
instrument panel for a Maruti car,
using HyperView. The analysis
carried out was for head form
impact in accordance with the
European safety standard ECE R -21. T Safety of the occupant is determined from a
HIC value which can be calculated using HyperView.
Head Form Impact Analysis Of Instrument Panel (Image source: Altair Calendar 2005,
Courtesy: Maruti Udyog Ltd.)
Most common standards in USA are FMVSS (Federal Motor Vehicle Safety Standards)
regulations and ECE (Economic Commission of Europe) regulations in Europe. In India,
ARAI testing standards are used now.
FMVSS standards can be obtained free of cost from the website and there are three 3
series listed in the Vol 49 Code of Federal regulations. The 100 series deals with active
safety or crash avoidance and 200 series forms the most important aspects of
crashworthiness tests. Listed in series 200 are several standards and an attempt has
been made to give the reader a brief survey of them. One can refer to the information
available online on NCAP (New car assessment program), IIHS (Insurance Institute for
Highway Safety) and NHTSA (National Highway Traffic safety association).

24
The interested reader is referred to the NHTSA website for a detailed discussion of
occupant safety related crash criteria description.

Consumer Goods Industry / Communications Industry

Crashworthiness also finds a lot of applications in drop test of components such as


television, plastic buckets and mobile phone. Points of interest here are to check the
structural integrity of the component and monitor any damage caused to the system.
In mobile industry any damage caused to the antenna and the LCD display are very
important as they make the device totally useless.

25
Applications In Other Industries

Although developed mainly for automotive applications, crash simulation software


have also found applications in train, ship and aircraft crashworthiness. The two main
standards associated with FAA (Federal Aviation Administration) requirements are
those of bird strike impact and engine blade containment. Other applications in
defence sector are simulating the explosive detonation process and design of
weapons. Computational Biomechanics also is continuously evolving with the
development of finite element models closely following the actual physical models.

3.5 Numerical Modeling Approaches


Based on the applications, there are various approached for numerical modeling. This
difference in the numerical modelling is due to the different approaches used for solid
mechanics and fluid mechanics. In solid mechanics, we always follow a particle
tracking or Lagrangian approach that is the Finite element mesh keeps on deforming
with the structure. Here the mesh continuously follows the deformed structure and
hence this approach has certain limitation with respect to the very large deformation
problems. Whereas for fluid mechanics the point of interest is put on a control volume
bound by control surfaces and the user always monitor the mass momentum and
energy exchange across it. This approach is also called the Eulerian approach. Thus,
the mesh remains fixed all the time in space and the user don’t have any severe or
stringent limitations on the mesh quality and fluid flow in a confined space or
controlled area is a typical example for Eulerian modelling.
The difference can be clearly observed in the following figure. In Lagrangian meshes,
element boundaries remain coincident with boundaries and material interfaces. In
Eulerian, they don’t and then you will have to use other approximate or tracking
methods for the treatment of moving boundaries.
The physics thus dictates the following requirements.

26
• Lagrangian meshes for slight to moderate deformation.
• Eulerian meshes for large deformation.
It is also possible to use another type of formulation called Arbitrary Lagrangian
Eulerian (ALE) formulation in which the nodes can be programmed to move arbitrarily.
The nodes on the boundaries are moved on the boundaries itself while the interior
nodes are moved to minimize the mesh deformation.

The above image clearly shows the movement of mesh in each case. In Lagrangian
mesh, the mesh deforms with the body. Both the nodes and the material points
change position as the body deforms. In Eulerian model, we can see that the body
flows through the mesh as it deforms. The nodes remain fixed and the materials
points move through the mesh.

Smooth Particle Hydrodynamics (SPH)

The principal idea of SPH is to treat hydrodynamics in a completely mesh-free fashion,


in terms of a set of sampling particles. This is a mesh free method and there is no
physical contact between the particles. The behavior of SPH particles after an impact

27
is shown in the image below. We can see that the object is splashing as soon as it
contacts with the floor. The typical examples of SPH formulations are bird strike
analysis, fuel sloshing…etc. Smoothed particle hydrodynamics is a remarkably
versatile and simple approach for numerical fluid dynamics

3.6 Time Step Calculation In RADIOSS


When solving dynamic problems with the Finite Element Method we use FEM only for
the spatial discretization and the temporal (Time) discretization is always by using the
Finite Difference Method. The total response time is divided into much smaller time
intervals called time steps or increments. The equilibrium equations are solved and
the value of unknowns are determined at (t + Δ t) based on the knowledge of their
values at time t.

What Is Time Step?


In the context of explicit analysis, time step is the time required for a shock wave
(governed by the speed of sound) to propagate across the smallest distance of the
element. A run is a single model simulation and it starts at time zero and continues till
the user defined endpoint is reached. The basic time unit in the run is the time step.

28
Calculating Time Step

From the above relation, it is clear that both material modulus and density control the
speed of sound. The solution time step is inversely proportional to the speed of sound.
So, if the material is stiffer (higher modulus of elasticity), the speed of sound increases
resulting in a lower solution time step. Thus, to keep a stable solution, the shock wave
through a medium must not pass across more than one element per time step. This
implies that the length of the element is the modeling parameter that affects the
overall solution time. In some cases, the run will be terminated due to very low time
step. There are various methods to control the time step which will be discussed in
the coming chapters.

29
Time Step Control Of Explicit Methods For Quasi-Static
Simulation

In order to increase the time step, we have following options:


• Increase the element length (Always preferred at least for initial runs)
• Increase the density of the component thereby increasing the mass of component.
• Decrease elasticity moduli of the material.
Option 2 is called as the “Mass Scaling“. Although this increases the computational
speed considerably, care must be taken to see that a bound is put on the percentage
of mass increase. Secondly it should not change the physics of the problem i.e. mass
increase is allowed for deformable bodies but if applied to rigid bodies then results
will be different than the original model. Option 3 is not used in practice.
Conventional way is to use implicit codes for quasi static simulation but the major
problem is issues of convergence and efficiency of the code when friction / sliding
contact is present. Usually it has been found that the contacts work better in an
explicit code rather than an implicit code.
Mass scaling is a much-accepted technique among the CAE Community to increase
the speed of the explicit software.

3.7 Hourglass
Hourglass is weird modes of deformation that occur in under integrated elements and
no stresses will be produced. This can affect solution accuracy by varying the
structure’s true response. This leads to inaccurate stress, strain and deflection results.
For an under-integrated element, displacements and forces may exist for each node
of the element. If sum of all these forces and displacements give null strain and stress
on the integration point of the element then that element fails to calculate stiffness
for certain mode of deformation. At this stage we say element has entered into
hourglass. Soon after this as element failed to calculate stiffness the element will
behave weirdly. There are ways to control hourglassing and in the subsequent units
we will discuss on it.

30
3.8 Some Points To Note While Running A
Dynamic Simulation
• Always mention the system of units which you have used for simulation.
• You must have a fairly good idea of what would be the total file size of the output
since a wrong input could generate a large output of ascii and binary files, and you
will run out of space very soon.
• A thumb rule to remember is: If T is the termination time, then your binary output file
generation time step should be T/10 (so that you have minimum 10 frames available)
and the ascii file size should be a further ten times lesser as you would require an
accurate representation of energy, contact forces and stresses data.
• There should be a zero penetration in the model so that there are no initial contact
related energies in the model. It is always a good practice to run a 0 msec run for the
model making you understand if there are gross modeling errors in the model.
• It is always better to set the global element length to 5 mm which corresponds to a
time step of 1 e-06 seconds. Try to make your initial runs with coarse mesh so that
you know how your model behaves.
• There are lot of element formulations and different material models, and a beginner
could naturally be very confused with understanding the software overall, it is hence
a good idea to understand the element formulations that perform best and why. This
comes through a lot of verification runs, reading the software manuals carefully and
experimental validation of the material models.

31
• Just defining the contacts and leaving to the software to do the rest of the things is
not enough. You must also see whether it is giving you a physically compatible
behavior or result.
• Usually a crash simulation software has different cards suitable for rigid bodies and
deformable bodies. Carefully study these options so that you get a full idea of the
capabilities of the software.
• Try to use restarts to the maximum extent as this is a powerful utility associated with
explicit dynamics code.
• Always check the CAE results with some hand calculations and make a note if it makes
some sense. Always check the hourglass energy if this is large enough than the internal
energy (strain energy) of the system, then it can mask the physical failure and can give
a totally erroneous result. A thumb rule is that the hourglass energy should be a
maximum of 10 % of the internal energy.
Note: By using reduced integration you save a lot of computational time which is
always necessary in transient dynamics. (Any philosophy which makes your run fast
should be used e.g. lumped mass matrix instead of consistent mass matrix). The risk
of using reduced integration is that it generates non –physical zero strain energy
modes which are to be compulsorily eliminated by using artificial stiffness / viscosity
thus generating additional strain energy which we call as hourglass energy. Another
strategy to make a fast run is to use mass scaling which increases the time step but
you should set reasonable limitations on the maximum mass increase.
• There can be many errors and mistakes by a beginner. Always remember that the
cause of these errors is very much present in the deck itself. Common are wrongly
formatted input, initial contact penetration, improper load curve definitions and
massless nodes due to temporary nodes created not deleted. Floating point
exceptions can be caused by parts with zero density or thickness, over constrained
nodes. Always take a closer look whether you have inputted the correct information
and whether the mesh is of a good quality. Normal run doesn’t necessarily guarantee
that you have solved the problem in a correct way as what is more important is always
checking the CAE simulation result with the test data or physical behavior of the
system.

32
• If the time step is too small due to the presence of a small element, it is always better
to coarsen the mesh thus allowing a considerable increase in the time step.
• Total energy of the system is always constant:
Total Energy = Internal Energy + Kinetic energy + External Work
For most of the impact problems, there will be initial kinetic energy and no external
work at time t=0.0. As the simulation advances Kinetic energy will decrease, Internal
energy will increase and the total energy of the system should always remain
Constant. If the total energy plot shows a very big jump then the model has an error
and you must check the definition of contacts. The momentum printout also shows
whether the bodies after hitting are going in the right direction. Always ensure that
parts are going in the right direction with right velocities at the right time.

References
1. Ted Belytschko, Wing Kam Liu and Brian Moran (2000): Nonlinear Finite Elements for
continua and structures, John Wiley and Sons
4. Y. C Fung (1965): Foundations of Solid Mechanics, Prentice Hall of India
5. J. N. Reddy (2006): An introduction to the Finite Element Method, Tata McGrawHill
6. T. J. R. Hughes (2000): The Finite Element Method, Linear Static and dynamic Finite
Element Analysis, Dover Publications
7. Jones A Zukas, Theodore Nicholas, Hallock F Swift, Longin B Greszczuk, Donald R
Curran (1992): Impact Dynamics, Krieger Publishing Company, USA
8. R. Cook, D. S. Malkus and M. E. Plesha (1989): Concepts and applications of Finite
Element Analysis, John Wiley and Sons
9. R. D. Cook (1995): Finite Element Modeling for Stress Analysis, John Wiley and Sons
10. Roy A Craig (1981): Structural Dynamics, An introduction to computer Methods, John
Wiley and Sons
11. Rao V Dukkipati, M A Rao, Rama Bhat (2000): Computer Aided Analysis and Design of
machine elements, New Age International Publishers

33
12. Code of Federal Regulations: Title 49 Transportation, Published by the office of the
Federal Register. National Archives and Records Administration, sale by US Govt.
Printing office Washington
13. Maurice Petyt (1990): Introduction to Finite element Vibration analysis, Cambridge
University Press
14. The following websites contain a lot of additional information
http://www.nhtsa.dot.gov (FMVSS Standards can be downloaded from this website)
http://crash-network.com/ (ECE regulations can be downloaded from this website)
http://www.faa.gov (The Federal Aviation Administration regulations can be obtained
from this website) http://www.euroncap.com http://www.iihs.org

34
4 RADIOSS Files
Preprocessing of RADIOSS models may be carried out using HyperCrash (i.e. you
typically import a meshed model into HyperCrash and then add definition of
materials, contacts, loads etc.)

Graphical User Interface Of HyperCrash


Alternatively, the entire model set up, starting from geometry cleanup-meshing, mesh
quality etc. may be completed in HyperMesh. Eventually, you decide in favor of your
preprocessor.

Graphical User Interface Of HyperMesh

35
Hence, the analysis file will be either exported from within HyperMesh or HyperCrash.

Which Files Do We Need At All?

Note, as depicted below, you need two input files: Starter Input Deck and Engine Input
File.

36
The Starter Input Deck

Starter Block Format

Each block defines one feature: a set of nodes or elements, materials, properties,
interfaces…
Each block starts with a “/” (slash) followed by a keyword and ends at the beginning
of the next block
The content of each block is entered as formatted lines of 100 characters divided in
10 fields of 10 characters.

37
Block Format Part Definition /PART

A Part in RADIOSS is a set of elements with the same property id (prop_ID) and
material id (mat_ID). In one part all elements are of the same type (e.g. all shells). The
subset_ID is optional. If omitted, the part belongs to the global model

Sample Card Image

38
RADIOSS Starter Listing File

Running A RADIOSS Simulation

Engine File
Generates animation & time history output files (Annn Tnnn). Outputs details of the
computation (runname_0001.out)
Generates runname_0001.rst file for restart(using runname_0002.rad engine file)

39
40
Typical Engine File Keywords

Running A Restart

Example: You wish to continue running a simulation that has finished to a new longer
termination time…

41
Running RADIOSS: Option 1 — From Run Manager

Note: This option is not available in the HyperWorks Student Edition

Running RADIOSS: Option 2 — From HyperMesh

42
Essential Run Options filename in form: filename_nnnn.rad
-starter Run just the starter
-nt nnn Specify number of CPUs in multiprocessor (SMP) run

-v ver Required version of solver (e.g. 13.0)

-sp Specify to run a single precision version of RADIOSS


-noh3d Do not generate h3d file from anim files
-screen Display results in the solver control window

RADIOSS File Structure (14 minutes)

(Video By Rahul Rajan, Altair India)

43
44
5 Frequently Asked Questions
(RADIOSS Files)

5.1 How To Set The Number Of Processors For A


Run

For RADIOSS SPMD

The number of processors has to be defined in command line argument (see


Command Line Arguments) by option –nspmd (or –np). It is used by RADIOSS Starter
for the “Domain Decomposition”.
The same number of processors must be input as a MPI parameter when you launch
RADIOSS Engine.
Example:
[user@machine]$ ./s_11.0_linux964 -nspmd 4
[user@machine]$ mpirun -np 4 ./e_11.0_linux64 plmpi -nt 2 -input
CRA2V51_0001.rad

For RADIOSS SMP

Define the number of SPMD domains to 1, and then define the number of threads by
domain to a value > 1 (see Help Documentation - RADIOSS > Reference Guide for
details).
Check if the number of processor is correctly set by looking at the RADIOSS Engine
display (standard output or listing file (Runname_0001.out, Runname_0001.lis or
RunnameLnn):

45
The number of processors has to be defined in Command Line Arguments by option:
-nthread [INTEGER] or -nt [INTEGER] at starter.
Example:
[user@machine]$ ./s_11.0_linux964 -nt 2 -i CRA2V51_0000.rad
The number of processors must be the same when launch engine and is defined in
command line argument by option:
-nthread [INTEGER] or -nt [INTEGER] at starter.
[user@machine]$ mpirun -np 1 ./e_11.0_linux64 -nt 2 -input CRA2V51_0001.rad

What Should I Do If The Engine Stops With The Message: **


ERROR IN OPENING INPUT FILE?

** ERROR IN OPENING INPUT FILE


This message appears, if:
• RADIOSS Engine input file in V44 (RunnameD01) contains the keyword /VERS/51
• RADIOSS Engine input file in V51 (or higher) (Runname_run#.rad) contains the
keyword /VERS/44.

46
5.2 Parallelism
What Is The Difference Between RADIOSS SMP And RADIOSS SPMD?
SMP Version
SMP stands for Shared Memory Parallelism.
This version is based on shared memory architecture of computers.
For this kind of computer, all processors can access to a common memory space. Each
process can access to all memory allocated by the program.

Shared Memory Architecture


Depending on the operating system, this version will be able to use up to 8, 12 or 16
processors with a reasonable efficiency. The efficiency can be considered as an
acceptable one if it is greater than 50% (that is to say, the speedup is greater than half
the number of processors).

SPMD Version

SPMD stands for “Single Program Multiple Data”.


RADIOSS Starter divides the model into several domains (equal to number of
processors asked for).
Up to RADIOSS Version 4:
The data is next sent to each process in RADIOSS Engine; this is the “Parallel
Initialization”. The first process reads the complete RESTART file, sending the data to
each other process.
Some additional memory is needed on first process for this operation in this version.
From RADIOSS Version 5.1:

47
RADIOSS Starter writes as many RESTART files as the number of processes which will
be used for the run. Each RESTART file is built from only the data that will be needed
by the relative process, and is read directly by this process when running RADIOSS
Engine. Additional memory is no longer needed for the “Parallel Initialization”.
Then each program executes on each domain. The management of computation at
the frontiers of the domains remains and the necessity to communicate some
information between the processes is implemented through the use of the MPI
‘Message Passing Interface’.
This SPMD version is free from the computer architecture. It can run on distributed
memory machines, shared memory machines, workstations cluster or high-
performance computation cluster.
With this version it is possible to use more processors than with the SMP version, with
a better efficiency. It becomes possible to use more than 128 cores with a good
speedup.

Distributed Memory Architecture For Computer

48
Version 100: Radioss File Extensions

49
50
6 Element Formulations
Discretization is the process of dividing a model body into an equivalent system of
many smaller bodies (finite elements) interconnected at points (nodes or nodal
points) common to two elements. Suppose we need to find the stresses and
displacement of an object upon loading. The object is irregular in shape. So in order
to model this, we need elements that have geometric shape similar to the real
structure or region. One geometric shape cannot represent all possible engineering
structural shapes. The user will model with elements that look like beam, truss, plate,
cylinder…etc. In FEA, by using basic element almost all structures are approximated.

6.1 1D Elements (Line elements)


An element created by connecting two nodes, as shown below, at its ends represents
a line element or 1D element. The cross-sectional area is assumed constant
throughout the element. Truss, beam and spring elements are examples of 1D
elements.
In the following we are going to have a look at Element Formulations available in
RADIOSS. This overview is based on the RADIOSS Standard Training Class Material and
the HyperWorks Help Documentation.
The different element formulations (1D, 2D, 3D) require different parameters to be
set. Hence, we strongly recommend to review the element definitions below before
you run your first model.

51
Typical Applications Of 1D Elements

Automotive Examples
1D Elements can be used to define a variety of connections between 2 nodes
• Beams
• Springs
• Pulleys
• Dashpots

52
6.2 Overview Of Element Library – 1D Elements

/PROP/TYPE2 or /PROP/TRUSS

Truss elements are a subset of beam elements which have no bending degrees of
freedom and they are modelled with 2 nodes. Truss elements can be oriented freely.
Truss elements transmit force axially only. That is, these elements allow translation
movement only and not rotation. Truss elements are normally used in conditions
where no bending is expected and are normally used to model towers, bridges…etc.
A constant cross section area is assumed and they are used for linear elastic structural
analysis. These elements are compatible with elastic and elastoplastic materials
(Material Law 1 and 2)

53
/PROP/TYPE3 or /PROP/BEAM

Beam elements are long and can be oriented anywhere in 3D space. Beam elements
have 6 degrees of freedom thereby allowing both translation and rotation at each end
node. This is the primary difference between beam and truss elements. The 1-2 node
define element geometry, the 3 node defines the cross sectional orientation. A
constant cross section area is assumed for beam elements.

54
/PROP/TYPE4 or /PROP/SPRING

Spring elements are also 1D elements which are modelled using 2 nodes. These
elements are commonly used to model connectors in the model.
• Simple physical spring with dashpot
• 1 DOF per node, tension-compression behavior:

• Linear elastic:

• Nonlinear elasto-plastic:

Linear and variety of non-linear springs can be defined for SPRING through function
reference (similar for SPR_PUL, SPR_BEAM and SPR_AXI)

55
Example Spring Property Card

/PROP/TYPE13 or /PROP/SPR_BEAM

• 2 or 3 nodes element.
• 6 degree of freedom per node according to the property
• For each degree of freedom (3 translation, 3 rotations), the behavior definition is
similar to the previous spring element.
• There is coupling between bending and shearing.
• Uncoupled failure parameters (min, max) for each direction or coupled failure model
• Failure defined function of elongation(strain),force or internal energy.
• Spring element can be activated of de-activated by sensor.

/PROP/TYPE45 or /PROP/KJOINT2

• Spring elements between node 1 and node 2.


• Node 1 and Node 2 nodes must be slave node of different rigid bodies.
• It is recommended to have initially node 1 and node 2 coincident.
• Local skew system is defined :
• Local X : between node 1 and node 3
• Local Y : between node 1 and node 4

56
Joints Types

x: denotes a blocked d.o.f. / 0: denotes a free d.o.f

Joints Types

x: denotes a blocked d.o.f. / 0: denotes a free d.o.f

The joint stiffness is automatically computed by default.

57
• Null for free degree of freedom
• Rigid body stiffness for the blocked degree of freedom Possibility to define specific
stiffness and viscous behavior for each degree of freedom

Spring Property Summary:

In the following videos, Rahul Rajan explains in some detail 1D Elements of type:
Beam, Spring Type 4, Spring Type 13, and Spring Type 8

1D Elements: BEAM & TRUSS (7 minutes)

58
1D Elements: Spring Type 4 (8 minutes)

1D Elements: Spring Type 12 (4 minutes)

1D Elements: Spring Type 13 (22 minutes)

59
1D Elements: Spring Type 8 (20 minutes)

Information Required For Spring (6 minutes)

Procedure to Create KJOINT2 in RADIOSS (9 minutes)

60
Some related questions:
Can we?
Is it possible to provide failure parameters for spring elements?
• Yes, it is possible to provide deformation, force and energy-based failure criteria for
spring elements.

How to?
How to model a revolute joint in RADIOSS?
• The easiest way to model joints in RADIOSS is to use a /PROP/TYPE8 spring. Provide
high stiffness to all degrees of freedom those are fixed and zero stiffness in the
direction which you want to be free and thus a revolute joint can be modelled.

6.3 2D Elements (Shell Elements)


Shell elements are 3 or 4 node 2D planar elements with constant thickness, and have
either a triangular or quadrilateral shape, that can be oriented in the space. They are
typically used to model structures such as pressure vessels, automobile bodies, ship
hulls, aircraft fuselages… etc. Shell elements support all translational degrees of
freedom as well as all rotational degrees of freedom, that is shell elements have 6
degrees of freedom. RADIOSS shells are based on Mindlin-Reissner theory that
includes transverse shear deformation valid for moderately thick and thin shells.
RADIOSS provides both reduced and fully integrated 4-noded shell elements.

61
Typical Applications Of Shell Elements

Overview Of Element Library – Shell Elements

• Geometry defined by a curved mid-surface defined with 3 or 4 nodes and a thickness


h—usually small compared to the two other dimensions
• Shells include membrane, bending and shear deformation while the normal stress
(through the thickness) is always zero
• RADIOSS shells are based on Mindlin-Reissner shell theory that includes transverse
shear deformation valid for moderately thick and thin shells
• Mass of the element is distributed to the nodes
• Large rotations with large strains with bilinear shape functions

62
RADIOSS provides both reduced and fully integrated 4-noded shell elements
A 3-noded shell element is also available

Shell Elements: Reduced Integration Schemes


• 1 integration point on the surface
• Widely used in programs with explicit time integration applied in crashworthiness
studies
• Dramatically decreases the computation time compared to 4-pt integration
• Very competitive if the hourglass modes (which result from the reduced integration
scheme) are “well” stabilized.

63
Hourglassing Due To Reduced Integration
• Reduced integration can lead to non-physical zero energy modes, called hourglass
modes
• Two methods used to control formation of hourglassing:
• Perturbation (Q4 Elements)
• Physical Stabilization (QEPH Elements)

How to?
How to control hourglass energy in shell elements?
• Normally hourglass occur in under integrated elements. The recommended parameters
for avoiding hourglass, for shell elements, are Ishell=24 (QEPH formulation) and keep N=5
(number of integration points).

Hourglass Perturbation Control Method

Perturbation method Ishell=1,2,3 or 4 (Q4)

• The hourglass deformation is detected with the relative motion (velocity) of the
nodes.
• When hourglassing is detected, a force is applied on the node to stabilize the
deformation.
• This force is defined according the element stiffness (Young modulus, yield stress,
plastic tangent..)
• This force introduces an artificial energy (named hourglass energy)

64
Perturbation method Ishell= 24 (QEPH)

• Hourglass yield depends on natural yield


• Hourglass loading stiffness depends on natural tangent modulus
• The hourglass treatment is done at the material level

Hourglass Control Methods Comparison

• The Q4 (Belytshko & Tsay) elements may have bad behavior in case of large rotation
and distortion and in case of elastic load-unloading cycles.
• The QEPH elements give accurate results in case of loading and unloading even if there
is hourglass mode deformation.

In this video Rahul Rajan talks about the meaning and effects of Hourglassing (7
minutes)

65
Key Parameters on/PROP/SHELL

Ishell Select element formulation for 4-noded shells


Ish3n Select element formulation for 3-noded shells

N Choose number of integration points through the thickness

Ismstr Specify small or large strain formulation

Iplas Choose method for calculating post-yield strains

Ithick Determine if thickness is constant under load, or if element thins

/PROP/TYPE1 or /PROP/SHELL: Q4

• Q4 is an under-integrated 4 node element (BT)


• 1 integration point on the surface
• Hourglass stabilization with perturbation
method
• Relatively inexpensive
• Default shell formulation is Ishell = 1

66
PROP/TYPE1 or /PROP/SHELL: QBAT

• QBAT is a fully integrated 4 node element (Batoz)


• 4 Gauss integration points on the surface
• No hourglass deformation
• 5 local degrees-of-freedoms per node—drilling (rotation around local z) can be
activated with Idrill
• No shear locking due to in-plane reduced integration for shear
• Transverse shear deformation is taken into account
• Most accurate element, but higher CPU cost (2 to 3 times that of a Q4)

/PROP/TYPE1 or /PROP/SHELL: QEPH

• QEPH is an improved under-integrated 4 node element with 1 integration point on


the surface
• Physical hourglass stabilization
• Best compromise between cost and quality
• Generally, QEPH elements costs only up to 15% more than a BT element and the
results are close to those of QBAT

/PROP/TYPE1 or /PROP/SHELL: C0

• C0 is a standard triangle (C0) with modification for large rotation


• 1 integration point on the surface
• No hourglass deformation
• Low CPU Cost

67
Shell Element Formulation: Summary

/PROP/TYPE1 or /PROP/SHELL

Integration Points:
• Up to 10 points through the thickness
• 1 integration point through the thickness give membrane behavior
• N = 5 integration points recommended

/PROP/TYPE1 or /PROP/SHELL

Small Strain
• Usually for problems with small deformations
• Strain computed based on the initial element shape
Large Strain
• For problems with large displacements and strains
• Strain computed based on the current element shape with the derivate of the shape
function

68
• Small strain activation possible in RADIOSS Engine in the event that time step reaches
D Tmin due to element deformation
Shell Stress Computation
From the deformation and the strain rate of the previous step, RADIOSS computes the
strain to each integration point
From this strain, the stress is computed with a elastic stress increment

With these computed stress values, the Von Mises stress is computed and compared
to the yield stress of the material

If the Von Mises stress is smaller than yield stress, the computed stress is the one used
in order to compute the internal force of the elements

/PROP/TYPE1 or /PROP/SHELL

If the Von Mises stress is bigger than yield stress, RADIOSS will :
• Update the strain through the thickness (ezz)
• Compute the new thickness of the shell element
• Compute the stresses using the material law plasticity surface, the strain and the
strain rate. The projection of the stress on the plastic surface is done with:
• Radial return (faster)
• Iterative projection with Newton’s method (more accurate)

69
/PROP/TYPE1 or /PROP/SHELL

• Initial thickness is used to compute strains and to integrate stresses, but the thickness
variation is still computed for post-processing reasons
• If a variable thickness is used, true thickness is computed not only for post-processing,
but also for strain computation and stresses integration.
• For accurate results, especially when necking or spring back, it is strongly
recommended to use iterative plastic projection and thickness variation.

Shell Elements: Recommendations For Crash Applications

70
Shell Elements: Comments On Shell Modeling Practices
• Mesh on neutral fiber
• Create homogeneous mesh size for each part
• Have minimum of 6 elements per buckle wave length
• Use 3 elements on flanges
• Suppress holes smaller than reference mesh size

71
Shell Elements: Shell Element Quality

In this video Rahul Rajan discusses different 2D elements (18 minutes)

72
6.4 3D Elements (Solid Elements)
The 3D solid elements are either a tetrahedral (4 faces) or hexahedral (6 faces)
element. Brick or tetrahedral elements may have different number of nodes and
support only translational DOF. They are normally used to model solid objects where
shell elements are not appropriate to model them

Typical Applications of 3D elements

• Geometry defined by nodal location


• Hexahedron and tetrahedron with linear and quadratic interpolation functions
• Linear elements are better in terms of time and memory consumption, especially due
to the low number of integration points and a larger time step

73
• Penta (6-noded) elements used in mesh transitions are degenerated hexahedrons
(handled by pre-processor)

Overview Of Element Library– Solid Elements

74
/PROP/TYPE14 or /PROP/SOLID
• Isolid Choose hex element formulation
• Itetra Choose tetra element formulation
• Ismstr Specify small or large strain formulation
• Iframe Select co-rotational formulation if required
(hex only)
• IHKT Choose tangent modulus estimation method

/PROP/TYPE14 or /PROP/SOLID: Brick8

• Brick8 is a reduced integration, 8 node element with linear interpolation


• 1 point integration
• Viscous (penalty) hourglass stabilization

• Characteristic length,
• Relatively inexpensive

75
/PROP/TYPE14 or /PROP/SOLID: HA8

• HA8 is a fully integrated, 8 node element


• Number of integration points, Inpts, can be specified, from 2x2x2 to 9x9x9
• No hourglassing
• Relatively expensive

/PROP/TYPE14 or /PROP/SOLID: H8C

• H8C is a fully integrated, 8 node element


• Number of integration points 2x2x2
• No hourglassing, Relatively expensive

/PROP/TYPE14 or /PROP/SOLID: HEPH

• HEPH is an improved reduced integration, 8 node element


• 1 point integration
• Physical hourglass stabilization
• Best compromise between cost and quality

/PROP/TYPE14 or /PROP/SOLID: Tetra4

• Tetra4 is a 4-noded solid


• 1 point integration

76
• No hourglassing
• Shear locking

• Characteristic length,
• Relatively inexpensive

/PROP/TYPE14 or /PROP/SOLID: Tetra4

• Tetra4 is a 4-noded solid with 6 DOF per node & 4 point integration
• Displacement of the dummy nodes is computed on the basis of rotational DOF
• High time step versus 10 nodes tetra element with same accuracy
• Shear locking effect low or negligible (it may appear in bending)
• Compatibility with shells
• Characteristic length,

77
Solid Element: Summary

/PROP/TYPE14 or /PROP/SOLID: Strain Formulation

/PROP/TYPE14 or /PROP/SOLID: Element Coordinate System

Co-rotational formulation is more accurate if large rotations are involved, at the


expense of higher computation cost
It is recommended in cases of elastic or visco-elastic problems with important shear
deformations or in the case of large rigid body motions

78
/PROP/TYPE14 or /PROP/SOLID: Hourglass Tangent Modulus

• Default is fine for most cases


• Recommend to use IHKT = 2 with HEPH and Lagrange type total strain, Ismstr = 10 for
foam or rubber materials (laws 38, 42, 69, and 82)

Solid Elements: Recommendations for Crash Applications

79
For Foam or Rubber Materials (Law 38, 42, 62, and 82):

Note: For Ismstr = 10 engineering stress-strain curves are required

For Foam Material (Law 70):

Note: For Ismstr = 1 engineering stress-strain curves are required

80
Note: The Isolid flag is ignored for /TETRA4 and /TETRA10 elements as these elements
have their own specific formulation

Overview of Element Library – 3D Thick Shell

Watch this video by Rahul Rajan to learn more about 3D element formulations (15
minutes)

81
Question
While running RADIOSS I’m getting this error as below
MESSAGE ID : 169
** STOP : ZERO OR NEGATIVE VOLUME
FOR 4 NODES TETRAHEDRON ID=121736, INTEGRATION POINT=1
Is there any options/ method to remove this error?
 The NEGATIVE VOLUME error happens when solid elements are very deformed
and their characteristic length goes to 0. For large strain formulation the time step
of an element goes to 0 when the element is compressed. In RADIOSS Starter
input file (Runname_0000.rad), use Ismstr =2 in the solid property and in
RADIOSS Engine file (Runname_0001.rad) use the option/DT/BRICK/CST which will
set the time step value tmin at which the solid elements will switch to small strain.
This means that the solid elements using Ismstr =2 will use large strain formulation
while their time step remains greater than tmin, and will then switch to small strain
formulation.
(/DT/BRICK/CST option is used to control the time step of the analysis. The time
step becomes constant after reaching Tmin specified by the user).

6.5 Frequently Asked Questions (Elements)


The Run Stopped With The Message: “Zero Or Negative Volume”, When Solid
Elements Are Used With Ismstr=2 and /DT/BRICK/CST; Is This Normal?
Ismstr =1, 2 and 3 are not available for the 8 integration points solids using
formulations Isolid =12 and 112.
This means that these solids continue to use large strain formulation, and hence the
following error message appears: “Zero or Negative Volume”
In order to use this small strain formulation with 8 integration points solid elements,
use the HA8 solid formulation which is available as of RADIOSS V44.
When using this formulation, set Isolid =14 with Inpts =222 (corresponding to Isolid
=222 in input format 44). Also set Icpre =1 for elastic or visco-elastic material law, and
Icpre =2 for elasto-plastic laws.

82
I Used Solid Elements And Several Integration Points, And Starter The Following
Error Message Appears, While The Element Seems To Be Well-Defined: ** ERROR:
ZERO OR NEGATIVE 3D SOLID VOLUME, is This Normal?
** ERROR: ZERO OR NEGATIVE 3D SOLID VOLUME

ZERO OR NEGATIVE VOLUME 3D-ELEMENT NB 1


In the figure above, the volume of the element is positive, but the sub-volume
associated to one integration point is negative. The solid is decomposed into sub-
volumes associated to each integration point. If the element is badly warped, one sub-
volume could be negative

How Many Integration Points Should Be Used In The Thickness Of Shell Elements?
If only one integration point is used, a membrane only behavior will be obtained
(except with law 1, up to V44). Some materials, such as fabric, can justify such a choice
(no bending strength).
In case of an elastic behavior, one gets the exact solution from three integration points
– that is to say that the bending moments are exactly integrated through the thickness
of the shell – and it is not necessary to use more integration points.
In case of a plastic behavior, the bending moments are not integrated exactly. Using
more integration points, the solution becomes more accurate; so it is recommended
to use five integration points.

83
I Am Using Material Law 1 With 1 Integration Point Through the Thickness And Still
Get Bending Moments; Why?
Up to V44, Material Law 1 for shells does not use integration points but switches to
the global formulation (corresponding to N=0); whatever the number of integration
points N is asked for. So up to V44, there is no way to use this material law with only
one integration point and membrane only behavior.
Workaround: use Material Law 2 with integration point and a huge value for the yield
stress.
As of V51, Material Law 1 for shells uses global formulation, except if one integration
point is asked for in the property, then a membrane only behavior occurs. Note that
this can explain some differences in the results between versions up to V44 and V51.

With Shell Elements Using The Same Material Law, But Different Types Of
Properties (While Keeping The Same Number Of Integration Points), I Do Not Get
The Same Results; Why?
The integration scheme which is used for property types 1 and 9 (relative to isotropic
shells through the thickness) sets the integration points and weights in order to
integrate exactly the bending moments in the elastic case (from three integration
points since for one integration point, no bending moments are computed).
The integration scheme which is used for property types 10 and 11 is a step-by-step
integration scheme and uses integration points at the center of each layer, and
weights which correspond to the relative thickness of each layer. So the integration
scheme is not the same one.
An error relatively important can occur in the elastic field, when there are a few layers
or large differences on the thicknesses of the layers. One way to work around this
problem is to subdivide the thicker layers. But it is generally not well-suited in case of
modeling the failure of the layers.

84
Elastic Case:
Stress is linear through the thickness; an integration of forces step-by-step is exact.
But the integration of moments, step-by-step is not exact since symbol(z) z is
quadratic.

When Is It Better To Use QEPH Shells Instead Of Belytschko Shells?


QEPH shells are more accurate for elastic or elasto-plastic loads, whatever the loading
type - quasi-static or dynamic; but they are not recommended with anisotropic and
orthotropic material laws.
QEPH shells will give better results if the mesh is fine enough. In case of a coarse mesh,
this formulation will be too stiff and some local buckling phenomena could be missed.
In case of a coarse mesh, the Belytschko shells often give better results.
QEPH/HEPH is not recommended for orthotropic materials because the stabilization
forces are computed based on isotropic assumptions

I Used Solid Elements And The Run Stopped Before The End Time, With The
Message: “Zero Or Negative Volume”: How Can This Problem Be Solved?
This happens when solid elements are very deformed and their characteristic length
goes to 0. You may notice in the output file that before getting this error message,
the time step of the element written into the message drops down.
In case of large strain formulation, the time step of an element goes to 0 when the
element is compressed. In a mathematical way, the element cannot reverse its
orientation since its stiffness increases to an infinite value; but due to numerical
accuracy, the element may go to reverse its orientation.

85
In order to solve the problem of both the drop in cycle time step and subsequent
termination of the run due to a negative volume, you might first check that the
material used is well-suited to the physics which is represented. Then switch the
elements to small strain formulation. This is done as follows:
In the RADIOSS Starter input file (Runname_0000.rad), use Ismstr =2 in the solid
property or in the option /DEF_SOLID; in RADIOSS Engine file (Runname_0001.rad)
use the option / DT/BRICK/CST which will set the time step value symbol_tritmin at
which the solid elements will switch to small strain.
This means that the solid elements using Ismstr =2 will use large strain formulation
while their time step remains greater than symbol_tritmin, and will then switch to
small strain formulation.
Their volume will then remain constant and the element can even reverse its
orientation. The drop of their time step normally stops except for some materials,
especially viscous materials.

86
6.6 Frequently Asked Questions (Property)
Question
I have a problem using PROP16 for fabric shell. THICK is defined in global thickness as
0.93, then three layers are set each 0.31 mm thickness.
What is the meaning of the following error?

Answer
This error means the number of layers in your model is higher than 1. As mentioned,
you defined three layers which is not allowed for /PROP/TYPE16. The following input
parameters are available in /PROP/TYPE16:
Thick: Shell thickness
N: Number of layers, with 1 < N < 100, Default set to 1 (Integer)
Ti: Thickness of layer i (currently only one is allowed)
Since this property is currently compatible with Elastic Anisotropic Fabric
(/MAT/LAW58 - FABR_A), only one layer is allowed. It is recommended:
Thick = 1*T1
/PROP/TYPE16 + /MAT/LAW58 will consider only membrane stress, no bending stress,
since fabric cannot resist bending.

87
Question
Which strain formulations are available for connection element with
/PROP/CONNECT?
/PROP/CONNECT is used only for /MAT/LAW59. This element formulation is based on
relative elongation of element in normal and shear direction. This allows avoiding the
dependence of element time step from element height. Element height can be even
equal to 0. Note that material stiffness is specified in normal and in shear direction.
Material has no stiffness in lateral direction. Therefore, it is always advised to couple
nodes of LAW59 elements to some Lagrange components either directly or through a
tied contact.
The parameter Ismstr in /PROP/CONNECT specifies how the nodal forces are
calculated from stresses.
For Ismstr=4 use the actual middle area of the element and for Ismstr=1 use initial
middle area. This is implemented in order to avoid time step drop, when LAW59
element separates from master surface (due to the failure of master surface
elements) and it starts to expand laterally.

88
7 Example: S - Beam Crash
(This example is taken from the online help documentation)

The purpose of this example is to study the behavior of a crashed S-beam using various
shell formulations and a number of different integration points. This test also
compares the initial velocity influence on results.
A sensitive study is performed on a crushed S-beam. The modeling includes a material
law using the elasto-plastic model of Johnson-Cook and a self-impacting interface
based on the Penalty method in order to model the buckling of the beam. An initial
velocity is applied on the left section via a kinematic condition: either a rigid body or
a rigid link. The impacting condition is sliding and is secured by specific boundary
conditions in the right section. Half of the structure is modeled.
The results are compared according to three different views:
• Shell element formulations (BATOZ, QEPH and BT hourglass type 3)
• Plasticity options (global and progressive plasticity)
• Influence of the initial velocity (5 and 10 ms-1) Several criteria are used to compare
the results:
• Deformation configuration
• Crushing force versus displacement (via momentum integration)
• Energy assessment
• Displacement of the left section
• Hourglass energy
• Kinetic energy
• Internal energy
• Maximum force
• Maximum plastic strain

89
BATOZ and QEPH element formulations provide accurate results. The BT hourglass
type 3 formulation is a low-cost method and the QEPH formulation provides a good
precision/ cost ratio (the cost is three times lower than the BATOZ formulation).
BATOZ and QEPH are element formulations which do not have hourglass energy.
The results show an over-estimation of the plastic strain in the case of the global
plasticity use. The input files used in this example are part of the HyperWorks
installation:
• QEPH:<install_directory>/demos/hwsolvers/radioss/03_S-
Beam/QEPH/Global_plasticity/QEPH*
• BATOZ:<install_directory>/demos/hwsolvers/radioss/03_S-
Beam/BATOZ/Global_plasticity/BATOZ*
• BT_type3_NiP0:<install_directory>/demos/hwsolvers/radioss/03_S-Beam/BT-type3/
Global_plasticity/Q4_NIP0*
• BT_type3_NiP5:<install_directory>/demos/hwsolvers/radioss/03_S-Beam/BT-type3/
NiP5/Q4_NIP5*

Physical Problem Description


An S–beam is crushed at an initial rate of 5 ms-1 against a rigid wall. The section is an
empty square-shaped tube (each side measuring 80 mm). The thickness is 1.5 mm.
The tube is made of steel, and plasticity is taken into account, but not failure. Using
symmetry, half of the crosssection is modeled

Problem Description And Beam Cross-Section.


The following system is used: mm, ms, g, N, MPa
The material used follows an isotropic elasto-plastic Johnson-Cook law.

90
Material properties:
Young’s modulus: 199355 MPa
Poisson’s ratio: 0.3
Density: 7.9x10-3 g/mm3
Yield stress: 185.4 MPa
Hardening parameter: 540 MPa
Hardening exponent: 0.32
Maximum stress: 336.6 MPa
All other properties are set to the default values.

Analysis, Assumptions and Modeling Description

Modeling Methodology

The mesh is a regular shell mesh. Each shell measures approximately 10 mm x 10 mm.
A sensitive study is performed on:
• Shell element formulations: BATOZ, QEPH and Belytschko hourglass type 3
• Plasticity options: global and progressive plasticity model
• Influence of the initial velocity: 5 and 10 ms-1

The rigid wall is modeled with boundary conditions on the right section of the beam
(X, Z translations and all rotations fixed).
The left section undergoes the following conditions:
• Fixed in the Z direction.
• Initial velocity of 5 m/s in the X direction.
• All nodes are rigidly connected in X, Y and Z directions.
• A 500 Kg mass is added on the left end.

91
Block format input specifications:
• Hierarchy organization: there is only one subset made up of three parts, one for each
side of the beam, and one for the top. The materials and properties are identical for
each part.
• Node groups: there are three node groups, one for each end of the beam, and one for
the symmetry plane. The boundary conditions are set on the left end.
• TH selection: DX is saved for node 1 (the node used to display displacement at the left
end).

RADIOSS Options Used

Taking account of symmetry, half of the structure is modeled. The symmetry plane
covers the y axis = 0 mm. Boundary conditions are also set at the right end to simulate
a rigid wall (slide).
Two equivalent possibilities are available for generating kinematic conditions
attached to the left extremity of the beam. The first consists of creating a rigid body
to connect all of the left section nodes to the gravity center of the beam cross-section,
with a mass being introduced on a master node. The second type of modeling retained
uses the rigid link option, which rigidly connects the left section nodes in the X, Y and
Z directions. A 500 kg mass is added to the master node.

Both models provide identical results; the rigid link will be used for this example.
• An initial velocity of 5 ms-1 is used for the master node of the rigid link or for the rigid
body.
• MODIF file:

92
• A MODIF file enables to add option(s) during a run. The MODIF files carry the name
ROOTNAMErun*.rad. Where, run# is the RADIOSS run number four digits from 0000
to 9999 and run# is the name of the last Restart file + 1.
• For example, to run a MODIF file after the first run (restart file ROOTNAME_0001.
rad), the run number for the MODIF file must be 2: ROOTNAME_0002.rad. MODIF files
use the same input format as the RADIOSS deck. Put all the input decks in one folder
and with Irun=2 RADIOSS will automatically recognize the MODIF file

After 20 ms, a self-impacting interface is required to deal with the buckling of the
beam. This is added using a MODIF file where the interface is defined and saved for
the TH. This type of interface corresponds to 7; all values are set to “default”. To
define the master side, a surface is defined using three parts of the model
(/SURF/PART). The safest and easiest method for defining the slave side of a self-

93
impacting interface consists of defining a node group with the master surface
(/GRNOD/SURF).
The MODIF file is ROOTNAME_0002.rad.
The next Engine file is ROOTNAME_0003.rad (final time = 30 ms).

Contact Force At The Start Of Self-Impacting.


The MODIF file options used in Engine file ROOTNAME_0002.rad are:

94
Simulation Results and Conclusions

Deformed Mesh For Belytschko


Hourglass Type 3 Formulation
(V=5 m.s-1)
The crushing force is obtained by time
derivation of the X-momentum. The
maximum displacement over a 20 ms long
computation corresponds to 96.4 mm.

Crushing Force (X-Direction) Versus Displacement For Different Element


Formulations (V=5 M.s-1)

95
Energy Assessment For Belytschko Hourglass Type 3 (V=5 M.s-1).

The structure does not absorb a lot of energy and that you should check the hourglass
energy, which may be relatively high compared with the total energy.
The following table shows the results obtained using different element formulations
and plasticity options:

Initial Velocity = 5 Ms-1 (Values In Brackets Are The Energy Percentages Compared
With The Initial Energy)

96
97
Initial Velocity Influence
The following table indicates the influence of the crushing velocity (5 ms-1 and 10 ms-
1
).

(Values In Brackets Refer To The Energy Percentages Compared With The Initial
Energy) BT hourglass type 3 formulation is used in this section.
The amount of internal energy stored in the beam during a crash is relatively higher
when the initial velocity is set to 10 ms-1, instead of 5 ms-1. The hourglass energy is
quite low with either initial velocity

Crushing Force Versus Displacement For The Different Initial Velocities

98
Kinetic Energy Normalized For The Different Initial Velocities

First self-contact:
Initial velocity = 5 ms-1: displacement = 120 mm;
Initial velocity = 10 ms-1: displacement = 94.15 mm.

99
100
8 Materials
Any good simulation of models depends on the inputs provided and an experienced
analyst will always evaluate and confirm the correctness of the inputs to get the best
results of the analysis. Good attention should be provided to the material type chosen
and the material properties provided. There could be different material properties
needed as input to fully characterize a material for analysis purposes. Any unclear or
illogical data provided during material characterization will lead to uncertain results.
Choosing the right material and providing the reliable material data in performing a
simulation is very important. RADIOSS provides one of the most comprehensive
material and rupture libraries in the industry. The material laws and rupture criteria
span across definitions for concrete, foam, rubber, steel, composites, biomaterials,
and more.
The overview of the materials in RADIOSS library can be classified as shown below

In this chapter, which again is based on the standard RADIOSS training course and the
HyperWorks help documentation, we will have a look at the following material laws:
• Law 1 – Elastic
• Law 2 – Johnson-Cook
• Law 27 – Elastic-Plastic Brittle
• Law 36 – Elastic-Plastic Tabulated
• Law 42 – Ogden (Visco Hyperelastic)
• Law 70 – Foam

101
The materials may be categorized as
• Isotropic Elasticity
• Isotropic Elasto-Plastic
• Composite and Anisotropic
• Viscous
• Hydrodynamic
• Explosives

8.1 Isotropic Elasticity

This keyword /MAT/ELAST defines an isotropic, linear elastic material using Hooke’s
law. This law represents a linear relationship between stress and strain. It is available
for truss, beam (type 3 only), shell and solid elements. Useful for elements attached
to rigid body.
This material law is used to model purely elastic materials. The material stiffness is
determined by only two values: the Young’s modulus (E), and Poisson’s ratio (n). The
shear modulus (G) can be computed using E and n, as shown below:

Note: The value of density is always used in explicit simulations.

102
8.2 Isotropic Elasto-Plastic

103
/MAT/LAW2 (PLAS_JOHNS) Johnson-Cook Material

This law represents an isotropic elasto-plastic material using the Johnson-Cook


material model. This model expresses material stress as a function of strain, strain
rate and temperature.
A built-in failure criterion based on the maximum plastic strain is available.

How To Determine Parameters

If material is not temperature or strain rate sensitive, then:


Calculate true stress and strain from engineering values from tensile test:
σy ≅ Sy σu = Su(1+ eu) εu = ln(1+ eu) Determine
Johnson-Cook parameters:

(Note: Mind the conversion of true stress vs. engineering stress. A brief “discussion”
about this topic is included in the free Study Guide “Practical Aspects of Finite Element
Simulation”)

104
/MAT/LAW2 or /MAT/JOHNS: Example Mild Steel
Mild steel material model with rupture defined on material card. Failure of the
element happens when 1 integration point reaches εpmax

/MAT/LAW2 or /MAT/JOHNS: Optional Failure Model


Using the Johnson-Cook failure model /FAIL/JOHNSON along with the Johnson-Cook
material model provides more user control of the failure Johnson-Cook failure strain
definition:

105
Include /FAIL/JOHNSON to define element failure criteria

Mild steel applied to tensile and compression specimen using failure on


/FAIL/JOHNSON

106
How Can I?
I only have yield strength, ultimate tensile strength, poisson’s ratio, and density data,
can I get value for a, b & n for Johnson-cook material?
For /MAT/LAW2, if user sets Iflag = 1, user will need yield stress ( SIG_y) , Ultimate
tensile engineering stress (UTS) and engineering strain ( EPS_UTS ) at necking point.
With this new input, RADIOSS automatically calculates the equivalent value for a, b
and n.

/MAT/LAW36 (PLAS_TAB) Elastic Plastic Piecewise Linear


Material

This law models an isotropic elasto-plastic material using user-defined functions for
the workhardening portion of the stress-strain curve (for example, plastic strain vs.
stress) for different strain rates.

107
Note: The first point of yield stress functions (plastic strain vs stress) should have a
plastic strain value of zero. If the last point of the first (static) function equals 0 in
stress, default value of epmax is set to the corresponding value of ep. (for more
information please see the Hyper-
Works Help Documentation)
• Isotropic elasto-plastic material
• User-defined function for the work-hardening portion of the stress-strain curve
• Available for brick and shell elements
• Elastic portion of material stress-strain curve defined by Young modulus and Poisson’s
ratio
• Material plasticity curves can be given for up an arbitrary number of strain rates
• Linear interpolation of strain-stress curve
• For a given strain rate
• For a given plastic strain

In order to avoid poor curve extrapolation


during the strain calculation, a dummy curve
defined for a high strain rate value can be added
in the model—typically the highest strain rate
curve is repeated.

108
The curves have to be defined such as they don’t intersect themselves in the strain
range being used. Otherwise the calculation may fail.

There are 4 damage and rupture parameters:


• εpmax is the maximum plastic strain for element deletion for any loading (tensile , shear
or compression).
• εt is the beginning of tensile damage. From this point the stress value defined in the
curve is reduced by the factor
where ε1 is the largest principal strain
• εm is the end of the tensile damage. The stress value is null but the element is not
deleted
• εf is the tensile failure strain value for element deletion
Note: The first point of yield stress functions (plastic strain vs stress) should have a
plastic strain value of zero.

109
Example /MAT/PLAS_TAB card image:

/MAT/PLAS_TAB or /MAT/LAW36 : Yield Stress Function

110
/MAT/PLAS_TAB or /MAT/LAW36 : Example Mild Steel
Results show failure in tension because of tensile failure criteria on material card

Example /MAT/PLAS_TAB card image with Johnson Cook failure model

Results using failure criteria on/FAIL/JOHNSON failure card

111
/MAT/PLAS_TAB or /MAT/LAW36 : Yield Stress vs Pressure
If the yield stress is pressure dependent, then a scale factor can be applied via a
function, fct_Idp
Function fct_Idp is used to distinguish the behavior in tension and compression for
certain materials (that is: pressure dependent yield)

Example of pressure dependent yield stress on crush tube

112
/MAT/PLAS_TAB or /MAT/LAW36 : Yield Stress Vs Pressure Example of pressure
dependent yield stress on crush tube—card image

Result of pressure dependent yield stress on crush tube:

113
Isotropic Elasto-Plastic (cont’d)

/MAT/LAW27 (PLAS_BRIT) Elasto-Plastic Material with Brittle


Failure

This law combines an isotropic elasto-plastic Johnson-Cook material model with an


orthotropic brittle failure model. Material damage is accounted for prior to failure.
Failure and damage occur only in tension. This law is applicable only for shells.
Yield surface definition same as Johnson-Cook (Law2)
Damage and rupture defined with 4 parameters for each principal direction:
et = Strain at the beginning of tensile failure
em = Maximum tensile strain at which the stress in the element is set to a value
dependent on dmax1
dmax = Maximum damage factor
ef = Maximum tensile strain for element deletion

114
Useful for modeling brittle failure of glass

Review the summary overview about Material Law2 and Law36 by Rahul Rajan (28
minutes)

/MAT/LAW27 or /MAT/PLAS_BRIT: Example Safety Glass

115
8.3 Composite and Anisotropic

116
8.4 Viscous

/MAT/OGDEN or /MAT/LAW42 : Ogden, Mooney-Rivlin

Visco Hyperelastic behavior (non-linear elastic). This law is generally used to model
incompressible rubbers, polymers, foams, and elastomers. This material can be used
with shell and solid elements
Generally used to model polymers and elastomers
Stress – strain relationship is based on the strain energy density, W:
The stress corresponding to
this strain energy is:
The Mooney Rivlin version
gives the strain energy as:

where:
Note: C10 and C01 are material parameters; I1 and I2 are first and second deviatoric
strain invariants:
Example Card Image of Hyper-Elastic Material

117
Best practice for corresponding property :
• HEPH solid element formulation (Isolid = 24)
• Lagrange type total strain (Ismstr = 10)
• Advanced tangent modulus estimation (IHKT = 2)

/MAT/LAW69

Incompressible Hyperelastic Material model specified using the Ogden, Mooney-


Rivlin material models. It is generally used to model incompressible rubbers,
polymers, foams, and elastomers. Material parameters are computed from
engineering stress-strain curve from uniaxial tension and compression tests. It is used
with shell and solid elements. Users can reproduce viscous behavior in this material
card using /VISC/PRONY card, which is an isotropic visco-elastic Maxwell model that
can be used to add visco-elasticity to certain shell and solid element material models.

118
The visco-elasticity is input using a Prony series. (the Card Image of Law69 is shown
below)

How to?
How to evaluate Mooney Rivlin constants from stress strain curve?
• You can use /MAT/LAW69 - Incompressible Hyperelastic Material with Tabulated
Input specified using the Mooney-Rivlin material models and this material model
is generally used to model rubbers. Material parameters are computed from
engineering stress-strain curve provided by the user.

/MAT/FOAM_TAB or /MAT/LAW70: Tabulated Foam

• Visco-elastic behavior
• Stress – strain relationship defined by functions according to strain rate
• Unloading behavior defined by:
• Curves
• Hysteresis parameters
Best practice for corresponding property:
• Belytschko solid element formulation (Isolid = 1, or 17 if hourglass issues)
• Small strain (Ismstr = 11)
• Element coordinate system non co-rotational (Iframe = 1)

119
8.5 Hydrodynamic

120
8.6 Failure Models
Failure models are independent and can be coupled with compatible material

8.7 Explosives

121
8.8 Generating Material Data
• Ask material supplier
• Do a web search (www.matweb.com)
• Use testing facilities to generate material parameters:
https://www.datapointlabs.com/ Catalog.asp

For the more advanced RADIOSS user:


You may also create & define your own materials in RADIOSS. For further information
please see

RADIOSS USER’S CODE INTERFACE


Also worthwhile to watch is the “lecture” by Paul du Bois about material models for
crashworthiness analysis”

122
8.9 Frequently Asked Questions
(Materials/Failure)
How To Get Law24 (Concrete Material) Data From Real Physical Tests?
You can get data for Law24 from the following tests:
• Get density by weighing a specimen
• Get Young modulus with a cylinder compression test
• Poisson coefficient (usually assumed to be near 0.2)
• Get ‘fc’ from Compression strength (using cube compression test)
With the following additional tests, material Law24 will be more accurate (in fact,
providing data fits the failure envelope):
• Get ‘ft’ from tensile test:
‘ft’ is direct tensile strength. This test provides ‘ft/fc’ value for Law24.
• Get couple of (ft,Ht) from Splitting tensile test:
This test is also called Brazilian test.
‘Fst’ is splitting tensile strength in Splitting tensile test. Assuming that Fst=ft. Then you
need to model this test to fit limit strength (ft) by validating ‘Ht’ value. This test should
be with very slow velocity, so use HA8 solid property for in Splitting tensile test.
• Get ‘fb’ from biaxial test and get (f2,s0) from confined test.

Which Curve Should Be Input In Material Law36 In Order To Characterize A


Material?
For elastic plastic laws in RADIOSS, one gives a true stress versus true strain curve, so
you have to convert the experimental curve. How to proceed is described in the online
documentation (Example 11: Tensile Test in the RADIOSS Tutorials and Examples
manual):
I Use Solid Elements And Epmax In The Corresponding Material After The Criteria
Epmax Is Reached, Elements Are Not Deleted; Why?
Depending on the material law, the solid elements are not deleted after the criteria
E1pmax is reached (they do not appear as “deleted elements” in post-processors).

123
For Material Laws2, 4 and 22, only the deviatoric part of the stress tensor is set to
zero, the internal pressure of the solid is still computed.
On the other hand, for Material Laws3, 23, 28 and 36 the solid elements are deleted
when E1pmax is reached.

124
9 Example: Law Characterization -
Tensile Test
(This example is taken from the Help documentation)

Aim of The Problem


It is not always easy to characterize a material law for transient analysis using the
experimental results of a tensile test. The purpose of this example is to introduce a
method for characterizing the most commonly used RADIOSS material laws for
modeling elasto-plastic material. The use of “engineering” or “true” stress-strain
curves is pointed out. Damage and failure models are also introduced to better fit the
experimental response.
Apart from the experimental results, the modeling of the strain rate effect on stress
will be considered at the end of this example using a sensitivity study on a set of
parameters for Johnson-Cook’s model.

Physical Problem Description


Traction is applied to an object. A quarter of the object is modeled using symmetrical
conditions. The material to be characterized is 6063 T7 Aluminum. A velocity is
imposed at the left-end.
Units: mm, ms, g, N, MPa.

125
The material undergoes isotropic elasto-plastic behavior which can be reproduced by
a JohnsonCook model with or without damage (/MAT/LAW27 and /MAT/LAW2,
respectively). The tabulated material law (/MAT/LAW36) is also studied.

Input File
• Law_2_Johnson_Cook:<install_directory>/demos/hwsolvers/radioss/11_Tensile_tes
t/ Law_2_Johnson-Cook/.../TENSIL2*
• Law_27_Damage:<install_directory>/demos/hwsolvers/radioss/11_Tensile_test/
Law_27_Damage/DAMAGE*
• Law_36_Tabulated:<install_directory>/demos/hwsolvers/radioss/11_Tensile_test/
Law_36_Tabulated/TENSI36*

Experimental Results Of The Tensile Test: Engineering Stress Vs. Engineering Strain.

126
Analysis, Assumptions and Modeling Description
Modeling Methodology
The average element size is about 2 mm in the mesh. There are 2014-node shells and
one 3-node shell.
The shell properties are:
• 5 integration points (progressive plastification).
• Belytschko elasto-plastic hourglass formulation (Ishell = 3).
• Iterative plasticity for plane stress (Newton-Raphson method; Iplas = 1).
• Thickness changes are taken into account in stress computation (Ithick = 1). •
Initial thickness is uniform, equal to 1.7 mm

Mesh Of The Object


• Node number 54 was renamed “Node 1” to be compliant with the Time History.
For node 54, only displacements in the x-direction (variable DX) are saved.

Sections Saved For Time History


For both sections, the variables FN and FTX, are saved; thus the following variables
will be available in /TH/SECTIO: FNX, FNY, FNZ (saved using “FN”), and FTX.
Engineering strains will be obtained by dividing the displacement of node 1 with the
distance up to the symmetry axis (75 mm). Engineering stresses will be obtained by
dividing the force through section 1 with its initial surface (10.5 mm2). Therefore, the
results shown correspond to the engineering stress as a function of the engineering
strain, equivalent to the force variation compared to displacement (similar curve
shape).

127
RADIOSS Options Used
An imposed velocity of -1.0 m/s in the x-direction is applied to the nodes, shown below
(abscissa less than or equal to 25 mm). The displacement is proportional to time.

Imposed Velocities

Variation Of Node 1 X-Displacement In Relation To Time


Only one quarter of the object is modeled to limit the model size and to eliminate the
rigid body motions. Symmetry planes are defined along axis x = 100 mm and axis y =
0. Two boundary conditions cannot be applied to the same node 13 (corner).

Boundary Conditions
The lower side is fixed in Y and Z translations and X, Y, and Z rotations.
The right side is fixed in X and Z translations and X, Y, and Z rotations; the node in the
corner is completely fixed.

128
Characterization Of The Material Law
There are two steps to characterize the material law:
• Transform the engineering stress versus engineering strain curve into a true stress
versus true strain curve (this step applies to any material law).
• Extract the main parameters from the true stress versus true strain curve, to define
the material law (Johnson-Cook law and material coefficients for /MAT/LAW2 or the
yield curve definition for /MAT/LAW36).

True Stress/True Strain Curve

Engineering strains are computed using the following relationship:

And true strains are computed with the relationship:

Both strains, therefore, are linked together by:

Engineering stresses are measured by dividing the force through one section with the
initial section. True stresses are measured by dividing the force with the true
deformed section:

Thus, to compute true stresses, the surface variation must be taken into account.
Assuming that Poisson’s coefficient is 0.5 during plastic deformation, the true surface
in mono-axial traction is:

Thus, the relationship between true and engineering stresses is:

129
Characterization Of The Material Law
The characterization will be made for /MAT/LAW2 (Johnson-Cook elasto-plastic),
/MAT/ LAW27 (elasto-plastic with damage model) and /MAT/LAW36 (tabulated
elasto-plastic). For each of the material laws, the yield stress and Young’s modulus are
determined from the curve.
The plastic strain can be defined as:

An important point to be characterized on the curve is the necking point, where the
slope of the force versus the displacement curve is equal to 0, and where the following
relationships apply:

Guidelines For Necking Point

Equations Used For Analysis

130
Simulation Results and Conclusions
Experimental Results
An experiment designed by the “Norwegian Institute of Technology” as part of an EC-
financed program, “Calibration of Impact Rigs for Dynamic Crash Testing” is used. The
following curve was obtained from the experiment:

Engineering Stress Versus Engineering Strain Curve (Experimental Data)


It is estimated that the necking point occurs between 6% and 8% (engineering strain).
After analysing the experimental data, the first point satisfying the necking condition
is at 6.68%.

Comparison Between Engineering And True Curves (From Experimental Data)


Engineering formulation is converted into true formulation using the relationship:

The true stress curve is higher than the engineering stress curve, as it takes into
account the decrease in the objects cross-section.

131
Law 2: Elasto-Plastic Material Law Using The Johnson-Cook
Model Johnson-Cook Material Coefficients

The stress versus plastic strain law is: (Johnson-Cook model)


where, a is the yield stress and is read from the experimental curve and then
converted into true stress.
To compute b and n, two states are needed. This leads to the following formulas for
b and n:

The first point is chosen at the necking point, then b and n are computed for each
other point of the curve and averaged out since the results tend to differ depending
on the point chosen.

Characterization Up To The Necking Point


The first stage when determining the material model is to obtain Johnson-Cook’s
coefficients. Neither the maximum stress, nor the failure plastic strain effects are
taken into account here (set at zero).
The values of coefficients are chosen so that the model adapts to the test data (see
following image)

Variation Of The Engineering Stress/Strain According To Johnson-Cook’s Model

132
Adapted To The Test The material coefficients used for Law2 are:
• Initial density: 2.7x10-3 g/mm3
• Yield stress: 90.27 MPa
• Poisson’s ratio: 0.33
• Hardening parameter: 223.14 MPa
• Young’s modulus: 60400 MPa
• Hardening exponent: 0.375
In the figure below the yield curve defined using the Johnson-Cook model with the
one extracted from experimental data is compared.

The true stress – true strain relationship can be described by:

The engineering stress deviations between experiment and simulation are described
in the table below:

Comparison is performed up to the necking point (engineering strain = 6.68%) because


after this state, a rapid decrease in the engineering stresses occurs in the object. The
rupture sequence is simulated in the following paragraphs. Results using Law 2 remain
within 8% of the experimental curve.

133
The curve could be improved by slightly adjusting some of the values. The purpose of
this test is to propose a method for deducing material law parameters using a tensile
test.

Beginning Of The Necking Point


Necking Point Simulation
The Johnson-Cook model previously defined
corresponds to the experimental results up to the
necking point. However, the slope
of the numerical response does not
enable the necking point to start at
the strain value observed
experimentally.
The necking point is characterized
by the slope value of the true stress
versus the true strain curve, which
must be approximately equal to the
true stress. The necking point numerically appears by continuing simulation until the
condition on the slope is observed.
The results are obtained using the Johnson-Cook model 1: Beginning of the necking
point using only the first coefficients of the Johnson-Cook model (a, b and n).

134
True Stress Versus True Strain Curve Up To The Beginning Of The Necking Point
The necking point can be simulated, either by adjusting the Johnson-Cook coefficients
to obtain an accurate slope, or by compelling curve with a maximum stress.

Simulation Of The Slope Near The Necking Point


By implementing an energy approach, the hardening curve can be modified to achieve
an engineering curve which resembles a horizontal asymptote near the necking point
with the purpose of simulating the behavior of the curve as observed in the test.
The Johnson-Cook coefficients used to describe the physical slope are:
• Yield stress: 79 MPa
• Hardening parameter: 133 MPa
• Hardening exponent: 0.17
For this model, the new true stress/true strain relationship is:

(Johnson-Cook model 2)
The results obtained with those
coefficients are provided below. In the
figure below the Johnson-Cook model 3
with the experiment is compared:
Adjusted Engineering Stress/Strain Curve
To Model The Beginning Of The Necking Point.
The shape of the yield curve versus the experimental data is depicted next

Yield Curves

135
The necking point is defined as
This condition is characterized by the intersection of the true stress versus the true
strain curve with its derivate.

Superposition Of Engineering Curve And True Curve With Its Derivate.

Beginning Of The Necking Point Using A Maximum Stress Limit,


lsa_smax
For this test, the Johnson-Cook coefficients input are those set in characterization up
to the necking point, the failure effect not being taken into account (the failure plastic
strain is set to zero). The beginning of the necking point is set using the choice of a
maximum stress value. In comparison to the experimental results, the necking point
is well defined for a maximum stress set at 175 MPa. The limit in stress appears on
the von Mises stress versus true strain curve on elements where the necking point
occurs.
The maximum true stress manages the beginning of the necking, as shown below:

Engineering Stress Versus Engineering Strain; Necking Point Characterization

136
Variation Of The Von Mises Stress With The True Strain From Shell 11

Maximum stress smax is reached for von Mises stress on shells where the necking
begins. To avoid overly-high stresses after the necking point, a maximum stress factor
must be set approximately equal to the true necking point stress.
The following curves show the evolution of the von Mises stress versus the true strain
shell at two characteristic locations of the object:

Von Mises Stress Curve With A Maximum Stress Limit


The beginning of the necking point is observed following the point where the stress is
equal to stress versus strain derivate.

137
Yield Curve With Maximum Stress

The yield curve is described by:

The derivate of the stress is very sensitive and strongly depends on the yield curve
definition. Thus, introducing the necking point into the simulation is very delicate (a
small change can result in many variations). The necking point should first begin on a
given element for numerical reasons. The preferred beginning of necking is addressed
below.

Preferred Beginning Of The Necking Point


Experimentally, the beginning of the necking point can appear anywhere on the
object. The beginning of the necking point should preferably be located on the right
end elements in order to propose a methodology for this quasi-static test. If the model
only uses a quarter part of the object, the necking point is found on elements 30, 125
and 78.
The beginning of the necking point is physically and numerically sensitive and can be
initiated on the right elements by changing a few of the coordinates along the Y-axis
of the node in the right corner (node 16) in order to decrease the cross-section and
privilege the necking point in this zone. Changing the node position by 0.01 mm is
enough for achieving the preferential beginning of the necking point.

138
Node 16 To Be Moved

A second approach also enables the necking point to be triggered on the right end
side by defining an extra part, including shells 3, 11 and 4 by using a maximum stress
slightly lower than the remaining part, in order to initiate the necking point locally
since the necking point stress is first reached in the elements having the lowest
maximum stress value, that is shells 3, 11 and 4. This method, based on material
properties, is quite appropriate for demonstrating the characterization of a material
law and will thus be used in the continuation of the example

Localization of the beginning of the necking point according to the models using smax.
The material is described as Johnson-Cook Model 1:

smax = 174 / 175 MPa


The following curves indicate the variation of the engineering stress versus the
engineering strain according to the beginning of the necking point zone and in
comparison, to the experiment.

139
Engineering Stress/Strain Curve For Each Starting Necking Point Location

There is a fast decrease in the engineering stress after the right-end necking point.
The necking point, due to the boundary conditions of the y-symmetry plane (y-
translation DOF released), becomes more pronounced.
The variations in the section where the necking point is found are quite similar up to
the necking point. After such point, there is a sharp surface decrease for the right-end
necking point, contrary to the second case where the surface decrease is more
moderate.

Variation Of Cross Section (Necking Point Zone).

Improvement Of The Elements’ Contribution During The Necking Point Sequence


In order to simulate physically the contribution of each element in the necking point,
it is advisable to adjust the curve by varying the Johnson-Cook coefficients in order to
increase the intensity of stress at the necking point. The main result is no longer the

140
variation of the stress/ strain curve but rather the surface under the curve which
characterizes the energy dissipated during the test. This energy-based approach is
relevant for crash tests since the final assessment is often more significant than how
it was achieved.

Engineering Stress/Strain Curve Obtained Using Adjusted Johnson-Cook Coefficients


The following graph compares the new yield curve with experimental data:

Yield Curves

141
Material as described in the Johnson-Cook coefficients are:
Johnson-Cook Model 3:

(true stress/strain)
• Yield stress = 50 MPa
• Hardening parameter = 350 MPa
• Hardening exponent = 0.38
• Maximum stress is set to 189 or 190 MPa (according to the parts)
The results of adjustment to the Johnson-Cook coefficients are depicted below:

Shell Contribution During The Necking Point Sequence (Von Mises Stress)
As the necking point progresses, more physical results are obtained due to the new
input data of the material law coefficients having a better element contribution.

Variation of The Von Mises Stress On Elements 110, 109, 108, 107, 11 And 106

142
Damage Modeling With Plastic Strain Failure
The elasto-plastic model of Johnson-Cook is used until failure, which is simulated using
a plastic strain failure option. The element is deleted if the plastic strain reaches a
user-defined value emax. This damage model shows good stability. A maximum plastic
strain is defined for each Johnson-Cook model:

emax = 75%; Yield Curve close to Experimental Data

emax = 47% ; Yield Curve Adjusted With Respect To Lower Stresses:

143
emax
Failure is reached for relatively high true strains.

Law 27: Elasto-Plastic Material Law With Model Damage


Law 27 is used to simulate material damage following a Johnson-Cook plasticity law.
Thus, model damage is associated with the previous law in order to take account of
failure.
The damage parameters are:
• Tensile rupture strain et1: damage starts if the highest principal strain reaches this
tension value.
• Maximum strain em1: the element is damaged if the highest principal strain is above
the tension value. The element is not deleted.
• Maximum damage factors symbolmax: this value should be kept at its default value
(0.999).
• Failure strain ef1: the element is deleted if the highest principal strain reaches the
tension value.

Stress/Strain Curve For Damage Affected Material.

144
The following graphs display the results obtained using the material coefficients of
two previous Johnson-Cook models. Damage parameters complete those models.
Damage Model A

145
Damage Model

146
Law36: Tabulated Elasto-Plastic Law
This is a tabulated law; therefore, the true stress versus plastic strain function can be
directly used. The rupture phase can be simulated by adding points to this hardening
function.

The hardening curve has to be defined with precision around the necking point while
the decrease of the curve is very sensitive to its adjustment. In order to improve the
modeling of the necking point, two points can be interpolated, one “just before” the
necking point, and one “just after” with the slope between those two points equal
to the necking point stress.

147
148
10 Interface / Contact Modeling
Normally the explicit analysis method solves highly complex nonlinear problems and
contact definition is an important part of many nonlinear problems. In order to get
the best results and real behavior of any simulations, including contacts, accurate
modeling of contact interfaces between bodies is a crucial part. Based on the
applications RADIOSS has a large number of contact types. Each of those interfaces
was developed for a specific application field, but this field is not the only choice
criteria. In RADIOSS an interface is created by defining the master segment and slave
set. The master segment can be parts, elements, contact surface, node sets and
similarly for slave the selection can be parts and node sets. Once these master slave
sets are identified the algorithm will check these locations for potential penetration
of a slave node through a master segment, and this check is done in every time step

In this chapter we will learn about:


• Interface Penalty Method
• TYPE7, variable stiffness, node to surface contact
• TYPE11, variable stiffness, edge to edge contact

149
The Penalty Method Used In Interfaces
The penalty method treats the behavior between slave nodes and master segments
as springs that generate resistive forces as a function of penetration

Spring elements are created between the slave node and its projection. The spring
stiffness is automatically computed according to the material and property of both
master segment and slave node

Variations Of Interfaces

150
Types Of Interfaces Used In RADIOSS

Out of the many interface types listed above, the most used ones for impact
simulation are:

Interface Type 7: Nodes To Surface Contact

• Automatic node to surface contact interface using penalty method with variable
stiffness
• A gap is defined around the master segment and the slave nodes
• The contact gap can be constant or variable

151
• Resistive forces are applied to keep slave nodes out of the contact gap
• Interface time step is computed to insure stability
Interface type 7 is a general-purpose interface and can simulate all types of impact
between a set of nodes and a master
surface. Contrary to interface types 3 and 5,
interface type 7 is non-oriented and slave
nodes can belong to the master surface.
Therefore, this interface can simulate auto-
impact, especially buckling during a high-
speed crash. Interface type 7 solves all
problems and limitations encountered with interface types 3 and 5. The search for the
closest segment is done via a direct search algorithm; therefore, there are no search
limitations and all possible contacts
are found. The energy jumps
induced by a node impacting from
the shell edges are removed by the
use of a cylindrical gap around the
edges.
Finally, the main advantage of
interface type 7 is that the stiffness
is not constant and increases with
the penetration preventing the
node from going through the shell midsurface. This solves many bad contact
treatments (common when using either interface type 3 or 5).
The gap used in interface type 7 is rather different than those of the previous
interfaces. The gap used is on both sides of the shell mid-surface, and a cylindrical gap
is added around the edges.

152
Important Interface Parameters

Igap Determines how the size of the gap is calculated


Gapmin Minimum gap for activation of interface
Inacti Action to take if initial penetrations exist
Istf Affects how the stiffness of the interface is calculated
Iform Friction formulation
Stmin Minimum stiffness to use in the interface
Idel What to do with slave nodes and master segments if element fails (deleted)
that they are attached to

The recommended contact parameters for interface Type 7:

• Istf=4
• Igap=2
• Fscale_Gap=0.8
• INACTI=6
• Gap_min=0.1
• Fric = 0.1
• Iform=2

Interface Type 7: Some Gap Definitions

• gm: master element gap:


• gm = t/2, t is the thickness of the master element for shell elements
• gm = 0 for brick elements
• gs: slave node gap:
• gs = 0 if the slave node is not connected to any element or is only connected to brick
or spring elements
• gs = t/2, t is the largest thickness of the shell elements connected to the slave node

• gs for truss and beam elements, with S being the cross section of the element.

153
• If the slave node is connected to multiple shells and/or beams or trusses, the largest
computed slave gap is used.

A Type 7 Interface uses a gap even if the user has not defined it
If Gapmin is not defined or set to zero, then Gapmin = min (tm, lmin/2)
where: tm is average thickness of master shell elements
lmin is smallest side length of all master segments (shell
or brick) If all master segments are brick, then Gapmin =
lmin/10 where: lmin is the small side length of all master
brick segments.

If not explicitly defined, the Constant Gap would equal Gapmin which is the master
segment thickness

Next, we look at a Variable Gap:

The Variable Gap = max [Gapmin, (gs + gm)] where:


Gapmin is defined as before gm = tm/2, where tm is thickness of master shell – or zero
for master brick elements gs = ts/2, where ts is largest thickness of shell connected
to slave node – or zero if slave node not connected to any element or only
connected to brick or springs

154
Caution With Variable Gap Igap = 1, 2, or 3
• If Gapmin is not specified or set to zero, you can get unexpected results:
• Gapmin = tm , (if Gapmin is not specified and tm < lmin/2)
• gm = tm/2
• gs = ts/2
So, Variable Gap = max [Gapmin, (gs + gm)] = max [tm, (tm/2 + ts/2)]

Solution: Specify Small Gapmin Specify small Gapmin so that:


Gapmin = small number (less than smallest shell thickness, like 0.1) gm = tm/2 gs = ts/2
So, Variable Gap = max [Gapmin, (gs + gm)] = max [small number, (ts/2 + tm/2)]

Interface Type 7: What About Solid Elements?


Recall that for lgap = 1, the Variable Gap = max [Gapmin, (gs + gm)] where:
Gapmin = min (tm, lmin/2) or Gapmin = lmin/10 if master segments are all brick
gm = tm/2, where tm is thickness of master shell – or zero for master brick elements gs
= ts/2, where ts is largest thickness of shell connected to slave node – or zero if slave

155
node not connected to any element or only connected to brick or springs Best to
always specify a small Gapmin

Note: Values for variable gap shown when Gapmin is not specified

Interface Type 7: Options For Defining The Gap

The Variable Gap = max {Gapmin, min [Fscalegap*(gs + gm), Gapmax]} where Gapmin is
defined as before gm = tm/2, where tm is thickness of master shell – or zero for
master brick elements gs = ts/2, where ts is largest thickness of shell connected to
slave node – or zero if slave node not connected to any
element or only connected to brick or springs
Gapmax is an upper limit on the gap (no upper limit if not specified)

156
Interface Type 7: Limitation Of Gap Definition
In case of self-impact definition, there is an initial penetration when the gap value is
larger than the mesh size

This situation can be solved with Igap= 3 or with the flag Irem_gap = 2

Interface Type 7: Options For Defining The Gap

The Variable Gap = max {Gapmin, min [Fscalegap*(gs + gm), %mesh_size*(gs_l+gm_l),


Gapmax]} where
Gapmin, gs, and gm are defined as before gm_l = length of smallest
edge of master element gs_l = length of smallest edge of element connected
to slave node
%mesh_size = percentage of mesh size (0.4 is default)

157
Interface Type 7: What About Initial Intersections?
• Intersections result if slave nodes are past the master segment (shells or solids)
• This condition can cause locking situations leading to computation failure
• If slave nodes lie on the master segment an ERROR will occur because of the Type 7
contact stiffness definition
• Initial intersections should be removed from the model (using HyperMesh or
HyperCrash).

• Penetrations result if a slave node is closer to the master segment than the calculated
gap
• Initial repulsion forces give an error in the energy balance due to the interfaces work

• There are options in the solver to remove initial penetrations

158
Interface Type 7: Options For Removing Initial Penetrations

The initial penetration, P0 should be smaller than 30% of the initial gap to avoid
contact stability issues

159
Interface Type 7: Resistive Forces
The resistive force is nonlinear with respect to the penetration of the slave node into
the gap and is a function of the penetration value and rate.
Resistive Force
where:

The instantaneous stiffness is then:

Interface Type 7: Initial Contact Stiffness, K0

The final value for the contact stiffness, K, can be bounded by lower and upper limits,
Stmin and Stmax:

160
Interface Type 7: Sliding Force Computation

Interface Type 7: Time Step Computation


Because the interface stiffness is nonlinear, the slave node cannot cross the master
segment. This gives good accuracy for the node location.
However, this means that the interface time step needs to be evaluated during the
solution. If there is large penetration, the interface stiffness becomes large, leading to
reduced time step.

Interface Type 7: Update For Deleted Elements Due To Failure

161
Interface Type 7: Recommendations For Crash Applications

Note: To get detailed information in the *out file regarding any initial penetrations,
set / IOFLAG parameter Ipri = 5

Non-Linear Type7
Advantages
• Physically correct representation
• High accuracy for describing self-contact and buckling
• Best used with shells to shells or solids to shells
Disadvantages

162
• Requires more stringent modeling practices
• Does not tolerate intersections (slave node passing master segment)
• Time step may reduce due to increasing interface stiffness
Limitation Of Nodes-To-Surface Interface
There are situations where the Type 7 node to surface interface will not “see” the
contact due to edges crossing each other at angles

Question?
Type 7 contact isn’t working in my simulation and the slave nodes are penetrating into
the elements. i already tried changing the various contact parameters and it didn’t
work. Any ideas to why this is happening?
Had a look in your model file and. Since it’s an edge to edge contact try defining the
same with
Type 11 contact which is used to define Edge to Edge or Line to Line Interface

Interface Type 11: Edge To Edge Contact

In the case of edges crossing, Type 11 interface can be added to capture the physics
correctly
Type 11 is defined by identifying master and slave lines (edges)
• Type 11 works similar to Type 7 in terms of:
• Penalty formulation
• Gap definition
• Search method

163
Question?
How to create line_id(s) set for the TYPE11 contact?
• Please create the line segments by finding the free edges for the component and you
can create a set for the same from Analysis>Contact Surfaces and choose card image
as LineSeg.

A General Remark About Master & Slave:


The coarser the mesh, the stiffer the structure will be provided the material and
structural properties are same. Hence, the Master is typically meshed coarser so that
the Slave does not penetrate into master. That is why we choose less deformable or
high stiffer structure as master and less stiff as slave.

Interface Post-Processing
Animation file
• contact forces : /ANIM/CONT
• contact pressure: /ANIM/PCONT
Time history file
• Total contact energy in the global parameter

164
• Detail of the forces and energy for each contact, forces in the global frame, normal
force and tangential force (friction): /TH/INTER
Interface Surface Definition
• Variety of ways to define interface surface
• Best way, if possible, is to use /SURF/PART
• This way, if the part is modified, the surface definition remains intact

Overview On Contact Interfaces In RADIOSS

by Rahul Rajan (3 minutes)

How Sliding Contact Works


by Rahul Rajan (9 minutes)

165
Summary On Contact Type 7 by Rahul Rajan (23 minutes)

10.1 Frequently Asked Questions (Contact


Interfaces)
When Defining An Interface Type 7 With Solid Element Parts On The Master Side, Is
It Necessary To Put Shells On The External Surface Of The Part?
No, only the external surface of the solid should be defined.
As of RADIOSS V9.0, it is possible to define the surface of an external skin of a solid
part (or a subset, a material, a property, and so on) with the option /SURF/PART/EXT.
With older RADIOSS versions, this surface should be built using segments (a surface
type / SURF/SEG will be created but no additional element added).
When Defining An Interface Type 7, Which Is Neither Self-Impacting Nor Symmetric,
How Are The Master Side And The Slave Side Chosen?
It is recommended to choose the coarse mesh on the master side for better contact
detection.
In the following example, contact is not detected if the fine mesh is chosen on the
master side.

166
In case both meshes are of the same size, the stiffer structure on the master side
should generally be chosen.

What Is the Meaning Of: WARNING ID 94?


WARNING ID: 94
** WARNING IN INTERFACE GAP
INPUT GAP xxx
HOWEVER GAP IS RECOMMENDED TO BE LESS THAN xxx
This message is especially important in case of a self-impacting interface. The
recommended value corresponds to the smallest side length of shell elements on the
master side, divided by 2.
Shown below is a self-impacting interface for which the GAP equals half the side
length of an element.

If this element is compressed more than 50%, the red node enters into the gap of the
neighboring element. A self-contact is then detected (which may not happen). This
leads to over stiffening of the structure.

167
So the message means that at least one element on the master side has a side length
less than twice the GAP and there is a risk of over stiffening.
If the side length L of an element is lower than the GAP, a self contact will be computed
from the beginning of the computation.
Such a situation can be accepted if it is local enough, but not if it is a frequent situation
over the self-impacting interface. It is possible to find the source of this message by
using a preprocessor and selecting the elements through size criteria.
In case the interface is not self-impacting, possible consequence of this message will
be low performance, but the model behavior will be correct.

What Is the Meaning Of: “Infinite Loop Detected”?


This message can be written by RADIOSS Starter or by RADIOSS Engine. It corresponds
to a contact interface (Type 7, Type 11, etc. ).
This message means that RADIOSS has found too many possible impacts2 for a given
interface, with respect to the memory it has allocated for the storage of these
impacts.
Increasing parameter Multimp (default value is 4 for SMP and 12 for SPMD) in the
interface will workaround this problem.
Due to memory decomposition, it might be necessary to increase Multimp when using
a large number of processors in SPMD.

168
11 Kinematic Conditions
In this chapter the following kinematic conditions will be discussed:
• Tied Contact /INTER/TYPE2
• Rigid Walls /RWALL
• Rigid Bodies /RBODY
• Boundary Conditions /BCS
• Imposed Displacement /IMPDISP
• Moving, fixed coordinate system /SKEW
• Imposed Velocity /IMPVEL
• Incompatible Kinematic Conditions

Tied Contact /INTER/TYPE2


This interface is used mainly for the connections. The interface is defined with slave
nodes which are linked to master segment. The slave node need be projected inside
the master segment, that is the nodes of the slave surface are assumed to be
connected to the faces of the master surface and the distance of the slave node from
the master segment should be small enough.
Tied Contacts define an interface that kinematically connects a set of slave nodes to a
master surface. It can be used to connect coarse and fine meshes, model spotwelds,
rivets, etc.

169
Tied Contact Interface /Inter/Type2 (Between Solid Spotweld And Sheet Metal)
A tied contact is defined by a group of slave nodes and a master surface. The slave
nodes are kinematically tied to the master surface within a defined search tolerance,
dsearch
Warnings are displayed in the starter listing file (_0000.out) if a slave node cannot find
a master segment with the search tolerance.
The slave nodes cannot be set in any other kinematic constraint otherwise there are
incompatible kinematic conditions.

Optional Penalty Method


RADIOSS has an option to define Type 2 interface with the penalty method which
removes potential incompatible kinematic conditions with Spotflag = 25
A spring element is defined between the slave node and its projection. The penalty
stiffness is constant, calculated as the mean nodal stiffness of master and slave side

Tied Contact /INTER/TYPE2 Formulations

170
Tied Contact /INTER/TYPE2 Ignore Option

The flag Ignore is available to automatically remove all slave nodes which cannot be
projected on the master surface

In case failure (rupture of the shell element) is defined on the master side of the
interface, it is necessary to update the interface in order to release kinematic
condition between the slave node and the deleted element.

Tied Contact /INTER/TYPE2 Card Image

171
Tied Contact by Rahul Rajan (12 minutes)

What do you think?


I’m getting mass error in the engine out file? No mass scaling was applied, however,
the mass error in the engine out file was not zero and it stayed the same all the
process.
Do you have any Type 2 interface in the model? And if Spotflag =1 is used in the
Interface Type 2, then there are added mass on master nodes.
If option Spotflag =1 is used in the Interface Type 2, then there are added mass on
master nodes at the first cycle of RADIOSS Engine. Indeed the slave node inertia with
respect to the center of the closest master segment – let Is + msd2 where, Is is the
slave node inertia; ms is its mass and d is its distance to the center of the segment - is
equilibrated by added mass on master nodes.
When creating a spot weld with spring and Interface Type 2, it is important to check
that the springs inertia (I =m*L2) are not too large - that is to say, that they correspond
to the physics – and that the slave nodes height with respect to the master segments
is as small as possible.
Indeed, the geometrical concordance is necessary for Interface Type 2 to work
normally, and will help getting a lower added mass due to term msd2.

172
Rigid Walls /RWALL
Rigid Walls allows the user an easy way to define an interface between a rigid surface
and nodes of a deformable body. A rigid wall is a non-yielding retaining wall which is
defined by the user Master Node and a group of Slave Nodes.
Rigid Walls allow for an easy way to define an interface between a rigid surface and
nodes of a deformable body

Four Types Of Rigid Walls Are Available

• Infinite plane Infinite cylinder of diameter Ф

173
• Sphere of diameter Ф Parallelogram

A Rigid Wall is defined by a Surface Type and a group of Slave Nodes


The slave nodes to a rigid wall can be defined as a group of nodes and/or as nodes
initially at a distance less than the distance (Dsearch) from the rigid wall
Each wall can be fixed (point M) or moving (node N) A Rigid Wall is a Kinematic
Condition
• If a slave node penetrates the rigid wall surface a new velocity is computed
• From this velocity, the impulse and rigid wall force are determined
Example of fixed infinite planar rigid wall:

Example of spherical rigid wall with mass and initial velocity:

174
Moving Rigid Wall
A moving rigid wall is defined by a node number, N, and a point, M1. This allows a
normal to be calculated. The motion of node N can be specified with fixed velocity,
or with an initial velocity. For simplification, an initial velocity and a mass may be
given at the wall definition level.
A moving wall is a master slave option. Master node defines the wall position at
each time step and imposes velocity on impacted slave nodes. Impacted slave node
forces are applied to the master node. The slave node forces are computed with
momentum conservation. The mass of the slave nodes is not transmitted to the
master node, assuming a large rigid wall mass compared to the impacted slave node
mass.

Moving Rigid Wall Defined With Node N And Coordinates Of Point M1

Some Limitations with Rigid Walls


Can only be fixed or have initial velocity – may want to move rigid object with an
imposed displacement function.
Can result in Incompatible Kinematic Condition – a result of nodes belonging to
another kinematic condition and the rigid wall Solution:
Best practice is to use rigid body connected to a simple mesh representing the rigid
object:
• Easy to define more arbitrary shapes
• Can impose displacements, velocities
according to function
• Removes issue with Incompatible Kinematic
Conditions

175
In the video below Rahul Rajan talks about Rigid Walls (6 minutes)

Rigid Body /RBODY


A Rigid Body (/RBODY) constrains a set of nodes to move together with no relative
motion and can be used to simplify the modelling of non-deformable objects. A
kinematic condition is applied on each node and for all directions. And by default,
the master node is moved to the centre of mass. A rigid body can be compared to a
part with an infinite stiffness. No relative displacement is allowed between slave
nodes, and the general motion of the rigid body manages the master node. An
example is shown below where the impactor and supports are made rigid.

3 Point Bending Test. Purple - Rigid Body, Yellow = Rigid Body

176
Salient Points
• A Rigid Body is defined by a set of slave nodes and a master node
• Mass and inertia may be added to the initial master node location
• By default the master node is moved to the center of mass of all nodes
• Best practice is to create a master node that is not attached to anything
• A Rigid Body is a kinematic condition
• No relative motion is allowed between the nodes
• A kinematic condition is applied on each slave node for all directions

Rigid Body /RBODY – Typical Applications


Rigid Body covering a part of a finite element model
The center of gravity is computed by taking into account the slave node mass; the
master node is moved to the computed center of gravity, added mass and inertia are
placed at the center of gravity

ICoG = 2: The center of gravity is only computed by taking into account the slave node
mass; the master node is moved to the computed center of gravity, added mass and
inertia are placed at the center of gravity.

177
Rigid Body Representing A Non-Modeled Component
ICoG= 3: The center of gravity is set at the master node coordinates; added mass and
inertia are placed on the master node coordinates; slave node mass and inertia are
transmitted to the center of gravity. The master node is not moved.

ICoG = 4: The center of gravity is set at the master node coordinates; added mass
and inertia are put on center of gravity. The slave node mass and inertia are ignored.
The master node is not moved.

What Do You Think?

I have created a block with Hexa elements and shell elements with rigid body at 1 mm
distance from block. Block is going to impact over rigid shell with 50 N-mm of torque.
But block is not rotating with respect to the master node. What went wrong?

Keep the ICOG of the rigid body to 3 and run the analysis. By keeping ICoG=3, the
master node is not moved

178
Some more information about Rigid Bodies provides Rahul Rajan in his video below
(18 minutes).

Boundary Conditions /BCS


Boundary Conditions or constraints define nodal translational and rotational degrees
of freedom that are fixed from moving. These parameters limit the motion of the
object. The user should properly apply the loads and boundary conditions to get
meaningful results and over or under-constraining the model will provide wrong
results.

179
There are six degrees of freedom at a node that can be constrained:

A DOF is free, if the code is set to 0 (default) and fixed, if set to 1.


Example: 101 111 means the x and z translations, as well as all rotations are fixed; the
y translation is free.
If skew_id is non-zero, then the boundary conditions are applied with respect to this
local skew (local coordinate system))
A Boundary Condition is a kinematic condition in which the velocity on a node’s degree
of freedom is set to zero

180
Skew and Frames
Skews and frames are used to define local directions. These directions can be used to
apply:
• boundary conditions
• concentrated load
• fixed velocity

To define:
• rigid link orientation
• rigid body added inertia frame
• general spring reference frame
• beam type spring initial reference frame
• nodal time history output frame
Two types of skew frames are available in RADIOSS:

Fixed Skew Frame Defined By Moving Skew Frame Defined By 3


Nodes N1, N2, And N3
Coordinates Of O’, Y’, And Z’
X’ = Y’ X Z’ And Y’ Recomputed As Z’ X X’

181
Skew Reference
It is a projection reference to define the local quantities with respect to the global
reference. In fact, the origin of the skew remains at the initial position during the
motion even though a moving skew is defined. In this case, a simple projection
matrix is used to compute the kinematic quantities in the reference
Example Of Skew:
In this example imposed velocity is applied in Y direction. In /IMPVEL skew is used,
then the imposed velocity is computed in the Y-axis of global coordinate system and
then projected onto the Y’ axis of local skew reference.

Frame Reference
It is a mobile or fixed reference. The quantities are computed with respect to the
origin of the frame which may be in motion or not depending on the kind of
reference frame. For a moving reference frame, the position and the orientation of
the reference vary in time during the motion. The origin of the frame defined by a
node position is tied to the node. Frame measures relative motion whereas skew
measures global motion and projects it to skew. Only few options use frame like
imposed velocity, TH/NODE and others use skew.

Example Of Frame:
In this example rotational velocity is applied around Z axis. In /LOAD/CENTRI frame is
used, then rotational velocity is around Z’ axis of frame reference not in Z axis of
global coordinate system anymore.

182
What Is The Difference?
What’s the difference between the two-system creation method, that is using node
reference and by axis direction?
Using axis direction method, the user defines the system based on the axis direction.
But for node reference, based on the node selected by the user the local system is
created. Moreover, using axis direction method user can select a number of nodes
and the local system will be created in every nodes selected whereas this is not
possible for node reference method. And the most important is for systems created
by node reference method the local system created is also updated when the node
is translated whereas for the other method it doesn’t.
A summary about Added Mass and Skew Frame is provided in the following video by
Rahul Rajan

183
Imposed Displacement & Velocity /IMPDISP, /IMPVEL
Imposed Displacements define nodal translational and/or rotational degrees of
freedom that are defined to move according to a function

A metal strip is passed through two rollers aiming at reducing its thickness. Both
rollers have an imposed angular velocity.

• The imposed condition will follow the displacement (or velocity) versus time curve
defined by fct_IDT
• The direction is defined by Dir and can be X, Y, or Z in translation, and XX, YY, or ZZ for
rotation
• /IMPDISP and /IMPVEL are kinematic condition

Can We?
Is it possible to assign the BCs to a cylindrical system?
Select the nodes and put imposed displacement (/IMPDISP, and keep icoor=1),for
each dof with zero magnitude function.

184
Guidelines to Apply Quasi-Static Load In Explicit Simulation

Why We Should Slowly Ramp the Load For Quasi Static Simulation In Explicit Solver

185
11.1 Frequently Asked Questions (Kinematic
Conditions)
What Is The Meaning Of: WARNING ID 147, Which Appears Often In The Listing File
Runname_0000out?
WARNING ID: 147
*** WARNING: INCOMPATIBLE KINEMATIC CONDITIONS
2 KINEMATIC CONDITIONS ON NODE

Conditions upon node velocity


Two kinematic conditions which apply at the same time along non-orthogonal
directions are considered incompatible.
For instance, a boundary constraint on a node in X direction of the global system and
an imposed velocity in the same direction are not compatible (such a behavior is not
possible). On the other hand, a boundary constraint on a node in X direction of the
global system and an imposed velocity in Y or Z direction of the global system are
compatible.
Since a rigid body controls the movement of its slave nodes, a rigid body slave node
cannot use another condition (specifically it cannot belong to another rigid body).
If flag Ikrem of the rigid body gets its default value =0, the slave nodes of the rigid
body are automatically suppressed of the rigid walls; the master node of the rigid
body will realize the contact with the rigid walls if it is declared as a slave node of the
rigid walls.
In case some conditions are incompatible and are applied at the same time on a node,
the model behavior is not guaranteed, since those conditions will not be correctly
taken into account. The /PARITH/ON option is not supported in this case, so it is
necessary to check the source of those WARNINGS and to adapt the model.

186
In case of two incompatible conditions that are never applied at the same time (for
instance, because their activation intervals [Tstart, Tstop] do not overlap); the
WARNINGS need not be taken into account.
Moreover, no error is written by RADIOSS Starter better than these WARNINGS, which
would prevent from running the computation, since the analysis of the compatibility
of the kinematic conditions is sometimes more complex than RADIOSS Starter can
manage. In such cases, the WARNINGS provide the only clue to error termination of
the run and it is the user’s responsibility to check the model.

Why Is It Not Recommended To Set A Node Of The Mesh As The Master Node Of A
Rigid Body?
It is not recommended to set a node of the mesh as the master node of a rigid body,
since:
• The local system of the master node is not compatible with options imposing rotations
(except boundary conditions and imposed velocities). If some elements having a
rotational stiffness are connected (such as beams, some springs and shells), an
incompatibility occurs.
• Depending on the flag ICoG, which is used in the rigid body definition, a lot of mass
and inertia can be added locally onto the node; as well as the node can be moved to
the center of mass of the rigid body.

What Is The Meaning Of: **ERROR: RIGID BODY 345886 UNSTABLE, Which Makes
RADIOSS Engine End
**ERROR: RIGID BODY 345886 UNSTABLE

This means that:


where, vr is the rotational velocity of the rigid body and dt is the time step.
The rigid body turns more than one radian per cycle. This situation is due to
divergence in the computation. The model and its behavior must be checked.

187
What Do WARNING Id: 345 And WARNING Id 86 Mean?
WARNING ID: 345
*** WARNING/CHECK TIED INTERFACE
93836629 1218 601311 601309 601310 601312 1.0218 -0.0095 1.339E-01
This message is written for Interface Type 2 using the search formulation for the
closest master segment Isearch =2 (the default value for this flag as of RADIOSS V4.1).

This message is not due to the height of the slave node with respect to the master
segment.
In case the Interface Type 2 with Spotflag =0, this can cause an error in RADIOSS Engine
at the first cycle, due to a negative mass or inertia of a master node of the segment.
This can also cause an error in RADIOSS Engine during the run due to a negative
stiffness on a master node of the segment:
** ERROR : NEGATIVE STIFFNESS NODE …
In case of the Interface Type 2 with Spotflag =1, the situation corresponding to this
message can cause a relatively important added mass at the first cycle on the master
nodes of the segment.
In case of the following message:
WARNING ID: 86
** WARNING: TIED INTERFACE TIED INTERFACE:
4533 5490 248 5430 5250 5490 5491 1.00734 1.14868
This message corresponds to search formulation for the closest master segment
Isearch =1.

188
The meaning and possible consequences of this message are the same as for
WARNING ID: 345.

I Get A Very Bad Energy Balance (The Energy Error Is Large And Positive: There Is
Energy Creation) The Message: WARNING Id 547, Was Written To The Listing File
Runname_0000out; Why?
WARNING ID: 547
** WARNING IN RIGID BODY DEFINITION
NODE ID= 167467 IS MASTER OF RIGID BODY NUMBER 25
AND SLAVE OF ANOTHER ONE
The standard formulation for rigid bodies is not working if a master node of a rigid
body appears as a slave node of another rigid body; except for rigid bodies using
Lagrange multipliers
(/RBODY/LAGMUL).
In the figure below, three rigid bodies A, B and C are defined. The master nodes of A
and B have an added mass so that the centers of mass of the rigid bodies A and B are
not centered.
For representing the rigid body motion of the whole system A+B, a rigid body C has
been defined using all slave nodes of A and B and also including the master nodes of
A and B in order to take into account the added mass on these master nodes.
Then the message is written by RADIOSS Starter and the rigid bodies A, B and C
absolutely can not be used at the same time (the results would be arbitrary).
On the other hand, we can stand from the same Runname_0000.rad to study the
whole system C=A+B rigid body motion by deactivating rigid bodies A and B
(/RBODY/OFF) and activating rigid body C (/RBODY/ON), whether the independent
behaviors of A and B by activating rigid bodies A and B (/RBODY/ON) and deactivating
rigid body C (/RBODY/OFF).
The activation/deactivation of rigid bodies can be realized by using options
/RBODY/ON and /RBODY/OFF in RADIOSS Engine or by defining SENSORS.

189
Incompatible Kinematic Conditions
A kinematic condition results from using a feature in RADIOSS that specifies a
kinematic constraint on nodal velocities
These constraints are mutually exclusive for each degree of freedom, and there can
be only one constraint per DOF
RADIOSS Starter does not check if the kinematic conditions are really incompatible. If
they are strictly orthogonal, or if they are not applied simultaneously, just ignore the
warning.

In the first case, a node is slave to a moving rigid wall and also has a fixed boundary
condition. If the wall is moving, it is impossible after impact to respect both
conditions.
In the second case, a node is defined as a slave for two parallel walls that are fixed,
there are no possible incompatible conditions, as the node cannot impact the two
walls at the same time.

190
An example of the above is shown below. All DOF’s of the sphere are
constrained/fixed. At the same time an imposed velocity is applied. As a result,
either of the Kinematic Condition might be ignored depending on which gets
processed first…

191
192
12 Load Definition
In this chapter you will be introduced to:
• Concentrated Load /CLOAD
• Gravity /GRAV
• Initial Velocity /INIVEL
• Initial Velocity about an Axis /INIVEL/AXIS
• Pressure Load /PLOAD

Concentrated Load /CLOAD


A concentrated load defines a
concentrated force applied to
each node of a prescribed nodal
group. The load is applied to the
nodes indirectly by means of a
nodal group /GRNOD/xxxxx
The load versus time is defined by
a function /FUNCT

The grnd_ID input is obligatory. The concentrated loads will only be applied to the
nodal group. The force is applied to each node individually, so:
total force = F(t) * number of nodes

193
Nodal Groups /GRNOD/xxxx

Most of the load definitions apply their function through the use of nodal groups.
Nodal groups are very flexible and can
reference nodes directly or via parts
(components), by box, material,
property, etc.
1. A node group is a set of nodes. It can be
defined explicitly by:
• a list of nodes
• a list of nodes with generation
• a list of submodels (all nodes belonging
to the listed submodel are included)
• a list of submodels (all nodes defined by /NODE or /CNODE in the listed submodel are
included)
• a list of subsets or parts (all nodes belonging to the listed subset/part are included)
• a list of property sets or materials (all nodes belonging to an element having those
Property Type’s/MID’s are included)
• a list of groups of nodes
• a list of element groups (all nodes connected to those elements are included)
• a box (all nodes within a defined box are included)
• an unsorted node list: same as list of nodes, but the nodes are not sorted (only needed
in interface type 8 and /XELEM option).
2. Node groups are used to define slave nodes of rigid walls, interfaces, rigid bodies or
nodes to which a load is applied, like a concentrated load or a fixed velocity.
3. Nodes are stored in this order. For all other input types, nodes are sorted.
Example of applying initial velocity to all nodes of three parts:

194
How To?

How to apply torque in RADIOSS?

Torque can be applied by selecting the concentrated load from the BCs Manager and
selecting the direction plane (XX, YY or ZZ).

Gravity /GRAV
Gravity defines an acceleration on a group of nodes. The gravity versus time is defined
by a function /FUNCT
/FUNCT - Defines a function (For example: stress (Y-axis) as a function of strain (X-
axis)) (see image below)
Note: The gravity loads are only
applied to the nodes defined in
grnd_ID, which must not be null.

195
Initial Velocity /INIVEL
Initial Velocity defines an initial velocity on a group of nodes. No function is required
for this option, just velocity in the x, y, and z-directions
Specifying the group of nodes by part is very convenient when making changes to the
model.

What is?
What is the difference between initial velocity and imposed velocity?
• An imposed velocity is a kinematic condition, whereas the initial velocity definition is
an initial condition assigned to the model.
Initial Velocity about an Axis /INIVEL/AXIS
Initial Velocity about an Axis initializes translational and/or rotational velocities on a
group of nodes in a given coordinate system

196
rotation is not around a global axis

Note: The unit of velocity in translation is [m/s] In rotation, the angle is in unit [rad]
and the velocity is in unit [1/s]

Procedure to Apply Initial Velocity Around An Axis

Pressure Load /PLOAD


A pressure load defines a distributed force applied to each element face of a
prescribed surface group.
The pressure is applied to the elements indirectly by means of a surface group
/SURF/xxxxx
The pressure versus time is defined by a function /FUNCT

197
Note: In 3D analysis, positive pressure acts in direction N = N1 N3, x N2 N4 with N1,
N2, N3, and N4 being the nodes of the segment in the surface definition.
In 2D analysis,
positive pressure acts in
direction N normal to N1 and
N2, obtained by a rotation in
the counterclockwise
direction.

Learn more about Loads & Boundary Conditions in the summary video by Rahul Rajan
(18 minutes)

198
12.1 Frequently Asked Questions (Initial Stresses)
Is It Possible to Begin A Run With An Initial Stress State?
It is possible to define an initial stress state by using option /INISTA. This option uses
an OUTPUT file (Runname_0000.sty) describing the structure of the model for which
a file Runname_nnnn.sty defines a stress state.
The files Runname_0000.sty and Runname_nnnn.sty can be built in a first
computation (see ASCII Output File (STY-File)). They can also be written using an
external program, with respect to their respective formats.
As of RADIOSS V4.4, it is also possible to enter the initial stress values element by
element in the Runname_0000.rad file, by using options /INISHE and /INISH3 for shells
and 3-node shells and /INIBRI for solid elements.
As of RADIOSS V9.0, it is also possible to write with RADIOSS Engine state files
Runname_ nnnn.sta for shells.
Such a file will contain the geometry of selected parts (nodes coordinates and element
connectivity), plus cards /INISHE and /INISH3 defining the stress, strain, and so on, of
the shells and 3-node shells belonging to these parts.
If the file is included in another input deck Runname_0000.rad, the state (geometry,
stress, strain, and so on) of these parts will thus be retrieved (refer to /STATE in
RADIOSS Engine Input).

199
200
13 Output Requests
This is the content of this chapter:
• Summary of Post-Processing Files
• Animation Output Requests & Best Practices
• Time History Output Requests & Best Practices

Summary Of Post-Processing Files

Animation Output (A01, A02, … Ann) Time History Output (T01)


Requested in Engine Defined in Starter
For viewing animations & contours For plotting X-Y data
Animation Output Requests
/ANIM/Keyword2/Keyword3/Keyword4
Request defined in Engine
Output is written to runnameA01, runnameA02…runnameAnn file at a frequency of
Tfreq
Define Tfreq with /ANIM/DT in Engine
One Ann file for each output step
A runname.h3d file is also generated at the end of a successful run

201
/ANIM/Keyword2/Keyword3/Keyword4

Creating Animation Output Request /ANIM

In HyperMesh

HyperMesh > View Solver Browser > Right click > Create > Engine Keywords > ANIM>
option

202
In HyperCrash

HyperCrash > Model > Control Card


Animation Output Best Practice
Around 20 to 50 animations are necessary to analyze correctly the model behavior:
e.g. Tfreq = Final time / 20.
Recommended requests in the animation files include:

Additional outputs can be set to get more information for shell elements

203
Additional outputs can be set to get more information for solid elements

* Information for each integration point in also available

Time History Output Requests


/TH/Keyword
Request defined in Starter
Data is written to runnameT01 file at a frequency of ΔThis
ΔThis defined with /TFILE in Engine
Needs usually around 1000 points for the full run to get enough information.
ΔThis = Final time / 1000
Restarts will generate additional time history files: runname_0002.rad =>
runnameT02runname_0003.rad => runnameT03etc.

204
Global Variables are always output by default

Example of typical energy balance, mass and time step display

205
Requesting output for
entities
/TH/Keyword/th_groupid
Var_IDn = DEF default output for that
type as defined in help.
Obj_IDn = Entity ID(s) to be output
HyperMesh: Create using Output
Blocks
HyperCrash: Create using Data
History >
Time History

Time History Output Requests: Example for Rigid Wall

Individual variables can be requested for output, or groups of variables


The variable group, DEF, can always be requested to output a set of default variables

206
Option For Time History File
If /ABF (Altair Binary File) option is used, then graphs will automatically have axes and
units plotted appropriately
Output is to *.abf

/TH/PART: request default output for all parts to monitor energies


/TH/INTER: request default output for all interfaces to monitor contact
forces Example:

Other outputs :
/TH/MONVOL: request default output for all monitored volumes
/TH/SECT: request default output for all section forces
/TH/ACCEL: request default output for all accelerometers
/TH/RBODY: request default output for the major rigid bodies where constraints are
set on the master nodes.
/TH/SPRING: request default output for the springs which important connection
between components (bolts).
/TH/NODE: request default output for the reference nodes (imposed displacement,
measurement points,…)
Additional Time History: /ATH, /BTH, … /ITH

207
In case outputs are defined for a large number of elements, they can be stored in
different time history files (T01a, T01b,…,T01i) with a smaller storage frequency in
order to avoid very large time history file.
/ATH/SPRING: request default output for the connection springs (spot-welds, seam-
weld).

Example

208
What Do You Think?
In explicit dynamic run using RADIOSS, I am trying to output energy plot. I have defined TH/
PART in a deck file. Whenever I am trying to load *.T01 file into HyperGraph, its showing
following error:
The time history file type was changed to Type 4, that is TFILE/4 was the change made (Binary
Format) in the engine file and the file is readable now (see image below).

13.1 Frequently Asked Questions (Post-processing)


Using Shell Elements, I Asked for Strain Time History Output (Runnamet01) And
Animation Files, But The Values Remain Equal To Zero; Why?
The strain tensor is not computed by default; it must be asked for in the RADIOSS
input file (Runname_0000.rad) by setting flag Istrain (flag to compute strains for post-
processing) to 1 in option /DEF_SHELL or in the shell property set.
On the contrary, the strain tensor is always computed and available for Law 25 and
Law 27.

209
What Are the Stresses SIGX, SIGY, And VONM In Animation Files If I Use
Integration Points For The Shells?
The stresses SIGX, SIGY… in Animation files represent the mean stresses through the
thickness of the shell element. The VONM stress represents the Mises criteria applied
to these mean stresses SIGX, SIGY… In the same way the stresses F1, F2, F12, Q1 and
Q2 given in Time History correspond to these mean stresses.
These mean stresses are computed by summation of the stresses at each integration
point, averaged by the integration weights (refer to the integration weights table in
the RADIOSS Theory Manual). They are used for the internal forces calculation.

Where Is the Plastic Strain Value Computed If I Ask For /ANIM/ELEM/EPSP for
Shell Elements?
This value corresponds to the plastic strain value on the neutral fibre.
Is It Possible to Get More (Or Less) Animation Files While A Computation Is
Running?
Yes, it is possible to write an Animation file by writing a control file in the data
directory.
For the run number nn (/RUN/Runname/nn in the RADIOSS Engine input file), you
have to write the file Runname_nn_0000_[C].rst with the process /ANIM in it.
RADIOSS Engine writes an Animation file at this time.
The other options available with control files are described in the Control File (C-File)
file.
In order to change the Animation files writing frequency, you have to stop your
RADIOSS computation while writing a RESTART file, by using a control file (option
/STOP). Then you can chain a second run with a different frequency for the Animation
files writing.
How Can I Plot Deleted Elements to Understand The Propagation Of A Fracture?
Select Preferences > Option > Visualization menu in HyperView to display Eroded
Elements. This will help you to understand the propagation of a fracture

210
Using Shell Elements With Qeph Formulation (Ishel=24), The Hourglass Energy Of
The Part And The Subset Are Not Equal To Zero In Time History; Why?
When looking to the SUBSET or the PART in Time History, the hourglass energy is not
zero.
This is because energy absorbed due to the numerical damping is output there. This
means, in output the place of hourglass energy has been used to present this viscous
energy.
The viscous energy is related to coefficient dn for shell property which using QEPH
(Ishel =24) and QBAT and DKT18 (Ishel =12 or Ish3n =30).
The energy corresponding to the physical stabilization of hourglass is counted as
internal energy for this formulation.

211
212
14 Example: Three Point
Bending (HyperMesh)
(this example is taken form the Online Help - Tutorial RD 3595)
This tutorial demonstrates how to set up 3-point bending model with symmetric
boundary conditions in Y direction.

Model Description
• UNITS: Length (mm), Time (s), Mass (ton), Force (N) and Stress (MPa)
• Simulation time: in Engine file [0 – 6.601e-002 s]
• Only one half of the model is modeled because it is symmetric.
• The supports are totally fixed. An imposed velocity of 1000 mm/s is applied on the
Impactor in the (–Z) direction
• Model size = 370mm x 46.5mm x 159mm
• Honeycomb Material /MAT/LAW28: HONEYCOMB
• [Rho_I] Initial density = 3.0e-10 ton/mm3
• [E11], [E22] and [E33] Young’s modulus (Eij) = 200 MPa
• [G11], [G22] and [G33] Shear modulus (Gij) = 150 MPa
• Elasto-Plastic Material /MAT/LAW36: Inner, Outer and Flat
• [Rho_I] Initial density = 7.85-9 ton/mm3
• [E] Young’s modulus = 210000 MPa
• [nu] Poisson’s ratio = 0.29

213
Strain Curve:

• Elastic Material /MAT/PLAS_JOHNS: Impactor


• [Rho_I] Initial density = 8e-9 ton/mm3
• [E] Young’s modulus = 208000 MPa
• [nu] Poisson’s ratio = 0.29

Step 1: Load The Radioss User Profile


1. Launch HyperMesh Desktop.
2. From the Preferences menu, select the User Profiles or click the userProfile-24 icon in
toolbar.
3. Select RADIOSS (Block140) and click OK.

Step 2: Retrieve The Radioss File


1. Click File > Import > Solver Deck or click fileImportGeometry-24.
2. Click the Select File icon files_panel to open the BENDING_0000.rad file you saved to
your working directory from the radioss.zip file. Refer to Accessing the Model Files.
3. Click Import.
4. Click Close to close the window.

Step 3: Create And Assign Material And Property For Hcfoam


1. In the Model browser, right-click and select Create > Material. The new material
appears in the Entity Editor.
2. For Name, enter Foam.
3. For Card Image, select M28_HONEYCOMB and click Yes to confirm.
4. Input values, as shown below:

214
5. In the Model browser, right-click and select Create > Property to create a new
property.
6. For Name, enter Foam and set the new property Card Image as P14_SOLID. Leave all
the settings as default, except for ISOLID which should be set to 24.
7. In the Model browser, right-click on the component HCFoam and select Assign. Assign
Foam as the Prop_Id and Foam as the Mat_Id.
8. Click Apply.

Step 4: Create And Assign Material And Property For The Component Inner
1. In the Model browser, right-click and select Create > Material. The new material
appears in the Entity Editor.
2. For Name, enter Inner.
3. For Card Image, select M36_PLAS_TAB and click Yes to confirm.
4. Input the values, as shown below:

215
5. In the Model browser, right-click and select Create > Property to create a new
property.
6. For Name, enter Inner and set Card Image as P1_SHELL. Leave all the settings as
default, except for Ishell which should be set to 4 and Thick which should be set to
9.119e-01.
7. In the Model browser, right-click on the component Inner and select Assign. Assign
Inner as the Prop_Id and Inner as the Mat_Id.

Step 5: Create And Assign Material And Property For The Component Outer
1.In the Model browser, right-click on the material Inner and select Duplicate. Name the
new material Outer. This creates a new material that is identical to the source
material.
2.In the Model browser, right-click on the property Inner and select Duplicate. Name
the new property Outer. This creates a new property that is identical to the source
property.
3.In the Model browser, right-click on the component Outer and select Assign. Assign
Outer as the Prop_Id and Outer as the Mat_Id.

216
Step 6: Create And Assign Material And Property For The Component Flat
1. Follow the procedure described in Step 5 with Outer replaced by Flat.

Step 7: Create And Assign Material And Property For Impactor


1. In the Model browser, right-click and select Create > Material. The new material
shows up in the Entity Editor.
2. For Name, enter Impactor.
3. For Card Image, select M1_ELAST.
4. Input the values, as shown below:

5. In the Model browser, right-click on the property Inner and select Duplicate. Name
the new property Impactor. This creates a new property that is identical to the source
property.
6. In the Model browser, right-click on the component Impactor and select Assign.
Assign Impactor as the Prop_Id and Impactor as the Mat_Id.

Step 8: Create And Assign Material And Property For Support


1. Follow the same procedures as in Step 5. Create a copy of Impactor property and
material with name support and assign it to component support.

Step 9: Create A Rigid Body To Make Impactor And Support Rigid 1. In the Model
browser, right-click and select Create > Component.
2. For Name, enter Impact rigid.

217
3. Select any color for easy visualization.
4. Set Card Image to None.
5. Go to the 1D page, select the rigids panel.
6. Verify that you are in the create subpanel.
7. For dependent switch to comps.
8. For primary node switch to calculate node.
9. Click comps.
10. Select Impactor, then click select.
11. Click create.
12. Click return to exit the panel.
13. Similarly, create rigid body for Support component in a collector with the name
“Support rigid” using Steps 9.1 to 9.12.

Step 10: Define Imposed Velocity And Boundary Condition For The Impactor
1. From the Utility page, start the BCs Manager.
2. For Name, enter IMPOSED_VELOCITY, set Select type to Imposed Velocity and set the
GRNOD to Nodes.
3. Click nodes and select the master node of the rigid body of the Impactor, as shown in
the following image.

218
4. Set the Direction as Z.
5. Set Scale Y to -1000.0 as the direction of velocity is opposite to the global Z-axis.
6. Set the Curve ID to Select curve.
7. Select the predefined curve to Func1.
8. Click create to create the imposed
velocity.
9. For Name, enter Impactor_constraints,
set Select type to Boundary Condition
and set the GRNOD to Nodes.
10. Click nodes and select the
master node of the rigid body.
11. Check all the degrees of
freedom to constrain, except Tz.
12. Click create to create the
boundary condition.

Step 11: Define Fixed Boundary Condition For The Support


1. From the Utility page, start the BCs Manager.
2. For Name, enter Support_fixed, set Select type to Boundary Condition and set the
GRNOD to Nodes.
3. Select the master node of the rigid body created on Supporter, as shown in the
following image.
4. Check all the degrees of freedom.

219
5. Click create to create the boundary condition

Step 12: Define Symmetry Boundary Condition For The Foam, Inner, Outer And
Flat
1. From the Utility page, start the BCs Manager.
2. For Name, enter SYMMETRY_XZ, set Select type to Boundary Condition and set the
GRNOD to Nodes.
3. Select the nodes of the foam, inner, outer and flat, as shown in the following image.
4. Check the degrees of translational degrees of freedom Y and rotational degrees of
freedom X and Z to constraint.
5. Click create to create the boundary condition.

220
6. Click close to exit the BC Manager.

Step 13: Define Contacts Between The Beam And The Support
1. Launch the HyperMesh Solver browser from View > Browsers > HyperMesh > Solver.
2. In the Solver browser, right-click and select Create > INTER > TYPE7.
3. Enter the values, as shown below:
4. Set the Surf_id (M) for the master selection to Components and select the Support
component.
5. Set the Grnod_id (S) for the slave selection to Components and select the Flat
component.
6. Similarly create the contact for Impactor with Outer, as shown below.

221
Step 14: Define The Self Contact Between The Beam
Components
1. Using the directions in Step 13, create a new Type 7
interface named Self with the components Outer, Inner,
and Flat as Master and the same components Outer,
Inner, and Flat as Slave. This will make the components
self-contact instead of self-penetrate. Verify that the
interface has a Fric of 0.1 and Gapmin of 0.2.

Step 15: Create Interface Time History


1. Right-click in the Solver browser and select Create > TH
> INTER.
2. For Name, enter IMPACTOR.
3. Switch the entity selector to groups.
4. Click groups and select the interfaces Impactor and Support from the list.
5. Click OK.
6. Set NUM_VARIABLES to 1 and Data: Var to DEF

222
Step 16: Creating Control Cards And Output Requests
1. Launch the HyperMesh Solver browser from View > Browsers > HyperMesh > Solver.
2. Right-click in the Solver browser general area to create the cards shown below with
the given values for each parameter:
Control Cards And Output Requests

223
Step 17: Export The Model
1. Click File > Export or click the Export icon fileExportSolver-24.
2. For File:, navigate to the destination directory where you want to export to.
3. For name, enter 3BENDING and click Save.
4. Click the downward-pointing arrows next to Export options to expand the panel.
5. Click Merge starter and engine file to export solver deck as one file (or export
separately).
6. Click on Export to export solver deck.

Step 18: Open RADIOSS Manager


1. Go to Start > Programs > Altair HyperWorks 2017 > RADIOSS.
2. For Input file, browse to the exercise folder and select the file 3PBENDING_0000.rad.

Step 19: Review The Listing Files For This Run And Verify On The Results
1. See if there are any warnings or errors in .out files.
2. Using HyperView, plot the displacement, strain contour and vectors

224
Verify The Results

Contact Force For Impactor Interface

225
226
15 Example: Three Point Bending
(HyperCrash)
In the course of this tutorial you will be introduced to the preprocessor Altair
HyperCrash as an alternative to Altair HyperMesh. After you have completed the
standard meshing task (i.e. geometry cleanup, 1D-2D-3D meshing) with HyperMesh,
it may be easier to continue the model set-up with Altair HyperCrash. HyperCrash
offers some very helpful crash/impact related routines which accelerate the model
building process enormously (of course, you can complete the entire model set-up
with HyperMesh, too).
Note: This tutorial corresponds pretty much to the previous example.

This exercise demonstrates how to set up 3-point bending model with symmetric
boundary conditions in Y direction (across the XZ plane).
Model Description
• UNITS: Length (mm), Time (s), Mass (ton), Force (N) and Stress (MPa)
• Simulation time: in Engine file [0 – 7.01e-002s]
• Only one half of the model is modeled because it is symmetric.
• The supports are totally fixed.
• An imposed velocity of 1000 mm/s is applied on the Impactor in the (–Z) direction
• Model size = 370mm x 46.5mm x 159mm

227
Honeycomb Material /MAT/LAW28: HONEYCOMB
• r= 3.0e-10 ton/mm3 ; [Rho_I] Initial density
• Eij = 200 MPa [E11], [E22] and [E33] Young’s modulus
• Gij = 150 MPa [G11], [G22] and [G33] Shear modulus

Elasto-Plastic Material /MAT/LAW36 Inner, Outer and Flat


• r = 7.85 -9 ton/mm3 [Rho_I] Initial density
• E = 210000 MPa [E] Young’s modulus
• n= 0.29 [nu] Poisson’s ratio
Strain Curve:

Elastic Material /MAT/PLAS_JOHNS: Impactor


• r= 8e-9 ton/mm3 [Rho_I] Initial density
• E = 208000 MPa [E] Young’s modulus
• n= 0.29 [nu] Poisson’s ratio
We first we’ll import the given RADIOSS model (BENDING_0000.rad; Starter File) into
HyperMesh, then define the elasto-plastic material used for the Inner, Outer and Flat
components. Subsequently, we will repeat similar working steps in HyperCrash (i.e.
we will create the Honeycomb material). In this way you’ll get an idea on how
HyperCrash works. Eventually, it will be up to your preference whether you will
continue with HyperMesh or whether you’ll switch over to HyperCrash after meshing
is completed in HyperMesh.
Note: A rather thorough introduction into HyperMesh e.g., meshing, organizing model
entities in collectors, creating materials and properties, assigning properties to
elements etc. is given in the free eBook “Practical Aspects Of Finite Element
Simulation”. In addition, you may also view the course material about HyperMesh
available in the Teaching & Certification Program.

228
Thus, general information about HyperMesh is kept deliberately short in here.

Import The RADIOSS Mesh Model In HyperMesh


Launch HyperMesh Desktop from the Start Menu
Load the RADIOSS User Profile
Retrieve the RADIOSS file; Click File > Import > Solver Deck > Select File and import
BENDING_0000.rad

Create And Assign Material And Property For The Component Inner (see image)
 Create > Material
 For name, enter: Steel.
 Select Card image as M36_PLAS_TAB.
 Elastic Plastic Piecewise Linear
Material; This law models an
isotropic elasto-plastic material using
user-defined functions for the
work- hardening portion of the stress-
strain curve (for example, plastic
strain vs. stress) for different strain
rates.

The stress-strain data (from above) are


already part of the imported model file.
The stress - strain data are defined as an
entity of type “curve” (named Steel- mild)).

229
In the Model Browser, right-click to select
Input the values, as shown below:

Create New Property


1. Right-click in the Model Browser to create a new property. Enter the name as
Shell-Skin and set the new property Card image P1_SHELL. Leave all the settings

230
as default except for Ishell which should be
set to 24 (QEPH element formulation) and
Thick which should be set to 9.119e-01.
N = Number of integration points through the
thickness with 0 < N < 10 (see Comment 16)
Ithick = Shell resultant stresses calculation flag
(Integer)
= 0: default set to value defined with /DEF_
SHELL
= 1: thickness change is taken into account
= 2: thickness is constant
Iplas = The default value for Iplas in case of Law 36 and global integration (N=0 in shell
property) is Iplas =1: iterative projection
Select the Inner, FLAT, and Outer components in the Model Browser and edit the
Prop_Id and Mat_Id
Assign Shell_Skin as the Prop_Id and Steel as the Mat-Id.

231
Create And Assign Material And Property For Impactor
• From the Model Browser, right-click > Create
> Material.
• For name, enter Impactor-Steel.
• For Card image, select M1_ELAST.
• Input the values, as shown (and provided at
the beginning)
• Next, you may even create a new Property
collector from scratch, or you duplicate and
edit the existing one.
Right-click on the property ShellSkin in the
Model Browser and select Duplicate.
Name the new property
Impactor. This creates a new
property that is identical to the
source property. Update the
element formulation to 4
(Belystcho-Tsay Q4 Improved
formulation). Number of integration points, N, can be set to 1 because these shells
will become part of a rigid body. Turn Ithick and Iplas off.
• Select the Impactor component in the Model Browser and edit the Prop_Id and
Mat_Id. Assign Impactor-Skin as the Prop_Id and Impactor-Steel as the Mat-Id.
• Save the model (HyperMesh file format, Engine file _0000.rad (here we named it:
3point_bending_continuation_0000.rad

Create And Assign Material And Property For Support


As mentioned in the intro to this tutorial, we are going to use HyperCrash as an
alternative Preprocessor as well. Please note, that the “cut” in here, i.e. the switch
from HyperMesh to HyperCrash at this particular point in the model building process
is chosen arbitrarily!

232
Import The RADIOSS Model in HyperCrash

Start HyperCrash from the Start Menu


• Click Run.
• From the menu bar, select File Import RADIOSS > Select the Engine file from the
respective directory (here we are going to import the model 3point_bending_
continuation_0000.rad.

233
In the following we will create and assign material and property for Support
• From the menu bar, select Model Material.
• In the Window (Model Browser), right mouse click and select Create New > Elastic >
Linear elastic

• For Title, enter Rigid Material; then enter all the material data as shown in the
following image

234
To assign the material to “Support” click the highlighted (yellow) icon “Include picked
parts”, then select the part “Support” in the graphics area. Confirm with a Yes-click in
the lower right corner.

• Finally, click Save.


In order to visually check whether the material was assigned correctly, select the
material of interest and then click on the “glasses” icon: all parts will be masked but
the once with the selected material.

You may also view the material named “Steel” which we previously defined in
HyperMesh.

235
(make sure that the “Scale factor for ordinate (stress) (Fscale_1] is set to 1.
• To view or to edit the stress - strain data, simply right click on the “Yield stress
function” (here Steel-Mild) and activate “Edit Function ...”

236
Create And Assign A New Material For HCFoam
• Create New > Other > Honeycomb orthotropic (28).

• For Title, enter Foam.


• Enter all the material data as in the following image:

• Right click on the Yield stress function 11 field

237
• “Select in Model…” to select a curve already present in the model: function with ID of
5 Honeycomb-Yfunc11

• Repeat this process for the Yield functions as shown in the following image.

• Click the Tree tab and select the part HCFoam (7) on the tree.
• Click to show only this part.
• Click on Material tab, then click “Include picked parts” and select the part HCFoam in
the graphics area.
• Confirm the selection with “Yes” (as before) and “save” the changes.

238
Review The Property Of Inner, Outer and Flat
The Property of Inner, Outer and Flat which was defined before (in HyperMesh) can
be reviewed in a similar fashion to the material check before.
From the menu bar, click Model > Property.

In the list of available Properties select the one of interest.


Here we review Shell-Skin (which was assigned to Inner, Outer and Flat)
In order to check to which Parts the Property “Shell-Skin” was assigned, select
Property in the tab, then select the property of interest (here “Shell-Skin”) and
eventually click on the glasses icon. Only parts with the respective property will be
shown.

239
In a similar way we review the Property named “Impactor”

Currently, this property is only assigned to the “Impactor” (as shown). However, the
“Support” will share the same property.
To assign the property to the part
“Support” we first display all parts

Then, while the property “Shell-Skin” is


still highlighted in the list of available
properties, click on the “Include picked
parts” (at the bottom), select the
“Support” from the display, confirm the
selection with “Yes” and finally save the
changes.
Cross-check which parts are assigned the
property “Shell-Skin” as illustrated above.

240
Create And Assign A Property For HCFoam
If you are not already inside the Property tab or the tab is not listed yet, just go to

Inside the Property tab, right mouse click


> Create New > Volume > General solid (14).
For Title, enter Foam.
In the Flag for solid elements formulation, select HEPH 8-node solid element.

In order to assign the property “Foam” to the respective part, you can follow the
process described earlier. Alternatively, you can check the following (other options
exist still):
• Click the Tree tab and select the part HCfoam
• Go back to Property tab, then

• click the icon “Add selected parts from Tree”. Finally, click on “save”.

241
Again, check whether the part named “HCfoam” is really linked to the property
“Foam” (same process as above; glasses icon)
To show all parts click the Tree tab and select the BENDING assembly on the tree.

242
Create Rigid Body For Impactor And Supports
In the menu bar, select Mesh Editing > Rigid Body

In the window right mouse click to select “Create New”, and enter the name
“Impactor”.
Use the “Include picked parts” (lower part of the window; you may need to toggle
[Gnod_Id
Slave nodes] option to select “Impactor” on graphic screen as shown)

the “slave nodes / dependent nodes” will be highlighted. Upon activating the “Save”
button the rigid body will be created. Its Master node /independent node will be
created automatically (as requested).

243
Create Rigid Body For Support
The process is like before. The rigid body is named “Support”. Once you have
completed the process with “Save” you can/should cross -check the result.
In the Menu bar go to Rigid Body (if it is not listed anymore, open the Rigid Body tab
in the
Menu bar, select Mesh Editing > Rigid Body

Define Boundary Conditions For The Model


From the menu bar, select LoadCase Boundary Condition In the window right mouse
click then Create New.
In the Boundary Condition creation field, enter “BC” for Title.

244
Use either “node by picking selection” (by box selection) to select master node of
“Impactor”. Constrain all DOF except translation in Z. To constrain the nodes, check
TX, TY, RX, RY and RZ.

245
BC as applied to the Master node of the part “Impactor” (the Rigid body is masked)
Repeat the same process to create “Support Fix” and “Symmetry” BCs as shown
below.
2. Constrain all DOF by selecting TX, TY, TZ, RX, RY and RZ as shown in the following
image. To create the constraint, don’t forget to activate the “save” button.

246
Create Symmetry Boundary Conditions
In the Boundary condition creation field, enter “Symmetry” for Title. Recall that the
plane of symmetry is the xz plane!
Before looking at the input data below, please ask yourself first: how should the
symmetry boundary conditions look like i.e., which DOF’s need to be fixed, which
DOF’s need to be released?

• Make sure and display the respective parts first.


• Click the Tree tab and select the parts Inner, Outer, HCfoam
and Flat.

• Click to show only these parts.


It may be helpful to press p (on the keyboard) to change the perspective visualization
• From visualization tool bar select the view YZ as shown in the image
• Click on Boundary Condition tab, and define the constraints of the nodes at the
symmetry plane as shown

247
• Select the nodes of the symmetry plane by using the “Add nodes by box selection”

For Time function, use the predefined curve


in the model “Funct 1”
For Direction, select Z (translation) and

• Eventually, click Save & Close.

Define Impactor Velocity


• From the menu bar, select LoadCase Imposed Velocity
• As with the creation of Material, Property, BC, ..., right mouse click in the window
(see below), then select “Create New”, for Title, enter IMPOSED_VELOCITY

248
Note: FscaleY defines the Ordinate scale factor for fct_IDT Default = 1.0 (Real). Its units
(default) are m/s, hence in this exercise we need to define Fscale as -1000 (to account
for 1000mm/s and the negative z-direction)
The imposed velocity is applied to the Master node of the “Impactor Rigid Body”.

As always, don’t forget to click on “save” to complete this working step. The model
should look like below.

249
Define Contacts for The Model
Display the parts Flat and Support only.
• Click the Tree tab and select the parts Flat and Support.
• Click to show only these parts.
From the menu bar, select LoadCase > Contact Interface
• Right mouse click in the window and choose Create New > Multi usage (Type 7)

250
• For Title, enter Support.
• Set Coulomb friction to 0.1.
• Set Min gap for impact active to 0.2.
Using select Flat component as Slave, and
Support as Master

Finally (for this step), click Save.

Repeat the same process to create contact between Outer and Impactor.

251
Display the parts Outer and Impactor only (in the Tree tab; same process as before)
Click the Contact Interface tab (if not listed /shown anymore you need to open it up
under LoadCase > Contact Interface)
• In Title, enter Imp_Outer
• Set Coulomb friction to 0.1.
• Set Min gap for impact active to 0.2
• Using select Outer Part as Slave and Impactor as Master
• Click Save

252
Repeat the same process for self-impact for Outer, Inner, Flat but as a self-impact

Start with displaying only the respective parts (process as before)


• In the Title, enter Self.
• Set the Coulomb friction to 0.1.
• Set the Min gap for impact active to 0.7
• Using select components Outer, Inner and Flat (selection process as before)
• Click Save.

253
Clean The Model
The intention behind the “cleaning” process is - like in real life - to get rid of unused
entities, such as properties, materials etc.

From the Menu bar, select Mesh Editing > Clean.


In the Mesh Editing > Clean panel activate the option ALL.
As stated before unused entities will be
removed.
In the context of this tutorial we added a
dummy Property (not referenced in the
model). This results in a Warning and
eventually in the message:
“1 PROP removed”
Export The Model, Write Starter
(_0000.rad) And Engine (_0001.rad) Files
Click Model > Control Card and select the
control cards in the images below.
(Be sure to save each control card before
moving on.)

254
/ANIM/DT - Frequency of Writing Animation Files
Write animation files (A-files) at a time frequency equal to Tfreq, the first file being
written at time Tstart. The animation file name will be “RunnameAnnn” where
Runname is the Run Name (see /RUN) and nnn is
the file number.
Tstart Start time for writing the animation files
Tfreq Frequency for writing the animation files
/ANIM/Eltyp/Restype - Animation of Element
Data for Specified Result
This generates animation files containing element
data for the specified result
Eltyp Type of element
ELEM - The variable is saved for all types of elements; except where not applicable.
BRICK - Brick elements
SHELL - Shell elements
BEAM - Beam elements
SPRING - Spring elements
TRUS - Truss elements
Restype (name of the variable to be saved)
EPSP - Plastic strain e1
VONM - von Mises stress

255
Eltyp Brick, Quad, Shell elements etc.; INTER-nterfaces (Type 7, 10, 11, 19,
and 20)
Keyword 3 Time step control

/DT/Eltyp/Keyword 3/Iflag Time Step with Entity Selection and Time Step Control
type: DEL-The element
which fixes the time step is removed. This option is the default for shell elements.
If “Keyword3”= DEL and “Eltyp” = INTER,
the impacted node which fixes the time
step is removed from the interface
DSTSCAScale factor on time step for the
option defined by “Eltyp”. Default = 0.9

/MON - CPU Time Estimation; Provides


an estimation of the CPU time spent for
each processor
Format /MON/Keyword2
Keyword2 Activation flag: ON (default) or OFF
For each process, the following information is written:
PROC Core number
CONT.SORT CPU time spent in sorting algorithms for contact interfaces

256
CONT.F - CPU time spent in computing interface forces Keyword2 (cont)
ELEMENT - CPU time spent in computing the elements, including material law
computation
KIN.COND. CPU time spent in computing kinematic conditions

INTEGR. CPU time spent in time integration

I/O CPU time spent in input/output subroutines

TASK0 CPU time spent in various non-parallel subroutines, including time


ASSEMB CPU time spent in forces assembly

RESOL Total CPU time, except time for reading first restart file

When /MON/ON is active, /PRINT will provide information about current and
remaining elapsed time.

/PRINT - Sets printout frequency for list file


Format /PRINT/Nprint
Nprint Printout every Nprint cycle on file “Runname_
nnnn.out”; Default = -1000
Note: If Nprint is negative, printout is made in the
standard output as well as in the list file. If it is
positive, the printout is only made in the list file

257
/RUN - Run Number for Restart - Identifies the run number
Format /RUN/Run Name/Run Number/Restart Letter
Run Name Character variable that identifies the problem solved (Character,
minimum 4; maximum 80 characters) Run Number 1, 2, ...
Restart Letter I, J, K , . . . if several restarts were saved in the previous run.
Tstop Final time for the run.

/TFILE - Frequency of Writing Time History File; Defines the frequency of writing
the time history file of T-file “RunnameTnn”.
Format /TFILE/Type
Type No value: Built-in format of current RADIOSS version
= 1: Binary (not readable by most RADIOSS post-processors)
= 2: Coded ASCII 32-bit
= 3: ASCII
= 4: Binary IEEE 32-bit (Default value)
Recommended options are default (no value) or four, which can be read by
HyperGraph
DThis Time frequency to write data on history plot file T-file
Finally, in the Menu bar, select File Export RADIOSS ….
In the output window that opens, enter the name 3PBENDING and click OK.
Leave the Header of RADIOSS file window empty and click Save Model.

258
The Starter/Engine file 3PBENDING_0000.rad is written.

Open RADIOSS Manager from windows Start menu, and run the model
3PBENDING_0000. rad
Review the listing files for this run and verify on the results
• Using HyperView, plot the Von Mises, and Plastic strain contour (*.h3d file)

Von Mises Stress (MPa)

259
Plastic Strain
Contact Force Between Impactor & Outer

260
16 Time Step Control
The time incrementation in RADIOSS is fully automatic and a priori requires no user
intervention. The step used for time integration (or moving forward in time) can be
calculated using two different methods. The method used depends on the type of
simulation being performed.

The two time step methods are:


• the element time step,
• the nodal time step.

The time step used by the solver is the largest possible time step, as determined by
the Courant condition that will maintain stability. If the default large strain
formulation is used, the time step is computed at each cycle.
Large element deformation can give a large time step decrease. If the deformation is
too large, negative volumes can result, which make it impossible to invert the Jacobian
matrix and to integrate the stress over the volume. If the small strain formulation is
used, assuming a constant Jacobian matrix during time and also a constant volume,
all spatial variables are defined at t = 0. This is either the beginning of the analysis or
the time at which the small strain formulation is initiated. If the sound speed is
constant, the time step thus becomes constant. Using this formulation, the time step
has no effect on the computation since the initial volume is used.

261
The following topics are addressed in this chapter:
• Analysis Stability Condition
• Time Step Calculation
• Elements
• Rigid Bodies
• Interfaces
• How to Speed up Computation:
• Time step control
• Parallel Processing
• Extended Single Precision
• Advanced Mass Scaling

Analysis Stability Condition


Critical Time Step – The time needed for information (stress, pressure wave) to cross
the smallest distance in an element
Stability – The model is stable if it is computed with a time step less than the smallest
critical time step of the model
Stability Condition depends on two factors:
Size of the smallest element > Numerical Sound propagation
speed > Physical

Element Time Step Control


The stable element time step is as:

where:
l is the element characteristic length, c is the speed of sound in the material.
This is the default setting.

262
The element time step is computed at the same time as the internal forces. The
characteristic length and the sound speed are computed for each element in every
cycle.
The computed time step is compared to a minimum time step value and a scale factor
is applied to ensure a conservative bound. Different minimum time step values can
be given to different element types by using the option: /DT/Keyword where Keyword
is defined in the user manual as the element type.
If deformation is large enough for the time step to reach the minimum defined value,
three options are possible under the user’s control:
• Stop the analysis when the minimum time step value is reached. This is the default for
brick and quadrilateral elements.
• Delete the element(s) defining the time step. This is the default for shell elements.
• Implement small strain formulation using a constant time step. This only works for
shell and brick elements
These options are defined using a third keyword: STOP, DEL, CST, AMS or SET

Shell Element Time Step


Critical length:

The critical length can vary with:


• Shell formulation
• Hourglass control parameters
• Numerical damping
• Membrane damping

263
Solid Element Time Step
Time Step

Critical Length

Nodal Time Step Control


The nodal time step is calculated after the computation of all the internal forces at
each node using the following equation:

where m is the nodal mass and k an equivalent nodal stiffness.


The nodal stiffness is one half of eigen value from element stiffness matrix; for a truss
element this value is equal to the diagonal term of the stiffness matrix. It is computed
from the accumulation of element and interface stiffness’. These stiffness’ are
obtained during internal force computation.
For a regular mesh, the element time step and nodal time step conditions are
identical.
To select the nodal time step when running RADIOSS the option /DT/NODA has to be
used. As for the element time step, minimum time step and scale factors are required.
The default value for the scale factor is 0.9. If the minimum time step is reached, the
analysis can either be stopped or a mass scaling formulation can be applied. In the
latter case, mass is added to the affected nodes so that the time step remains constant
at the minimum value. This option can be enabled using the same third keyword as
used in the element time step control. It must be checked that added masses do not

264
affect the accuracy of results. Note: If one uses the nodal time step, the element time
step is ignored.

Time Step Control Limitations


Many of the time step control options influence the solution results. The solution of
the nonlinear dynamic response of a finite element system accurate is the numerical
model correctly represents the physical model. The critical time step given for finite
element system is determined by a theoretical approach in which the highest
frequency of the discretized system controls this value. Therefore, the time step
limitations are related to the model and cannot be changed without incidence on the
quality of results.
Using the “DEL” option can significantly alter the model, since elements and nodes
are removed without replacement. In fact, mass and/or volume is lost. Using either
/DT/NODA/CST or / DT/INTER/CST will add mass to the model to allow mathematical
solution. The added mass will increase the kinetic energy. This should be checked by
the user to see if there is a significant effect. Switching to a small strain option using
brick or shell elements also introduces errors.
Generally, in the study of the nonlinear dynamic response of a system, three physical
laws have to be respected:
• Conservation of mass,
• Conservation of energy,
• Conservation of momentum dynamic equilibrium.
The last one is generally respected as the equation of motion is resolved at each
resolution cycle. However, in the case of adding masses especially when using
/DT/NODA/CST option, it is useful to verify the momentum variation. The two other
conservation laws are not explicitly satisfied. They should be checked a posteriori
after computation to ensure the validity of the numerical model with respect to the
physical problem.

265
Time Step Checking In Pre-Processor
The element time step can be checked with HyperMesh (or HyperCrash). Example
with the quality panel of HyperMesh for shell elements:

Time Step Checking In RADIOSS Starter Listing File


The estimated smallest time step are printed in the starter output file
(runname_0000.out) for each type of
element (solid, shell, spring,…) and also
the smallest estimated nodal time step
The time step for the computation is
reduced by the scale factor (default=0.9)
and can be smaller than the estimated
one according the loading.

Why is it so?
When the RADIOSS engine is getting started the time step is decreasing and thereby
increasing the computation time? What may be the possible cause of it?

266
From the engine out file we can see that the model is having very low time step. It can
be an element, a node or an interface.
• For an element, check the related material (especially its Young modulus and density
in case of an elastic-plastic material; and its viscosity in case of a visco-elastic
material). There must not be an error in the units system that this data is given in.
Check the size of the element, since elemental length is proportional to time step.
• For a node, check the characteristics of connected elements. If the node is on the
master side or the slave side of an interface, this interface must be verified.
• For an interface, the gap of the interface must be verified if some failure happens on
the master or the slave side of the interface.
Once this is cleared you can run the simulation without time step being dropped. You
can also impose a time step in order to control the time step drop.

Scale Factor
The time step scale factor value is set to 0.9 by default to ensure the stability, taking
into account the approximations made in its evaluation. However, in some special
cases, a smaller value is required. For example, if the structure becomes stiff suddenly
(Fig. 1), due to the material behavior, the critical time step computed for the next step
will be smaller than the current one. Therefore, the critical value computed by
RADIOSS may be greater than the theoretical critical time step.

Typical Behavior Of Foam


To avoid this problem, the time step scale factor is reduced. The numerical experience
shows that a scale factor around 0.66 may be small enough to avoid time
discretization stabilities.
Another typical case is when only one element is used in the model. The shock wave
propagation for a mesh with only one element cannot be studied correctly, the time

267
step approximation by RADIOSS may return a higher value, which can be found by
theory. In this case, the use of DTsca= 0.1 is recommended.

Comments
1. For foams: DTsca = 0.66
2. Model with one element: DTsca = 0.1
3. Model with two finite elements: DTsca = 0.2
4. Model with more than three finite elements: DTsca = 0.9 (good in general)
5. Never use a scale factor greater than 1.0

How To Speed Up Computation: Time Step Control For Impact


Theoretical best solution:
• No time step control in the model (element time step)
• The computation can take lot of time because the time step will drop due to the
element deformation
Element Time Step Control:
/DT/Eltyp/CST : switch to small strain when the element is deformed
Where Eltyp is BRICK or SHELL
/DT/INTER/CST : add mass on the contact when time step drops
E.g. The shell elements will switch to small strain formulation when their time step
become smaller than target time step (1E-06) / scale factor (1.11E-06 in the example
below)

Nodal Time Step Control:


Most often used option
/DT/NODA/CST : add mass on the node when time step drops
Need to check amount of mass added and where it is added
• Total mass of the model from time history file (T01)
• Animation : /ANIM/NODA/DMAS and /ANIM/MASS

268
E.g. Mass will start to be added on the node when its time step becomes smaller than
target time step (1E-06) / scale factor (0.67) = 1.5E-06.

How To Speed Up Computation: Multi-Threading SMP


Multi-threading on a single socket (or physical CPU) : SMP
• Use option to define the number of thread for the computation : -nt #threads
• For efficiency, the number of threads should equal or less than the number of cores
in the socket (CPU)

Typical Usage For Desktop Or Laptop Machine

How To Speed Up Computation: Extended Single Precision


By default RADIOSS is computed with double precision
Extended single precision allows reduction of the CPU cost from 10 to 40% depending
on the load-case and the model.

Example of executable names


Results may be different than running with double precision, but care is taken still use
double precision on critical RADIOSS calculations

269
How To Speed Up Computation: Advanced Mass Scaling
Typical speed-up from 3 to 5 times natural time step
Additional terms on the mass matrix to increase the time step without changing the
physical mass and inertia
The mass matrix can no longer be inverted, thus the equation of motion is solved using
the conjugated gradient
Applications:
Quasi-static load-case
• Roof Crush Computed With AMS
• Roof crash
• Seatbelt anchorage test (ECE-R14)
• Loading with imposed displacement
• Solid element model with small
mesh size
Method not recommended for dynamic load-cases (crash, high speed impact,…)

Is It Possible To Use /AMS For Crash Applications?


• AMS is not recommended for crash simulations. The conventional mass scaling, that
is /DT/NODA/CST (mass error preferably less than 2 percent) is advisable for crash
simulations whereas AMS is recommended for quasi static simulations. AMS modifies
high frequencies in the system and in a dynamic event this can change results. In
quasi-static simulation this isn’t an issue.
Starter deck (_0000.rad)
Add /AMS to model along with a group of parts to which AMS will be applied
If no part group is defined, AMS will be applied to the full model

270
Engine file (_0001.rad)
Add /DT/AMS with:
Recommended scale factor: 0.67
Target minimum time step value: 10 to 20 times
the natural element time step
Comments:
The AMS treatment has an associated CPU cost, so
this method becomes efficient when the target
time is > 10 times the elementary time step
AMS is compatible with Single Precision
This option is not compatible with all RADIOSS options. In these cases, the target time
step will not be honored
In the case that contact interface TYPE7 and TYPE11 are defined in the model, the
target time step should not be larger than 1.0E-05 seconds to avoid numerical
issues. However, this is case dependent.

Time Step Control (32 Minutes)

(Video by Rahul Rajan, Altair India)

271
16.1 Frequently Asked Questions (Time Step)
Why Is The Cycle Time Step Very Small?
First, check which object is responsible for the low time step in the listing file
(Runname_nnnn. out) of RADIOSS Engine. It can be an element, a node or an
interface.
• In case of an element, check the related material (especially its Young modulus and
density in case of an elastic-plastic material; and its viscosity in case of a visco-elastic
material). There must not be an error in the units system that this data is given in.

Check the size of the element, since where, lc is the characteristic


length of the element; c is the material sound speed; and dtsca is the scale factor.
• In case of a node, check the characteristics of connected elements. If the node is on
the master side or the slave side of an interface, this interface must be verified.
If a master node of a rigid body gives the time step, check the rigid body inertia
reported in the listing file (Runname_0000.out), at:
Rigid Body Initialization (for a rigid body made of only a few nodes, it may be necessary
to set a spherical inertia).
• In case of an interface, the gap of the interface must be verified if some failure
happens on the master or the slave side of the interface.

How Is It Possible That The Time Step Is Controlled By A Node Even While Using
Elementary Time Step (ie Did Not Specify /DT/NODA)?

The input file probably includes an Interface Type 7, Type 10 or Type 11.

272
In case of an impact between a slave node and a master segment, the stiffness of the
penalty spring for this impact is applied to both slave and master nodes, in order to
ensure stability.
More precisely, Interface Type 7 associates a penalty spring and damping, so that an
equivalent stiffness to the whole system is considered. The nodal stability time step

is computed as follows: where, M is the nodal mass and K is the nodal


stiffness

is the sum of the stiffness applied to the node by all elements and contact interfaces.
In case the deck includes Interface Type 7, Type 10 or Type 11, this time step is always
computed at all nodes of the deck, in order to ensure stability even if an elementary
time step is used.
But when the time step is taken onto a node, this node is probably impacting

(otherwise nodal time step computed with is greater than elementary


time step).

What Is The Difference Between The Messages “Delete Slave Node” And “Remove
Slave Node” Concerning Interfaces Type 7 In The Listing File (Runname_0001out)?
These messages are related to the option /DT/INTER/DEL. They both consist in the
slave node suppression from the interface. This node will no longer impact any master
segment for this interface.

273
The message “Remove Slave Node” is written when the interface kinematic time step
is lower than the minimum time step; which was specified in option /DT/INTER/DEL.
For the slave/ master couple which is considered, the kinematic time step is:

If it is such that Dtsca Dt ≤ Dtmin, then the slave node is suppressed from the interface.
The message “Delete Slave Node” is written when the stability time step of the
interface for the slave/master couple; which is considered is lower than the minimum
time step which was specified in option /DT/INTER/DEL.
The stability time step for the slave/master couple is computed as follows:

where, M is the minimum nodal mass among the slave node and the nodes of the
master segment, and Ki is the interface stiffness for the slave/master couple which is
considered (Ki is an equivalent stiffness to the penalty spring and the damping which
are applied between the slave node and the master segment).
When the stability time step for the slave/master couple is such that Dtsca Dt ≤ Dtmin,
then the slave node is suppressed from the interface.
It is possible that those messages would be written several times within the same
cycle for the same node and the same interface (if the node impacts several master
segments within the same cycle and several contacts get a lower time step than the
minimum time step which was specified in option /DT/INTER/DEL).

274
17 Helpful Tools
In this chapter you will learn more about the following helpful tools:
• Functions /FUNCT
• Table of Functions /TABLE
• Added Mass /ADMAS
• Group of Nodes /GRNOD
• Group of Elements /GRoption
• Surface Definition /SURF
• Sections /SECT
• Sensors /SENSOR
Functions /FUNCT
Functions are used to define X – Y curves. Examples are material curves, failure
models, and imposed motion.

For values that are outside the domain defined, the function is linearly extrapolated
with a slope defined by the last two points of the function.
This behavior is different than some other explicit solvers.

Table Of Functions /TABLE


Tables are used for some specific materials that may have temperature and/or strain
rate dependencies. Up to 4 dimensions are available for some table types.

275
Added Mass /ADMAS
Added Masses are used to model non-structural mass such as connectors, wires,
carpet, etc.
A number of options are available for how the mass is distributed, defined by the
ADMAS Type:

Note: HyperMesh only supports Type 0 and 1


Example: Added mass representing instrumentation is added to selected nodes in
trunk

276
Group Of Nodes /GRNOD/Keyword
Node Groups are used to define a list of nodes needed for any option that requires a
group of nodes. An example is a group of nodes to which an initial velocity is applied.
Some most often used keyword options are shown below.

277
Example Card Image:

Group Of Elements /GRoption


Element Groups are used to define a list of elements needed for any option that
requires a group of elements. An example is a group of elements which is needed for
a section definition.
Typical options:

Surface Definition /SURF


A surface is a set of 3-node or 4-node segments. Surfaces are used to define interfaces
and pressure loads. All nodes must belong to a 3- or 4-node shell or to a brick element.
*Some keywords have additional options:
/SURF/Keyword/EXT – Take the outer surface of the solid blocks (even it is made with
several parts)
/SURF/Keyword/ALL – Take all the facets of the
shell and solid elements.
Sections /SECT

278
A section is a cut in the structure where forces, moments or energy will be computed
and stored in output files (using /TH/SECTIO).
In RADIOSS /SECT, /SECT/PARAL and /SECT/CIRCLE can be used to define a section.

Sensors /SENSOR
A sensor can be used to activate or deactivate a variety of features within RADIOSS.
A variety of sensor types are available. Some common ones are shown in the table
below.

Examples of sensor usage:


• To deactivate the global rigid body defined on the car during the “flying” phase just
before the impact on the ground
• To define fire time of safety devices (seatbelt pretentioner, airbag,..)

279
How Can I Do It?
I need the simulation to terminate when the contact force reaches 100kN.Can I
perform it in
RADIOSS?
You can do this using /STOP/LSENSOR option, which stops the engine due to sensor
activation when the first sensor in the defined list is activated. Define /SENSOR/INTER
in the model and provide the normal force (100 kN) to activate this sensor. Recall this
sensor in /STOP/ LSENSOR which will stop the engine when the /SENSOR/INTER
sensor is activated

Restarting A Run In RADIOSS


If you need to continue the calculation without modifying any definitions in Starter,
run the (i+1)th Engine file (Runname_000(i+1).rad).
For Example:
run Starter: TEST_0000.rad run the first Engine: TEST_0001.rad
run the second Engine: TEST_0002.rad > to continue the run after TEST_0001.rad
RADIOSS could restart the calculation from the last restart file which was created by
running TEST_0001.rad

280
Procedure To Restart A RADIOSS Run

281
282
18 Checks During & After
Analysis
Numerical validation includes the following steps which will be discussed below:
a. Starter output file diagnostics
b. Engine output file diagnostics
c. Review of simulation energies
d. Typical debugging scenarios
In general recall:
• Many modeling errors can be checked in the pre-processor
• The RADIOSS Starter will check the model completely
• All errors need to be corrected for the simulation to run
• All warnings must be reviewed to ensure a valid simulation
• Review penetrations
• Fix real incompatible kinematic conditions
• Check time step (element and nodal)
• Check mass and center of gravity calculated by starter

A Review The Diagnostics In Starter Output File: runname_0000.out


Check Mass And Inertia

• Total mass and inertia should conform to expected values


• Mass of the parts (if necessary)

283
Review Element And Nodal Time
Minimum time step value should conform to what is expected
Otherwise check element shape, material and/or property

Examine any Incompatible Kinematic Condition (IKC) warnings


The RADIOSS starter does not check if the IKCs are really incompatible
For example, if they are strictly orthogonal or if they are not applied simultaneously,
ignore the warning
If the IKC involves nodes that are slave to a rigid wall which are also part of a Type 2
interface, the nodes will be automatically removed from the rigid wall

Initial Run Checks For RADIOSS

284
B Review Diagnostics In Engine Output File: runname_0001out
Once the RADIOSS starter file check is completed without errors, the model can be
computed and the Engine output file should be reviewed…

Example Engine File Output:


In this model there is mass scaling & only gravity loading (no initial energy)

Time Step: Notice at t=0 the natural time step is reported, then at cycle=100, the time
step shown is that set on /DT/NODA/CST in engine file (0.100E-05). Note: If the time
step decreases and then it goes up quickly, there is not a problem. If it varies greatly
from one cycle to another, it is maybe due to the interface stiffness. If the time step
remains low, a problem has occurred. In this case, you have to find the node (or
element) controlling the time step and try to understand why the decrease occurred.
Mass Error: Non-zero value printed because of mass added to meet specified time
step
Energy Error: In the first cycle of a run for which the initial energy of the system is null,
a false energy balance introduces large relative error which can be ignored

Check at the end of the computation in the engine output file:


• Time step
• Energy error

285
• Mass error

“NORMAL TERMINATION” message should be printed.

C Energy Error

• If the error is negative, it means that some energy has been dissipated.
• In case of under integrated elements (Belytschko shells, solids with 1 integration
point), the Hourglass Energy can explain a negative Energy Error since it is not counted
in the energy balance. The normal amount of Hourglass energy is about 10% to 15%.
• If the error is positive, there is an energy creation.
• In case of using QEPH shell formulation or fully integrated elements, the Energy Error
can be slightly positive since there is no Hourglass energy and the computation is
much more accurate. An error of 1% or 2% will be acceptable.
• In case of a larger positive energy, the source of this energy has to be identified.
Incompatible kinematic conditions can lead to such a situation.
• If the error reaches ±99%, the computation has diverged; except in the first cycle of a
run for which the initial energy of the system was null, so that a relatively false energy
balance often introduces large relative error.
• It is recommended that

• If using mass scaling to control time-step, the variation should be smaller than 1%

286
(Above 1%, check whether the nodes with added mass are moving), use
/ANIM/NODA/ DMAS, /ANIM/MASS
• Check Energies of parts critical to your analysis (Specify critical parts in Time History)
/TH/PART, /TH/SUBS
• Review animations closely along with stress and strain contours
• If simulating a quasi-static loading then, ideally:

Check Simulation’s Energy Balance

Energy balance (from time history file - T01)


• Plot all available energies on the same graph to evaluate globally the results

In Case Of Imposed Motion Or Loading

• Energy is added during the computation --> External work.


• The external work is transformed to kinematic and internal energy.
• Total energy will follow the external works
• Numerical energies should be as small as possible.
Remark : Due to its formulation, the energy error in the engine listing file can show
±99% if there is no initial energy in the model.

287
If The Energy Is Initially Set In The Model:
• Translational Total Energy (TTE) should be constant or slightly decreasing (loss of
energy)
• TTE = IE + KE + RKE + CE + HE
• If energy absorption is the goal then most of the kinematic energy should be
transformed to internal energy.
• Numerical energies should be as small as possible:
• Contact energy, CE
• Hourglass energy, HE
• Spring energy
• Rotational kinematic energy, RKE (except rotary machine application)

288
What If Contact Energy Is High?
Normal value for contact energy:
• 0% to 5% of the total energy
Contact internal = penetration value * contact force
High contact energy in case:
• Too soft contact stiffness
• Too stiff contact stiffness
To avoid high contact energy:
• Set physical material parameters (young modulus, thickness)
• Set recommended parameter for the contact definition
• Set “physical” contact interface gap value.

What About Hourglass Energy?


Normal value for hourglass energy:
• 10% to 15% of the total energy when using under integrated elements (Belytschko
shells, solids with 1 integration point)
• Hourglass Energy is not counted in Energy Error, so negative energy error can occur
High hourglass energy in case :
• Under integrated elements without optimized parameters.
To avoid high hourglass energy:
• Set recommended parameter for the property definition
• Fully integrated elements or
• HEPH for solid elements
• QEPH for shell elements (HEPH & QEPH can result in slightly positive energy error, 1
to 2%)
Some more “in depth” information:

289
Energy Balance

Taking into account the external works, the total energy must remain constant or
decrease slightly. The total energy can increase at the end of the computation, during
the spring back or at the beginning during the first cycles.
Internal energy + Kinematic energy + Hourglass energy + Contact Energy + … =
Variation of the External Work
If under-integrated elements are used, the total hourglass energy must remain lower
than 10% of the total energy. If this is not the case, the mesh should be reworked or
elements with physical stabilization method should be used. The contact energy is not
really physical. For each subset and for each part the following limitation is
recommended:

Energy error should be 0.0 <Energy Error< -5.0%

The energy error is computed by RADIOSS as follows:


Ek: Translational kinetic energy
Ekr: Rotational kinetic energy
Ei: Internal energy
Ek1: Initial translational kinetic energy
Ekr1: Initial rotational kinetic energy
Ei1: Initial internal energy
Ewk: External work
Ewk1:Initial external work
The Energy Error computed by RADIOSS is a percentage. If the error is negative, it
means that some energy has been dissipated.
In case of under integrated elements (Belytschko shells, solids with 1 integration
point), the Hourglass energy can also explain a negative Energy Error since it is not

290
counted in the energy balance. The normal amount of Hourglass energy is about
10% to 15%.
If the error is positive, there is an energy creation. In case of using QEPH shell
formulation or fully integrated elements, the Energy Error can be slightly positive
since there is no Hourglass energy and the computation is much more accurate. An
error of 1% or 2% will be acceptable.
In case of a more important positive energy, the source of this energy has to be
identified. Incompatible kinematic conditions can lead to such a situation.
If the error reaches +/-99%, the computation has diverged; except in the first cycle
of a run for which the initial energy of the system was null, so that a relatively false
energy balance often introduces large relative error.

Mass Balance

If the mass increases, its variation must remain smaller than 1% for each subset and
for each part (dM/M < 1%). If the mass variation is between 1% to 3%, you have to
check if the nodes with the added mass are moving or not. If this is the case, the added
mass results in an increase in kinetic energy. For more than 3% of variation, the results
are probably bad.
• The added mass can be due to Interface Type 2, Spotflag =1. In this case, the added
mass is totally made at time t=0.
• In case of added mass in the model, it is necessary to check if it is not too important
with respect to the total mass of the model (see the DM/M value in the last column
in the RADIOSS Engine listing file (Runname_nnnn.out)).
• It is also important to post-process this added mass in order to check that it is not
too large locally, since this could mean false results (for checking this, the
Animations written with /ANIM/NODA/DMAS have to be visualized).

291
Momentum Balance

The dynamic equilibrium of each node is satisfied by the Newton law at the end of
each cycle. As RADIOSS resolves the equilibrium equations at each cycle, normally the
momentum balance is satisfied. However in case of a problem, a cross-check between
nodal accelerations and the impactor forces (interface, rigid wall, barrier, etc.) can
help to better understand problems.

Visual inspection of Deformed Parts

After getting some animation files, the overall deformation of the structure can be
compared to the physical behavior. You can inspect the deformed shape to ensure:
• There are no flying nodes (parts).
• The deformed shape is smooth.
• The chord angle between adjacent elements is sufficiently small.
• There are no intersections.
Good physical behavior is obtained when:

• The plastic strain is less than 30%,•

292
Meaningless Results

Contact forces, von Mises stress, nodal velocities and accelerations have to be
checked carefully. If the values are meaningless (for example: von Mises stress = 1
GPa), the first check may be the unit system consistency. Refer to the Unit Systems
conversion table below for more information.
Basic Relations

293
Unit System

294
D Typical Debugging Scenarios
Energy Error

• Check the last animation file (written before the computation stops with option /STOP
in the engine file).
Usual issue :
• Imposed time step too high
• Material badly defined
• Elastic material with large deformation
• Contact interface (stiffness, gap) not correctly defined

Mass Error

• Check the last animation file (written before the computation stops with option /STOP
in the engine file).
Usual issue :
• Imposed time step too high
• Material and property (spring) badly defined
• Elastic material with large deformation
• Contact interface (stiffness, gap) not correctly defined
Negative Volume Usual issue :
• The material is not correctly defined for large deformation (densification phase)
Countermeasure
• Improve the material definition
• Add a self impact interface in the solid block with a gap value corresponding to the
maximum compression value of the solid element.
• Add failure in the material model

295
Time Step is Negative

• Issue in the tied interface (/INTER/TYPE2). The slave node is not correctly projected
to the master surface
• Improve the projection or the interface master and slave group definition
• Set the flag Spotflag=1 in the contact interface TYPE2
In the following video series Rahul Rajan describes some typical debugging techniques

Starter Crash (14 minutes)

Engine Crash Before First Cycle (5 minutes)

296
Engine Crash After Several Cycles (8 minutes)

Troubles With Energy and Time Step (7 minutes)

297
Troubles With Material Definition (4 minutes)

Troubles with Interfaces (7 minutes)

298
Results Checking
How Is The Energy Error In The Listing File (Runname_0001Out) Computed?

Ek is the translational kinetic energy at current time


Ekr is the rotational kinetic energy at current time
Ei is the internal energy
Ewk is the work of external forces (energy brought to the system)
E,1 is the energy at beginning of the RUN (not at time t=0)
The Hourglass energy is not counted in this energy balance, so that a negative energy
error generally occurs (except if using QEPH or BATOZ shells, and fully integrated
solids, for which there is no Hourglass).
It is bounded to symbol_plus-minus99%. The energy error is reset after each RESTART.
It is possible to stop a job and rerun it by using a control file containing /CHKPT. In this
case, error and energy values will restart from their last value.

How To Appreciate The Energy Error, And Which Values Are Reasonable Values?
The Energy Error computed by RADIOSS is a percentage.
If the error is negative, it means that some energy has been dissipated.
In case of under integrated elements (Belytschko shells, solids with 1 integration
point), the Hourglass energy can also explain a negative Energy Error since it is not
counted in the energy balance. The normal amount of Hourglass energy is about 10%
to 15%.
If the error is positive, there is an energy creation.
In case of using QEPH shell formulation or fully integrated elements, the Energy Error
can be slightly positive since there is no Hourglass energy and the computation is
much more accurate. An error of 1% or 2% will be acceptable.
In case of a more important positive energy, the source of this energy has to be
identified. Incompatible kinematic conditions can lead to such a situation.

299
If the error reaches symbol_plus-minus99%, the computation has diverged; except in
the first cycle of a run for which the initial energy of the system was null, so that a
relatively false energy balance often introduces large relative error.

300
19 Best Practice
Meshing Recommendations

Warping
• Warping = h/L
• Avoid zones with majority of warped elements
• Warped elements introduce weak spots in the structure and may lead to Buckling at
unphysical locations

Aspect Ratio
• Aspect ratio = L/l
• Bad aspect ratio can generate membrane instability with nodal time step control

Number Of Elements Per Wave Length


• Minimum recommended size for shell elements = 5 mm

• Average size for shell elements = 10 mm

301
Minimum of 6 elements per wave length in the energy absorbing parts, preferably 10.
If less, you may not present buckling.

Element Size to Target a Time Step of 1 Microsecond. If the


element is too small time step can be low

Percentage Of Triangles
• Triangular elements have no hourglass modes
• Membrane stiffness is too high
• Recommended percentage < 10%
• Used for : mesh transition, replace highly warped 4 node shells

Number Of Elements Per Side Length


• On flanges, parts with holes
• Minimum of 3 elements per side of any section
• If less, bad behavior in bending

Note The Different Numbers Of Elements

302
Meshing On Neutral Plane
• Keep the true geometry of the part (section, inertia …)
• Keep the true distance between the different parts (contacts)
Keep All Flanges
• Need for welding and keep the true geometry

Suppress The Small Details


• Delete holes smaller than the reference mesh size
• Remove the shape details which are not necessary
Details deleted and number of elements

303
Flanges Transitions

Quality Model

• Follow the meshing recommendations


• Average mesh size: 10mm
Shell
• Q4 Belytchko & Tsai (BT)
• 5 IP through the thickness
• Ishell = 4 (improved hourglass formulation)
• Ismstr = 0 (default value)
• Ithick = 1 (thickness changes)
• Iplast = 1 (Iterative plasticity)

304
Solid
• Isolid = 0 (default with 1 IP)
• Ismstr = 0 (default value)

1 – Pay Attention To The Mesh Size (For Time Step Issues)

2 – Simplify The Model (Geometry, Symmetry)


Example: Remove or simplify details (fillet, edges..) which are << than mesh size

3 – Adjust Your Meshing Style As A Function Of Part Dimension


if mesh size > 4* thickness then use shell mesh
We recommend using QEPH elements for shells on the critical parts with the following
options
• IPLAST set to 1
• N = 5 (number of integration point through the thickness)
if mesh size < 4* thickness then use solid mesh If you have:
• 1 element through the thickness (not recommended) then use fully integrated solids
(HA8) or thick shells
• more than 1 element through the thickness then use HEPH

305
4 – Choose A Unit System And Be Consistent With It
Most commonly used unit system :
[mm][ms][kg] => [Gpa][KN] | g=0.00981 mm/ms2
[mm][ms][g] => [Mpa][N] | g=0.00981 mm/ms2
[mm][s][Mg] => [Mpa][N] | g=9810 mm/s2

5 – Always Check The Model Before Running It (Listing File) Here Is List Of Check
That You Should Be Perform:
• Check the mass and center of gravity
• Time Step of different entities
• Check for free parts | the connectivity of the model
• Intersections / Penetrations
• Use ModelChecker in HyperCrash if possible
• Common mistakes :
• loose nodes of spring,
• incompatible kinematics conditions,
• nodes in multiple rigid bodies,
• master node of rigid body on elements,
• null thickness…

6 – Always Check Your Energy Balance And Added Mass PercentageFor common
structural problems (frontal crash, Impact,…):
• Added mass < 2% ; Energy error negative and < 5% > Watch for sudden changes in
time step and energy spikes

7 – Energy Balance In Initial Velocity Problems:


• Energy Balance is not respected in Imposed Conditions

306
8 – Check Energies Of Parts Critical To Your Analysis
• Specify parts critical to your analysis in Time History

9 – Review Animations Closely


• Review Stress and Strain contours and perhaps other contours relevant to your
analysis
• Some CAE Engineers have a bad habit of just reviewing the TH curves and making a
decision or changing the model – Not Recommended

10 – Narrowing Technique Helps


• If you cannot identify problem with input deck, keep deleting entities to rule out
possibilities

11 – Track Changes And Document Results


• You could loose track of changes very easily – Recommended to document in the
Header Card of the RADIOSS Deck
• Programs like Xdiff (linux/unix) and Compare (Windows) sometimes help

307
308
20 Postprocessing -
Crashworthiness
This chapter will cover the necessary tools for post processing a crash analysis
(transient event).
The topics covered include:
· Measures – Distance Between and Position
· Section cuts · Vector plots
· Tracking
· Tracing
· Exploded views
· Synchronizing data
· Video/Image overlay
· Export deformed shape
· Crash Tools for Plotting
· Collision detection

20.1 Measures – Distance Between And Position


This section will focus on how to use the Measure panel to measure the Distance
Between two entities and also to measure an entities Position.

309
Measure Panel. Click on “Add” to add a new measure entity. Alternatively, right
mouse click in the Result Browser.

Distance Between
By setting the Measure Type to “Distance Between”, we can measure the distance
between two entities over a period of time.
You can measure the distance between either two Nodes or two Systems. Create a
measure using Results browser > Model View + Entity Editor.
In the example below:
Standard > Entity IDs > Nodes, the right front tire and the right rear tire nodes are
selected.

Then the Standard >Items > Select … > Create Curves option was selected.

This opened a new dialog box to set Axis and Live Link option > OK and a curve was
created.

310
This curve has the magnitude as the Y axis and time as the X axis.
Once the model is animated, the curve is created and plotted:

Position
Setting the Measure Type to “Position” allows you to measure the X, Y and Z position
of a Node or System.

For example, the position of 1 or multiple nodes over a crash event


Set the Measure Group type to “Position” and select 1 or multiple nodes.

311
Using the same model and nodes as in the image above, a Measure Group was created
at each node. Notice how in the Graphics Area the X, Y, and Z positions are reported.
Curves are created for each Measure Group as well using the X position as the Y axis
and time as the X axis:

312
20.2 Section Cuts
The Section Cut panel allows you to cut planar or deformable sections through a
model. This is useful if you want to see details inside a model.
Create a section cut using Results browser > Model View + Entity Editor.
To access this panel, click the Section Cut icon , on the Display toolbar to display
the Section Cut panel.
You can also use the Results Browser > Model View > Create > Section Cut to set the
section cut in the Entity Editor.

313
The section cuts defined for the current model are listed in the Section Cuts list.

You can add and delete section cuts in the list using the Create and Delete. The
browser allows you to display one or more section cuts.
To create a section cut, you must first define the plane of the section that you want
to view.
A plane is defined by using one of the following methods Standard > Type > Section
Plane dialog box:
• X, Y, or Z Axis: Cuts the model along the designated axis. The base node is optional if
you want to specify a certain location to cut the model.
• Vectors: Uses a predefined vector.
• N1 N2 N3: Defines the normal direction of the cutting plane when nodes are picked
on the model for N1, N2, and N3 using the mouse.
• Normal to screen: Specify a plane normal to the screen by picking a base node on the
model.
• Base: Allows you to check the current position of the section cut. Double-click on Base
to enter a Base Node ID or enter user-defined X, Y, and Z coordinates.

Once a plane is defined, Standard > Position and Manipulator (drag it with your left
mouse button) in the graphic area can be used to examine your model by controlling
the position of the section cut.

314
You can also choose how you want the section cut to be applied to the model. By
default, it displays a section cut of the model when the model passes through the
plane.
Display > Deformable allows the sections cut to be deformed, along with the
components, during animation.

There are also various display options available for a section cut. The section cut can
be displayed as:

315
Cross Section Only

Cross Section Width Widths

Clip Elements

316
This will clip the elements that pass through the section cut so that a smooth section
cut is shown. Unchecking this option will not clip the elements so that the whole
element is shown.
The Section color option allows you to assign a color to the cross section. The Grid
option simply adds a gridline to the section cut.
There are also options for displaying the clipped geometry as:

Feature Lines

Transparency Mode

Cross-Section Measurements And Path Plots


Measures can be defined on the nodes of a cross-section enabling a deeper analysis
of the deformation of the section.
Typical applications of measuring values on the cross-section include plotting the
position or distance between nodes on the cross-section (as measuring the intrusion
shown in the image below).

317
A distance between measures on the opposing sides of the vehicle B-pillar section
shows how close the structure is intruding into the occupant compartment.
The measurement on a cross-section can be helpful in identifying the minimum
distance separating the two sides that may not be possible by considering only the
mesh nodes.

In addition to making position based measurements on the nodes of a cross section,


you can also graph the section or any portion of the continuous segments by defining
a node path on the cross-section cuts and cross plotting node X, Y, or Z locations (or
the distance between them).
Result variations along a section cut can also be studied by defining a node path along
the cut. This helps engineers to assess the variation of a result along a user specified
path on a cut section.
For stamping, it is often necessary to validate simulation against test data.
One of the common ways to obtain test data is to cut the stamped part into two pieces
and then measure the thickness along the length of the cut.
To be able to compare FEA results against such data, a path needs to be defined along
the cut section of the FEA model and must be able to plot a selected FEA results type
as illustrated below:

318
The nodes on the section cuts are also accessible from the Notes panel to attach any
text based annotations to the cross sections.

20.3 Vector Plots


The Vector panel allows you to create vector plots that can be used to display any
vector data associated to nodes. Examples include displacement, velocity, and
acceleration.
To access the Vector panel:
• Results browser > Result View
• Click the Vector button, , on the Results toolbar
• From the main menu select Results > Plot > Vector

319
The Vector panel works in the same way that the Contour or Tensor panels work.

First you must define the Result Type.


For the Vector panel, all results are nodal vector results, and when available the Layers
can also be specified.
You can also select the Components that should be shown. Once the Result type and
Components are selected, the Selection can be made as Nodes, Elements, or
Components.
If a cross section has been created and is active, there will be a Section selection
available. When entities are not selected for the Selection, the entire model has the
vector plot applied to it.

The Resolved in system also works the same as in the Contour panel.
Once you have created a vector plot, you can change how the vectors are displayed
using the Display tab options.

The Size scaling sets how the vectors are scaled or sized. The three options are Auto,
By Magnitude, and Uniform.
• Auto: Automatically scales the vector length at approximately 2% of the model size.

320
• By Magnitude: The vector size is displayed relative to the value of the vectors. The
Scale value text box allows you to increase or decrease the size according to a scaling
value.
• Uniform: This displays all vectors in a uniform size. You can increase or decrease the
size using the Scale value text box.

The Color by option either shows the vector colored by the Value of the vector or by
the Direction. When Direction is selected, the X component is colored red, the Y
component is colored green, and the Z component is colored blue.

The next Display option sets how the vector is displayed. The Vector heads can either
be set to None or Arrow and the vector head can either be placed at the vector Tail,
Tip, or Center.

321
The options included in the Plot tab.

The available settings are Resultant or Components. This shows the vectors as either
a single resultant vector or the components.
The final option, Show values, allows you to display the vector values in the graphics
area.
The options included in the Legend tab and Result display control options, behave the
same as in the Contour panel.
This includes the Plot Styles and accessing the plot styles from the Results Browser.

Once the Selection and Display options have been selected, there are additional
options available under the Section tab, where you can select options for projecting
the vectors to the cross section plane.
A section cut has been applied to the model: the Selection input collector is set to
Sections and the Resolved in system is set to Global

322
There are two additional options available in the Section tab.

These are Projected and Evenly distributed.


Projected will project the vectors onto the cross section. Below is an example of
velocity vectors before and after the projection:

No Projection Projected to the Cross Section


Evenly distributed will evenly distribute the vectors on the cross section by the specific
Number of rows and Number of columns.
Below is an example of the default distribution of vectors and an evenly distributed
set of vectors.

Default Distribution Even Distribution

323
20.4 Tracking Systems
The Tracking panel allows you to track any entity during animation. When an entity is
tracked, the view remains constant with respect to the selected displacements and
rotations of that entity. This enables you to view the dynamics of a model as the
selected entity remains in the same location.
To access the Tracking System:
· Results browser > Model View > Create > Tracking System

Click the Tracking icon , on the Results toolbar


Right click in the graphic area

Main menu Results > Create > Tracking System

324
HyperView now supports both model based and window based tracking systems.
Tracking systems still belong to a particular model and are still defined using entities
from a particular model, however the mode for a tracking system can now be set to
either model based or window based.
By default, a newly created tracking system will be model based. To change the mode
to window based, activate the Window track option (located on the right side of the
panel).
Tracking Systems
The tracking systems defined for the current model are listed in the Tracking Systems
list or in the Results > Model View browser.
Tracking systems can be added to the Tracking Systems list (using the Add button) and
then defined using a node, plane, component, or system.
Keyboard shortcuts are also available for selecting items within the list:
• Press the SHIFT or CTRL key on the keyboard + the left mouse button to select multiple
items from the list.
• Press the CTRL + A keys on the keyboard to select all items in the list.
• Press the DELETE key on the keyboard to delete the selected tracking system(s).
A tracking system is activated, or deactivated, using the radio button.
In addition, you can access the following context menu options by right-clicking
anywhere within the Tracking Systems list:

Delete - Deletes the selected tracking system(s) from the list.


Rename - Displays the Rename dialog, which allows you to rename the selected
tracking system.

325
Make Current - Allows you to select which tracking system will be activated and
displayed in the graphics window. This option performs the same action as clicking
the radio button “on” next to a tracking system.
You can sort the items in the Tracking Systems list by clicking on the list heading.
If Clear Tracking System is selected, no tracking system will be applied to the model.
Activating the Clear Tracking System option clears the contour, and deactivates all of
the tracking systems listed.
A tracking system is associated with a model. When overlaying two models in the
same window, you need to define a tracking system for each model for model shape
comparison.
Tracking systems from other models in the window will be displayed in the Tracking
Systems list provided those tracking systems are “window” based.
For example, the first tracking system (T1) from the first model (M1) will be displayed
as M1:T1 in another model’s tracking systems list.

Track
The Track options allow you to set the parameters that define a tracking system.

Select an option from the drop-down menu and then define a Node, Plane and Plane
type, Component,
The model is animated with respect to the selected entity. A new node, plane,
component, or coordinate system can be defined and tracked at any time during
animation by picking different nodes or components in the window.
Plane type - This option is activated for Plane tracking systems only!

326
For tracking systems defined by a Plane, the local axis of the tracking systems can be
aligned with any specific global direction by specifying the Plane type (available
options are: OXY,
OYZ, and OZX). This is useful in situations where the global coordinate system
alignment differs from the model coordinate system.
This option allows you to specify the local plane defined by three nodes. The normal
vector of the plane is computed by a cross product of planar vectors formed by N1,
N2, and N3:
• The first node (N1) always defines the Origin of the local coordinate system (the
tracking system).
• The direction from the first node (N1) to the second node (N2) defines the local axis.
Use the Plane type drop-down menu to select one of the following:
OXY - X-axis
OYZ - Y-axis
OZX - Z-axis
• The third node (N3) lies in the plane of interest.
Input Collector - The input collector allows you to select, or change, the entities which
will be used to define the tracking system. You can do this by picking entities directly
from the screen or by clicking the Nodes, Components input collector.

Displacements

Use the check boxes to select which displacements remain fixed within the window.
• Global X - locks the view to the X position of the selected entity as it moves. The entity
appears to remain stationary along the global X direction during animation.

327
• Global Y - locks the view to the Y position of the selected entity as it moves. The entity
appears to remain stationary along the global Y direction during animation.
• Global Z - locks the view to the Z position of the selected entity as it moves. The entity
appears to remain stationary along the global Z direction during animation.
For example, selecting Global X locks the view to the X position of the selected entity
as it moves. The entity appears to remain stationary along the global X direction
during animation. Entities can also be tracked with respect to rotation. This is done by
activating the Lock rotations check box to lock the rotation movement of the selected
entity.

Rotations
Entities can be tracked with respect to rotation. Activate the Lock rotations check box
to lock the rotational movement of the selected entity.

Window Track
Activate this option to change the tracking system mode from “model” based to
“window” based. Window based tracking systems from other models in the window
will be displayed in the current model’s tracking systems list (provided the current
model also contains window based tracking systems).
For example, the first tracking system (T1) from the first model (M1) will be displayed
as M1:T1 in another model’s tracking systems list. The entity definition for other
models cannot be changed, however other properties of the tracking system can be
changed (even though they are not part of the current model).
A tracking system will be removed from another model’s tracking systems list
whenever the Window track option is deactivated (changing the mode from
“window” based back to “model” based).
Once a tracking system is created, the Use Tracking System option will become
available on other panels. For example, this option is available in the Contour panel.
This allows you to create a contour plot with respect to this tracking system. Below is
an example with and without a tracking system in a contour plot:

328
20.5 Tracing
System, component, and node entities can be traced during animation. Line tracing
for metal forming can also be performed by selecting nodes directly on the model or
by importing a trimmed line definition. When an object is traced, an outline of the
object is left behind at certain intervals along the path of motion.
To access the Tracing:
• Click the Tracing icon , on the Annotations toolbar
• Main menu Annotations > Tracing

You can trace nodes, components, systems and line by selecting one of the following
Trace options: Node, Component, System or Line. Select the entity you want to trace
By ID using the input collector or by picking them from the screen.
The selected entities are displayed in the list box below the input collector.
You can remove highlighted entities from the list by clicking on the Delete button.
You can also trace lines by activating the Line option, under Trace, and selecting either
nodes or a trimmed line definition file as the line source.

329
The image below shows a couple of examples.

Once the entities are selected, the Tracing mode is defined.


When From first step is selected, the selected entities are traced starting at the first
step to the current step as the model is animated.
For components, an outline of the selected entities is left behind at each time step.
The traced path is erased when the animation starts over.
If you are tracing nodes or systems, a line is drawn to illustrate the traced entities.
When All steps is selected, all steps of the selected entities are traced.
The final option, Last, restricts the traced path to the last user-specified time steps
before the current step.
This gives the appearance that the path is following the selected entities.
Finally, the Display options are set. This allows you to see the Color and Thickness of
the trace (except when component is selected).
Move with tracking system allows the trace to move along with entities relative to an
active tracking system.

330
20.6 Exploded View
The Exploded View panel allows you to explode a model or selected components for
improved visualization.
You can then continue to select entities and perform post-processing procedures on
the exploded model.
Multiple explosions can be saved as part of a session file (*.mvw file) to be used later
for design review or presentations. To access the Exploded View:
Click the Exploded View icon , on the Visualization toolbar

By default there are two options available in the Explosions listing; Clear Explosion
and Automatic Explosion.
When Clear Explosion is selected, no explosion is applied on the model.
When Automatic Explosion is selected, the explosion is applied to the entire model
using pre- defined settings for center of gravity, direction, and magnitude.
Below is an example of the Automatic Explosion:

The explosions defined for the current model are listed in the Explosions list.
Customized explosions can be added and deleted using the Add and Delete buttons.

331
The radio buttons allow you to select which explosion will currently be displayed in
the graphics window.
You can also select multiple explosion(s) and click the Delete button, in order to delete
the selected explosions from the Explosions list.
To create a custom explosion, the Components to be included in the explosion need
to be selected under Selection.
Then the Mode needs to be selected to determine the origin and orientation of the
explosion:
• Translate allows you to define the direction and distance (measured in model units)
of the model translation.
• Explode from model center sets the origin of the explosion at the model center
• Explode from selection center sets the origin at the center of the selected component.
For both, the orientation of the explosion is defined using the Scale factor.

The Scale factor can either be set to Uniform (all 3 directions translated uniformly) or
Component (each direction can be translated a specified distance). Below is an
example where a single component is exploded:

332
Finally, to translate the components in the explosion, is selected.
This moves the model explosion in (-) and out (+) based on the mode and distance
(measured in model units) or the designated scale factor.
There are two final options; Auto fit and Reset.

When Auto fit is activated, each time an explosion is added or deleted, the view is
adjusted so that the exploded model fits in the window. If this option is not selected,
the view is not adjusted and parts of the model could move out of the window.
Reset simply returns the components to their default locations.

20.7 Synchronizing Data And Windows


The data obtained from a solver is generally on a different time scale than that of
videos captured while prototype testing.
In order to correlate the two, it is essential to bring them both to the same time scale
and then apply the necessary changes so that they can be synchronized and studied
together.
This chapter outlines the procedure on how to synchronize data obtained from a
solver to that of video captured while actual testing.

333
Synchronizing Data
To access the Synchronize Animation:
• Tools menu > Synchronize Animation to display the dialog

• Visualization toolbar > Animation Controls icon > Time Scales… button
• Click the Image Planes icon , on the Annotations toolbar

The Synchronize Window Animation dialog allows you to synchronize the animation,
plot, and video windows on a page by applying appropriate scaling and offset values.
This tool is useful for synchronizing windows that contain animation results in
different units and/or begin at different time steps.
Each window in the current page is listed in the Synchronize Window Animation
dialog.

334
The slider bar represents the global time line for all windows. It reflects all window
time frames regardless of whether they are active or inactive. Each time frame is a
subset of the range represented on the slider bar.
Current time displays the current slider bar setting.

Synchronize Window Animation Slider Bar


There are two methods you can use to synchronize multiple animations:

Using the slider bar


Use the slider bar to specify a time for all active Time A and Time B fields.
When the desired value is displayed in a specific field, click the lock button to set the
value. After all fields have been defined, click Apply to automatically calculate and
synchronize the animations.
When you click Apply, the scale factors for the time values for each window are
automatically calculated. The scale factor is calculated by subtracting Time A from
Time B for each window. The end result is all the frames specified for Time A will
appear at the same time during the animation. Similarly, all frames for Time B will
appear at the same time. All other frames will be scaled accordingly.
In this mode, you cannot enter data in the Scaling and Offset fields. They are used for
reference only.

335
Editing Scaling And Offset
• Scaling - Multiplies the total running time of the animation by the specified amount.
• The animation is then scaled to run within the new time boundaries. An animation
file that runs from 0 to 10 seconds scaled by a factor of 2 runs from 0 to 20 seconds.
Likewise, if the scale factor is 0.5, the animation is adjusted to run from 0 to 5 seconds.
• Offset - Specifies how long to wait before animating the corresponding window.
• When the page reaches the specified time value, animation starts in the delayed
window.
• T1-T0 - Allows you view the time scale before the scale factor is applied (for reference
only).
• Scaled T1-T0 - Allows you view the time scale after the scale factor is applied (for
reference only).

The check boxes in the animate column allow you to select the animations to
be included in the synchronization. You can select all or select none of the animations,
as well as reverse the current selection.

Master Animation
One of the animations listed in the Synchronize Window Animation dialog can be
selected as the master animation with which all other animations will be
synchronized. The time scale assigned to the master frame skips steps that are not
relevant.

336
An activated radio button indicates which window contains the master time frame
and the Master check box designates that the master window should be used to
synchronize all the windows in the list.

Synchronizing Windows
The Synchronize Windows feature allows you to manipulate the view matrix of
multiple models in multiple windows at the same time. To access the Synchronize
Animation:

Start/Set Synchronization
The Synchronize Windows dialog allows you to set the synchronization settings and
begin the synchronization.
The number of available windows on the current page is represented by colored
rectangles (which are arranged in the same page/window layout), and the active
window is also displayed.
By default, all windows are displayed in the current graphics window background
color, indicating that these windows are active for view synchronization. You can click
on each of these rectangles (which represent a page window) to exclude it from the
synchronization. The background color of the window will change to the same color
as the dialog for all windows that are excluded.

337
Synchronized Windows Dialog With Three Excluded Windows
To include a currently excluded window, simply click on the rectangle again (and the
background color will change back to the current graphics window background color).
The Synchronize Windows dialog remains on top until you click OK, or close it. All
activated windows will synchronize the view with the active one. You can now use
keyboard shortcuts and also the 3D view controls to move (rotate, zoom, etc.) the
model. Models in the linked windows will rotate or obey the same view upon applying
view actions. If a selected window is a plot, text, or video window, then it is ignored.
Changes to the window, page, and window layout (for example, adding new windows
or pages) can still be made while the view synchronization is running. If a new page
is added using the add page or cut/copy/paste options, and the synchronized view
mode was “on” prior to the addition, then the new page (and all the windows in that
page) will be linked by default - until the mode is changed or stopped.

Stop Synchronization
Stops the synchronization across all windows/pages.
To stop the view synchronization, click on the active Synchronize Windows icon
on the Page Controls toolbar.
Note: No pop-up dialog is associated with this option.

338
20.8 Overlaying Images And Videos
The Image Plane tab allows you to select and overlay images (such as company logos,
watermarks, etc.) or videos in the graphics area for correlation and presentation
purposes. Static positioning of these images with the model is available, as well as
synchronization of the videos with HyperView model animations.
Multiple images and videos can be loaded into HyperView, and each of these can be
edited in the Image Plane tab.
It is possible to load static images and video files at the same time (AVI, PNG, BMP,
etc.) interchangeably on the image planes.
The image planes can be 2D or 3D:
• A 2D image plane places the media in the graphics area as a foreground or background
to the models that are loaded.
• A 3D image plane embeds the media in the model space so that when the model view
is rotated or panned, the related media will also change correspondingly. This is
particularly useful for correlating simulation with multiple videos or static images that
are captured with various views.

Example Model With 2D And 3D Image Planes

339
To access the Image Planes:
• Click the Image Planes icon , on the Annotations toolbar
• Results Browser > Model View > Create > Image Plane

Image Plane Tab


The width of the entire Image Plane tab can be resized to adjust for the image/video
resolution. To resize the width of the browser, click on the vertical line which
separates the tab from the graphic display area and drag and release the mouse
button when the width is the desired size.
You can save, and later restore, the settings of image planes (position and orientation
with a given model) by saving your work session as a Session File.
To enable comparison of multiple sets of simulation results with each other, or against
various physical test data, Report Templates will accept the image/video files as
parameters for every instance of comparison.
This is useful for comparing the test and simulation results, as well as comparing
across variants of a design.
In addition, you can use the Apply Style To option to apply a 2D or 3D Image Plane to
multiple windows.

340
This has the effect of copying the image planes from the source window to all the
target windows that require the image planes.
A background image like a company logo or a test video positioned in one animation
window can be quickly copied to all relevant animation windows using the Apply Style
option.
H3D Export is also supported for static images.
All image planes that contain background or foreground, as well as 3D planes, will be
embedded in the H3D file for visual display in HyperView Player.
Only static images are supported at this time.

The main Image Plane tab options and sub-tabs that are available are dependent on
the type of image plane selected:
• 2D Image Planes
• 3D Image Planes

341
2D Image Planes
The 2D Image Plane tab is primarily used for displaying static images and videos in the
foreground or background of the graphics window. The image planes can be stacked
on top of each other, or against the model. To accommodate multiple media files in
the graphics to be clearly visible, the 2D image planes can be resized and positioned
in the graphics area. Various image processing filters and transparency settings allow
you to adjust the media visibility to achieve a desired effect.
Select the 2D option from the Type drop-down menu to display the 2D Image Plane
tab.

2D Image Plane Tab

342
A Label is then entered which is used to reference the image plane in the Results
Browser. Labels can also be updated using the Rename option on the Results Browser
context menu. Then the file containing the image is selected in the File field.
Many file types are supported including JPEG, TIF, BMP, AVI, and AMF (please refer to
the online help for a full listing of the supported file types).
Once the file is loaded, the image is shown in both the Image Plane tab > Viewing Area
and the graphics area.

The Crop tool allows you to crop an image so that only a desired portion of the
media can be seen in the graphics window. The content that is displayed in the image
plane is controlled by what is cropped. Click the Crop icon to activate the crop tool
inside the Image Plane tab viewing area, and use your mouse to select the crop area
(a rectangle which defines what part of the image you want to keep). The image in
the graphics area will be updated immediately upon release of the mouse button,
however the image in the Image Plane tab viewing area will remain unchanged. Click
the Crop icon again to deactivate the Crop tool. Right-clicking on the Crop icon resets
the crop tool and draws the full picture in the image plane.
The other tools in the viewing area include the Zoom tool . This tool allows you
to zoom the image in the Image Plane tab viewing area in or out. Left-click on the
Zoom icon to zoom in on the image, or right-click on the image to zoom out. The
current zoom percentage is displayed to the Zoom field. You can also manually enter

343
a zoom percentage directly into the Zoom field. The Zoom tool does not have an effect
on the image in the graphics area. The size of the image in the graphics area is
controlled by the graphical manipulator. To resize the image in the graphics area,
simply left click on one of the corners of the image and drag the corner to resize the
image.
There is also a Time Scales tool. Click this button to display the Synchronize Window
Animation dialog, which allows you to synchronize windows that contain animation
results in different units and/or begin at different time steps. Both model and video
animations will be displayed in this dialog. The Time Scales option is only valid for use
in the Transient animation mode.
Sync Mode - The drop-down menu allows you to adjust the scale used to synchronize
the video time scale with the animation by selecting one of the following options:
• Manual - The default setting in which no changes or adjustments are made. The mode
will automatically return to Manual whenever you manually update the scale.
• Length - Adjusts the scale so that the video starts and finishes at the same time as the
animating model.
• Framerate - Adjusts the scale so that the video and model play at the same speed. In
other words, during animation the index of the video frame will correspond with the
index of the animation frame.
Length and Framerate adjustments can be seen in the row of the image plane file and
the column labeled Scaling in the Synchronize Window Animation dialog. The time
scales for the video are automatically adjusted whenever the Sync Mode is changed
from Manual to Length or Framerate, or when switching between the Length and
Framerate modes. The Manual mode simply retains whatever is the last setting of the
time scales.
To undo any Sync Mode adjustments, you can manually change the scale or click the
Reset button in the Synchronize Window Animation dialog.
Adjusting the scale of any image plane will have an effect on the overall speed of the
animation.
Sync Model - The automatic sync modes of Length and Framerate are dependent on
the animating frames of the model to which they are referenced. If more than one

344
model is present in the animation window, you can choose the reference model from
the Sync Model dropdown menu. The scales are adjusted based on the number of
frames in the selected model to synchronize the video.
Each image plane has a mode set for synchronization and depending on the mode, a
reference model to sync with. No model reference is available when the Sync Mode
is set to Manual. If a model that is referenced by an image plane is deleted, the Sync
Mode will be set back to Manual.
The Sync Mode and the model to which an image plane is referenced, Sync Model, are
rightly propagated when the information is copied to other windows in a session using
the Apply Style or Copy/Paste operations.
At the bottom of the Image Plane tab there are additional sub-tabs available.
For 2D Image Planes, these tools include Anchor, Display Options, Z-Stack, and Filter.

Anchor allows you to select or change the location of the image plane relative to the
graphics area. The options available are listed below:

The Reset option will return the image/video in the graphics area to the state when it
was originally loaded.
In addition to the Anchor tools, the graphical manipulator can be used to reposition
the image/ video. In the lower left corner of the image there is a 2D manipulator
(shown above in the Altair image). Each arrow can be freely dragged along a straight
line. The entire rectangular area will be translated accordingly, while its size and shape

345
will remain the same. The square handle can also be dragged on the plane of the area
to reposition the rectangular area for any in-plane movements.
The Display Options sub-tab allows you to control how the image or video is displayed
on the image plane and in the graphics area. There are four options available.

When Maintain aspect ratio is activated the aspect ratio of the image will be
maintained while the image is being manipulated by the graphical manipulator. When
Stretch is activated, the image will be stretched to fill the graphics area. Backface
visibility is only available for 3D Image Planes. Transparency allows you to control the
transparency of the image plane using the slider bar. This is particularly useful for
applying watermarks onto the graphics area, as well as controlling the visibility of the
stacked images.
The Z-Stack sub-tab is used to control the stacking order of all of the image planes and
the model. The arrows located on the right side of the dialog can be used to move the
items up and down the list. The item at the top of the list will be placed on the
foreground with the other items are placed in the background (based on their order
in the list). The analysis model loaded is represented by <<MODEL>> in the list. Other
image planes are listed with their Label.

The final sub-tab is Filter. This sub-tab contains tools allowing for certain colors can
be filtered out, and the media can be enhanced by applying the appropriate filters
and pixel mapping algorithms. Please refer to the online help for a complete listing of
the available features in this sub-tab.

346
Above, the options have been set so that the color white is removed from the image.
This makes the Altair logo appear to have a transparent background. In order to
remove all of the color white, the Type has been set to ColorAlpha, the Pixel mapping
has been set to Linear2, the Mode has been set to Greater Than, and the Color white
has been selected.
This just leaves the lighter grey surrounding the letters and removes the white from
the image:

3D Image Planes
The 3D Image Plane panel is primarily used for positioning the static images and videos
in the model space.
A 3D image plane embeds the media in the model space so that when the model view
is rotated or panned, the related media will also change correspondingly.
This is particularly useful for correlating simulation with multiple videos or static
images that are captured with various views, as in correlating a crash test video or a
drop test captured from different angles.
Each of the views can be positioned with respect to the model in 3D space such that
the view of the model is tied to an image plane for best correlation.
The 3D image planes that are embedded in model space can be incrementally
adjusted to match the model dimensions.

347
In addition to using the 2 Point and 3 Point alignment modes for scaling the image
planes to the model size, the image planes now allow minor adjustments to the
rectangular plane by the sizing manipulators available on the periphery.
This helps in making small size adjustments without necessarily going through picking
alignment and scaling points again.
Small changes to plane orientations can of course be adjusted by the tripod
manipulator attached to the image planes.
Also, once the scaling and view alignment of the 3D image plane is accomplished, the
plane can be moved in or out of the model to match with any section cut views.
As the 3D image plane is scaled in the model space, the measurements applied on the
model are valid for the media as well, and are thus no longer dependent on the
accuracy of picking points based on the pixel resolution.
This provides a close to accurate measurement in film analysis tasks.
With appropriate image processing filters, a good correlation (or the lack thereof) can
be determined easily with 3D planes.
In addition, an image plane can be tied to a model, so that any tracking system based
deformations are considered on the media as well.
The tracking mode is extended to reference points on the video as well, so that the
image plane’s relative position to the model is maintained.
See the Model tab section (below) for additional details.
Select the 3D option from the Type drop-down menu to display the 3D Image Plane
tab. Many of the options available when for 2D image planes are also available for 3D
image planes.
When the Type is set to 3D, additional options are added.
There is a new sub-tab named Alignment that is added. Here the inputs required for
scaling and aligning of the image plane with respect to the model are entered. The
available options allow you to align the image plane in 3D space with respect to the
models. The alignment uses a simple two or three point input scheme for scaling the
image plane to the model size and orienting it in the appropriate view.

348
3D Image Plane Tab
For example, to align a video file and a model file of a crash event, you might choose
to select the center point of two tires to align the two.

First the Alignment mode would be set to 2 Point, and then the appropriate Plane
would be selected to align the model and video in. Then the two reference points for
the model (N1 and N2) would be selected and then the two reference points for the
video file (IN1 and IN2) would be selected.
Within the Display Options sub-tab, the Backface visibility option is available. When
this option is activated the backside of image plane regardless of the image plane
orientation is always shown. Also, the Cull angle option is added. This option allows
you to adjust the fading of the image (using the slider bar). The combination of
settings within the Transparency, Backface visibility, and Cull angle options allow you

349
to effectively control the manner in which the media is displayed on the image plane
(with model rotations and pan) without obscuring important details.

The final sub-tab added is the Model sub-tab. This allows you to associate an image
plane with a tracking system. This option is only applicable if a tracking system has
already been added to a model. When a model is selected the tracking behavior of
any associated model based tracking systems will be applied to the image plane. This
is useful for studying the deformation motion with the tracking system applied
(provided that there is no out of plane rotation in the video).

Also with the 3D image plane is the 3D graphical manipulator.

The 3D manipulator allows you to make incremental adjustments to the rotating plane
or move the image plane in 3-D space.

350
This manipulator is created based upon the nodes specified using the 2 Point node
input collector (N1N2 and IN1IN2) or the 3 Point node input collector (N1N2N3 and
IN1IN2IN3). Each arrow can be freely dragged along a straight line. The arc handle
connecting arrows can also be grabbed to allow rotation of the arrows and the image
plane.
In addition, the square handles connecting arrows can be freely dragged in space for
planar translations to modify the location of the image plane.

20.9 Exporting A Deformed Shape


The Export Deformed Shape feature allows you to export the deformed shape as
follows:
• Abaqus
• DynaKey
• OptiStruct Analysis/Nastran
• RADIOSS
• STL
These files can be read into a pre-processor, such as HyperMesh, or into a CAD tool to
view the deformed geometry.
To access the Export Deformed Shape:

• Menu bar File > Export > Solver Deck

351
This utility writes out unscaled coordinated only. If scaled coordinates are needed, use
the Linear Superposition utility to create a subcase, and then use this utility to export
the deformed shape (with scaled coordinates).

20.10 Crash Tools For Plotting


This section will cover the Vehicle Safety Tools available in HyperGraph 2D as well as
some of the basic information for creating a bar chart. These are two common plotting
tools that are used for post processing crash results.

Vehicle Safety Tools


The Vehicle Safety Tools (VST) is a collection of over 200 macros, Templex functions,
and external functions for use by a vehicle safety engineer. The toolkit contains many
basic curve operations, current injury assessment mechanisms, and several additional
data transformation features. After opening the VST file, seven new menus that
contain the macros and functions are displayed: Math, Filter, Units, Injury, Misc,
Transforms, and Safety Tools.
To load the Vehicle Safety Tools
• File > Load > Preference File from the menu bar. Then within the Preferences dialog,
select Vehicle Safety Tools and click Load.

352
The VST macros and functions can also be accessed from the Macro list on the Plot
Macros panel in addition to the menu bar.
The parameters associated with the selected macro or function are displayed in the
Parameters list. Before executing some macros and functions, you may need to enter
additional parameter information. When you execute a plot macro or function by
clicking Apply, the macro or function uses the parameters to create curves and notes
within the session. Undo reverses the implementation of the plot macro.
See Plot Macros for more information on working with macros.

Below is a brief description of each menu. For a full description of the menu items,
please refer to the online help.

Math Menu
The vehicle safety tools’ Math menu allows you to perform basic math operations on
one, two, or three curves at a time. The operations include curve addition and
multiplication and resultants, for example. The start, end, or sample frequency of the
various curves is not an issue - the VST automatically crops and samples the input
curves as needed. The curves must be using the same units, in most cases.
The Math menu contains three secondary menus: one curve, two curves, and three
curves.

Filter Menu
The vehicle safety tools’ Filter menu provides curve filtering functions, including
Butterworth, FIR, and a new SAE 1995 that are not available to regular HyperGraph

353
2D users. Filtering requires the time unit of seconds. You must convert your data to
seconds before applying a filter.

Units Menu
The vehicle safety tools’ Units menu allows you to convert to/from common unit sets.
It is often necessary to convert data when comparing disparate sources (for example
test vs CAE). This menu provides a simple and fast way to convert one or many data
curves simultaneously. All conversion constants are displayed for you to review
before applying the conversion.
The Units menu contains secondary menus. From each menu, you can convert various
units of measure.

Injury Menu
The vehicle safety tools’ Injury menu offers calculations for all known injury
assessment tools. Each of the assessments may be used to generate either a note or
a curve at your discretion. The assessments clearly indicate the required units for
input data. If your data does not match the required input, you can use the Units
menu to convert the data or, optionally, supply conversion factors directly to the
input.

Misc Menu
The vehicle safety tools’ Misc menu provides some functions that are not commonly
used. For example, you can simulate a bar chart or histogram using this menu.

Transforms Menu
The vehicle safety tools’ Transforms menu provides several functions to calculate
running statistics, regressions, FFTs, etc.

354
Bar Charts
The Bar Chart toolbar is displayed when you select HyperGraph 2D from the
application menu and the active window is set to the Bar Chart mode, .

Select An Application From The Menu Plot Type Options


Most of the panels available for plotting bar charts are similar to the ones available
for XY plotting, although there are a few differences.

Define Curves
The Define Curves panel allows you to edit existing bar charts individually and
add new bar charts to the current plot. The panel also provides access to the
application’s curve calculator.

355
Define Bars Panel With Values As A Data Source Enabled
The curve list displays the names of all the bar charts in the active window.
The Cut, Copy, Paste, and Add buttons under the bar chart list are used to maintain
the bar chart list for the active window. Bar charts can be copied and pasted into
other bar chart windows or within the current bar chart window. Bar charts can be
also be added and deleted from the list, and they can be renamed.
Individual bar charts can be turned on or off. To turn a bar chart off, select the bar
chart from the window and deactivate the Series check box next to the bar chart
name. To turn a bar chart on that was turned off, activate the Series check box.

Data Sources
Bar charts are comprised of data and categories. Data can be read from an external
file, defined as a mathematical expression, or entered as values. A label identifies
categories and tic marks are used to separate them from other categories. If there are
more data points than there are categories, the application adds blank categories to
accommodate the extra data points.
If Values is selected as the source (as shown above), the Values table is displayed,
allowing you to directly enter data point values.
The bar chart display is divided into Categories. Categories can be renamed by
selecting a cell, typing a new name, and clicking Apply.
You can have more than one bar series in each category. Data point values can be
entered directly into the Curve column in the Values table. Enter numbers in the Curve
columns and click Apply.
After adding a curve to the Curve list, a Curve column is displayed next to the
Categories column.

356
By default, a new curve is named Curve 1, Curve 2 etc. If you change the name of the
bar chart curve in the Curve list, the heading for the Curve column in the Values table
also changes.

Collision Detection
The Collision Detection panel can be used to perform collision interference checking.
It allows you to define a collision set, select components (parts) in the collision set,
and then detect penetration between two groups. The ability to define multiple
collision sets allows you to quickly perform design reviews. To access the Collision
Detection panel:

• Tools toolbar > Collision Detection


• Tools > Collision Detection from the menu bar.

The options within the panel allow you to add collision sets and define the properties
of each set individually.

Collision Sets
The collision sets defined for the current model are listed in the Collision Sets list. You
can add or delete collision sets using the Add and Delete buttons. You can also select
multiple collision sets and click the Delete button, in order to delete the selected
collision sets from the list.
Keyboard shortcuts are also available for selecting items within the list:
Press the SHIFT or CTRL key on the keyboard + the left mouse button to select multiple
items from the list.
Press the CTRL + A keys on the keyboard to select all items in the list.
Press the DELETE key on the keyboard to delete the selected collision set(s). A collision
set is activated, or deactivated, using the radio button.

357
In addition, you can access the following options by right-clicking anywhere within the
Collision Sets list:
• Delete - Deletes the selected collision set(s) from the list.
• Rename - Displays the Rename dialog, which allows you to rename the selected
collision set.
• Make Current - Allows you to select which collision set will be activated and displayed
in the graphics window. This option performs the same action as clicking the radio
button “on” next to a collision set.
You can sort the items in the Collision Sets list by clicking on the list heading.

358
21 Exercise: Bird Strike Analysis
On An Aircraft Leading Edge
(Using SPH Bird Model)

The front facing edges of an airplane, which include the wings, nose cone and engine,
are more prone to bird strikes. In this tutorial, we simulate the bird strike at the
leading edge of an airplane wing structure using smoothed particles hydrodynamics
(SPH) method. The preprocessing of the model is done in Altair HyperMesh 2017 and
is solved in RADIOSS, which is the highly nonlinear explicit solver of Altair.
Pre-processing of the Model:
1. Launch HyperMesh 2017.
2. Choose User Profile RADIOSS 2017
From File>>Open>>Select file Aircraft_Wing.hm

3. Observe the model file loaded. The model has four components which are the parts
of an aircraft wing. The parts are modeled with shell elements.

359
Creating and Updating Materials
1. Please notice that the unit system followed for this analysis is Tons, mm, seconds and
Newton.
2. We are using Johnson-Cook Material (/MAT/LAW2) in this model, with card image
M2_PLAS_JOHNS_ZERIL. The material parameters are provided as shown below.

3. Now we will create four different shell properties for the components. Right click on
model browser > Create > Property. The property name is renamed as skin in the
entity editor. Card image is P1_Shell. The thickness is assigned as 5mm. And also
updated with other parameters as shown below.

360
4. Similarly, for the other three components the above property is duplicated (right Click
skin property, select duplicate and rename) with change in thickness. That is, for ribs-
3mm, spar-1mm, ribs_outer-5mm.
5. The four components are updated with the material created above and with the
respective properties created in the above step.
6. Now we are creating a rigid at one side of the aircraft wing, thereby we are assuming
that this end is fixed on to the aircraft fuselage. A new component collector called
Rigids is created. From 1D > Rigids > Create, now we will create a rigid body. Multiple
nodes selection is invoked and the nodes of the ribs_outer component is selected as
shown below and the rigid is created.

361
Selecting Nodes for Creating Rigids

Wing After Creation Of Rigids


7. Please note that for Rigids component no card image is assigned and no materials or
property is assigned.
8. Now we will create the bird. The bird shape used in this tutorial is a cylindrical model
with hemispherical ends. As per the regulations, for a leading edge bird strike a 4lb
(1.81kg) bird is used.
With this mass, the bird density, length and diameter is calculated (Note: Many open
web resources are available on this calculation). Density: 9.42e-10 tons/mm3
Diameter: 98mm
Length: 189mm

362
Steps for Creating Bird
1. Create a new component collector for Bird. The bird impact location is selected and
two nodes are created as shown for creating the cylindrical surface.

2. The cylindrical body is created with the dimensions as below

3. Now two circles are created at the cylinder ends with 49mm radius. After that, two
spheres of 49mm radius are created at the ends of the cylinder as shown below.

363
4. Next, the sphere is trimmed (Geom>Surface edit) with the lines created therefore
hemispherical ends are created for the cylinder. The inner surface of the sphere is
deleted.

5. Now the bird geometry has been created as per the required dimensions.

6. With the bird component collector active, go to 1d>SPH.


Select comps from the selector and select the bird component created. The pitch
between particles provided is 5 and enter the density value and click Create. A
property collector, sph_1 is also automatically created now.

364
7. Now a material is created for the bird. /MAT/LAW6 material law is used and card
image is MLAW6 and the input is provided as below. Assign the material and sph_1
property to the bird component.

8. The loads and boundary conditions for the model are created now. From Utility Menu,
select BCs Manager. Provide a name for the load collector and select type as Boundary
Condition. From the selector choose nodes and select the master node of the rigid
body created. Fix all the degrees of freedom and click create.

365
9. From BC’s Manager choose Initial Velocity as the type. Parts selection is activated and
select the Bird component. A velocity of 130000mm/s is assigned in X-direction. For
that, the value is entered in Vx field.

10. From Tools> mass calc option, select the Bird component and check the mass and it
shows it is equal to 4lbs, that is 1.8kgs or 1.8 e-3 tons.

366
11. In the next step we will define contacts between the components.
From Analysis>Interfaces option is selected. The interface is renamed as Inter_1 and
type is selected as Type 7. From the Add option, for master, components option is
selected and ribs, ribs_outer and skin components are selected as master. For slave,
nodes option is selected and the bird particles are selected as slave and updated. A
Gapmin value of 10 is provided and Inacti is set to 5.

12. From Analysis Menu > Output Block is selected, from selector comps is selected > all
the parts are selected > then click Create > Edit. For NUM_VARIABLES option 1 is given
and for Var option DEF is provided. It means, all default outputs for this part will be
provided in the time history file.

367
13. Now the engine file is created which contains the run controls. We will start with run
time, /RUN. Right click in Solver browser > Create >Engine Keywords > RUN. This
specifies the run time for this model. We will provide .005 seconds, that is 5
milliseconds.

A good practice is to create 25 - 30 animation files for the simulation. Here we will
create 25 animation files. For that, right click in Solver browser > Create >Engine
Keywords >ANIM> / ANIM/DT is selected.

For frequency of writing time history file TFILE is selected from solver browser and
frequency is given as .00002.

From Solver browser > Create >Engine Keywords >ANIM> ANIM/VERS is selected and
input is given as 44. When the deck includes SPH, one may use /ANIM/VERS/44.

368
Output results like energy density, von mises stress, plastic strain..etc. are requested
now. As shown above, right click in Solver browser > Create >Engine
Keywords>ANIM>ANIM/ELEM is selected and required results are checked.

14. Now all the required output requests and engine control cards are provided. The file
can now be exported into RADIOSS file type. The model file or deck file will be saved
as xxxx_0000.rad, which is also called as starter file. If the merge starter and engine
file option is not checked a separate file xxxx_0001.rad will also be created which is
the engine file containing the engine controls and output requests created.
RADIOSS solver is launched and the deck file is uploaded for running the simulation.
Starter out file will be generated (xxxx_0000.out) which gives the summary of the FE
model created (property, materials, parts, contacts, time step...) and summary of load
cases in the form of warnings and/or error messages. Once the FE model runs fine,
the engine file will run according to control card parameters. The engine file will also
create an engine out (xxxx_0001.out) file. The engine out file gives the cycle summary
and termination notification. The user should continuously monitor the time step in
the engine out file. Check for the abnormal changes in time step.

369
15. Now the run has been completed and the animation sequence is as shown:

16. User can check for the various results from the contour panel.
17. The global energy graph for the simulation is plotted:

370
22 Example: Crash Tube
In the following we look rather “blindly” (because we don’t know much about
RADIOSS) at a ready to solve RADIOSS model. This exercise is meant to give you an
idea on how to carry out a RADIOSS run. Details about parameters / settings are
provided in the subsequent chapters. Thus, we don’t worry about element quality,
size, element type, predefined contacts, type of contact formulation etc. at this
point.
Note: In reality you must check a given model for possible modeling errors first.
Model Description
The model simulates an impact of a thin-walled, closed hat section on a rigid wall.

• Simulation time: [0 – 50] ms


• Uniform thickness = 1.5 mm
• Initial velocity of node 81500 (Master node) Vx= -15.6 m/s
• Input files for this tutorial: FIRST_RUN_0000.rad and FIRST_RUN_0001.rad
• Units: Length (mm), Time (ms), Mass (g), Force (N) and Stress (MPa)
• Results requested:
/ANIM/VECT/DISP Displacements
/ANIM/VECT/VEL Velocities
/ANIM/ELEM/EPSP Plastic Strain
/ANIM/ELEM/VONM VonMises Stress
/ANIM/ELEM/ENER Energies
/ANIM/ELEM/HOURG Hourglass Energy
• Johnson-Cook Elasto-Plastic Material /MAT/LAW2

371
ρ = 7.8e-3 g/mm3 [Rho_I] Initial density
E = 210,000 MPa [E] Young’s modulus

ν= 0.29 [nu] Poisson’s ratio

sy= 180 MPa [a] Yield stress

Κ= 450.0 MPa [b] Plasticity hardening parameter

n = 0.5 [nu] Plasticity hardening exponent

sMAX= 350 MPa [SIG_max] Max Stress

Objectives
1. Become familiar with format of RADIOSS Engine and Starter files
2. Know how to use HyperWorks Solver Run Manager to execute a RADIOSS simulation
3. Analyze starter and engine output files
4. Be able to post-process key results from an explicit analysis with HyperView

Step 1: Scroll through the RADIOSS Engine (FIRST_RUN_0001.rad ) and Starter Files
(FIRST_ RUN_0000.rad) in a text editor and observe Keyword cards denoted with a
“/”.
With a modulus of 210,000 MPa and a density of 0.0078 g/mm 3, calculate the speed
of sound in the material.

Speed of sound in steel rail


=_______________________________________________
The length of the rail is about 1,000 mm. Calculate the time it takes for a shock wave
to travel from one end of the rail to the other.
Time for sound to travel length of rail =
______________________________________

372
Edit the Engine file so that the stress wave can be monitored, moving from one end
of the rail to the other during impact—this will require a termination time equal to
the time it takes for the sound to travel the length of the rail (set on /RUN card).
Set the frequency of animation output to a time that will give 20 animation steps
(/ANIM/DT).

Engine File

Step 2: Open RADIOSS Run Manager

373
Step 3: Input the STARTER file (*0000rad) and Run with –starter option

Notes: Option “-nt 2” denotes that there are 2 CPU threads on the machine for
parallel processing
Option “-starter” tells RADIOSS to only run the starter and stop (if not specified, then
the engine will run automatically if the starter finds no fatal errors)

Step 4: Review the listing file FIRST_RUN_0000out and verify the model:
1. Check if the starter finished with no errors and/or warnings.
2. Find the minimum shell, spring and nodal time steps and write them down.
Shell time step = _____________________
Spring time step = ____________________
Nodal time step = _____________________

Step 5: Run the ENGINE file using the RADIOSS Run Manager

374
Step 6: Review the starter output file FIRST_RUN_0001out with a Text Editor
1. Check the % (Energy) Error and Mass Addition at the end of the computation

Note: In one of the following chapters we discuss “errors” in some detail. The
following information is meant to give you a first idea:
The error must be negative and decreasing (except for the first cycles or in spring back
stage). The error must be less than 15% at the end of computation (5% for a good
model). If the error increases a little, then decreases may be normal. If the error
increases, this means that a problem has occurred. Then, the error may grow to 99.9%
with computation failure message.”
So ideally the error should be at most -15% and better to have -5%. The closer to 0%
error the better. For QEPH shell formulation or fully integrated elements, the Energy
Error can be slightly positive since there is no Hourglass energy and the computation
is much more accurate. An error of 1 or 2% will be acceptable. Otherwise any positive
energy error means energy has been created which indicates a model issue.

Step 7: Review the results file (FIRST_RUNh3d) with HyperView


Open HyperView, load the file FIRST_RUN.h3d, make the Rigid Wall transparent in the
Result Browser and contour the velocities at Time=0 to verify initial conditions.

375
Initial velocity (T=0) 15.6 mm/ms
Contour VonMises stress and step through the animation of the Stress Contour to see
the stress wave propagating from the front of the rail to the rear of the rail. Check the
time it takes the stress wave to reach the back end of the rail. How does it compare
to your calculation?
Note that in the first part of this exercise, the end time of the run was set to the
estimated time the stress wave needs to travel from one end of the rail to the other.

376
To view the individual steps activate the Animation Controls (highlighted icon below).
then browse through the simulation by clicking the “forward” button

Create a copy of the RADIOSS engine file with the same name, but with the _0002.rad
extension. Update the engine file with a new end time of 50 msec and an animation
output of 2.5 msec. Re-run the model with this new engine file.
Notice that the animation files continue numbering from where the first run left off.
Also notice that a second time history file (T02) is created.

377
To view the new results, reload the *.h3d file and repeat the steps from above.
time: 22.7 ms

time: 50 ms

378
Step 8: Review the time history files (T01, T02) with HyperGraph

Open HyperGraph (e.g. Start Menu > HyperGraph) or if you are using HyperWorks
Desktop by splitting the Page in multiple windows (for the course of this exercise we
would need 4 windows)
Plot the Rigid Wall Force as a function of time.
• load the Time History File FIRST_RUN_T01 (note that the second file FIRST_ RUN_T02
will also be loaded).
• Select the Y-type, Y-Request, and Y-Component as Rigid wall > RIGID WALL 1 >
Resultant and click Apply.

Plot energy curves using the Build plot panel:


• select the Global Variables, then Internal Energy, Kinetic Energy, Hourglass Energy,
and TTE-Total Translational Energy
Plot Time Step versus Time
• select the Global Variables, then Time Step and Apply.

379
Step 9: Save the session to a file
Once the post-processing done, it can be saved by writing a session file:
• A *.mvw file reloads exactly your post-processing; Click on File > Save > Session enter
the name Rail_Results and Save
• A *.tpl file reloads your post-processing, with the possibility to modify some
parameters (input file, results file, etc.); Click on File > Save as > Report Template enter
the name Rail_Results and click Save
Note: A comprehensive introduction into post-processing with
HyperView/HyperGraph with a special focus on crashworthiness is provided in the
chapter- Postprocessing Crashworthiness

380
23 Example: Bumper Analysis
In the video below we will show you:
• Model set-up with HyperCrash
• Analysis run with RADIOSS
• Postprocessing with HyperView and HyperGraph We do hope that you’ll find
it helpful.

381
382
24 Example: Honeycomb Crash-
Box
In this tutorial Moritz Guenther (a former student at Altair Engineering GmbH)
explains how to simulate a honeycomb crash-box as quite often used in student race
cars.

We start this project by launching HyperMesh in order to mesh the given geometry.
The geometry is available as a HyperMesh binary file: Model_start.hm

1. Start HyperMesh and select the User Profile RADIOSS: Block130 and open the file
Model_start.hm from file > open.

383
2. Check the Model Browser for details about the loaded model. Notice: We have three
surfaces:
Anti-Intrusion Plate, Barrier, and Bulkhead. The crashbox is a solid geometry.

3. Next, we mesh the crashbox.


Click Mesh > Create > Solid Map Mesh, then select the option: “one volume”; select
the Crashbox (which is modeled as a solid). Change “source shells” to “quads” and set
the elem size to 10mm. Finish your operation with a final click on mesh and return.

384
4. Now we have to mesh the surfaces. Click Mesh > Create > 2D AutoMesh; element size
19mm and quads as mesh type. Choose the surface of the Barrier and start meshing.

385
5. The Anti-Intrusion-Plate needs to be meshed with an element size of 6mm; mesh type:
mixed.
The Bulkhead surface needs to be meshed with an element size of 10mm.
Eventually, the meshed model should like the one displayed:
6. Click on File > Export > Solver Deck, pick the file type RADIOSS and export the model
to your desired folder.

Note: we export the Starter File *_0000.dat

We continue the model set up with HyperCrash.


1. Start HyperCrash, set the working directory to a suitable location, choose RADIOSS
V13, set the Unit system to kN mm ms kg, switch the User Interface Style to “new”,
then hit “Run”.

386
2. Import your RADIOSS mesh and ignore the transformation of the unit system.

Materials
1. Now we have to define two different materials. On the menu panel click on Model
> Material and select Create a new object > Elasto-plastic > Piecewise linear (36).
Type the name EN AW 7075 in and insert the values in the image.

387
We next need to define the function for the material, for this, right click on the field
next to “Yield stress” > “Yield Stress function” and click on “Define Function”. Now
manually insert the data points/pairs of the function from the image below. At the
end click on “Save”

Now assign this material to the relevant components in the model, namely Anti-
IntrusionPlate, Bulkhead and Barrier. Hereto, switch to the “Model Tree”

and pick the respective components.

388
Next, go back to the “Material” browser, activate the “Support” option (at the bottom
of the panel) and click on “Add selected parts of Tree. Don’t forget to save the
changes.

2. Now we repeat the whole procedure of step 3 defining a Honeycomb structure


for the crashbox. Go to Create a new object > Other > (28) and copy the entries
from the image below.
Notice there are six different functions.
Important: The Initial Density is 9.13E-8!

You can define the yield- and shear stress functions by rightclicking in the respective
field and selecting “Define Function…” (as shown above). The function window opens

389
where you can define a function by entering the coordinates (as before) or import an
existing curve.
Note: In this tutorial all functions are imported.

Yield Stress Function 11

Within the function window click on Import curve and choose the file HoneycombT.txt
from your working directory. The function should look like below.

Also provide a function name. In case you need to edit the data later, just make use
of the option “Edit Function”

Define the other five functions by importing the following curves:


Yield stress function 22 HoneycombLW.txt
Yield stress function 33 HoneycombLW.txt

Shear yield stress function 12 HoneycombF3.txt

Shear yield stress function 23 HoneycombF4.txt

Shear yield stress function 31 HoneycombF3.txt


Now to assign this material to the crashbox component, first pick the crashbox
component from the model tree, go back to the material browser, click on add
selected parts from tree, finish with save (same working step as with the material EN
AW 7075).

390
To review whether the material was correctly assigned, you may go to the Material
Browser, select the material of interest, and then hit the “glasses icon”. The
component the material was/is assigned to will be displayed. In parallel, you can
access /edit the data inside the material definition.

Properties
1. To define the properties follow the path Model > Property and select Creates a
new object

Use Include picked parts to choose the Anti-Intrusion-Plate to apply the property
settings from below to that component.

391
Don’t forget to save the changes.
For safety sake, check whether the property was correctly assigned: Clicking on the
glass icon should automatically display the Anti-Intrusion-Plate in the model window.

Some General Questions And Explanations About


The Before Used Settings
Why QEPH Shell Formulation?
• QEPH (Quadrilateral ElastoPlastic Physical Hourglass Control) element
• QEPH is the best compromise between cost and quality. Generally, it costs no more
than 15% of a BT element and the results obtained by this element are close to those
of QBAT.
• With one-point integration formulation, if the non-constant part follows exactly the
state of constant part for the case of elasto-plastic calculation, the plasticity will be
under-estimated due to the fact that the constant equivalent stress is often the
smallest one in the element and element will be stiffer. Therefore, QEPH, defining a
yield criterion for the non-constant part seems to be a good ideal to overcome this
drawback.
• QEPH shells are more accurate for elastic or elasto-plastic loads, whatever the loading
type - quasi-static or dynamic
• QEPH shells will give better results if the mesh is fine enough. It is not recommended
for coarse mesh.

Why Use 5 Integration Points In The Thickness Of Shell Elements?


• In case of an elastic behavior, one gets the exact solution from three integration points
– that is to say that the bending moments are exactly integrated through the thickness
of the shell – and it is not necessary to use more integration points.

392
• In case of a plastic behavior, the bending moments are not integrated exactly. Using
more integration points, the solution becomes more accurate; so it is recommended
to use five integration points.

Why Ithick =1 ?
• By default, shell thickness is supposed constant during shell deformation. Initial
thickness is used to compute strains and to integrate stresses, but the thickness
variation is still computed for post-processing reasons. If a variable thickness (Ithick
=1) is used, true thickness is computed not only for post-processing, but also for strain
computation and stresses integration

Why Iplas =1?


• In plasticity computation, two fundamental assumptions must be satisfied. First, the
stress in the plastic region must verify the plasticity criteria (for example, von Mises
criteria). Second, in the principal stress space the direction (Ds), due to work
hardening, is normal to the yield surface.
• An iterative orthogonal projection is used in this tutorial, this method enables to
satisfy the second assumption using an iterative algorithm, the main reason we use
iterative projection is to achieve a high accuracy on stress distribution after
projection.

Let’s continue with the model building - Definition of Properties of Barrier and
Bulkhead
2. For the components Barrier and Bulkhead, similar properties are needed. Follow
the same steps as in (1) but with the entries as shown in the image below,

393
3. Quite obviously, the Crashbox needs a different property as it is a solid
component, therefore use Volume > Orthotropic solid (6) and these settings:

After selecting the Crashbox (in the Tree), save to finish creating the new property.

394
Contact Definition
1. Click on LoadCase > Contact Interface and pick the Multi usage (Type 7). At first we
define the contact between Crashbox and Barrier. Make sure the details are as below,

After clicking selecting “Include picked parts”, pick first the Crashbox and then the
Barrier.

Approve with Save.


2. Do the same procedure for the contact between Bulkhead and Anti-Intrusion-Plate:
Again, recall that you can review the contact definition by clicking on the glasses icon
which will not only list the contact parameters but also display the components
(elements) of this particular contact definition.
3. The third contact is a Kinematic condition. So choose type 2.

395
To pick the location of the Slave nodes, first turn off the perspective view by pressing
‘P’. Then change the view to the 2D view of the x- and z-axis.

Now choose with Add nodes by box selection the row of nodes of the Crashbox next
to the Anti-Intrusion-Plate and for Master surface pick the Anti-Intrusion- Plate. Finish
with Save.

396
Why Type 7? What Is Interface Gap?
• This interface simulates the most general type of contacts and impacts
• Impacts occur between a master surface and a set of slave nodes
• It is a fast search algorithm without limitations. Interfaces have a gap that determines
when contact between two segments occurs.

What Is Coulomb Friction?


• Type 7 interface allows sliding between contact surfaces. Coulomb friction between
the surfaces is modelled
• Coulomb’s friction law is a classic friction law which states that FT < m* FN with m =
Fric

397
Rigid Body And Boundary Condition
1. Select Barrier in the Tree and go to Mesh Editing > Rigid Body
> Create. Name it and click on Add selected parts of Tree.
Switch to Properties, Type in 305kg at Mass and choose at
center of gravity flag the third option: Mast. Fixed, mass &
inertia on COG.
Finish with Save.

2. Create a new Rigid Body, click on Include picked parts and select the Bulkhead.
Additionally, select every four nodes around the boreholes of the Anti-Intrusion-
Plate using “Add nodes by picking selection”. Confirm this operation with Save.

3. Go to Loadcase > Boundary Condition > Create and select “Add nodes by picking
selection”.

398
Pick the Master
node of the first
Rigid Body.
Please select the
translation to y
and z and every
rotation. Confirm with Save. Now the Barrier
can only move towards the Crashbox.

4. Do the same with the second Rigid Body but select every translation and rotation.

Some additional information about Rigid Bodies:

What Is A Rigid Body?


a. A rigid body is an idealization of a solid body in which deformation is neglected. In
other words, the distance between any two given points of a rigid body remains
constant in time regardless of the external forces exerted on it.
b. Rigid body elements are used to:
i. impose equal displacement to a set of nodes;

399
ii. model rigid connections and pin-joints;
iii. enforce symmetry;
iv. model transitions, connections, spot-welds, seam-welds between non-matching
(dissimilar) meshes;
v. distribute concentrated loads/masses to a set of nodes

Initial Velocity
Go to Loadcase > Initial Velocity and click on “Creates a new object”.
With “Add nodes by picking” selection choose the masternode of the Rigid Body of
the Barrier. Apply a velocity of 7 mm/ms.

Accelerometer
To be able to measure the acceleration you have to fix an accelerometer. For this
purpose, click on Data History > Accelerometer > Create and give it a name.

Then pick the same master node where the velocity is located and save it.

400
Mass Check
To check if you entered correct values in the Material-, Property- and Rigid Body
definition (i.e. density, shell thickness, rigid body mass) click on Mass > Part Details.

Here you see a summary of the model masses. Control the values and the decimal
powers for physical correctness.

Control Cards
Within the control card (Model > Control Card) you define
important and necessary settings like the value of the time
step.
Copy the settings from the image below to your model.
Make sure you save every setting before you move to the
next Control Card.
Start at 0 ms
Write one animation
file per 1 ms

401
Generates animation files containing element data for the selected result type. These
types of data are typically contoured in HyperView.
ELEM - the variable is saved for all types of elements, except where not applicable.

Generates animation files containing nodal masses

Generates animation files containing nodal scalar data. For Nodal data type DT - Nodal
time step is written out

402
Generates animation files containing vectorial data for the specified variable.
Here the variables are velocities, contact forces, forces and moments for rigid bodies,
rigid walls and sections
Provides an estimation of the CPU time spent for each processor

Specify the Printout Frequency


Printout every “Nprint” cycle on file “Runname_nnnn.out. If Nprint is negative, printout
is made in the standard output as well as in the list file

Number of Processors used for this run

403
Run Number

Time when the run stops


Defines the frequency of writing the time history file of T-file, option 4 is used which
can be read by HyperGraph

Start Of Simulation
Now you can export your model. For this go to File > Export > RADIOSS and save your
model.
Then open RADIOSS Solver Run Manager (from the Start Menu), select the Starter File
(location of the file determines the working directory for this simulation), make use
of relevant options like number of CPU cores, .. More options are shown in the image
below.
To check how many CPU’s are available on your laptop/PC - start the Windows Task
Manager.
Here we see there are 4 CPU available.

404
Eventually, start the analysis.
Once the run is completed, we draw our attention to Postprocessing.

Postprocessing
Before we are going to review some results, please note that we explain some basic
principles about HyperView (contouring) and HyperGraph (plotting) in the Appendix
“PostProcessing With HyperWorks Desktop - Crashworthiness”. In addition, please
view the course material available in the Learning & Certification Program
(http://certification.altairuniversity.com/)
Click on the results button to open HyperView, and load the *h3d file from the
working directory.
Note, that the initial distance between Barrier and Bulkhead is 202,2 mm.

405
The displacement plot of the system

The distance between the barrier and the bulkhead is shown in upper right graph. It
starts with its initial value of 202,2 mm. After 55 ms the distance is reduced down to
59 mm. The graph indicates, that the crashbox “relaxes” (jumps back) at 44 ms
(minimum value of 57,1 mm)

In the graph below, the displacement of the Master node of the Barrier is shown.
The “relaxation” is somewhat more obvious in this representation (of course, we
could also zoom in the other graph)

406
The energy balance of the system is shown in the lower right corner in the image
above:
The total energy remains constant (no loss or generation of energy), kinetic energy
decreases while internal energy increases. Hourglass energy and Contact energy are
negligible small. Note that because of the initial velocity condition applied to the
Master node of the Barrier, the maximum kinetic energy of the system is at t=o ms
(no deformation at that time).
The velocity “profile” of the Master node of the Barrier is shown below. Again, the
“jump back” of the system at 44 ms is more apparent to see in the representation.

407
408
25 Capstone Project: Cell Phone
Drop
Objectives
1. Demonstrate ability to set up proper initial conditions for a drop analysis
2. Show appropriate use of interface(s) to model contacting parts
This capstone project is of a cell phone drop. Many parts of the phone are already
connected. Often, this type of analysis would be done with second order tetrahedral
elements. The model supplied is first order to save computation time.
Your task is to review the model for completeness and finish the model setup using
the information below to get a valid simulation.
• Input file for this capstone: phonedrop_start_0000.rad

409
Typical phone model
Front cover (plastic)
Window* with Touch sensor (galss)
Dus gasket (foam)
Buttons (plastic, rubber)
LCD* module (plastic, glass)
Carrier (steel plate w. plastic parts)
PCB* (equivalent elastic material)
Battery (equiv. elastic material)
Frame (plastic)
Screws (solid steel)
Back cover (plastic)
with membrane cover for accurate stress post-postprocessing

Elements
Second order tetras (plastic parts) 110 000
Hex8 full integration (PCB, Window, LCD etc.) 18 000
Hex8 red. integr. (gasket, battery, buttons etc.) 19 000
Shell elements (PCB shield lids) 2 000
Membrane elements (stress cover) 23 000

Already defined: Materials


• Parts Elastic (/MAT/ELAST)

• Properties Elasto-plastic (/MAT/PLAST_TAB)

• Material Hyperelsatic (/MAT/OGDEN)

• Tie interfaces (part-to-part) Interfaces

Tied interface (/INTER/TYPE2)

Contact with friction (/INTER/TYPE24)

410
To define:
1. Contact Loads
 Contact Surface Impact with floor (/RBODY, /BCS)
 Contact Interface Initial Velocity (/INVEL)

2. Loads and Boundary Conditions Gravity (/GRAV)

3. Engine File

Complete Connections And Contact


The screws need to be connected to the front cover. Assume the parts interact with
each other with a friction value of 0.2.

Loading Conditions
The phone drops from a height of 1.5 meters in the z-direction to a flat floor on the Z
plane. The forces between the phone and the floor need to be recovered.
The total simulation time is 1 millisecond. Animation output requests should include
plastic strain for solids and shells, vonMises stress, tensor of shell membrane stress,
and nodal added mass.

Run Analysis And Post-Process


Run the simulation and postprocess the results, comparing to the results on the
following page. Look at the force between the phone and the floor as a function of
time. Compare with the weight of the phone at rest.

411
Expected Results

412
26 Capstone Project: Roof Crush

Objective
Create a simulation of the FMVSS 216 roof-crush test

Initial Models:
• Car
• Impactor

TAURUS_A00.inc FMVSS_ROOF_IMPACTOR.inc

413
Step #1- Import and Create Submodels
• Import the Car.
• Import the impactor as an include file
• Create a //Submodel for the impactor
• Create offsets in the impactor //Submodel that no values are repeated

Step #2 - Transforms
This will require (in order)
• A 180° rotation about the global z-axis
• A 5° rotation about the axis through axis AB
• A 25° rotation about the axis through axis AC
• A translation to put point A at global ( -2145.26, 584.822, 1343.06)

Step #3 - Interfaces
Create a Type 7 contact interface between the impactor and the car
• choose an appropriate stiffness definition
• choose an appropriate minimum gap
• change the coefficient of friction to 0.2
• choose the appropriate friction penalty formulation Save the data from the
interface to the time history file.

414
Step #4 - Gravity
Create a constant gravity load

Step #5 – Impactor Boundary Conditions


The impactor assembly contains a spring attached for stability. Create a BCS collector
to fix the free end of the spring.
Create a moving SKEW to define the direction normal to the impactor’s face.
Create a BCS collector to guide the master node of the impactor rigid body so that it
is free to translate normal to the face of the impactor, but is fixed in all other DOFs.

Step #6 – Imposed Displacement


Impose the velocity of the impactor starting from 0 mm/s at t = 0. The displacement
of the impactor should be 200 mm @ 200 ms.
Avoid abrupt changes in the acceleration of the impactor.

Step #7 – Control Cards


Add
• a /ANALY card
• a /DEF_SHELL card,
• a /DEF_SOLID card, and
• a /IOFLAG card
• a /SPMD card
• a /SPHGLO card
• a /ALE/DONEA card
• a /TITLE card

415
Step #8 - Engine File
Create the engine file for the model:
• Request solver version 13.0
• Use a constant nodal time step
• Print time history every 0.0001 seconds
• Create restart file every 1,000,000 cycles
• Suppress printing of the list file
• Solve for 200 ms
• Create animation file every 0.005 seconds
• Include selected output in animation files:
• Elemental energy
• Elemental equivalent plastic strain
• Elemental hourglass energy
• Elemental von Mises stress
• Nodal added mass
• Request output of CPU time for each processor
• Turn on parallel arithmetic

Step #9 - Solution And Post-Processing


• Run the model
• Plot force vs. displacement. Check that the FMVSS 216 target load of 47,000 N (= 3 *
GVW) has been met.
• Plot the energy vs. time curves.

416
27 Appendix: Recommended
Videos
For your convenience we compiled the introductory videos (by Rahul Rajan)
embedded in the eBook below:

Why & When To Run An Explicit Analysis (14


minutes)

Effects of Hourglassing (7 minutes)

417
RADIOSS File Structure (14 minutes)

1D Elements: BEAM & TRUSS (7 minutes)

418
1D Elements: Spring Type 4 (8 minutes)

1D Elements: Spring Type 12 (4 minutes)

419
1D Elements: Spring Type 13 (22 minutes)

1D Elements: Spring Type 8 (20 minutes)

420
Information Required For Spring Elements (6
minutes)

2D Elements (18 minutes)

421
3D Elements (15 minutes)

Material Law2 and Law36

422
Contact Interfaces (3 minutes)

Tied Contact (12 minutes)

423
How Sliding Contact Works (9 minutes)

Contact Type7 (23 minutes)

424
Rigid Walls (6 minutes)

Rigid Bodies (18 minutes).

425
Added Mass and Skew Frame

Loads & Boundary Conditions (18 minutes)

426
Time Step Control (32 Minutes)

Debugging Techniques - Starter Crash (14 minutes)

427
Debugging - Engine Crash Before First Cycle (5
minutes)

Debugging - Engine Crash After Several Cycles (8


minutes)

428
Troubles With Energy and Time Step (7 minutes)

Troubles With Material Definition (4 minutes)

429
Troubles with Interfaces (7 minutes)

Altair India released a new RADIOSS Introduction


video series recently:

Mesh Criteria, Materials and Properties (by Prashanth A.R.)

430
Interface 2 Connections (by Prashant A.R.)

Creating Loads and Boundary Conditions (by Prashant A.R.)

Contacts Type 7 and Type 11 (by Prashant A.R.)

431
Creating Monitoring Volumes (by Prashant A.R.)

Creating Output Block (by Prashant A.R.)

Model Check and Engine File(by Prashant A.R.)

432
28 Appendix: Conversion Tables
And Constants
The following overview is part of the RADIOSS Theory (see HyperWorks Help
Documentation for more details)

433
434
435
436
437
438

You might also like