1 #@ MODIF macro_elas_mult_ops Macro DATE 05/07/2005 AUTEUR DURAND C.DURAND
2 # -*- coding: iso-8859-1 -*-
3 # CONFIGURATION MANAGEMENT OF EDF VERSION
4 # ======================================================================
5 # COPYRIGHT (C) 1991 - 2002 EDF R&D WWW.CODE-ASTER.ORG
6 # THIS PROGRAM IS FREE SOFTWARE; YOU CAN REDISTRIBUTE IT AND/OR MODIFY
7 # IT UNDER THE TERMS OF THE GNU GENERAL PUBLIC LICENSE AS PUBLISHED BY
8 # THE FREE SOFTWARE FOUNDATION; EITHER VERSION 2 OF THE LICENSE, OR
9 # (AT YOUR OPTION) ANY LATER VERSION.
11 # THIS PROGRAM IS DISTRIBUTED IN THE HOPE THAT IT WILL BE USEFUL, BUT
12 # WITHOUT ANY WARRANTY; WITHOUT EVEN THE IMPLIED WARRANTY OF
13 # MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. SEE THE GNU
14 # GENERAL PUBLIC LICENSE FOR MORE DETAILS.
16 # YOU SHOULD HAVE RECEIVED A COPY OF THE GNU GENERAL PUBLIC LICENSE
17 # ALONG WITH THIS PROGRAM; IF NOT, WRITE TO EDF R&D CODE_ASTER,
18 # 1 AVENUE DU GENERAL DE GAULLE, 92141 CLAMART CEDEX, FRANCE.
19 # ======================================================================
23 def macro_elas_mult_ops(self,MODELE,CHAM_MATER,CARA_ELEM,NUME_DDL,
24 CHAR_MECA_GLOBAL,CHAR_CINE_GLOBAL,LIAISON_DISCRET,
25 CAS_CHARGE,SOLVEUR,**args):
27 Ecriture de la macro MACRO_ELAS_MULT
33 # On met le mot cle NUME_DDL dans une variable locale pour le proteger
35 # On importe les definitions des commandes a utiliser dans la macro
36 CALC_MATR_ELEM =self.get_cmd('CALC_MATR_ELEM')
37 NUME_DDL =self.get_cmd('NUME_DDL')
38 ASSE_MATRICE =self.get_cmd('ASSE_MATRICE')
39 FACT_LDLT =self.get_cmd('FACT_LDLT')
40 CALC_VECT_ELEM =self.get_cmd('CALC_VECT_ELEM')
41 ASSE_VECTEUR =self.get_cmd('ASSE_VECTEUR')
42 RESO_LDLT =self.get_cmd('RESO_LDLT')
43 CREA_RESU =self.get_cmd('CREA_RESU')
44 CALC_ELEM =self.get_cmd('CALC_ELEM')
45 CALC_NO =self.get_cmd('CALC_NO')
46 # La macro compte pour 1 dans la numerotation des commandes
49 # Le concept sortant (de type mult_elas ou fourier_elas) est nommé
50 # 'nomres' dans le contexte de la macro
52 self.DeclareOut('nomres',self.sd)
58 ielas=1 # mot clé NOM_CAS présent sous CAS_CHARGE
61 ifour=1 # mot clé MODE_FOURIER présent sous CAS_CHARGE
62 tyresu = 'FOURIER_ELAS'
63 if ielas==1 and ifour==1:
65 self.cr.fatal("""<F> <MACRO_ELAS_MULT> On ne peut avoir a la fois NOM_CAS et MODE_FOURIER""")
68 if (numeddl in self.sdprods) or (numeddl==None):
69 # Si le concept numeddl est dans self.sdprods ou n est pas nommé
70 # il doit etre produit par la macro
71 # il faudra donc appeler la commande NUME_DDL
78 if CHAR_MECA_GLOBAL: motscles['CHARGE'] =CHAR_MECA_GLOBAL
79 elif CHAR_CINE_GLOBAL: motscles['CHARGE'] =CHAR_CINE_GLOBAL
80 if CHAM_MATER : motscles['CHAM_MATER']=CHAM_MATER
81 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
82 __nomrig=CALC_MATR_ELEM(OPTION='RIGI_MECA',MODELE=MODELE,**motscles)
85 # On peut passer des mots cles egaux a None. Ils sont ignores
88 motscles['METHODE'] =SOLVEUR['METHODE']
89 motscles['RENUM'] =SOLVEUR['RENUM']
91 motscles['METHODE'] ='MULT_FRONT'
92 motscles['RENUM'] ='METIS'
94 self.DeclareOut('num',numeddl)
95 num=NUME_DDL(MATR_RIGI=__nomrig,**motscles)
97 _num=NUME_DDL(MATR_RIGI=__nomrig,**motscles)
102 __nomras=ASSE_MATRICE(MATR_ELEM=__nomrig,NUME_DDL=num)
104 __nomraf=FACT_LDLT(MATR_ASSE=__nomras,NPREC=SOLVEUR['NPREC'],STOP_SINGULIER=SOLVEUR['STOP_SINGULIER'])
106 #####################################################################
107 # boucle sur les items de CAS_CHARGE
116 if CHAR_MECA_GLOBAL: motscles['CHARGE'] =CHAR_MECA_GLOBAL
117 elif CHAR_CINE_GLOBAL: motscles['CHARGE'] =CHAR_CINE_GLOBAL
118 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
119 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
120 motscles['MODE_FOURIER'] =m['MODE_FOURIER']
121 __nomrig=CALC_MATR_ELEM(OPTION='RIGI_MECA',MODELE=MODELE,**motscles)
124 _num=NUME_DDL(MATR_RIGI=__nomrig,METHODE=SOLVEUR['METHODE'],RENUM=SOLVEUR['RENUM'])
128 __nomras=ASSE_MATRICE(MATR_ELEM=__nomrig,NUME_DDL=num)
130 __nomraf=FACT_LDLT(MATR_ASSE=__nomras,NPREC=SOLVEUR['NPREC'],STOP_SINGULIER=SOLVEUR['STOP_SINGULIER'])
133 if m['VECT_ASSE']==None :
135 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
136 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
137 if ifour : motscles['MODE_FOURIER'] =m['MODE_FOURIER']
138 if m['CHAR_MECA'] : motscles['CHARGE'] =m['CHAR_MECA']
139 elif m['CHAR_CINE'] : motscles['CHARGE'] =m['CHAR_CINE']
140 __nomvel=CALC_VECT_ELEM(OPTION='CHAR_MECA',**motscles)
141 __nomasv=ASSE_VECTEUR(VECT_ELEM=__nomvel,NUME_DDL=num)
143 __nomasv=m['VECT_ASSE']
146 __nomchn=RESO_LDLT(MATR_FACT=__nomraf,CHAM_NO=__nomasv,TITRE=m['SOUS_TITRE'])
147 nomchn.append(__nomchn)
149 # fin de la boucle sur les items de CAS_CHARGE
150 #####################################################################
155 if CHAM_MATER : motscle2['CHAM_MATER']=CHAM_MATER
156 if CARA_ELEM : motscle2['CARA_ELEM']=CARA_ELEM
157 if CHAM_MATER or CARA_ELEM :
161 motscles['AFFE'].append(_F(MODELE=MODELE,
162 CHAM_GD=nomchn[iocc],
163 NOM_CAS=m['NOM_CAS'],
169 motscles['AFFE'].append(_F(MODELE=MODELE,
170 CHAM_GD=nomchn[iocc],
171 NUME_MODE=m['MODE_FOURIER'],
172 TYPE_MODE=m['TYPE_MODE'],
179 motscles['AFFE'].append(_F(MODELE=MODELE,
180 CHAM_GD=nomchn[iocc],
181 NOM_CAS=m['NOM_CAS'],) )
186 motscles['AFFE'].append(_F(MODELE=MODELE,
187 CHAM_GD=nomchn[iocc],
188 NUME_MODE=m['MODE_FOURIER'],
189 TYPE_MODE=m['TYPE_MODE'],) )
193 nomres=CREA_RESU(OPERATION='AFFE',TYPE_RESU=tyresu,NOM_CHAM='DEPL',**motscles)
195 #####################################################################
196 # boucle sur les items de CAS_CHARGE pour CALC_ELEM ete CALC_NO
207 if type(m['OPTION'])==types.StringType:
208 if m['OPTION'] in ('FORC_NODA','REAC_NODA',
209 'EPSI_NOEU_DEPL','SIGM_NOEU_DEPL','EFGE_NOEU_DEPL',
210 'EQUI_NOEU_SIGM','EQUI_NOEU_EPSI','FLUX_NOEU_TEMP',):
212 liste_no.append(m['OPTION'])
215 liste_el.append(m['OPTION'])
217 for opt in m['OPTION']:
218 if opt in ('FORC_NODA','REAC_NODA',
219 'EPSI_NOEU_DEPL','SIGM_NOEU_DEPL','EFGE_NOEU_DEPL',
220 'EQUI_NOEU_SIGM','EQUI_NOEU_EPSI','FLUX_NOEU_TEMP',):
230 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
231 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
233 motscles['NOM_CAS']=m['NOM_CAS']
235 motscles['NUME_MODE']=m['MODE_FOURIER']
238 for chargt in m['CHAR_MECA'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
239 elif m['CHAR_CINE'] :
240 for chargt in m['CHAR_CINE'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
241 if CHAR_MECA_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_MECA_GLOBAL))
242 elif CHAR_CINE_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_CINE_GLOBAL))
243 CALC_ELEM(reuse=nomres,
246 NIVE_COUCHE=m['NIVE_COUCHE'],
247 NUME_COUCHE=m['NUME_COUCHE'],
248 OPTION=tuple(liste_el),
252 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
253 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
255 motscles['NOM_CAS']=m['NOM_CAS']
257 motscles['NUME_MODE']=m['MODE_FOURIER']
260 for chargt in m['CHAR_MECA'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
261 elif m['CHAR_CINE'] :
262 for chargt in m['CHAR_CINE'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
263 if CHAR_MECA_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_MECA_GLOBAL))
264 elif CHAR_CINE_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_CINE_GLOBAL))
265 CALC_NO(reuse=nomres,
268 OPTION=tuple(liste_no),
271 # fin de la boucle sur les items de CAS_CHARGE
272 #####################################################################