test/data/meshio_meshes/med/README.md

   1 `cylinder.med` is first generated by salome 9.2.2. The mesh version is then modified by HDFView by changing the mesh version from 4.0.0 to 3.0.0 so that it can also be read in gmsh.
   2
   3 `box.med` is generated by code_aster 13.6 using the following command file. A specific displacement field is prescribed to an orthotropic hexahedral element, and we verify if meshio is able to read the current stress/strain/energy data.
   4
   5 ```
   6 DEBUT()
   7
   8 mesh = LIRE_MAILLAGE(FORMAT='MED', INFO=2, UNITE=20)
   9
  10 mesh = DEFI_GROUP(
  11     reuse=mesh, MAILLAGE=mesh,
  12     CREA_GROUP_NO=(_F(NOM='Left', NOEUD=('N1', 'N2', 'N3', 'N4')),
  13                    _F(NOM='Right', NOEUD=('N5', 'N6', 'N7', 'N8'))))
  14
  15 model = AFFE_MODELE(
  16     AFFE=_F(MODELISATION=('3D', ), PHENOMENE='MECANIQUE', TOUT='OUI'),
  17     MAILLAGE=mesh
  18 )
  19
  20 mat = DEFI_MATERIAU(
  21     ELAS_ORTH=_F(RHO=1.37486e-09, E_L=5329.16, E_T=2418.93, E_N=1746.66,
  22                  G_LT=1202.6, G_TN=515.638, G_LN=569.752,
  23                  NU_LT=0.450006, NU_TN=0.564862, NU_LN=0.397544,
  24                  ALPHA_L=1.0163e-05, ALPHA_T=2.67311e-05, ALPHA_N=6.86122e-05))
  25
  26 load = AFFE_CHAR_MECA(
  27     MODELE=model, DDL_IMPO=(_F(GROUP_NO='Left', DX=0, DY=0, DZ=0),
  28                             _F(GROUP_NO='Right', DX=1, DY=1, DZ=1))
  29 )
  30
  31 fieldmat = AFFE_MATERIAU(
  32     AFFE=_F(MATER=(mat, ), TOUT='OUI'),
  33     MODELE=model
  34 )
  35
  36 resu = MECA_STATIQUE(
  37     CHAM_MATER=fieldmat,
  38     EXCIT=_F(CHARGE=load),
  39     MODELE=model,
  40     OPTION='SANS',
  41     SOLVEUR=_F(METHODE='MUMPS')
  42 )
  43
  44 resu = CALC_CHAMP(
  45     reuse=resu, RESULTAT=resu, CONTRAINTE='SIEF_ELNO', DEFORMATION='EPSI_ELNO', ENERGIE=('ENEL_ELEM', 'ENEL_ELNO'))
  46
  47 IMPR_RESU(
  48     FORMAT='MED',
  49     RESU=_F(RESULTAT=resu),
  50     UNITE=80
  51 )
  52
  53 FIN()
  54 ```