1 #@ MODIF macro_elas_mult_ops Macro DATE 14/09/2004 AUTEUR MCOURTOI M.COURTOIS
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
50 # Le concept sortant (de type mult_elas ou fourier_elas) est nommé
51 # 'nomres' dans le contexte de la macro
53 self.DeclareOut('nomres',self.sd)
59 ielas=1 # mot clé NOM_CAS présent sous CAS_CHARGE
62 ifour=1 # mot clé MODE_FOURIER présent sous CAS_CHARGE
63 tyresu = 'FOURIER_ELAS'
64 if ielas==1 and ifour==1:
66 self.cr.fatal("""<F> <MACRO_ELAS_MULT> On ne peut avoir a la fois NOM_CAS et MODE_FOURIER""")
69 if (numeddl in self.sdprods) or (numeddl==None):
70 # Si le concept numeddl est dans self.sdprods ou n est pas nommé
71 # il doit etre produit par la macro
72 # il faudra donc appeler la commande NUME_DDL
79 if CHAR_MECA_GLOBAL: motscles['CHARGE'] =CHAR_MECA_GLOBAL
80 elif CHAR_CINE_GLOBAL: motscles['CHARGE'] =CHAR_CINE_GLOBAL
81 if CHAM_MATER : motscles['CHAM_MATER']=CHAM_MATER
82 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
83 __nomrig=CALC_MATR_ELEM(OPTION='RIGI_MECA',MODELE=MODELE,**motscles)
86 # On peut passer des mots cles egaux a None. Ils sont ignores
89 motscles['METHODE'] =SOLVEUR['METHODE']
90 motscles['RENUM'] =SOLVEUR['RENUM']
92 motscles['METHODE'] ='MULT_FRONT'
93 motscles['RENUM'] ='METIS'
95 self.DeclareOut('num',numeddl)
96 num=NUME_DDL(MATR_RIGI=__nomrig,**motscles)
98 _num=NUME_DDL(MATR_RIGI=__nomrig,**motscles)
103 __nomras=ASSE_MATRICE(MATR_ELEM=__nomrig,NUME_DDL=num)
105 __nomraf=FACT_LDLT(MATR_ASSE=__nomras,NPREC=SOLVEUR['NPREC'],STOP_SINGULIER=SOLVEUR['STOP_SINGULIER'])
107 #####################################################################
108 # boucle sur les items de CAS_CHARGE
117 if CHAR_MECA_GLOBAL: motscles['CHARGE'] =CHAR_MECA_GLOBAL
118 elif CHAR_CINE_GLOBAL: motscles['CHARGE'] =CHAR_CINE_GLOBAL
119 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
120 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
121 motscles['MODE_FOURIER'] =m['MODE_FOURIER']
122 __nomrig=CALC_MATR_ELEM(OPTION='RIGI_MECA',MODELE=MODELE,**motscles)
125 _num=NUME_DDL(MATR_RIGI=__nomrig,METHODE=SOLVEUR['METHODE'],RENUM=SOLVEUR['RENUM'])
129 __nomras=ASSE_MATRICE(MATR_ELEM=__nomrig,NUME_DDL=num)
131 __nomraf=FACT_LDLT(MATR_ASSE=__nomras,NPREC=SOLVEUR['NPREC'],STOP_SINGULIER=SOLVEUR['STOP_SINGULIER'])
134 if m['VECT_ASSE']==None :
136 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
137 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
138 if ifour : motscles['MODE_FOURIER'] =m['MODE_FOURIER']
139 if m['CHAR_MECA'] : motscles['CHARGE'] =m['CHAR_MECA']
140 elif m['CHAR_CINE'] : motscles['CHARGE'] =m['CHAR_CINE']
141 __nomvel=CALC_VECT_ELEM(OPTION='CHAR_MECA',**motscles)
142 __nomasv=ASSE_VECTEUR(VECT_ELEM=__nomvel,NUME_DDL=num)
144 __nomasv=m['VECT_ASSE']
147 __nomchn=RESO_LDLT(MATR_FACT=__nomraf,CHAM_NO=__nomasv,TITRE=m['SOUS_TITRE'])
148 nomchn.append(__nomchn)
150 # fin de la boucle sur les items de CAS_CHARGE
151 #####################################################################
156 if CHAM_MATER : motscle2['CHAM_MATER']=CHAM_MATER
157 if CARA_ELEM : motscle2['CARA_ELEM']=CARA_ELEM
158 if CHAM_MATER or CARA_ELEM :
162 motscles['AFFE'].append(_F(MODELE=MODELE,
163 CHAM_GD=nomchn[iocc],
164 NOM_CAS=m['NOM_CAS'],
170 motscles['AFFE'].append(_F(MODELE=MODELE,
171 CHAM_GD=nomchn[iocc],
172 NUME_MODE=m['MODE_FOURIER'],
173 TYPE_MODE=m['TYPE_MODE'],
180 motscles['AFFE'].append(_F(MODELE=MODELE,
181 CHAM_GD=nomchn[iocc],
182 NOM_CAS=m['NOM_CAS'],) )
187 motscles['AFFE'].append(_F(MODELE=MODELE,
188 CHAM_GD=nomchn[iocc],
189 NUME_MODE=m['MODE_FOURIER'],
190 TYPE_MODE=m['TYPE_MODE'],) )
194 nomres=CREA_RESU(OPERATION='AFFE',TYPE_RESU=tyresu,NOM_CHAM='DEPL',**motscles)
196 #####################################################################
197 # boucle sur les items de CAS_CHARGE pour CALC_ELEM ete CALC_NO
208 if type(m['OPTION'])==types.StringType:
209 if m['OPTION'] in ('FORC_NODA','REAC_NODA',
210 'EPSI_NOEU_DEPL','SIGM_NOEU_DEPL','EFGE_NOEU_DEPL',
211 'EQUI_NOEU_SIGM','EQUI_NOEU_EPSI','FLUX_NOEU_TEMP',):
213 liste_no.append(m['OPTION'])
216 liste_el.append(m['OPTION'])
218 for opt in m['OPTION']:
219 if opt in ('FORC_NODA','REAC_NODA',
220 'EPSI_NOEU_DEPL','SIGM_NOEU_DEPL','EFGE_NOEU_DEPL',
221 'EQUI_NOEU_SIGM','EQUI_NOEU_EPSI','FLUX_NOEU_TEMP',):
231 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
232 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
234 motscles['NOM_CAS']=m['NOM_CAS']
236 motscles['NUME_MODE']=m['MODE_FOURIER']
239 for chargt in m['CHAR_MECA'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
240 elif m['CHAR_CINE'] :
241 for chargt in m['CHAR_CINE'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
242 if CHAR_MECA_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_MECA_GLOBAL))
243 elif CHAR_CINE_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_CINE_GLOBAL))
244 CALC_ELEM(reuse=nomres,
247 NIVE_COUCHE=m['NIVE_COUCHE'],
248 NUME_COUCHE=m['NUME_COUCHE'],
249 OPTION=tuple(liste_el),
253 if CHAM_MATER : motscles['CHAM_MATER'] =CHAM_MATER
254 if CARA_ELEM : motscles['CARA_ELEM'] =CARA_ELEM
256 motscles['NOM_CAS']=m['NOM_CAS']
258 motscles['NUME_MODE']=m['MODE_FOURIER']
261 for chargt in m['CHAR_MECA'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
262 elif m['CHAR_CINE'] :
263 for chargt in m['CHAR_CINE'] : motscles['EXCIT'].append(_F(CHARGE=chargt))
264 if CHAR_MECA_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_MECA_GLOBAL))
265 elif CHAR_CINE_GLOBAL: motscles['EXCIT'].append(_F(CHARGE=CHAR_CINE_GLOBAL))
266 CALC_NO(reuse=nomres,
269 OPTION=tuple(liste_no),
272 # fin de la boucle sur les items de CAS_CHARGE
273 #####################################################################