1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
|
# CONFIGURATION MANAGEMENT OF EDF VERSION
# ======================================================================
# COPYRIGHT (C) 1991 - 2012 EDF R&D WWW.CODE-ASTER.ORG
# THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
# IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
# THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
# (AT YOUR OPTION) ANY LATER VERSION.
#
# THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
# WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
# MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
# GENERAL PUBLIC LICENSE FOR MORE DETAILS.
#
# YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
# ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
# 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
# ======================================================================
# person_in_charge: albert.alarcon at edf.fr
def assemblage_prod(self,NUME_DDL,MATR_ASSE,VECT_ASSE,**args):
if ((not MATR_ASSE) and (not VECT_ASSE)): raise AsException("Aucun concept a assembler")
if not NUME_DDL : raise AsException("Impossible de typer les concepts resultats")
if NUME_DDL.is_typco():
self.type_sdprod(NUME_DDL,nume_ddl_sdaster)
if MATR_ASSE !=None:
for m in MATR_ASSE:
opti=m['OPTION']
if opti in ( "RIGI_MECA","RIGI_FLUI_STRU",
"MASS_MECA" , "MASS_FLUI_STRU" ,"RIGI_GEOM" ,"RIGI_ROTA",
"AMOR_MECA","IMPE_MECA","ONDE_FLUI","MASS_MECA_DIAG",
"MECA_GYRO","RIGI_GYRO" ) : t=matr_asse_depl_r
if opti in ( "RIGI_ACOU","MASS_ACOU","AMOR_ACOU",) : t=matr_asse_pres_c
if opti in ( "RIGI_THER","MASS_THER","RIGI_THER_CONV" ,
"RIGI_THER_CONV_D",) : t=matr_asse_temp_r
if opti == "RIGI_MECA_HYST" : t= matr_asse_depl_c
self.type_sdprod(m['MATRICE'],t)
if VECT_ASSE !=None:
for v in VECT_ASSE:
self.type_sdprod(v['VECTEUR'],cham_no_sdaster)
return None
ASSEMBLAGE=MACRO(nom="ASSEMBLAGE",
op=OPS('Macro.assemblage_ops.assemblage_ops'),
UIinfo={"groupes":("Matrices et vecteurs",)},
sd_prod=assemblage_prod,
regles=(AU_MOINS_UN('MATR_ASSE','VECT_ASSE'),),
fr="Calcul des matrices et vecteurs assemblés ",
MODELE =SIMP(statut='o',typ=modele_sdaster),
CHAM_MATER =SIMP(statut='f',typ=cham_mater),
INST =SIMP(statut='f',typ='R',defaut=0.),
CARA_ELEM =SIMP(statut='f',typ=cara_elem),
CHARGE =SIMP(statut='f',typ=(char_meca,char_ther,char_acou),validators=NoRepeat(),max='**'),
CHAR_CINE =SIMP(statut='f',typ=(char_cine_meca,char_cine_ther,char_cine_acou) ),
NUME_DDL =SIMP(statut='o',typ=(nume_ddl_sdaster,CO)),
SOLVEUR =FACT(statut='d',
METHODE =SIMP(statut='f',typ='TXM',defaut="MULT_FRONT",into=("MULT_FRONT","LDLT","GCPC","MUMPS","PETSC") ),
b_mult_front =BLOC(condition="METHODE=='MULT_FRONT'",fr="Paramètres associés à la méthode multifrontale",
RENUM =SIMP(statut='f',typ='TXM',into=("MD","MDA","METIS"),defaut="METIS" ),
),
b_ldlt =BLOC(condition="METHODE=='LDLT'",fr="Paramètres associés à la méthode LDLT",
RENUM =SIMP(statut='f',typ='TXM',into=("RCMK","SANS"),defaut="RCMK" ),
),
b_mumps =BLOC(condition = "METHODE == 'MUMPS' ",fr="Paramètres de la méthode MUMPS",
RENUM =SIMP(statut='f',typ='TXM',defaut="AUTO",into=("AMD","AMF","PORD","METIS","QAMD","SCOTCH","AUTO")),
),
b_gcpc =BLOC(condition="METHODE=='GCPC'",fr="Paramètres associés à la méthode gradient conjugué",
RENUM =SIMP(statut='f',typ='TXM',into=("RCMK","SANS"),defaut="RCMK" ),
),
b_petsc =BLOC(condition = "METHODE == 'PETSC'",fr="Paramètres de la méthode PETSC",
RENUM =SIMP(statut='f',typ='TXM',defaut="RCMK",into=("SANS","RCMK") ),
),
),
MATR_ASSE =FACT(statut='f',max='**',
MATRICE =SIMP(statut='o',typ=CO),
OPTION =SIMP(statut='o',typ='TXM',
into=("RIGI_MECA","MASS_MECA","MASS_MECA_DIAG",
"AMOR_MECA","RIGI_MECA_HYST","IMPE_MECA",
"ONDE_FLUI","RIGI_FLUI_STRU","MASS_FLUI_STRU",
"RIGI_ROTA","RIGI_GEOM","MECA_GYRO","RIGI_GYRO",
"RIGI_THER","MASS_THER",
"RIGI_ACOU","MASS_ACOU","AMOR_ACOU",)
),
b_rigi_meca = BLOC( condition = "OPTION=='RIGI_MECA'",
MODE_FOURIER =SIMP(statut='f',typ='I',defaut= 0),
),
b_rigi_geom = BLOC( condition = "OPTION=='RIGI_GEOM'",
SIEF_ELGA =SIMP(statut='o',typ=cham_elem),
MODE_FOURIER =SIMP(statut='f',typ='I',defaut= 0),
),
b_rigi_ther = BLOC( condition = "OPTION=='RIGI_THER'",
MODE_FOURIER =SIMP(statut='f',typ='I',defaut= 0),
),
#
), # fin MATR_ASSE
#
VECT_ASSE =FACT(statut='f',max='**',
VECTEUR =SIMP(statut='o',typ=CO),
OPTION =SIMP(statut='o',typ='TXM',into=("CHAR_MECA","CHAR_ACOU","CHAR_THER") ),
b_char_meca =BLOC(condition = "OPTION == 'CHAR_MECA'", fr="chargement mécanique",
CHARGE =SIMP(statut='f',typ=char_meca,validators=NoRepeat(),max='**'),
MODE_FOURIER =SIMP(statut='f',typ='I',defaut= 0 ),
),
b_char_ther =BLOC(condition = "OPTION=='CHAR_THER'", fr="chargement thermique",
CHARGE =SIMP(statut='f',typ=char_ther,validators=NoRepeat(),max='**'),
),
b_char_acou =BLOC(condition = "OPTION=='CHAR_ACOU'", fr="chargement acoustique",
CHARGE =SIMP(statut='f',typ=char_acou,validators=NoRepeat(),max='**'),
),
#
), # fin VECT_ASSE
#
TITRE =SIMP(statut='f',typ='TXM',max='**'),
INFO =SIMP(statut='f',typ='I',defaut=1,into=(1,2)),
);
|
