From 2b72f4636335ac3c3bd8c5eadb9733c29c7b3bd6 Mon Sep 17 00:00:00 2001 From: eficas <> Date: Tue, 15 Nov 2005 16:24:36 +0000 Subject: [PATCH] PN : exemples de .comm --- Editeur/Patrons/modes_3D.comm | 57 +++++++++++++++++++++ Editeur/Patrons/modes_pou.comm | 63 +++++++++++++++++++++++ Editeur/Patrons/statlin_3D.comm | 57 +++++++++++++++++++++ Editeur/Patrons/statlin_pou.comm | 70 ++++++++++++++++++++++++++ Editeur/Patrons/tradir_3D.comm | 72 ++++++++++++++++++++++++++ Editeur/Patrons/tramod_3D.comm | 86 ++++++++++++++++++++++++++++++++ Editeur/Patrons/vide.comm | 13 +++++ 7 files changed, 418 insertions(+) create mode 100644 Editeur/Patrons/modes_3D.comm create mode 100644 Editeur/Patrons/modes_pou.comm create mode 100644 Editeur/Patrons/statlin_3D.comm create mode 100644 Editeur/Patrons/statlin_pou.comm create mode 100644 Editeur/Patrons/tradir_3D.comm create mode 100644 Editeur/Patrons/tramod_3D.comm create mode 100644 Editeur/Patrons/vide.comm diff --git a/Editeur/Patrons/modes_3D.comm b/Editeur/Patrons/modes_3D.comm new file mode 100644 index 00000000..aafe1f2a --- /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 index 00000000..fc72e239 --- /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 index 00000000..39f7057f --- /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 index 00000000..0dd06eb4 --- /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 index 00000000..41d17e57 --- /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 index 00000000..9a68bc5d --- /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 index 00000000..8340546a --- /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(); +# -- 2.39.2