PN : exemples de .comm

diff --git a/Editeur/Patrons/modes_3D.comm b/Editeur/Patrons/modes_3D.comm
new file mode 100644 (file)
index 0000000..aafe1f2
--- /dev/null
+++ b/Editeur/Patrons/modes_3D.comm
@@ -0,0 +1,57 @@
+
+# PATRON
+# ANALYSE MODALE
+# 3D
+
+DEBUT();
+# MISE EN DONNEES
+
+MAIL=LIRE_MAILLAGE();
+
+MODELE=AFFE_MODELE(MAILLAGE=MAIL,
+                   AFFE=_F(GROUP_MA='OBJET',
+                           PHENOMENE='MECANIQUE',
+                           MODELISATION='3D',),);
+# MATERIAU
+
+ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
+                            NU=0.3,
+                            RHO=7800.0,),);
+
+CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
+                    AFFE=_F(GROUP_MA='OBJET',
+                            MATER=ACIER,),);
+
+# CONDITIONS AUX LIMITES ET CHARGEMENT
+# CONDITIONS AUX LIMITES 
+
+BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
+                        DDL_IMPO=_F(GROUP_MA='ENCAST',
+                                    LIAISON='ENCASTRE',),);
+# CALCUL DES MATRICES ASSEMBLEES K ET M
+
+MACRO_MATR_ASSE(MODELE=MODELE,
+                CHAM_MATER=CHMAT,
+                CHARGE=BLOQUAGE,
+                NUME_DDL=CO('NUMEDDL'),
+                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
+                              OPTION='RIGI_MECA',),
+                           _F(MATRICE=CO('MASSE'),
+                              OPTION='MASS_MECA',),),);
+# CALCUL DES CINQ PREMIERS MODES
+
+MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
+                       MATR_B=MASSE,
+                       CALC_FREQ=_F(
+                       OPTION='PLUS_PETITE',
+                       NMAX_FREQ=5,),);
+# ECRITURE DES MODES AU FORMAT GMSH
+
+IMPR_RESU(MODELE=MODELE,
+          FORMAT='GMSH',
+          UNITE=19,
+          RESU=_F(MAILLAGE=MAIL,
+                  RESULTAT=MODES,
+                  NOM_CHAM='DEPL',),);
+
+FIN();

diff --git a/Editeur/Patrons/modes_pou.comm b/Editeur/Patrons/modes_pou.comm
new file mode 100644 (file)
index 0000000..fc72e23
--- /dev/null
+++ b/Editeur/Patrons/modes_pou.comm
@@ -0,0 +1,63 @@
+
+# PATRON
+# ANALYSE MODALE
+# POUTRE
+
+DEBUT();
+# MISE EN DONNEES
+
+MAIL=LIRE_MAILLAGE();
+
+MODELE=AFFE_MODELE(MAILLAGE=MAIL,
+                   AFFE=_F(GROUP_MA='OBJET',
+                           PHENOMENE='MECANIQUE',
+                           MODELISATION='3D',),);
+# MATERIAU
+
+ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
+                            NU=0.3,
+                            RHO=7800.0,),);
+
+CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
+                    AFFE=_F(GROUP_MA='OBJET',
+                            MATER=ACIER,),);
+
+CARELEM=AFFE_CARA_ELEM(MODELE=MODELE,
+                       POUTRE=_F(GROUP_MA='OBJET',
+                                 SECTION='CERCLE',
+                                 CARA=('R','EP',),
+                                 VALE=(1.0,0.01,),),);
+# CONDITIONS AUX LIMITES 
+
+BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
+                        DDL_IMPO=_F(GROUP_MA='ENCAST',
+                                    LIAISON='ENCASTRE',),);
+#MATRICES ASSEMBLEES K, M
+#
+
+MACRO_MATR_ASSE(MODELE=MODELE,
+                CHAM_MATER=CHMAT,
+                CARA_ELEM=CARELEM,
+                CHARGE=BLOQUAGE,
+                NUME_DDL=CO('NUMEDDL'),
+                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
+                              OPTION='RIGI_MECA',),
+                           _F(MATRICE=CO('MASSE'),
+                              OPTION='MASS_MECA',),),);
+#5 PREMIERS MODES
+
+MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
+                       MATR_B=MASSE,
+                       CALC_FREQ=_F(
+                       OPTION='PLUS_PETITE',
+                       NMAX_FREQ=5,),);
+# MODES AU FORMAT GMSH
+
+IMPR_RESU(MODELE=MODELE,
+          FORMAT='GMSH',
+          UNITE=19,
+          RESU=_F(MAILLAGE=MAIL,
+                  RESULTAT=MODES,
+                  NOM_CHAM='DEPL',),);
+
+FIN();

