commit 8a75033fbfe7a5fe3fd3bef6f9f3a2aaea3b6270 Author: Denis V. Dedkov Date: Sat Jun 7 19:45:53 2014 +0600 Initial commit diff --git a/d1.comm b/d1.comm new file mode 100644 index 0000000..4361441 --- /dev/null +++ b/d1.comm @@ -0,0 +1,151 @@ +#No fiber + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + _F(GROUP_MA='Fiber', + MATER=empty,), + # _F(GROUP_MA='fs1', + # MATER=empty,), + # _F(GROUP_MA='fs2', + # MATER=mat,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_ELEM(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + OPTION=('SIGM_ELNO','EPSI_ELNO','SIEF_ELNO','SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_NO(reuse =RESU, + RESULTAT=RESU, + OPTION=('EPSI_NOEU','SIEF_NOEU','SIGM_NOEU','SIEQ_NOEU',),); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=28, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fiber', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=29, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fiber', + IMPR_COOR='OUI',),); + +FIN(); \ No newline at end of file diff --git a/d2.comm b/d2.comm new file mode 100644 index 0000000..34bdb41 --- /dev/null +++ b/d2.comm @@ -0,0 +1,156 @@ +#No fiber + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + _F(GROUP_MA='Fiber', + MATER=mat,), + # _F(GROUP_MA='fs1', + # MATER=empty,), + # _F(GROUP_MA='fs2', + # MATER=mat,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_CHAMP(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + CONTRAINTE=('SIGM_ELNO','SIEF_ELNO',), + DEFORMATION=('EPME_ELNO',), + CRITERES=('SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_CHAMP(reuse =RESU, + RESULTAT=RESU, + CONTRAINTE=('SIEF_NOEU','SIGM_NOEU',), + DEFORMATION=('EPME_NOEU',), + CRITERES=('SIEQ_NOEU',), + ); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=28, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fiber', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=29, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fiber', + IMPR_COOR='OUI',),); + +FIN(); \ No newline at end of file diff --git a/d3.comm b/d3.comm new file mode 100644 index 0000000..def1c5d --- /dev/null +++ b/d3.comm @@ -0,0 +1,151 @@ +#Fiber break + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + # _F(GROUP_MA='Fiber', + # MATER=empty,), + _F(GROUP_MA='fs1', + MATER=empty,), + # _F(GROUP_MA='fs2', + # MATER=mat,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_ELEM(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + OPTION=('SIGM_ELNO','EPSI_ELNO','SIEF_ELNO','SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_NO(reuse =RESU, + RESULTAT=RESU, + OPTION=('EPSI_NOEU','SIEF_NOEU','SIGM_NOEU','SIEQ_NOEU',),); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=28, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=29, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +FIN(); \ No newline at end of file diff --git a/d4.comm b/d4.comm new file mode 100644 index 0000000..2c0ad33 --- /dev/null +++ b/d4.comm @@ -0,0 +1,151 @@ +#Fiber break + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + # _F(GROUP_MA='Fiber', + # MATER=empty,), + _F(GROUP_MA='fs1', + MATER=mat,), + # _F(GROUP_MA='fs2', + # MATER=mat,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_ELEM(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + OPTION=('SIGM_ELNO','EPSI_ELNO','SIEF_ELNO','SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_NO(reuse =RESU, + RESULTAT=RESU, + OPTION=('EPSI_NOEU','SIEF_NOEU','SIGM_NOEU','SIEQ_NOEU',),); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=28, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=29, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +FIN(); \ No newline at end of file diff --git a/d5.comm b/d5.comm new file mode 100644 index 0000000..bcd0f2c --- /dev/null +++ b/d5.comm @@ -0,0 +1,168 @@ +#Two fibers break + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + # _F(GROUP_MA='Fiber', + # MATER=empty,), + _F(GROUP_MA='fs1', + MATER=empty,), + _F(GROUP_MA='fs2', + MATER=empty,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_ELEM(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + OPTION=('SIGM_ELNO','EPSI_ELNO','SIEF_ELNO','SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_NO(reuse =RESU, + RESULTAT=RESU, + OPTION=('EPSI_NOEU','SIEF_NOEU','SIGM_NOEU','SIEQ_NOEU',),); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=28, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=29, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs2', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=30, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=31, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs2', + IMPR_COOR='OUI',),); + + +FIN(); \ No newline at end of file diff --git a/d6.comm b/d6.comm new file mode 100644 index 0000000..dd83b3b --- /dev/null +++ b/d6.comm @@ -0,0 +1,166 @@ +#Two fibers break + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + # _F(GROUP_MA='Fiber', + # MATER=empty,), + _F(GROUP_MA='fs1', + MATER=mat,), + _F(GROUP_MA='fs2', + MATER=mat,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_ELEM(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + OPTION=('SIGM_ELNO','EPSI_ELNO','SIEF_ELNO','SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_NO(reuse =RESU, + RESULTAT=RESU, + OPTION=('EPSI_NOEU','SIEF_NOEU','SIGM_NOEU','SIEQ_NOEU',),); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=28, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=29, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs2', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=30, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs1', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=31, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='fs2', + IMPR_COOR='OUI',),); + +FIN(); \ No newline at end of file diff --git a/regular.comm b/regular.comm new file mode 100644 index 0000000..ac0d6ce --- /dev/null +++ b/regular.comm @@ -0,0 +1,133 @@ + +DEBUT(); + +import sys +sys.path.append('/home/denis/Documents/disser/model/v4_w_friction') +from materials import * + +empty=DEFI_MATERIAU(ELAS=_F(E=m_empty['E'], + NU=m_empty['NU'],),); + +mat=DEFI_MATERIAU(ELAS=_F(E=m_matrix['E'], + NU=m_matrix['NU'],),); + +fib=DEFI_MATERIAU(ELAS=_F(E=m_fibers['E'], + NU=m_fibers['NU'],),); + +MAIL=LIRE_MAILLAGE(FORMAT='MED',); + +MAIL=MODI_MAILLAGE(reuse =MAIL, + MAILLAGE=MAIL, + ORIE_PEAU_3D=_F(GROUP_MA=('ZY_Side','XY_Side',),),); + +MODE=AFFE_MODELE(MAILLAGE=MAIL, + AFFE=_F(TOUT='OUI', + PHENOMENE='MECANIQUE', + MODELISATION='3D',),); + +#Cont=DEFI_CONTACT(MODELE=MODE, +# FORMULATION='DISCRETE', +# ZONE=_F(GROUP_MA_MAIT='Matrix', +# GROUP_MA_ESCL='fShell', +# ALGO_CONT='CONTRAINTE',),); + +MATE=AFFE_MATERIAU(MAILLAGE=MAIL, + AFFE=(_F(GROUP_MA='Matrix', + MATER=mat,), + _F(GROUP_MA='Fibers', + MATER=fib,), + # _F(GROUP_MA='Fiber', + # MATER=empty,), + # _F(GROUP_MA='fs1', + # MATER=empty,), + # _F(GROUP_MA='fs2', + # MATER=mat,), + ), + ); + +CHAR=AFFE_CHAR_MECA(MODELE=MODE, + DDL_IMPO=( + _F(GROUP_MA='YX_Side', + DZ=0.0,), + _F(GROUP_MA='YZ_Side', + DX=0.0,), + _F(GROUP_MA='ZX_Bott', + DY=0,), + _F(GROUP_MA='ZY_Side', + DX=-1.0), + _F(GROUP_MA='XY_Side', + DZ=1.0), + ), + LIAISON_MAIL=_F(GROUP_MA_MAIT='Matrix', + GROUP_MA_ESCL='fShell', + TYPE_RACCORD='MASSIF',), + ); + +RESU=MECA_STATIQUE(MODELE=MODE, + CHAM_MATER=MATE, + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_ELEM(reuse =RESU, + MODELE=MODE, + CHAM_MATER=MATE, + RESULTAT=RESU, + OPTION=('SIGM_ELNO','EPSI_ELNO','SIEF_ELNO','SIEQ_ELNO',), + EXCIT=_F(CHARGE=CHAR,),); + +RESU=CALC_NO(reuse =RESU, + RESULTAT=RESU, + OPTION=('EPSI_NOEU','SIEF_NOEU','SIGM_NOEU','SIEQ_NOEU',),); + +IMPR_RESU(FORMAT='MED', + UNITE=80, + RESU=_F(MAILLAGE=MAIL, + RESULTAT=RESU, + NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','SIGM_ELNO','DEPL',),),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=22, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=23, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Matrix', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=24, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=25, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + GROUP_MA='Fibers', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=26, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIGM_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +IMPR_RESU(FORMAT='RESULTAT', + UNITE=27, + RESU=_F(RESULTAT=RESU, + NOM_CHAM='SIEQ_NOEU', + FORM_TABL='OUI', + IMPR_COOR='OUI',),); + +FIN(); \ No newline at end of file