Nonlinear finite element analyses confront users with many choices. An understanding of the fundamental concepts of
nonlinear finite element analysis is necessary if you do not want to use the finite element program as a black box.
The purpose of this manual is to describe the numerical methods included in Radioss.
Kinematic constraints are boundary conditions that are placed on nodal velocities. They are mutually exclusive for each degree
of freedom (DOF), and there can only be one constraint per DOF.
The stability of solution concerns the evolution of a process subjected to small perturbations. A process is considered
to be stable if small perturbations of initial data result in small changes in the solution. The theory of stability
can be applied to a variety of computational problems.
A large variety of materials is used in the structural components and must be modeled in stress analysis problems.
For any kind of these materials a range of constitutive laws is available to describe by a mathematical approach the
behavior of the material.
Explicit scheme is generally used for time integration in Radioss, in which velocities and displacements are obtained by direct integration of nodal accelerations.
The performance criterion in the computation was always an essential point in the architectural conception of Radioss. At first, the program has been largely optimized for the vectored super-calculators like CRAY. Then, a first parallel
version SMP made possible the exploration of shared memory on processors.
A large variety of materials is used in the structural components and must be modeled in stress analysis problems.
For any kind of these materials a range of constitutive laws is available to describe by a mathematical approach the
behavior of the material.
Nonlinear Pseudo-plastic Orthotropic
Solids (LAWS28, 50 and 68)
Conventional Nonlinear
Pseudo-plastic Orthotropic Solids (LAW28 and LAW50)
These laws are generally used to model honeycomb material structures as crushable foams.
The microscopic behavior of this kind of materials can be considered as a system of three
independent orthogonal springs. The nonlinear behavior in orthogonal directions can then be
determined by experimental tests. The behavior curves are injected directly in the
definition of law. Therefore, the physical behavior of the material can be obtained by a
simple law. However, the microscopic elasto-plastic behavior of a material point cannot be
represented by decoupled unidirectional curves. This is the major drawback of the
constitutive laws based on this approach. The cell direction is defined for each element by
a local frame in the orthotropic solid property. If no property set is given, the global
frame is used.
The Hooke matrix defining the relation between the stress and strain tensors is diagonal,
as there is no Poisson's effect:(1)
An isotropic material may be obtained if:(2)
Plasticity may be defined by a volumic strain or strain dependent yield curve (Figure 2). The input yield stress function is
always positive. If the material undergoes plastic deformation, its behavior is always
orthotropic, as all curves are independent to each other.
The failure plastic strain may be input for each direction. If the failure plastic strain
is reached in one direction, the element is deleted. The material law may include strain
rate effects (LAW50) or may not (LAW28).
Cosserat Medium for Nonlinear
Pseudo-plastic Orthotropic Solids (LAW68)
Conventional continuum mechanics approaches cannot incorporate any material component
length scale. However, a number of important length scales as grains, particles, fibers, and
cellular structures must be taken into account in a realistic model of some kinds of
materials. To this end, the study of a microstructure material having translational and
rotational degrees-of-freedom is underlying. The idea of introducing couple stresses in the
continuum modelling of solids is known as Cosserat theory which returns back to the works of
brothers Cosserat in the beginning of 20th century. 1 A recent renewal of Cosserat mechanics is
presented in several works of Forest 2345 A short summary of these publications is
presented in this section.
Cosserat effects can arise only if the material is subjected to non-homogeneous straining
conditions. A Cosserat medium is a continuous collection of particles that behave like rigid
bodies. It is assumed that the transfer of the interaction between two volume elements
through surface element dS occurs not only by means of a traction and shear forces, but also
by moment vector as shown in Figure 3.
Surface forces and couples are then represented by the generally non-symmetrical
force-stress and couple-stress tensors and (units MPA and MPa-m):(3)
The force and couple stress tensors must satisfy the equilibrium of momentums:
(4)
Where,
Volume forces
Volume couples
Mass density
Isotropic rotational inertia
Signature of the perturbation (i,k,l)
In the often used couple-stress, the Cosserat micro-rotation is constrained to follow the
material rotation given by the skew-symmetric part of the deformation
gradient:(5)
The associated torsion-curvature and couple stress tensors are then traceless. If a
Timoshenko beam is regarded as a one-dimensional Cosserat medium, constraint Equation 5 is then the
counterpart of the Euler-Bernoulli conditions.
The resolution of the previous boundary value problem requires constitutive relations
linking the deformation and torsion-curvature tensors to the force- and couple-stresses. In
the case of linear isotropic elasticity, you have:(6)
Where, and are respectively the symmetric and skew-symmetric part of the Cosserat
deformation tensor. Four additional elasticity moduli appear in addition to the classical
Lamé constants.
Cosserat elastoplasticity theory is also well-established. von Mises classical plasticity
can be extended to micropolar continua in a straightforward manner. The yield criterion
depends on both force- and couple-stresses:(7)
Where,
Stress deviator
and
Material constants
Cosserat continuum theory can be applied to several classes of materials with
microstructures as honeycombs, liquid crystals, rocks and granular media, cellular solids
and dislocated crystals.
1Cosserat E. and Cosserat F., “Theory of Deformable Bodies”,
Hermann, Paris, 1909.
2Forest S. and Sab K., “Cosserat overall modeling of heterogeneous
materials”, Mechanics Research Communications, Vol. 25, No. 4, pp.
449-454, 1998.
3Forest S., Cailletaud G. and Sievert R., “A Cosserat Theory for
Elastoviscoplastic Single Crystals at Finite Deformation”, Archives of
Mechanics, Vol. 49, pp. 705-736, 1997.
4Besson J., Bultel F., and Forest S., “Plasticity Cosserat media.
Application to particular composites particles”, 4th Symposium
Calculation of Structures, CSMA/Teksea, Toulouse, pp. 759-764, 1999.
5Forest S., “Cosserat Media”, ed. by K.H.J. Buschow, R.W. Cahn,
M.C. Flemings, B. Ilschner, E.J. Kramer and S. Mahajan, Encyclopedia of
Materials, Science and Technology, Elsevier, 2001.