diff --git a/Editeur/Patrons/statlin_3D.comm b/Editeur/Patrons/statlin_3D.comm
new file mode 100644 (file)
index 0000000..39f7057
--- /dev/null
+++ b/Editeur/Patrons/statlin_3D.comm
@@ -0,0 +1,57 @@
+
+# PATRON 
+# ANALYSE STATIQUE
+# 3D
+
+DEBUT();
+# MISE EN DONNEES
+
+MAIL=LIRE_MAILLAGE();
+
+MAIL=DEFI_GROUP(reuse =MAIL,
+                MAILLAGE=MAIL,
+                CREA_GROUP_MA=_F(NOM='TOUT',
+                                 TOUT='OUI',),);
+
+MODMECA=AFFE_MODELE(MAILLAGE=MAIL,
+                    AFFE=_F(GROUP_MA='TOUT',
+                            PHENOMENE='MECANIQUE',
+                            MODELISATION='3D',),);
+
+ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E11,
+                            NU=0.3,),);
+
+CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
+                    AFFE=_F(TOUT='OUI',
+                            MATER=ACIER,),);
+# CONDITIONS AUX LIMITES ET CHARGEMENT
+#POINT A : ENCASTREMENT
+
+BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODMECA,
+                        DDL_IMPO=_F(GROUP_MA='ENCAST',
+                                    LIAISON='ENCASTRE',),);
+#POINT B : FORCE IMPOSEE
+
+CHARGT=AFFE_CHAR_MECA(MODELE=MODMECA,
+                      FORCE_NODALE=_F(GROUP_NO='PB',
+                                      FY=1000.0,),);
+# RESOLUTION DU PROBLEME ELASTIQUE
+
+RESU1=MECA_STATIQUE(MODELE=MODMECA,
+                    CHAM_MATER=CHMAT,
+                    EXCIT=(_F(CHARGE=BLOQUAGE,),
+                           _F(CHARGE=CHARGT,),),);
+# CALCUL DES EFFORTS
+
+RESU1=CALC_ELEM(reuse =RESU1,
+                RESULTAT=RESU1,
+                OPTION='EFGE_ELNO_DEPL',);
+# RESULTATS FORMAT GMSH
+
+IMPR_RESU(MODELE=MODMECA,
+          FORMAT='GMSH',
+          RESU=_F(MAILLAGE=MAIL,
+                  RESULTAT=RESU1,
+                  NOM_CHAM=('DEPL','EFGE_ELNO_DEPL',),),);
+
+FIN();

diff --git a/Editeur/Patrons/statlin_pou.comm b/Editeur/Patrons/statlin_pou.comm
new file mode 100644 (file)
index 0000000..0dd06eb
--- /dev/null
+++ b/Editeur/Patrons/statlin_pou.comm
@@ -0,0 +1,70 @@
+
+# PATRON 
+# ANALYSE STATIQUE
+# POUTRE
+
+DEBUT();
+# MISE EN DONNEES
+
+MAIL=LIRE_MAILLAGE();
+
+MAIL=DEFI_GROUP(reuse =MAIL,
+                MAILLAGE=MAIL,
+                CREA_GROUP_MA=_F(NOM='TOUT',
+                                 TOUT='OUI',),);
+
+MODMECA=AFFE_MODELE(MAILLAGE=MAIL,
+                    AFFE=_F(GROUP_MA='TOUT',
+                            PHENOMENE='MECANIQUE',
+                            MODELISATION='POU_D_T',),);
+
+ACIER=DEFI_MATERIAU(ELAS=_F(E=2.E10,
+                            NU=0.3,),);
+
+CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
+                    AFFE=_F(TOUT='OUI',
+                            MATER=ACIER,),);
+R1 = 1.0;
+
+EP1 = 0.001;
+
+
+CARA_POU=AFFE_CARA_ELEM(MODELE=MODMECA,
+                        POUTRE=_F(GROUP_MA='TOUT',
+                                  SECTION='CERCLE',
+                                  CARA=('R','EP',),
+                                  VALE=(R1,EP1,),),);
+
+# CONDITIONS AUX LIMITES ET CHARGEMENT
+
+#POINT A : ENCASTREMENT
+
+BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODMECA,
+                        DDL_IMPO=_F(GROUP_NO='PA',
+                                    LIAISON='ENCASTRE',),);
+#POINT B : FORCE IMPOSEE
+
+CHARGT=AFFE_CHAR_MECA(MODELE=MODMECA,
+                      FORCE_NODALE=_F(GROUP_NO='PB',
+                                      FY=1000.0,),);
+# RESOLUTION DU PROBLEME ELASTIQUE
+
+RESU1=MECA_STATIQUE(MODELE=MODMECA,
+                    CHAM_MATER=CHMAT,
+                    CARA_ELEM=CARA_POU,
+                    EXCIT=(_F(CHARGE=BLOQUAGE,),
+                           _F(CHARGE=CHARGT,),),);
+# CALCUL DES EFFORTS
+
+RESU1=CALC_ELEM(reuse =RESU1,
+                RESULTAT=RESU1,
+                OPTION='EFGE_ELNO_DEPL',);
+# RESULTATS FORMAT GMSH
+
+IMPR_RESU(MODELE=MODMECA,
+          FORMAT='GMSH',
+          RESU=_F(MAILLAGE=MAIL,
+                  RESULTAT=RESU1,
+                  NOM_CHAM=('DEPL','EFGE_ELNO_DEPL',),),);
+
+FIN();

