1
!-------------------------------------------------------------------------------
3
! This file is part of Code_Saturne, a general-purpose CFD tool.
5
! Copyright (C) 1998-2011 EDF S.A.
7
! This program is free software; you can redistribute it and/or modify it under
8
! the terms of the GNU General Public License as published by the Free Software
9
! Foundation; either version 2 of the License, or (at your option) any later
12
! This program is distributed in the hope that it will be useful, but WITHOUT
13
! ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
14
! FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
17
! You should have received a copy of the GNU General Public License along with
18
! this program; if not, write to the Free Software Foundation, Inc., 51 Franklin
19
! Street, Fifth Floor, Boston, MA 02110-1301, USA.
21
!-------------------------------------------------------------------------------
23
! Module for ALE structure movement with internal coupling
27
!=============================================================================
31
!=============================================================================
33
! Methode ale - mouvement de structures en couplage interne
35
! nbstru : nombre de structures mobiles
37
! xmstru : matrice de masse des structures (kgl)
38
! xcstru : matrice de friction des structures (kg/s)
39
! xkstru : matrice de raideur des structures (kg/s2 = N/m)
40
! xstreq : vecteur ecart de la position des structure dans le maillage
41
! initial par rapport a leur position d'equilibre (m)
42
! xstr : vecteur deplacement des structures par rapport a leur position
43
! dans le maillage initial (m)
44
! xpstr : vecteur vitesse des structures (m/s)
45
! xppstr : vecteur acceleration des structures (m/s2)
46
! forstr : vecteur force exerce sur les structures (N)
47
! xsta : valeur de xstr au pas de temps precedent
48
! xpsta : valeur de xpstr au pas de temps precedent
49
! xppsta : valeur de xppstr au pas de temps precedent
50
! forsta : valeur de forstr au pas de temps precedent
51
! xstp : valeur predite de xstr
52
! forstp : valeur predite de forstr
53
! dtstr : pas de temps associe au mouvement des structures
54
! aexxst : coefficient de prediction du deplacement (sur xpstr)
55
! bexxst : coefficient de prediction du deplacement (sur xpstr-xpsta)
56
! cfopre : coefficient de prediction des efforts
58
integer, save :: nbstru
60
double precision, save :: xmstru(3,3,nstrmx)
61
double precision, save :: xcstru(3,3,nstrmx)
62
double precision, save :: xkstru(3,3,nstrmx)
63
double precision, save :: xstr(3,nstrmx) ,xsta(3,nstrmx)
64
double precision, save :: xstp(3,nstrmx) ,xstreq(3,nstrmx)
65
double precision, save :: xpstr(3,nstrmx) ,xpsta(3,nstrmx)
66
double precision, save :: xppstr(3,nstrmx),xppsta(3,nstrmx)
67
double precision, save :: forstr(3,nstrmx),forsta(3,nstrmx)
68
double precision, save :: forstp(3,nstrmx)
69
double precision, save :: dtstr(nstrmx)
70
double precision, save :: aexxst, bexxst, cfopre
72
! methode de Newmark hht
74
! alpnmk : coefficient alpha de la methode de Newmark hht
75
! betnmk : coefficient beta de la methode de Newmark hht
76
! gamnmk : coefficient gamma de la methode de Newmark hht
78
double precision, save :: alpnmk, gamnmk, betnmk
80
!=============================================================================