A domain decomposition method allowing the combination of nonlinear sub-domains with linear modal
sub-domains has been proposed. 1 With this technique, the displacement field in the linear
sub-domains is projected on a local basis of reduction modes calculated on the detailed
geometry and the kinematic continuity relations are written at the interface in order to
recombine the physical kinematic quantities of reduced sub-domains locally. The method
yields promising save of computing time in industrial applications. However, the use of
modal projection is limited to linear sub-domains. In the case of overall rigid-body motion
with small local vibrations, the geometrical nonlinearity of sub-domains must be taken into
account. Therefore, the projection cannot be used directly even thought the global
displacements may still be described by a small number of unknowns; for example six
variables to express motion of local frame plus a set of modal coordinates in this frame.
This approach is used in the case of implicit framework. 2 In the case of direct integration with an explicit scheme
an efficient approach is presented. 3 One of the main problems is to determine the stability
conditions for the explicit integration scheme when the classical rotation parameters as
Euler angles or spin vectors are used. A new set of parameters, based on the so-called
frame-mass concept is introduced to describe the global rigid body motion.
The position and the orientation of the local frame are given by four points where the
distances between the points are kept constant during the motion. In this way, only the
displacement type DOF is dealt and the equations of motion are derived to satisfy perfectly
the stability conditions. This approach, which was integrated in Radioss V5, will be presented briefly here.
Linear Modal Reduction
A modal reduction basis is defined on one or more sub-domains of the decomposition. The definition of this basis is completely arbitrary. Any combination of eigen modes and static corrections can be used. All these modes are orthogonalized with respect to the finite element mass matrix in order for the projected mass matrix to be diagonal and suitable for an explicit solver.
Considering the case of a structure divided into two sub-domains, assume that modal reduction is
used for linear Sub-domain 1. Thus, the displacement field of this
sub-domain is projected onto the reduction vectors:
(1)
with
vector of discretized displacements in Sub-domain 1,
vector of modal participations. This naturally
yields:
(2)
The number of modal unknowns
chosen is much smaller than the original number of degrees of freedom of Sub-domain 1.
In order to obtain the new coupled system, the dynamic equilibrium of sub-domain 1 must be
projected onto the reduction basis and the velocities involved in the
kinematic relations must be expressed in terms of the modal coordinates.
Thus, write the new matrix system for a single time scale
as:
(3)
Where,
(4)
The structure of this system is strictly identical to that which existed before reduction. Therefore, use exactly the same resolution process and apply the multi-time-step algorithm.
The time step for a reduced sub-domain is deduced from the highest eigen frequency of the projected system in order to preserve the stability of the explicit time integration. This time step is often larger than that given by the Courant condition with the finite element model before reduction.
Model Reduction with Finite Overall
Rotation
Since large rotations are highly nonlinear 4, the displacement field in a sub-domain
undergoing finite rotations cannot be expressed as a linear combination of
constant modes. However, the rigid contribution to the displacement field
creates no strain. In the case of small strains and linear behavior, the
local vibrating system can still be projected onto a basis of local
reduction modes. Then, take into account the rotation matrix from the
initial global coordinate system to the local coordinate system and its time
derivatives. A classical parameterization of this rotation, for example,
Euler angles, would introduce nonlinear terms involving velocities. Since
these quantities, in the central difference scheme, are implicit, this would
require internal iterations in order to solve the equilibrium problem, a
situation you clearly want to reduce the computation time due to the
reduction.
Classically, the displacement field of a rotating and vibrating sub-domain is decomposed into a
finite rigid-body contribution and a small-amplitude vibratory contribution measured in a
local frame. The large rigid motion is represented using the so-called
four-mass approach. 3
Four points
in space are arbitrarily chosen to represent the position of a
local frame attached to the sub-domain. In order to simplify the local equations, choose
these points so that they constitute an ortho-normal frame.
注: The four points
defining the local frame do not have to coincide with nodes of the mesh.
Displacement Field Decomposition
The global displacements of these four points are the unknowns describing the rigid motion of the
sub-domain. These are classical displacement-type parameters. The relative distances between
these points are kept constant during the dynamic calculation through external links. This
enables us to express the total continuous displacement field of the sub-domain
as:
(5)
Where,
-
- Coordinates in the local frame
-
- Rotation matrix expressing the transformation from the local to the global
coordinates: since
are unit vectors,
.
In order to express the dynamic equilibrium,
式 4 must be derived with
respect to time to yield velocities and accelerations.
(6)
The time derivatives of the rotation matrix are given by:
(7)
Thus, the fields in
式 6 have the following
expression:
(8)
Where,
are the components of the local
displacement in the local frame. The assumption of small perturbations in
the local frame enables us to consider that the rigid and the deformed
configurations are the same, that is, you can neglect
compared to the local coordinates,
. Thus, you get simplified expressions of
the velocity and acceleration fields:
(9)
To express the weak form of the dynamic equilibrium, you also need the variation
of the displacement field:
(10)
Where,
Again, the same assumption as above allows us to simplify this expression:
(11)
Local Reduced Dynamic System
The local dynamic equilibrium of the sub-domain is given by:
(12)
The principle of virtual work yields a weak form of this equilibrium, taking into account,
Dirichlet-type boundary conditions:
(13)
Where,
-
- Verifies the kinematic boundary conditions
-
- Volume of the sub-domain
To introduce 式 4
into this weak form of the equilibrium, you must express with the new
parameterization the virtual works of both the internal and external forces
and the virtual work, due to the rigid links among the points defining the
local frame.
The internal forces can be calculated from the local part of the displacement, using 式 5 and taking into
account the rigid links, for example, the fact that displacement
creates no strain.
(14)
Where, index
expresses that the coordinates and
the spatial derivatives are taken in the local frame.
The virtual work of the internal forces is then:
(15)
The integration by parts in the local frame introduces external surface forces
:
(16)
Where,
-
- The boundary of
-
- The normal to
expressed in the local frame
To compute forces associated to the rigid links, first, new Lagrange multipliers are introduced
to express the energy of a link:
(17)
Where,
and
are the initial coordinates of point
and the rigid link between points
and
is given by:
.
Then, the differentiation of this energy is used to obtain the virtual work to be introduced into
the weak form of the equilibrium:
(18)
(19)
注: The quantity
can be viewed as the resisting force applied to
point
to preserve the distances from this point to the other
points of the local frame.
Weak Form of Equilibrium
Now, express the displacement field using
式 5 and the local field
projected on a Ritz basis:
(20)
where
,
,
,
,
is a basis of the global frame,
is a basis of local Ritz vectors obtained,
for example, by finite element discretization or by modal analysis,
is the vector of the discrete
unknowns:
, with
and
,
is the projection basis:
.
式 9
and
式 11
yield:
(21)
Where,
is the gyroscopic contribution to the acceleration, given
by:
The final expression of the complete weak form of the dynamic equilibrium is obtained
as:
(22)
Where,
-
- Vector formed by the 6 relations preserving the relative distances of points
-
- Vector of the Lagrange multipliers corresponding to each rigid link
-
- Vector of the link forces given by 式 19
式 22
can be rewritten using classical matrix and vector operators obtained by
finite element discretization:
(23)
Where,
-
-
- Classical mass matrix of sub-domain
-
- Projection matrix consisting of vectors of
discretized on the nodes of the mesh.(24)
, with
being the sub-domain's local stiffness matrix and
and
deduced (as was
) from
,
and the mesh,
, with
being the classical vector of the external forces assembled on
the sub-domain.
Now, you are able to reduce the number of unknowns on the sub-domain drastically by choosing as
the Ritz vectors, instead of classical finite element shape functions, an
appropriate (and small) family of local reduction vectors. The modal
vibration problem is purely local and guidelines found in the literature for
the proper choice of the projection basis apply here.
注: As far as inertia
coupling with local vibration and overall large motion is concerned,
two separate contributions must be considered. The first contribution
appears in the projected mass matrix, which as now the following
form:
(25)
Where,
is the constant mass matrix corresponding only to the global
displacement field given by
,
is the constant mass matrix corresponding to the local
vibration given by
,
is a coupling matrix, variable with overall rotation, arising
from the interaction between the local vibratory acceleration field expressed in the global
frame
and the overall virtual displacement field
;
naturally comes from the symmetric interaction between virtual local
displacement field and the overall acceleration field.
The second contribution to the inertia coupling is the gyroscopic forces.
注: In
Radioss a special procedure is used to linearize the rigid links.
3
The rigid body motion component of the displacement increment is computed in unconditionally stable way by the use of Lagrange Multiplier to impose the rigid links. The deforming part is generated by the local vibration modes retained in the reduction basis. Therefore, you can conclude that the stability condition is the same as that given by the local vibrating system. The critical time step is constant throughout the calculation and can be derived from the highest eigen frequency of the local reduced stiffness matrix with respect to the local reduced mass matrix.
The highest eigen frequency
is given by the system:
(26)
Where,
and
.
Having determined
, the maximum time step which can be used
on the reduced sub-domain while ensuring the stability of the time
integration is:
(27)