1
!-------------------------------------------------------------------------------
3
! This file is part of the Code_Saturne Kernel, element of the
4
! Code_Saturne CFD tool.
6
! Copyright (C) 1998-2009 EDF S.A., France
8
! contact: saturne-support@edf.fr
10
! The Code_Saturne Kernel is free software; you can redistribute it
11
! and/or modify it under the terms of the GNU General Public License
12
! as published by the Free Software Foundation; either version 2 of
13
! the License, or (at your option) any later version.
15
! The Code_Saturne Kernel is distributed in the hope that it will be
16
! useful, but WITHOUT ANY WARRANTY; without even the implied warranty
17
! of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
18
! GNU General Public License for more details.
20
! You should have received a copy of the GNU General Public License
21
! along with the Code_Saturne Kernel; if not, write to the
22
! Free Software Foundation, Inc.,
23
! 51 Franklin St, Fifth Floor,
24
! Boston, MA 02110-1301 USA
26
!-------------------------------------------------------------------------------
29
!===============================================================================
31
! METHODE ALE - MOUVEMENT DE STRUCTURES EN COUPLAGE INTERNE
33
! NBSTRU : NOMBRE DE STRUCTURES MOBILES
35
! XMSTRU : MATRICE DE MASSE DES STRUCTURES (kg)
36
! XCSTRU : MATRICE DE FRICTION DES STRUCTURES (kg/s)
37
! XKSTRU : MATRICE DE RAIDEUR DES STRUCTURES (kg/s2 = N/m)
38
! XSTREQ : VECTEUR ECART DE LA POSITION DES STRUCTURE DANS LE MAILLAGE
39
! INITIAL PAR RAPPORT A LEUR POSITION D'EQUILIBRE (m)
40
! XSTR : VECTEUR DEPLACEMENT DES STRUCTURES PAR RAPPORT A LEUR POSITION
41
! DANS LE MAILLAGE INITIAL (m)
42
! XPSTR : VECTEUR VITESSE DES STRUCTURES (m/s)
43
! XPPSTR : VECTEUR ACCELERATION DES STRUCTURES (m/s2)
44
! FORSTR : VECTEUR FORCE EXERCE SUR LES STRUCTURES (N)
45
! XSTA : VALEUR DE XSTR AU PAS DE TEMPS PRECEDENT
46
! XPSTA : VALEUR DE XPSTR AU PAS DE TEMPS PRECEDENT
47
! XPPSTA : VALEUR DE XPPSTR AU PAS DE TEMPS PRECEDENT
48
! FORSTA : VALEUR DE FORSTR AU PAS DE TEMPS PRECEDENT
49
! XSTP : VALEUR PREDITE DE XSTR
50
! FORSTP : VALEUR PREDITE DE FORSTR
51
! DTSTR : PAS DE TEMPS ASSOCIE AU MOUVEMENT DES STRUCTURES
52
! AEXXST : COEFFICIENT DE PREDICTION DU DEPLACEMENT (SUR XPSTR)
53
! BEXXST : COEFFICIENT DE PREDICTION DU DEPLACEMENT (SUR XPSTR-XPSTA)
54
! CFOPRE : COEFFICIENT DE PREDICTION DES EFFORTS
59
common / istruc / nbstru
61
double precision xmstru(3,3,nstrmx)
62
double precision xcstru(3,3,nstrmx)
63
double precision xkstru(3,3,nstrmx)
64
double precision xstr(3,nstrmx) ,xsta(3,nstrmx)
65
double precision xstp(3,nstrmx) ,xstreq(3,nstrmx)
66
double precision xpstr(3,nstrmx) ,xpsta(3,nstrmx)
67
double precision xppstr(3,nstrmx),xppsta(3,nstrmx)
68
double precision forstr(3,nstrmx),forsta(3,nstrmx)
69
double precision forstp(3,nstrmx)
70
double precision dtstr(nstrmx)
71
double precision aexxst, bexxst, cfopre
73
common / rstruc / xmstru, xcstru, xkstru, xstr , xsta , xstp, &
74
xstreq, xpstr , xpsta , xppstr, xppsta, &
75
forstr, forsta, forstp, dtstr, &
76
aexxst, bexxst, cfopre
78
! METHODE DE NEWMARK HHT
80
! ALPNMK : COEFFICIENT ALPHA DE LA METHODE DE NEWMARK HHT
81
! BETNMK : COEFFICIENT BETA DE LA METHODE DE NEWMARK HHT
82
! GAMNMK : COEFFICIENT GAMMA DE LA METHODE DE NEWMARK HHT
83
double precision alpnmk, gamnmk, betnmk
84
common / rnemrk / alpnmk, gamnmk, betnmk