diff --git a/Editeur/Patrons/tradir_3D.comm b/Editeur/Patrons/tradir_3D.comm
new file mode 100644 (file)
index 0000000..41d17e5
--- /dev/null
+++ b/Editeur/Patrons/tradir_3D.comm
@@ -0,0 +1,72 @@
+
+# PATRON
+# ANALYSE MODALE
+# 3D
+
+DEBUT();
+# MISE EN DONNEES
+
+MAIL=LIRE_MAILLAGE();
+
+MODELE=AFFE_MODELE(MAILLAGE=MAIL,
+                   AFFE=_F(GROUP_MA='OBJET',
+                           PHENOMENE='MECANIQUE',
+                           MODELISATION='3D',),);
+# MATERIAU
+
+ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
+                            NU=0.3,
+                            RHO=7800.0,),);
+
+CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
+                    AFFE=_F(GROUP_MA='OBJET',
+                            MATER=ACIER,),);
+# CONDITIONS AUX LIMITES ET CHARGEMENT
+# CONDITIONS AUX LIMITES 
+
+BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
+                        DDL_IMPO=_F(GROUP_MA='ENCAST',
+                                    LIAISON='ENCASTRE',),);
+# CALCUL DES MATRICES ASSEMBLEES K ET M
+
+MACRO_MATR_ASSE(MODELE=MODELE,
+                CHAM_MATER=CHMAT,
+                CHARGE=BLOQUAGE,
+                NUME_DDL=CO('NUMEDDL'),
+                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
+                              OPTION='RIGI_MECA',),
+                           _F(MATRICE=CO('MASSE'),
+                              OPTION='MASS_MECA',),),);
+# CALCUL DES CINQ PREMIERS MODES
+
+MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
+                       MATR_B=MASSE,
+                       CALC_FREQ=_F(
+                       OPTION='PLUS_PETITE',
+                       NMAX_FREQ=5,),);
+#CHARGEMENT
+
+EXCIT=AFFE_CHAR_MECA(MODELE=MODELE,
+                     PRES_REP=_F(GROUP_MA='PRESSE',
+                                 PRES=100.e5,),);
+
+EXCEL=CALC_VECT_ELEM(OPTION='CHAR_MECA',
+                     CHARGE=EXCIT,);
+
+EXCASS=ASSE_VECTEUR(VECT_ELEM=EXCEL,
+                    NUME_DDL=NUMEDDL,);
+OMEGA = 10.0;
+
+
+SINU = FORMULE(VALE='sin(OMEGA*INST)',
+               NOM_PARA='INST',);
+#CALCUL TRANSITOIRE
+
+TRADIR=DYNA_LINE_TRAN(MATR_MASS=MASSE,
+                      MATR_RIGI=RIGIDITE,
+                      NEWMARK=_F(),
+                      EXCIT=_F(VECT_ASSE=EXCASS,),
+                      INCREMENT=_F(PAS=.01,
+                                   INST_FIN=10.,),);
+
+FIN();

