Modelling 2D Beam

DEBUT();


#######################################################################
## The following lines read in the different meshes of the flap model
Beam=LIRE_MAILLAGE(UNITE=21,
					FORMAT='MED',
					NOM_MED='beam_160x8_2D',
					INFO_MED=1,);
					
#########################################################################
## Material properties
## The material properties are adjusted to the actual physical model

stiff=DEFI_MATERIAU(ELAS=_F(E=207000, #MPa
                            NU=0,
                            RHO=2.713e-9,),); #T/mm3

######################################################################
## Physical phenomenon & type of FE
hexaQ=AFFE_MODELE(MAILLAGE=Beam,
                  INFO=1,
                  AFFE=(_F(TOUT='OUI',
							PHENOMENE='MECANIQUE',
							MODELISATION='COQUE_3D',),
                        ),);
##########################################################################                  
cara1=AFFE_CARA_ELEM(	MODELE=hexaQ,
						COQUE=_F(GROUP_MA='beam',
								 EPAIS=5,
								 ),
								 );                        
#######################################################################						
## Allocate material properties to the mesh
MATE1=AFFE_MATERIAU(MAILLAGE=Beam,
                    AFFE=(_F(TOUT='OUI',
                             MATER=stiff,),
                         ),);  								
###########################################################################
## Define boundary conditions
char1=AFFE_CHAR_MECA(MODELE=hexaQ,
						DDL_IMPO=(_F(GROUP_MA='clamp',			
									DX=0.0,
                                    DY=0.0,
                                    DZ=0.0,
                                    DRX=0.0,
                                    DRY=0.0,),),
								
					INFO=1,);


#########################################################################                  
################################Definition of load : force imposed


load1=AFFE_CHAR_MECA(MODELE=hexaQ,
						FORCE_NODALE=(_F(GROUP_NO='force_n',
									 	FX=6468.75),
									 ),
					INFO=1,);

############################################################################
## Ramp

LREEL=DEFI_LIST_REEL(DEBUT=0.0,
                     INTERVALLE=_F(JUSQU_A=1.0,
                                   NOMBRE=100,),);

LINST=DEFI_LIST_INST(DEFI_LIST=_F(METHODE='MANUEL',
                                  LIST_INST=LREEL,
                                  SUBD_PAS=10,),
                     INFO=2,);

RAMPE=DEFI_FONCTION(
                    NOM_PARA='INST',
                    VALE=(0.0,0.0,
                          1.0,1.0,
                          ),
                    INTERPOL='LIN',);

# Solve system of equations

RESU=STAT_NON_LINE(MODELE=hexaQ,
                   CHAM_MATER=MATE1,
                   CARA_ELEM=cara1,
                   EXCIT=(_F(CHARGE=load1,
                          FONC_MULT=RAMPE,),
                          _F(CHARGE=char1,),),
                   COMPORTEMENT=_F(RELATION='ELAS',
                                   DEFORMATION='GROT_GDEP',
                                   GROUP_MA='beam',),
                   INCREMENT=_F(LIST_INST=LINST,),
                   NEWTON=_F(REAC_INCR=1,
                             PREDICTION='TANGENTE',
                             MATRICE='TANGENTE',
                             REAC_ITER=1,),
                   CONVERGENCE=_F(RESI_GLOB_MAXI=1e-6,
                                  ITER_GLOB_MAXI=30,),
                   SOLVEUR=_F(SYME='OUI',),);
###Write results on a file
IMPR_RESU(FORMAT='MED',
		  UNITE=81,
          RESU=_F(MAILLAGE=Beam,
                  RESULTAT=RESU,
                  NOM_CHAM=('DEPL'),
                  ),);

FIN(FORMAT_HDF='OUI',);