Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
1/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
Organization (S):
EDF-R & D/AMA
Manual of Reference
R4.03 booklet: Analyze sensitivity
Document: R4.03.04
Calculation of sensitivity in dynamics
Summary
Obtaining the sensitivities is important in dynamics where small disturbances of the parameters can
to have a great influence on the sizes of interest. The objective of the document is to clarify the calculation of
sensitivities (i.e derived) of the main sizes: modes and Eigen frequencies, dynamic responses in
frequential and temporal, into linear and nonlinear). For each size, the “derived problem”,
i.e. the problem whose derivative is solution, is written (in variational and/or matric form).
numerical algorithms of calculation of these derivative are then detailed and discussed. All establishments
corresponding data processing in Code_Aster were not carried out. A table thus gives the current state
developments.
Derivation compared to the parameters of the matrices of mass, stiffness, damping, the loading
and also the particular case of derivation compared to a variation of field is presented.
calculations by the finite element method, which is specific to the derivative, are clarified on the level of
elementary calculations.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
2/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
Contents
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
3/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
1 Introduction
In a mechanical context of calculation on a model of structure, one understands by “sensitivities” them
derived from a size result of the model compared to parameters of model input of the aforesaid. In
dynamics, of small disturbances of the parameters can have a great influence on
dynamic and vibratory sizes [bib1]. The calculation of the sensitivities proves thus primordial, because it
allows to carry out analyzes of influences of the parameters, analyzes probabilistic, of
identification of parameters (or retiming), optimization…
This documentation develops obtaining the sensitivities of the dynamic sizes compared to
mechanical parameters, in the boundary conditions, with the loading, and a variation of field.
The dynamic sizes treated are the dynamic responses in displacement, the frequencies
clean and clean modes. Sensitivities of the other dynamic sizes (forced,
efforts….) which is in general calculated by postprocessing can result some a posteriori by
ad hoc calculations.
The general method of calculation retained is a method of exact direct differentiation (by
opposition to the methods of the associated state, the methods of finished differences or other methods
semi-analytical (the calculation of derived by an semi-analytical method consists in calculating them
derived from the elementary matrices by finished differences and then to calculate them analytically
derived from the sizes results. This method is easier to establish than the analytical method
complete and was retained in Nastran or Abaqus for example (cf [bib2]))).
This method has as main advantage of being exact. One avoids the problems thus of
convergence of the algorithms of optimization based on the gradient or the problems of precise details of
finished differences (when the results result from an iterative algorithm, the finished differences can
to be particularly vague). It is also more general than the method of the assistant state where
a functional calculus to be derived must be specified before calculation (it thus makes it possible for example to obtain
directly fields derived).
This method can be very economic in calculating times if the algorithms of resolution of
derived problems - i.e. the problems whose derivative are solutions - are well chosen and if
the programming is neat. In the most favorable cases, the calculating time of a derivative
represent that a very small percentage of the calculating time of the standard size.
Its main disadvantage is the extent of the data-processing developments necessary to sound
establishment (more especially as the architecture of Code_Aster does not make it possible to use the differentiation
automatic).
The document attempts to present the various “derived problems”. They are written in form
variational and/or matric according to the need. These derived problems often have the very maid
properties. For example, in linear dynamics, the first members of the systems to be solved are
preserved. In nonlinear dynamics, the derived problem is a linear problem.
The methods of resolution which meet a priori more the needs for the users of Code_Aster
are presented. In order to minimize the calculating time, the calculation algorithm of derived must exploit
as much as possible intermediate results of the resolution of the standard problem. In addition,
the effectiveness of certain algorithms depends on the properties of the matrices of the problems to solve, of
a number of parameters per report/ratio for which it is necessary to derive, and of the required precision.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
4/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
2
Derived from dynamic responses
2.1 Problem
standard
In the frequential field, equations of linearized dynamics and the dualisation of the conditions
to the limits lead to the following matric equation (it is the discretization by the method of
finite elements of the variational formulation which makes it possible to obtain this equation. During
document, the variational formulations and with the method of the elements are clarified when
necessary):
(
)
-
+
+
=
2
M
C
K
X
F
()
()
()
(,)
(,)
p
I
p
p
p
p
éq
2.1-1
where:
·
X (,)
p
R
N
(
= +
N
N2
) is, for the pulsation
, the vector made up of N ddl physical
structures discretized and
2
multipliers of Lagrange [R3.03.01],
·
p is a scalar parameter of mechanical property, loading or of geometry which
intervenes directly explicitly in the matrices of mass, stiffness or
of damping, or in the loading. Let us notice that that can be a parameter
of geometry for the elements of structure (for example a thickness of hull) but not
a geometrical parameter of the mesh for which the field of integration is modified and
for which thus it is necessary to return to a variational formulation of the problem and to establish one
parameter setting of the field (cf [§4]),
·
M ()
p
,
C ()
p
and
K ()
p
are the real symmetrical matrices of mass, damping and
of stiffness “generalized” (thus not necessarily definite positive),
·
F (,)
p
is the vector of the external forces.
Code_Aster solves a more general equation for the systems fluid-structures [R5.05.03].
The writing of the derived problem and its resolution are not basically modified by it. For more
of clearness and concision, one thus limits oneself in the continuation of this Doc. to the equation [éq 2.1-1].
In Code_Aster, the resolution of [éq 2.1-1] is carried out by the direct methods multifrontale and
“LDL
T
“[R6.02.02]. These two methods proceed in three stages:
·
renumerotation of the unknown factors,
·
factorization of the matrix,
·
gone up descent/(resolution of two triangular systems).
If several linear systems, of the same matrix, are to be solved, only them
gone up descents/are to be carried out several times. The same if several of the same matrices
structure are to be factorized, the renumerotation of the unknown factors will not be to remake (this phase
preliminary has a considerable cost, even if its relative weight in calculating times decreases with
cut matrices to be factorized). These two properties are very important and will allow
to calculate the derivative at numerical lower cost.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
5/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
2.2
Derived from a harmonic dynamic response
2.2.1 Method
direct
The derivation of [éq 2.1-1] gives (after handling for exhiber a “good” second member):
(
)
-
+
+
=
- -
+
+
2
2
M
C K X
F
M
C
K X
I
p
p
p
I
p
p
()
()
()
éq
2.2.1-1
The resolution of the linear system [éq 2.2.1-1] gives the size
X ()
p
sought.
One can observe that the system [éq 2.2.1-1] has the same matrix coefficient as [éq 2.1-1], and
a second different member but whom one can obtain without inversion starting from elementary calculations.
This second member is called in the literature “pseudo-load” giving to the direct method it
name also of “method of the pseudo-loads”.
There are as many resolutions of [éq 2.2.1-1] that parameters by report/ratio to which one derives. For
to carry out these resolutions by minimizing the calculating time, one uses the properties of the methods
direct multifrontale and “LDL
T
“recalled in [§2.1]. The renumerotation of the unknown factors and
factorization of the matrix are thus carried out only once for all the resolutions.
Also let us notice that the majority of the terms of the second member are null. Indeed, if p is one
nodal force only remainder the term
F ()
p
, if p is a density only remainder the term
2
MR. X
p
()
, etc… This remark subsequently remains valid of Doc., for all them
terms containing of the derivative.
2.2.2 Simple modal method
To accelerate the resolutions of [éq 2.1-1] frequency per frequency, it is usual to project
the operator on a truncated modal basis
[] [
]
,
L
L
=
1
2
,
L
where
I is the i-ième vector
clean (cf [§3.1]) and with
L N
<
[R5.06.01]. The new system to be solved becomes:
[]
(
)
[]
(
)
[]
LT
L
LT
I
-
+
+
=
2
M
C K
F
()
()
éq
2.2.2-1
with the projection of
X ()
on the incomplete modal basis which is written
[]
X
L
L
= ()
.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
6/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
The derivation of [éq 2.2.2-1] gives:
[]
(
)
[]
(
)
()
[]
()
[]
()
[]
[]
[]
(
)
[] []
(
)
[]
()
p
I
I
p
p
p
I
p
p
p
p
I
L
T
L
L
T
L
L
T
L
T
L
T
L
L
T
L
+
+
-
+
+
+
-
+
+
+
-
-
+
=
+
+
-
K
C
M
K
C
M
K
C
M
F
F
K
C
M
2
2
2
2
éq
2.2.2-2
The resolution of [éq 2.2.2-2] provides
()
p
. One obtains then the derivative of projection on the basis
modal of
X ()
, i.e.
X
L
p
by:
[]
[]
X
L
L
L
p
p
p
=
+
()
()
éq
2.2.2-3
The remark made for the direct method is valid for the simple modal method: the pre ones
processing of the matrix coefficient in the linear system [éq 2.2.2-1] (factorization
primarily) can advantageously be re-used for the resolution of [éq 2.2.2-2].
In this method, only truncation carried out is that of the projection of
X ()
on the basis
modal incomplete. The calculation of the derivative is then carried out without another approximation. In
counterpart, the resolution of [éq 2.2.2-2] and the calculation of [éq 2.2.2-3] require calculation as a preliminary
of
[]
L
p
, for example by one of the methods presented at the chapter [§3.3].
2.2.3 Double modal method
In the same way that
X ()
is projected on truncated modal basis,
X ()
p
can be projected on
this same base. Applying this projection to [éq 2.2.1-1], one obtains:
[]
(
)
[]
(
)
()
[]
()
[]
()
X
K
C
M
F
K
C
M
+
+
-
-
=
=
+
+
-
p
p
I
p
p
I
T
L
T
L
L
T
L
2
2
µ
éq
2.2.3-1
with projection on the basis of modal
X
()
p
who is written
[]
µ
L
()
.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
7/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
In practice, in [éq 2.2.3-1],
X
()
is known only after the first modal projection. Method
by double projection thus consists in carrying out two projections successively, and thus writing
(substitution of
X
()
by
X
L
in [éq 2.2.3-1]):
[]
(
)
[]
(
)
()
[]
()
[]
()
L
T
L
T
L
L
T
L
p
p
I
p
p
I
X
K
C
M
F
K
C
M
+
+
-
-
=
=
+
+
-
2
2
^
µ
éq 2.2.3-2
It should well be seen that the resolution of [éq 2.2.3-2] gives an approximation of
X ()
p
after one
double truncation. The first truncation is usual truncation carried out to obtain
projection of
X ()
. The second truncation is of comparable nature but concerns
X ()
p
.
Techniques of acceleration of the convergence of the series giving
X ()
p
were proposed in
literature with [bib3]. The user of Code_Aster will use the simple modal method preferentially
to obtain more precision in the results.
2.2.4 Establishment by the finite element method
For the calculation of the second members of [éq 2.2.1-1], [éq 2.2.2-2] and [éq 2.2.3-2], several choices are
possible according to whether the matric products are calculated on an elementary level or not.
For reasons of performance, it was selected to carry out the matric products of the second member
of [éq 2.2.1-1] on an elementary level. In Code_Aster, the second member of [éq 2.2.1-1] is thus
calculated directly by the assembly of the following elementary matrices:
()
(
)
(
)
E
G
E
G
E
T
G
E
G
E
G
E
T
G
E
G
has
G
E
G
E
T
G
E
G
G
G
p
C
p
I
p
p
F
X
H
X
X
)
(
))
(
(
)
(
)
(
)
(
))
(
(
)
(
)
(
))
(
(
2
+
+
-
éq 2.2.4-1
where:
·
same diagrams of integration (same functions of forms, points of Gauss,…) who have
allowed to obtain the standard linear system are used,
·
G
are the parametric co-ordinates of the point of Gauss of weight
G,
·
NR
E
G
(
)
is the matrix of the functions of form of the element
E
,
·
B
E
G
(
)
is the matrix connecting the deformations to the nodal variables for the element
E
.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
8/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
Thus, the assembled vector
-
+
+
2
M
C
K X
p
I
p
p
()
is obtained without calculation and
the assembly preconditions of the matrices
M
p
,
K
p
and
C
p
. This procedure is thus economic in
assembly and storage of assembled matrices. It has on the other hand the disadvantage of
to require the calculation of the elementary matrices derived at each frequency since the product by
X ()
is carried out on an elementary level and thus prevents the storage of these matrices
elementary derived (at least with the architecture of Code_Aster).
The second members of [éq 2.2.2-2] and [éq 2.2.3-2] will result them from a shandy between operations
elementary and matric operations.
Once again, one can notice that many terms of [éq 2.2.4-1] are null. For example if
p is Young the modulus of certain elements, only remainder the term
(
)
(
)
G
E
G
T
E
G
E
G
G
p
(
)
(
(
))
(
) (
)
H
X
for these elements and [éq 2.2.4-1] is null everywhere
elsewhere.
2.3
Derived from a temporal dynamic response into linear
One limits oneself in this chapter and the following if one does not derive compared to
parameters intervening under the initial conditions. This case raises theoretical difficulties and
practical out of the field of this Doc.
2.3.1 Method
direct
Let us leave the differential equation of dynamics in matric form:
M
X
C
X
K
X
F
()
& & (,)
()
& (,)
(,) (,)
(,)
p
T p
p
T p
T p
T p
T p
+
+
=
éq
2.3.1-1
where:
·
X ()
T
R
N
(
= +
N
N2
) is, at the moment T, the vector made up of N ddl physical of
structures discretized and of
2
multipliers of Lagrange [R3.03.01],
·
F (,)
T p
is the vector of the external forces, with, to calculate the response in transient
linear in Code_Aster [U4.53.02], the loading which must be written in the form
F
F X
=
I
I
I
T
() ()
.
The differential equation governing the derivative compared to p of the vector displacement
X
,
Y
X
= p
, is
obtained while deriving directly [éq 2.3.1-1] compared to
p
. After change about
derivation compared to time and with
p
(the sufficient regularity is supposed), and after rearrangement
terms, one obtains:
() () () ()
() ()
()
() () () ()
() ()
p
T
p
p
p
T
p
p
p
T
p
p
p
p
T
p
T
p
K
p
T
p
p
T
p
,
,
,
,
,
,
,
X
K
X
C
X
M
F
Y
Y
C
Y
M
-
-
-
=
+
+
&
&&
&
&&
éq 2.3.1-2
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
9/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
The vector of displacements and its derivatives are obtained while solving successively [éq 2.3.1-1]
then [éq 2.3.1-2]. For that, one uses a diagram of numerical integration (implicit or clarifies in
Code_Aster). On the same principle as that of the derivation of a harmonic answer, one exploits
intermediate results of the resolution of the standard problem in order to minimize calculations related to
resolution of the derived problems.
Let us specify this in the case of an implicit scheme, the identical remaining principle in the case of one
explicit diagram: For the standard problem, the vectors accelerations and speeds, with the pitch
of integration
T
N
+1, is obtained by the linear approximations of the following form:
&&
(
,
,
&, & &)
X
X
X X X
N
N
N
N
N
L
+
+
=
1
1
1
éq
2.3.1-3
&
(
,
,
&, & &)
X
X
X X X
N
N
N
N
N
L
+
+
=
1
2
1
éq
2.3.1-4
One solves the standard problem then, i.e one obtains
X
N
+1, by the substitution of [éq 2.3.1-3] and
[éq 2.3.1-4] in [éq 2.3.1-1]. This gives a linear system to solve.
The idea for the problems derived is to use the same diagram of integration:
&&
(
,
,
&, & &)
Y
Y
Y Y Y
N
N
N
N
N
L
+
+
=
1
1
1
éq
2.3.1-5
&
(
,
,
&, & &)
Y
Y
Y Y Y
N
N
N
N
N
L
+
+
=
1
2
1
éq
2.3.1-6
Thus, factorization used to solve the system can be also used to obtain
Y
N
+1
.
Indeed, since the matrices coefficient of [éq 2.3.1-1] and [éq 2.3.1-2] are identical, and if one uses
the same diagram of integration for the calculation of
X
N
+1
and
Y
N
+1
, then the matrix coefficient is
even in the two linear systems to solve after substitution of [éq 2.3.1-3] and [éq 2.3.1-4]
in [éq 2.3.1-1] and of [éq 2.3.1-5] and [éq 2.3.1-6] in [éq 2.3.1-2].
The calculating time of derived can then be negligible compared to the calculating time associated with
standard problem. It should however be noted that, in the case of a resolution by explicit diagram,
each calculation of derived can be as expensive as standard calculation. Indeed, in this case,
the main part of time CPU for the standard problem is associated at the stage of calculation of the second member
with each pitch of time; and it is also necessary to carry out this stage for the derivative (calculation of the second
member of [éq 2.3.1-2]).
2.3.2 Method
modal
Two modal methods simple and double for calculation of derived from a dynamic response
frequential ([§2.2.2] and [§2.2.3]) are easily transposable with the calculation of the dynamic response in
temporal. One obtains a linear system of the same type as [éq 2.3.1-2] with second members
different according to whether one uses the simple or double modal method. In the case, method
modal double, one obtains:
[]
[]
() ()
[]
()
[]
()
[]
()
[]
()
[]
()
[]
() ()
[]
() ()
[]
() ()
p
T
p
p
p
T
p
p
p
T
p
p
p
p
T
p
T
p
p
T
p
p
T
p
L
T
L
L
T
L
L
T
L
T
L
L
L
L
T
L
L
T
L
,
,
,
,
,
^
,
^
,
^
X
K
X
C
X
M
F
Y
K
Y
C
Y
M
-
-
-
=
+
+
&
&&
&
&&
éq 2.3.2-1
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
10/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
The response in generalized co-ordinates is obtained by solving this system by a diagram
of numerical integration of the same type as [éq 2.3.1-3] and [éq 2.3.1-4]. And in the same way, them
factorizations used for the calculation of the answer can be used for the calculation of the derivative.
2.3.3 Establishment by the finite element method
In a way similar to the preceding cases [§2.2.4], the second member of [éq 2.3.1-2] can be calculated
after assembly of the matrices
M
p
,
K
p
and
C
p
, or in a more powerful way directly to one
elementary level.
The temporal integration of
X
with the pitch
T
N
+1
is carried out the temporal integration before of
its derivative
Y
. With the pitch of integration
T
N
+1
of
Y
, terms
X
(
,)
T
p
N
+1
,
& (
,)
X
T
p
N
+1
, and
& & (
,)
X
T
p
N
+1
are thus known (the last two terms are calculable by [éq 2.3.1-3] and [éq 2.3.1-4] by
example).
With the pitch of integration
T
N
+1
of
Y
,
M
X
C
X
K
X
() && (,)
() & (,)
() (,)
p
p
T
p
p
p
T
p
p
p
T
p
N
N
N
+
+
+
+
+
1
1
1
result thus from the assembly of the following elementary matrices (cf [§3.2.3] for the others
notations):
)
(
)
(
))
(
(
)
(
)
,
(
)
(
)
(
))
(
(
)
(
)
(
)
(
))
(
(
1
1
1
+
+
+
+
+
N
E
T
G
E
G
E
T
G
E
G
N
G
G
T
G
E
G
has
G
N
E
G
T
G
E
G
G
T
p
T
C
p
T
p
X
B
H
B
X
NR
NR
X
NR
NR
&
&
&&
éq
2.3.3-1
2.4
Derived from a temporal dynamic response into nonlinear
Let us consider the class of the problems of nonlinear dynamics of which the equation discretized of
movement is form:
[
]
M
X
C
X
G X
F
()
& & (,)
()
& (,)
(,),
(,)
p
T p
p
T p
T p p
T p
+
+
=
éq
2.4-1
where
[
]
G X (,),
T p p
is a function (nonlinear) of
X
(,)
T p
and of
p
. It represents the forces
interns of the system and all the other forces which are dependant on
X ()
T
.
Noting the tangent matrix
[
]
K X
G X
X
((,),)
(,),
T p p
T p p
=
, the equation of derived compared to
p
vector displacement is obtained:
()
() () () ()
()
(
)
fixed
X
X
G
X
C
X
M
F
Y
X
K
Y
C
Y
M
p
p
p
T
p
T
p
p
p
T
p
p
p
p
T
p
T
p
p
T
p
T
p
p
T
p
-
-
-
=
+
+
,
,
,
,
,
)
,
(
)
),
,
(
(
)
,
(
)
(
)
,
(
)
(
&
&&
&
&&
éq
2.4-2
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
11/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
With share the term
G
X
()
p
p
fixed
, all the terms derived from [éq 2.4-2] are calculated in a conventional way
(i.e by assembly of the matrices or elementary vectors derived). The term
G
X
()
p
p
fixed
can
to require a particular processing. Indeed, [éq 2.4-1], the dependence enters
G
and the value of
X
(,)
T p
is explicit at the moment T; however
G
also depends on the history in plasticity. This
dependence complicates calculations clearly of
G
X
()
p
p
fixed
that we will not approach here.
What is remarkable, it is that the derived problems [éq 2.4-2] are linear problems then
that the standard problem is nonlinear. This property allows a calculation of really derived
very economic compared to the cost calculation of the standard problem. This is all the more true as,
for the resolution of [éq 2.4-1] and [éq 2.4-2], factorizations of matrices used for the calculation of
vector displacement can be re-used for the calculation of the derived vector; as it is the case
in the case of a linear standard problem. Indeed, the resolution of [éq 2.4-1] by an algorithm of
newton type requires the resolution of a linear system whose matrix coefficient is the same one as
in [éq 2.4-2]. Once the algorithm of the Newton type converged to obtain
X
(
)
T
N
+1, the vector
gradient
Y
(
)
T
N
+1 can thus be obtained directly not a simple resolution of two systems
triangular without convergence of an iterative algorithm.
3
Derived from the frequencies and clean modes from vibration
3.1 Problem
standard
The problem with the eigenvalues (problem known as quadratic) associated [éq 2.1-1] or [éq 2.3.1-1]
is written:
(
)
()
()
() ()
()
()
p
p
p
p
p
p
2
0
M
C
K
+
+
=
éq
3.1-1
where
()
p
and
()
p
are complex values and clean vectors.
In the continuation, to gain in compactness of the equations, we will omit the dependence of the sizes
compared to the parameter
p
and/or compared to time
T
, except when the clarification of this
dependence will be able to avoid ambiguities.
The reduction of the quadratic problem [éq 3.1-1] in a problem generalized equivalent [R5.01.02]
give the system of size
NR NR
N
N
× =
×
(
) (
)
2
2
according to:
-
=
K O
O
M
C M
MR. O
éq
3.1-2
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
12/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
One can standardize the clean vectors so that they check:
J
J
T
J
J
J
-
=
K
0
0
M
éq
3.1-3
J
J
T
J
J
=
C M
M
0
1
éq
3.1-4
3.2
Derived from the Eigen frequencies
The multiple eigenvalues are not differentiable with the direction running of derived (i.e within the meaning of
Fréchet). We must thus distinguish the case from the simple eigenvalues of the case of the eigenvalues
multiples for which some precautions must be taken.
3.2.1 The case of the simple eigenvalues
The derivation of the equation [éq 3.1-1] written for the J
ième
mode compared to a scalar parameter
p
give:
(
)
J
J
J
J
J
J
J
J
J
p
p
p
p
p
2
2
2
M
C K
MR. C
M
C
K
+
+
+
-
+
+
= -
(
)
éq 3.2.1-1
While multiplying on the left by
jT two members of [éq 3.2.1-1] and while using [éq 3.1-1] and
[éq 3.1-4], the derivative of the eigenvalue can be obtained from:
J
jT
J
J
J
p
p
p
p
= -
2
M
C
K
+
+
éq
3.2.1-2
The equation [éq 3.2.1-2] above requires the calculation of the clean vector
J associated with
J. One cannot
thus not to exempt calculation of the clean vectors for the calculation of derived from eigenvalues.
3.2.2 The case of the multiple eigenvalues
In a clean subspace whose associated eigenvalue is multiple of command
m
, any combination
linear of vectors is clean vector. One cannot however apply the formula directly
[éq 3.2.1-2] with any of these vectors, because the clean subspace can be divided into
m
distinct clean subspaces when a parameter
p
vary.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
13/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
On the theoretical level, the multiple eigenvalues are not differentiable as a function with
several variables, only of the derivative directional exist (this can pose problems of
convergence with many algorithms of optimization which are based on the differentiability of
function objectifies as a function with several variables). As example, let us consider the values
clean of the following matrix:
With
=
+
4 2
4
y X
X
,
1
2
2
2
2
2
4
4
= + +
+
= + -
+
y
X
y
y
X
y
For
()
()
X y
,
,
= 0 0
, the eigenvalue of
With
is multiple. And, one a:
1
0
1
X
(
, 0)
-
= -
,
1
0 0)
1
X
(
,
+
=
,
2
0
1
X
(
, 0)
-
=
,
2
0
1
X
(
, 0)
+
= -
.
It is possible “to connect”,
1
for
X
< 0
with
2
for
X
> 0
, and “to connect”
2
for
X
< 0
with
1
for
X
> 0
(let us take for example
0
=
y
. Then, “the smallest eigenvalue of
With
“
(
)
X
-
2
is not differentiable in 0, and “the largest either”
(
)
X
+
2
. On the other hand, “the eigenvalue of
With
corresponding to the clean vector
()
T
1
,
1
“
(
)
X
+
2
and “the eigenvalue of
With
corresponding to
clean vector
(
)
T
1
,
1
-
“
(
)
X
-
2
are differentiable into 0). If one does not consider independently
two Eigen frequencies of
With
, one could thus write
/X
= ±1
. (One can proceed of
even with the derivative in
y
and to write
/y
= 0
2
or
). However, eigenvalues of
With
are thus not differentiable as functions with two variables (i.e one cannot
to write
D
D
D
=
+
X X
y y
), and have only derivative directional. For of
to convince, it is enough to notice that:
(cos, sin)
lim (cos, sin)
(,)
sin
=
=
-
=
±
0
0
0 0
1
On the plan practices, for the calculation of these derivative directional, the choice of the clean vector
J
associated
J in [éq 3.2.1-2] is ambiguous when the eigenvalue is multiple. This choice thus must
to be specified.
In the case of an eigenvalue
m of command of multiplicity
m
, the problem with the eigenvalues
[éq 3.1-1] can rewrite itself in the following matric form:
[]
[]
[]
M
C
K
m
m
m
m
m
2
0
+
+
=
éq
3.2.2-1
with
m
m
m
I
=
and
[] [
]
,
m
I
I
I m
=
+
+
+
1
2
,
L
where them
I K
+ are the clean vectors
associated
Mr.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
14/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
Note:
In any rigor and in the general case, the matrices
M
and
K
not being definite positive (with
cause dualisation of the boundary conditions), it does not exist inevitably
m
vectors
clean generating a clean subspace of dimension equalizes with the multiplicity of the value
clean. This point is never approached in the literature. By convenience, we will not approach
either this pathological point which becomes delicate when it is about calculation of derivative.
The conditions of orthogonality of the clean vectors are written with the matric notations
above:
[]
(
)
[]
m T
m
m
m
2
MR. C
I
+
=
éq
3.2.2-2
Let us note
[]
X
m
clean vectors of the same multiple subspace clean which become
m
clean vectors
()
[
]
p
X
of
m
subspaces two to two distinct when when a parameter
p
vary. When
p
tends towards its face value
()
[
]
p
p
O
X
,
tends towards
[]
X
and them
m
eigenvalues
distinct tend towards the multiple eigenvalue
m
. In the literature, such clean vectors are
called clean vectors “adjacent”. It should be noted that these vectors depend on the parameter which varies
i.e. parameter by report/ratio to which one derives. In other words, if parameter is changed, it
is then necessary to change adjacent vectors.
Adjacent vectors
[]
X
can be connected to
[]
m by an orthogonal transformation
T
(
[]
X
and
[]
m
two bases of the same clean subspace constitute):
[]
[]
X
T
=
m
éq
3.2.2-3
with
T T
I
T
m
=
.
The procedure consists in finding the transformation
T
for then deducing some
[]
X
and
m
p
. It is it
that we develop below.
The adjacent vectors check the problem with the eigenvalues [éq 3.2.2-1]:
[]
[]
[]
MR. X
C X
K X
m
m
2
0
+
+
=
éq
3.2.2-4
The derivation of [éq 3.2.2-4] compared to
p
give:
(
)
[]
(
)
[]
[]
m
m
m
m
m
m
p
p
p
p
p
2
2
2
M
C K
X
MR. C X
M
C
K X
+
+
= -
+
+
+
-
éq 3.2.2-5
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
15/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
While multiplying [éq 3.2.2-5] on the left by
[]
m T, using [éq 3.2.2-2] and substituent
[]
[]
X
T
=
m
, one
the new problem with the eigenvalues obtains:
DT T
=
m
p
éq
3.2.2-6
with
[]
[]
D
M
C
K
=
+
+
m T
m
m
m
p
p
p
2
.
The resolution of the problem to the eigenvalue [éq 3.2.2-6] gives the derivative
m
p
eigenvalues
and stamps it orthogonal transformation
T
.
It should be stressed that there will be a subproblem with the eigenvalues of size
m
to solve for each
eigenvalue of multiplicity
m
to derive and for each parameter.
3.2.3 Establishment by the finite element method
In a conventional way, several choices are possible for the calculation of the second member of [éq 3.2.1-2] and
matrix
D
of [éq 3.2.2-6] according to whether the matric products are calculated on an elementary level
or not.
The calculation of the second member of [éq 3.2.1-2] on an elementary level is a priori most powerful.
T
p
p
p
2
M
C
K
+
+
can be written directly on an elementary level by the scalar:
(
) (
)
(
) (
)
(
)
(
)
G
G
G
T
G
G
E has
G
G
T
G
G
E
G
G
T
E
G
E
G
G
p
PC
p
2
((
))
(
)
(
)
(
(
))
(
)
(
)
(
) (
)
(
(
))
(
) (
)
+
+
B
H
B
éq
3.2.3-1
In particular, a procedure not to start the calculation of each elementary term of
[éq 3.2.3-1] that if quantities
p ()
,
p
E
G
(
(
))
H
and
PC
has
G
((
))
are nonnull allows
to minimize time CPU.
In a way similar to the case of the simple eigenvalues, the matrix
D
system [éq 3.2.3-1] can
to be calculated directly on an elementary level or after effective assembly of the matrices
M
p
,
K
p
and
C
p
.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
16/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
3.3
Derived from the clean modes of vibration
3.3.1 The case of the simple eigenvalues
One cannot directly obtain the derivative of the clean vectors by solving the linear system
[éq 3.2.1-1] because it is a noninvertible system since
J
J
2
M
C K
+
+
is singular (by
definition of
J
).
Many methods were proposed in the literature to solve or circumvent this
difficulty. Among them, two types of methods emergent: methods of projection on basis
modal incomplete [bib4], [bib5] and algebraic methods [bib6], each one of them with theirs
alternatives. These methods are presented in [bib1]. The algebraic methods are methods
exact which, moreover, does not require the other clean vectors to calculate the derivative of one
clean vector given. One presents only one algebraic method here consisting in adding one
equation with [éq 3.2.1-1] to obtain a regular and symmetrical system to solve [bib6].
The rewriting of the condition of orthonormality [éq 3.1-4] gives:
(
)
jT
J
J
2
1
MR. C
+
=
éq
3.3.1-1
The derivation of this equation compared to the parameter
p
leads to:
(
)
jT
J
J
jT
J
J
J
p
p
p
p
2
1
2
2
0
MR. C
M
M
C
+
+
+
+
=
éq
3.3.1-2
The equations [éq 3.2.1-1] and [éq 3.3.1-2] can be rewritten together:
(
)
(
)
(
)
(
)
+
+
-
+
+
-
+
-
=
+
+
+
+
J
J
J
T
J
J
J
J
J
J
J
J
J
T
J
J
J
J
J
p
p
p
p
p
p
p
p
C
M
M
K
C
M
C
M
C
M
C
M
K
C
M
2
2
1
2
0
0
2
2
2
2
éq 3.3.1-3
The resolution of the system [éq 3.3.1-3] is the key idea of the algebraic method. It can be shown that
this system is always regular. The derivative
J
p
clean vectors can be obtained
directly by solving this system of dimensions
(
) (
)
+ × +
N
N
1
1
with a factorization
LDL
T
for example.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
17/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
The algebraic method preserves the structure in tape and the symmetry of the matrices. Moreover, for
to calculate the derivative of a clean vector, only the knowledge of this clean vector and value
clean corresponding is necessary. Indeed, [éq 3.3.1-3] does not utilize the other vectors
clean that that to derive.
However, the algebraic method requires the resolution of a system of size
(
) (
)
+ × +
N
N
1
1
different for each vector suitable to derive and each parameter.
3.3.2 The case of the multiple eigenvalues
Taking again the notations of the paragraph [§3.2.2], the derivative
[]
X
p
clean vectors “adjacent”
[]
[]
X
T
=
m
are solutions of the equation [éq 3.2.2-5].
The matrix
(
)
m
m
2
M
C K
+
+
being singular, one cannot solve [éq 3.2.2-5] directly.
The idea, still, is to add an additional equation for
[]
X
p
.
The condition of standardization of
[]
X
is written:
(
)
jT
J
J
2
1
MR. C
+
=
éq
3.3.2-1
The derivative of [éq 3.3.2-1] gives the sought additional equation:
[]
(
)
[]
[]
[]
[]
[]
[]
[]
X
MR. C
X
X
MR. X
X MR. X
X
C X
T
m
T
m
T
m
T
p
p
p
p
2
1
2
+
= -
-
-
éq 3.3.2-2
[éq 3.2.2-5] and [éq 3.3.2-2] form the following system to solve:
(
)
(
)
[]
[]
(
)
[]
(
)
[]
[]
[]
[]
[] []
[]
[]
-
-
-
+
+
-
+
-
=
+
+
+
+
X
C
X
X
M
X
X
M
X
X
K
C
M
X
C
M
X
C
M
X
X
C
M
K
C
M
p
p
p
p
p
p
p
p
T
m
T
m
T
m
m
m
m
m
T
m
m
m
2
1
2
0
0
2
2
2
2
éq 3.3.2-3
The system [éq 3.3.2-3] to solve is of dimension
(
) (
)
+
× +
N
m
N
m
. To calculate the derivative of one
clean vector, the knowledge of the clean vectors associated the eigenvalue is necessary.
As for the case of the nonmultiple eigenvalues, the structure in tape and the symmetry of
matrices are preserved and it can be shown that the system to be solved is always regular. Of
more, the system must be solved for each vector suitable to derive and each parameter.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
18/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
3.3.3 Establishment by the finite element method
In the equations [éq 3.3.1-3] and [éq 3.3.2-3], the second members can be calculated either in
assembling the matrices initially
M
p
,
K
p
and
C
p
then by carrying out matric operations, that is to say
by directly calculating the various terms on an elementary level.
A direct calculation at the elementary level being more powerful,
2
M
C
K
p
p
p
+
+
, by
(
)
(
)
G
G
E
G T
G
G
has
G
E
G T
G
G
E
G
T
E
G
E
G
G
p
PC
p
2
((
))
(
)
(
)
((
))
(
)
(
)
(
)
(
(
))
(
) (
)
+
+
H
éq
3.3.3-1
4
Derived compared to a variation from field
The matric equation [éq 3.1-1] uses an explicit dependence of the matrices of mass,
of damping, and stiffness compared to the variable parameter
p
. When it is a question of deriving by
report/ratio with a variation of field, one cannot in the general case write this dependence
directly. In this case, obtaining derived from the dynamic sizes requires a work
precondition. It is necessary on the one hand to describe the variation of field in a mathematical form which one can
to handle, and in addition to return to the variational formulation of the equations of dynamics that
one derives directly.
This chapter discusses these various items. One gives in particular the derivative of the various operators
of mass, stiffness necessary to the calculations of derived from all the dynamic sizes. One
limit then with the calculation of derived from the Eigen frequencies and clean vectors. One limits oneself to
conservative problem and thus with the real clean modes.
4.1
Variational formulation of the standard problem
One considers an elastic solid on a field
, subjected to surface static loadings
F
and voluminal
F
. In variational form, the problem of balance is written:
() ()
0
0
V
V
+
=
v
F.v
f.v
v
:
U
With
U
F
D
D
D
:
that
such
To find
éq
4.1-1
where
With
is the tensor of elasticity,
()
U
the conventional tensor of the deformations linearized and
V
0
space
of Sobolev of
(
)
H
1
()
3
displacements kinematically acceptable.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
19/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
The clean modes around this balance are obtained by the resolution of the following spectral problem:
() ()
()
0
2
0
0
V
=
-
+
>
v
w.v
v
:
U
W
v
:
W
With
W
D
D
D
:
that
such
no one
not
mode
and
real
clean
frequency
To find
V
éq
4.1-2
(notation: if
B
With
and
are two tensors of command 2,
B
With
is the product contracted of
B
With
and
)
where
is a parameter of possible amplification of the solution of [éq 4.1-1] for the loadings
statics above,
is the tensor of the stresses given by
()
()
U
With U
=
and
mass
voluminal of material considered.
This problem can be rewritten by using the following conventional multilinear operators:
·
elastic operator of stiffness
() ()
has
D
(,)
W v
With W: v
=
·
geometrical operator of stiffness
()
B
D
(,)
U W v
W
U: v
=
()
=
ik
L I
kl
D
U W
v
,
()
·
operator of mass
C
D
(,)
W v
W. v
=
in the form:
()
(
)
()
(
)
=
=
-
+
>
1
0
,
0
0
0
2
W
W
v
v
W
v
W
U
v
W
W
,
C
,
C
,
B
,
has
V
V
:
that
such
no one
not
and
reality
To find
éq
4.1-3
4.2
Description of the variation of field and form of derivation
It is considered now that the field
is likely to change form, by a variation
of its edge (not included/understood the edge where are applied imposed displacements). To describe this
variation, one uses a bijective transformation making correspond the area of reference (that
that one nets) with the field modified (representation known as Lagrangienne of variation of field,
cf [R7.02.04]).
The edge of the field
()
is controlled by the scalar
(such as
()
0
=
), according to the field of
vectors
.
The derivative of a field of vector
X
compared to a variation of field can be now written
like the Lagrangian derivative, which one will note with a point (cf [R4.03.01]):
&X
X
X
X
=
=
+
=
=
D
D
0
0
éq
4.2-1
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
20/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
And, the derivative of the gradient of
X
is written:
(
)
D
D
D
D
= -
X
X
X
éq
4.2-2
. These integral sizes are
form:
()
()
()
()
(
)
()
I
U
U
=
E
D
éq
4.2-3
where
()
U
is the solution of the problem of balance [éq 4.1-1] on the field
()
.
This integral for derived (theorem of Reynolds):
()
()
(
)
()
()
(
)
D
D
E
of
D
D
I
=
=
=
+
0
0
U
U
U
div
éq
4.2-4
The derivative intervening in the intégrande is expressed:
()
()
(
)
()
()
of
D
E D
D
U
.
U
=
with
()
()
()
D
D
D
D
T
T
U
U
U
U
=
- -
éq
4.2-5
Note:
Implicitly, in the continuation the derivative will be taken in
= 0
.
4.3
Derivation of the variational formulation
4.3.1 Derivation of the elastic operator of stiffness
Using the symmetry of the tensor of elasticity,
has (,)
W v
is also written:
() ()
has
D
(,)
W v
With W:
v
=
éq
4.3.1-1
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
21/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
The derivation of this expression and the use of [éq 4.2-4] give:
(
)
() ()
() ()
()
()
D
D has
D
D
D
D
D
(,)
div
W v
With W:
v
With W: v A W: v
=
+
+
éq
4.3.1-2
By [éq 4.2-5], and knowing
()
()
With W
With
W
=
and
v = 0
, [éq 4.3.1-2] becomes:
(
)
() ()
()
()
() (
))
-
-
+
=
D
D
D
has
D
D
v
:
W
With
v
:
W
With
v
:
W
With
v
:
W
With
v
W
)
(
div
)
,
(
éq
4.3.1-3
that one can rewrite, by using amongst other things again the symmetry of the tensor of elasticity:
(
)
() ()
(
)
-
-
=
+
=
D
has
has
D
D
has
has
D
D
v
:
W
v
:
W
v
:
W
With
v
W
v
W
v
W
v
W
)
(
)
(
div
)
,
(
^
)
,
(
^
)
,
(
)
,
(
with
éq
4.3.1.- 4
4.3.2 Derivation of the geometrical operator of stiffness
In the same way, one obtains:
-
-
+
-
+
=
D
D
D
D
D
B
D
D
v
:
U
W
v
:
U
With
W
v
:
U
With
W
v
:
U
W
v
:
U
W
v
:
U
W
v
W
U
)
(
)
(
)
(
div
)
)
(
(
)
,
,
(
éq 4.3.2-1
that one can rewrite:
D
D B
B dd
B
D
D
B
(,)
(
,)
(,
,)
$ (,)
U W v
U W v
U W v
U W v
=
+
+
éq
4.3.2-2
with:
(
)
-
-
-
=
D
With
B
v
:
U
W
v
:
U
W
v
:
U
W
v
:
U
W
v
W
U
)
(
)
(
div
)
)
(
(
)
,
,
(
^
éq 4.3.2-3
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
22/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
In [éq 4.3.2-2], one sees that the calculation of derived from the geometrical operator of stiffness requires
au préalable the derivative
D
D
U
solution of the problem of prestressed [éq 4.1-1].
This problem of prestressing, which is thus to derive, is written:
To find
U
V
such as:
()
()
has,
D
D
F
F
U v
F. v
F. v
v
v
=
+
=
V
0
éq
4.3.2-4
The derivation of this variational equation gives:
()
()
dd has,
D
D F
has,
U v
v
U v
=
- $
éq
4.3.2-5
with, by supposing that forces
F F
,
are independent of
and
:
(
) (
)
(
)
(
)
(
)
(
)
D
D F
D
dS
F
()
. div
. div
,
v
f.v
F. .v
F.v
.n
F. .v
N
=
+
+
-
+
éq 4.3.2-6
(The formula of derivation of the surface integral is given by proposal 4 of appendix 2 of
[R4.03.01]).
The resolution of the linear system associated [éq 4.3.2-5] giving
D
D
U
have the same matrix of
rigidity that the problem of standard prestressing, only the second member changes. Code_Aster solves
this problem in deformations plane and 2D-axisymmetric, cf [R4.03.01].
4.3.3 Derivation of the operator of mass
The derivation of the operator of mass gives:
()
(
)
D
D C,
D
D
D
v W
v.w
. v.w
v. W
=
+
+
div
that one can rewrite:
(
)
()
()
(
)
+
=
+
=
D
,
C
C
D
D
C
C
D
D
v.w
.
v.w
W
v
v
W
v
W
v
W
div
^
)
,
(
^
)
,
(
)
,
(
with
éq
4.3.3-1
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
23/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
4.3.4 Derivation of the variational formulation of the problem of modal analysis
Derivation in
= 0
problem [éq 4.1-3] of calculation of the clean modes on
()
, brought back on
()
0
=
, the parameter of initial loading
being fixed, gives now:
() ()
(
) (
) (
)
(
)
()
() ()
(
)
(
) (
)
=
+
=
+
+
-
+
+
+
+
0
^
2
0
^
(
2
,
^
,
,
^
W
W
W
W
v
v
W
v
W
v
W
v
W
U
v
W
U
v
W
U
v
W
v
W
W
,
C
,
C
,
C
,
C
,
C
,
B
,
B
,
B
,
has
,
has
&
&
&
&
&
&
&
&
V
V
:
that
such
and
reality
To find
éq
4.3.4-1
4.3.5 Obtaining derived from the Eigen frequencies and the clean modes
In a preoccupation with a simplification, we suppose here who the eigenvalue of the mode which one drift is
simple. The case of a multiple eigenvalue results from the adaptation of the developments from
paragraphs [§3.2.2] and [§3.3.2].
While using
(
)
(
)
(
)
(
)
has
,
B
,
C,
C,
W W
U W W
W W
W W
&
,
&
&
+
-
=
=
2
0
1
and while taking
v
W
=
in [éq 4.3.4-1], one
obtains the following expression of derived from the own pulsation:
(
)
(
) (
)
(
)
(
)
&
$
&,
$,
$
=
+
+
-
has
,
B
,
B
,
C
,
W W
U W W
U W W
W W
2
2
éq
4.3.5-1
Obtaining [éq 4.3.5-1] is to be brought closer that of [éq 3.2.1-2], except for standardization. In
cases of eigenvalues multiple, the choice
v
W
=
condition the calculated directional derivative.
[éq 4.3.4-1] gives the following linear system:
()
(
)
()
()
(
) (
)
(
)
() ()
(
)
has
,
B
,
C,
has
,
B
,
B
,
C,
C
,
&
, &
&
$
&,
$,
&
$
W v
U W v
W v
W v
U W v
U W v
W v
W v
+
-
= -
-
+
+
+
2
2
éq
4.3.5-2
This system is same form as the system [éq 3.2.1-1]. It must be solved same manner,
for example by using the addition of the equality
(
)
(
)
C,
C
,
&
$
W W
W W
= -
1
2
.
4.3.6 Establishment by the finite element method
In a way similar to the preceding chapters, terms of the second members [éq 4.3.5-1] and
[éq 4.3.5-2] can be calculated, either directly on an elementary level, or after assembly
matrices of mass, stiffness and damping.
In addition, of new procedures of elementary calculations must be developed for
calculations of
()
$a,
W v
(this term A already was developed in deformations plane and 2D-axisymmetric,
cf [R4.03.01]),
(
)
$,
B
,
U W v
,
()
$c,
W v
, and
D
D F
()
v
.
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
24/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
5
Effective establishment in Code_Aster
The numerical calculation of derived from the clean modes or the dynamic responses proceeds in
two main stages:
·
stage 1: calculation of the elementary terms constituting the new terms (often in
second member) associated the calculation of the derivative,
·
stage 2: resolution of the derived problem. The calculation of derived from the clean modes requires
the use of specific algorithms different from those used to solve the problem
standard. The calculation of derived from the other sizes, i.e. Eigen frequencies and
the dynamic responses, re-uses the standard algorithm by modifying some terms (in
General second members of the systems to be solved).
Sensitivities of the other dynamic sizes (forced, efforts…) require a third
stage which consists of the derivation of calculations of “postprocessing”. We do not treat these
sizes in this document.
The following table recapitulates the derivative available in Code_Aster. It is updated at the fur and at
measure developments of the derivatives.
Sizes and methods
Operators, and operands private individuals (i.e.
data necessary to the method of
calculation)
Established
Derived from a dynamic response
harmonic
Direct method
DYNA_LINE_HARM
YES
Simple modal method
DYNA_LINE_HARM
clean vectors and derived clean vectors
NOT
Double modal method
DYNA_LINE_HARM
clean vectors
NOT
Derived from a dynamic response
temporal into linear
Direct method
DYNA_LINE_TRAN
YES
Simple modal method
DYNA_TRAN_MODAL
clean vectors and derived clean vectors
NOT
Double modal method
DYNA_TRAN_MODAL
clean vectors
NOT
Derived from a dynamic response
temporal into nonlinear
standard loading Dirichlet
DYNA_NON_LINE
for some laws of behavior
NOT
YES
Derived from an Eigen frequency
simple modes
MODE_ITER_SIMULT, MODE_ITER_INV
multiple modes
YES
NOT
Derived from a clean mode
simple modes
MODE_ITER_SIMULT, MODE_ITER_INV
multiple modes
YES
NOT
Derived compared to a variation
of field
* NOT
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
25/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
6 Bibliography
[1]
R.A. IBRAHIM: Structural dynamics with parameter uncertainties, Applied Mechanics
Reviews, vol. 40, No 3, pp.309-328, 1987.
[2]
S. CAMBIER: Vibratory sensitivities of sizes. Theory and algorithms for one
establishment in Code_Aster. Note HT-62/01/011.
[3]
Z.Q. Qu: Accurate methods for frequency responses and to their sensitivities off proportionally
damped systems, Computers & Structures, vol. 79, pp.87-96, 2001.
[4]
Mr. YU, Z.S. LIU, D.J WANG: Modal Comparison off several approximate methods for
computing shape mode derivative, Computers & Structures, vol. 62, No.2, pp.381-393,
1996.
[5]
S. ADHIKARI, Mr. I. FRISWELL: Not-conservative Eigenderivative analysis off asymmetric
systems, Int. J. Numer. Meth. Engng; Vol. 51, pp. 709-733, 2001.
[6]
I.W. LEE, D.-O. KIM, G. - H. JUNG: Natural frequency and mode shape sensitivities off
damped systems: Leaves I, distinct natural frequencies, Journal off Sound and Vibration 223 (3),
pp. 399-412 (1999). Leaves II, multiple natural frequencies, Journal off Sound and Vibration
223 (3), pp. 313-424, (1999).
[7]
I. VAUTIER: Isoparametric elements [R3.01.00].
[8]
J. PELLET: Dualisation of the boundary conditions [R3.03.01].
[9]
G NICOLAS: Sensitivity of the mechanical thermo fields to a variation of the field
[R4.03.01].
[10]
D. SELIGMANN, R. MICHEL: Algorithms of resolution for the quadratic problem with
eigenvalues [R5.01.02].
[11]
O. BOITEAU: Algorithm of resolution for the generalized problem [R5.01.01].
[12]
G. JACQUART: Methods of RITZ in linear and non-linear dynamics [R5.06.01].
[13]
G. DEBRUYNE: Lagrangian representation of variation of field [R7.02.04].
[14]
D. GIRARDOT: Operator DYNA_LINE_TRAN [U4.53.02].
Code_Aster
®
Version
7.4
Titrate:
Calculation of sensitivity in dynamics
Date:
02/05/05
Author (S):
H. ANDRIAMBOLOLONA,
S. CAMBIER
Key
:
R4.03.04-B
Page
:
26/26
Manual of Reference
R4.03 booklet: Analyze sensitivity
HT-66/05/002/A
Intentionally white left page.