
DEBUT();

#直交異方性材料の設定
MA=DEFI_MATERIAU(ELAS_ORTH=_F(E_L=300.0,
                              E_T=300.0,
                              E_N=7500.0,
                              NU_LT=0.016,
                              NU_LN=0.016,
                              NU_TN=0.016,
                              G_LT=500.0,
                              G_LN=500.0,
                              G_TN=500.0,),);

MAIL=LIRE_MAILLAGE(FORMAT='MED',);

#MAIL=MODI_MAILLAGE(reuse =MAIL,
#                   MAILLAGE=MAIL,
#                   ORIE_PEAU_3D=_F(GROUP_MA='saika_m',),);

MODE=AFFE_MODELE(MAILLAGE=MAIL,
                 AFFE=_F(TOUT='OUI',
                         PHENOMENE='MECANIQUE',
                         MODELISATION='3D',),);

#局所座標系（オイラー角）の設定
CARA_EL=AFFE_CARA_ELEM(MODELE=MODE,
                       MASSIF=(_F(GROUP_MA='moto',
                                  ANGL_EULER=(0.0,0.0,0.0,),),
                               _F(GROUP_MA='alp330',
                                  ANGL_EULER=(0.0,330.0,270.0,),),
                               _F(GROUP_MA='bt30',
                                  ANGL_EULER=(90.0,30.0,270.0,),),),);

MATE=AFFE_MATERIAU(MAILLAGE=MAIL,
                   AFFE=_F(TOUT='OUI',
                           MATER=MA,),);

CHAR=AFFE_CHAR_MECA(MODELE=MODE,
                    DDL_IMPO=_F(GROUP_MA=('kotei_m','kotei_a','kotei_b',),
                                DX=0.0,
                                DY=0.0,
                                DZ=0.0,),
                    FORCE_FACE=(_F(GROUP_MA='saika_m',
                                   FY=0.1,),
                                _F(GROUP_MA='saika_a',
                                   FX=0.1,),
                                _F(GROUP_MA='saika_b',
                                   FY=0.1,),),);

RESU=MECA_STATIQUE(MODELE=MODE,
                   CHAM_MATER=MATE,
                   CARA_ELEM=CARA_EL,  #<--ここを追加
                   EXCIT=_F(CHARGE=CHAR,),);

RESU=CALC_CHAMP(reuse =RESU,
                RESULTAT=RESU,
                CONTRAINTE=('SIGM_ELNO','SIGM_NOEU',),
                CRITERES=('SIEQ_ELNO','SIEQ_NOEU',),);

IMPR_RESU(FORMAT='MED',
          UNITE=80,
          RESU=_F(RESULTAT=RESU,
                  NOM_CHAM=('SIGM_NOEU','SIEQ_NOEU','DEPL',),),);

FIN();