diff --git a/Editeur/Patrons/tramod_3D.comm b/Editeur/Patrons/tramod_3D.comm
new file mode 100644 (file)
index 0000000..9a68bc5
--- /dev/null
+++ b/Editeur/Patrons/tramod_3D.comm
@@ -0,0 +1,86 @@
+
+# PATRON
+# ANALYSE MODALE
+# 3D
+
+DEBUT();
+# MISE EN DONNEES
+
+MAIL=LIRE_MAILLAGE();
+
+MODELE=AFFE_MODELE(MAILLAGE=MAIL,
+                   AFFE=_F(GROUP_MA='OBJET',
+                           PHENOMENE='MECANIQUE',
+                           MODELISATION='3D',),);
+# MATERIAU
+
+ACIER=DEFI_MATERIAU(ELAS=_F(E=2.1E11,
+                            NU=0.3,
+                            RHO=7800.0,),);
+
+CHMAT=AFFE_MATERIAU(MAILLAGE=MAIL,
+                    AFFE=_F(GROUP_MA='OBJET',
+                            MATER=ACIER,),);
+# CONDITIONS AUX LIMITES ET CHARGEMENT
+# CONDITIONS AUX LIMITES 
+
+BLOQUAGE=AFFE_CHAR_MECA(MODELE=MODELE,
+                        DDL_IMPO=_F(GROUP_MA='ENCAST',
+                                    LIAISON='ENCASTRE',),);
+# CALCUL DES MATRICES ASSEMBLEES K ET M
+
+MACRO_MATR_ASSE(MODELE=MODELE,
+                CHAM_MATER=CHMAT,
+                CHARGE=BLOQUAGE,
+                NUME_DDL=CO('NUMEDDL'),
+                MATR_ASSE=(_F(MATRICE=CO('RIGIDITE'),
+                              OPTION='RIGI_MECA',),
+                           _F(MATRICE=CO('MASSE'),
+                              OPTION='MASS_MECA',),),);
+
+P=POST_ELEM(MASS_INER=_F(TOUT='OUI',
+                         GROUP_MA='M',),);
+# CALCUL DES CINQ PREMIERS MODES
+
+MODES=MODE_ITER_SIMULT(MATR_A=RIGIDITE,
+                       MATR_B=MASSE,
+                       CALC_FREQ=_F(
+                       OPTION='PLUS_PETITE',
+                       NMAX_FREQ=5,),);
+#CHARGEMENT
+
+EXCIT=AFFE_CHAR_MECA(MODELE=MODELE,
+                     PRES_REP=_F(GROUP_MA='PRESSE',
+                                 PRES=100.e5,),);
+
+EXCEL=CALC_VECT_ELEM(OPTION='CHAR_MECA',
+                     CHARGE=EXCIT,);
+
+EXCASS=ASSE_VECTEUR(VECT_ELEM=EXCEL,
+                    NUME_DDL=NUMEDDL,);
+OMEGA = 10.0;
+
+
+SINU = FORMULE(VALE='sin(OMEGA*INST)',
+               NOM_PARA='INST',);
+#PROJECTION
+
+MACRO_PROJ_BASE(BASE=MODES,
+                MATR_ASSE_GENE=(_F(MATRICE=CO('MASSPROJ'),
+                                   MATR_ASSE=MASSE,),
+                                _F(MATRICE=CO('RIGIPROJ'),
+                                   MATR_ASSE=RIGIDITE,),),
+                VECT_ASSE_GENE=_F(VECTEUR=CO('EXCPROJ'),
+                                  VECT_ASSE=EXCASS,),);
+#CALCUL TRANSITOIRE
+
+TRAMOD=DYNA_TRAN_MODAL(MASS_GENE=MASSPROJ,
+                       RIGI_GENE=RIGIPROJ,
+                       AMOR_REDUIT=(.02,.01,),
+                       INCREMENT=_F(INST_INIT=0.0,
+                                    INST_FIN=10.0,
+                                    PAS=.001,),
+                       EXCIT=_F(VECT_GENE=EXCPROJ,
+                                FONC_MULT=SINU,),);
+
+FIN();

diff --git a/Editeur/Patrons/vide.comm b/Editeur/Patrons/vide.comm
new file mode 100644 (file)
index 0000000..8340546
--- /dev/null
+++ b/Editeur/Patrons/vide.comm
@@ -0,0 +1,13 @@
+
+
+DEBUT(CODE=_F(NOM='VIDE',
+              NIV_PUB_WEB='INTERNET',),);
+#
+
+MA=LIRE_MAILLAGE();
+
+MO=AFFE_MODELE(MAILLAGE=MA,
+               AFFE=_F(PHENOMENE=None,),);
+
+FIN();
+#