forma03aをためしてみました

DOCUMENTによると２次元の要素C_PLANで線形材料を用いた解析です。

MESHは次のとおり。穴の開いた板の1/4の部分をモデル化していました。

commの内容です。

DEBUT(identifier='1:1',
      CODE=_F(NIV_PUB_WEB='INTERNET'),
      DEBUG=_F(SDVERI='OUI'))

# LECTURE DU MAILLAGE
MAIL = LIRE_MAILLAGE(identifier='3:1',
                     FORMAT='MED',
                     UNITE=20)

MAIL = MODI_MAILLAGE(identifier='4:1',
                     reuse=MAIL,
                     MAILLAGE=MAIL,
                     ORIE_PEAU=_F(GROUP_MA_PEAU='haut'))

# AFFECTATION DU MODELE (CONTRAINTES PLANES)
MODE = AFFE_MODELE(identifier='6:1',
                   AFFE=_F(MODELISATION='C_PLAN',
                           PHENOMENE='MECANIQUE',
                           TOUT='OUI'),
                   MAILLAGE=MAIL)

# DEFINITION DU MATERIAU
ACIER = DEFI_MATERIAU(identifier='8:1',
                      ELAS=_F(E=200000.0,
                              NU=0.3))

CHMATE = AFFE_MATERIAU(identifier='9:1',
                       AFFE=_F(MATER=ACIER,
                               TOUT='OUI'),
                       MAILLAGE=MAIL)

# RAMPE APPLICATION CHARGEMENT
RAMPE = DEFI_FONCTION(identifier='11:1',
                      NOM_PARA='INST',
                      VALE=(0.0, 0.0, 1000.0, 1000.0))

# DISCRETISATION EN TEMPS
L_INST = DEFI_LIST_REEL(identifier='13:1',
                        DEBUT=0.0,
                        INTERVALLE=_F(JUSQU_A=10.0,
                                      NOMBRE=1))

DEFLIST = DEFI_LIST_INST(identifier='14:1',
                         DEFI_LIST=_F(LIST_INST=L_INST))

# CONDITIONS DE SYMETRIE
CHSYME = AFFE_CHAR_CINE(identifier='16:1',
                        MECA_IMPO=(_F(DY=0.0,
                                      GROUP_MA='bas'),
                                   _F(DX=0.0,
                                      GROUP_MA='gauche')),
                        MODELE=MODE)

# CHARGEMENT EN FORCE
CHFORC = AFFE_CHAR_MECA(identifier='18:1',
                        FORCE_CONTOUR=_F(FY=1.0,
                                         GROUP_MA='haut'),
                        MODELE=MODE)

# CALCUL
RESU = STAT_NON_LINE(identifier='20:1',
                     CHAM_MATER=CHMATE,
                     COMPORTEMENT=_F(RELATION='ELAS'),
                     EXCIT=(_F(CHARGE=CHSYME),
                            _F(CHARGE=CHFORC,
                               FONC_MULT=RAMPE)),
                     INCREMENT=_F(INST_FIN=10.0,
                                  LIST_INST=DEFLIST),
                     MODELE=MODE)

# INTERPOLATION DES CONTRAINTES AUX NOEUDS
RESU = CALC_CHAMP(identifier='22:1',
                  reuse=RESU,
                  CONTRAINTE='SIGM_ELNO',
                  RESULTAT=RESU)

RESU = CALC_CHAMP(identifier='23:1',
                  reuse=RESU,
                  CONTRAINTE='SIGM_NOEU',
                  RESULTAT=RESU)

# INTERPOLATION DES CONTRAINTES EQUIVALENTES AUX NOEUDS
RESU = CALC_CHAMP(identifier='25:1',
                  reuse=RESU,
                  CRITERES=('SIEQ_ELGA', 'SIEQ_ELNO'),
                  RESULTAT=RESU)

RESU = CALC_CHAMP(identifier='26:1',
                  reuse=RESU,
                  CRITERES='SIEQ_NOEU',
                  RESULTAT=RESU)

# IMPRESSION DES RESULTATS
IMPR_RESU(identifier='28:1',
          FORMAT='MED',
          RESU=_F(RESULTAT=RESU),
          UNITE=80)

TEST_RESU(identifier='29:1',
          RESU=_F(CRITERE='RELATIF',
                  GROUP_NO='B',
                  NOM_CHAM='SIGM_NOEU',
                  NOM_CMP='SIYY',
                  NUME_ORDRE=1,
                  PRECISION=0.021,
                  REFERENCE='ANALYTIQUE',
                  RESULTAT=RESU,
                  VALE_CALC=30.6077958737,
                  VALE_REFE=30.0))

TEST_RESU(identifier='30:1',
          RESU=_F(CRITERE='RELATIF',
                  GROUP_NO='A',
                  NOM_CHAM='SIGM_NOEU',
                  NOM_CMP='SIXX',
                  NUME_ORDRE=1,
                  PRECISION=0.02,
                  REFERENCE='ANALYTIQUE',
                  RESULTAT=RESU,
                  VALE_CALC=-9.9967257536,
                  VALE_REFE=-10.0))

TEST_RESU(identifier='31:1',
          RESU=_F(CRITERE='RELATIF',
                  GROUP_NO='B',
                  NOM_CHAM='SIGM_NOEU',
                  NOM_CMP='SIYY',
                  NUME_ORDRE=1,
                  RESULTAT=RESU,
                  VALE_CALC=30.6077958737))

TEST_RESU(identifier='32:1',
          RESU=_F(CRITERE='RELATIF',
                  GROUP_NO='A',
                  NOM_CHAM='SIGM_NOEU',
                  NOM_CMP='SIXX',
                  NUME_ORDRE=1,
                  RESULTAT=RESU,
                  VALE_CALC=-9.9967257536))

FIN(identifier='33:1',
    )